cropCSM: designing safe and potent herbicides with graph-based signatures.

Herbicides have revolutionised weed management, increased crop yields and improved profitability allowing for an increase in worldwide food security. Their widespread use, however, has also led to a rise in resistance and concerns about their environmental impact. Despite the need for potent and safe herbicidal molecules, no herbicide with a new mode of action has reached the market in 30 years. Although development of computational approaches has proven invaluable to guide rational drug discovery pipelines, leading to higher hit rates and lower attrition due to poor toxicity, little has been done in contrast for herbicide design. To fill this gap, we have developed cropCSM, a computational platform to help identify new, potent, nontoxic and environmentally safe herbicides. By using a knowledge-based approach, we identified physicochemical properties and substructures enriched in safe herbicides. By representing the small molecules as a graph, we leveraged these insights to guide the development of predictive models trained and tested on the largest collected data set of molecules with experimentally characterised herbicidal profiles to date (over 4500 compounds). In addition, we developed six new environmental and human toxicity predictors, spanning five different species to assist in molecule prioritisation. cropCSM was able to correctly identify 97% of herbicides currently available commercially, while predicting toxicity profiles with accuracies of up to 92%. We believe cropCSM will be an essential tool for the enrichment of screening libraries and to guide the development of potent and safe herbicides. We have made the method freely available through a user-friendly webserver at http://biosig.unimelb.edu.au/crop_csm.

Structural and biochemical analyses of concanavalin A circular permutation by jack bean asparaginyl endopeptidase.

Over 30 years ago, an intriguing posttranslational modification was found responsible for creating concanavalin A (conA), a carbohydrate-binding protein from jack bean (Canavalia ensiformis) seeds and a common carbohydrate chromatography reagent. ConA biosynthesis involves what was then an unprecedented rearrangement in amino-acid sequence, whereby the N-terminal half of the gene-encoded conA precursor (pro-conA) is swapped to become the C-terminal half of conA. Asparaginyl endopeptidase (AEP) was shown to be involved, but its mechanism was not fully elucidated. To understand the structural basis and consequences of circular permutation, we generated recombinant jack bean pro-conA plus jack bean AEP (CeAEP1) and solved crystal structures for each to 2.1 and 2.7 Å, respectively. By reconstituting conA biosynthesis in vitro, we prove CeAEP1 alone can perform both cleavage and cleavage-coupled transpeptidation to form conA. CeAEP1 structural analysis reveals how it is capable of carrying out both reactions. Biophysical assays illustrated that pro-conA is less stable than conA. This observation was explained by fewer intermolecular interactions between subunits in the pro-conA crystal structure and consistent with a difference in the prevalence for tetramerization in solution. These findings elucidate the consequences of circular permutation in the only posttranslation example known to occur in nature.

The genetic origin of evolidine, the first cyclopeptide discovered in plants, and related orbitides.

Cyclic peptides are reported to have antibacterial, antifungal, and other bioactivities. Orbitides are a class of cyclic peptides that are small, head-to-tail cyclized, composed of proteinogenic amino acids and lack disulfide bonds; they are also known in several genera of the plant family Rutaceae. Melicope xanthoxyloides is the Australian rain forest tree of the Rutaceae family in which evolidine, the first plant cyclic peptide, was discovered. Evolidine (cyclo-SFLPVNL) has subsequently been all but forgotten in the academic literature, so to redress this we used tandem MS and de novo transcriptomics to rediscover evolidine and decipher its biosynthetic origin from a short precursor just 48 residues in length. We also identified another six M. xanthoxyloides orbitides using the same techniques. These peptides have atypically diverse C termini consisting of residues not recognized by either of the known proteases plants use to macrocyclize peptides, suggesting new cyclizing enzymes await discovery. We examined the structure of two of the novel orbitides by NMR, finding one had a definable structure, whereas the other did not. Mining RNA-seq and whole genome sequencing data from other species of the Rutaceae family revealed that a large and diverse family of peptides is encoded by similar sequences across the family and demonstrates how powerful de novo transcriptomics can be at accelerating the discovery of new peptide families.

Plant asparaginyl endopeptidases and their structural determinants of function.

Asparaginyl endopeptidases (AEPs) are versatile enzymes that in biological systems are involved in producing three different catalytic outcomes for proteins, namely (i) routine cleavage by bond hydrolysis, (ii) peptide maturation, including macrocyclisation by a cleavage-coupled intramolecular transpeptidation and (iii) circular permutation involving separate cleavage and transpeptidation reactions resulting in a major reshuffling of protein sequence. AEPs differ in their preference for cleavage or transpeptidation reactions, catalytic efficiency, and preference for asparagine or aspartate target residues. We look at structural analyses of various AEPs that have laid the groundwork for identifying important determinants of AEP function in recent years, with much of the research impetus arising from the potential biotechnological and pharmaceutical applications.

Structural Characterization of the PawL-Derived Peptide Family, an Ancient Subfamily of Orbitides.

Plants are an excellent source of bioactive peptides, often with disulfide bonds and/or a cyclic backbone. While focus has predominantly been directed at disulfide-rich peptides, a large family of small, cyclic plant peptides lacking disulfide bonds, known as orbitides, has been relatively ignored. A recently discovered subfamily of orbitides is the PawL-derived peptides (PLPs), produced during the maturation of precursors for seed storage albumins. Although their evolutionary origins have been dated, in-depth exploration of the family's structural characteristics and potential bioactivities remains to be conducted. Here we present an extensive and systematic characterization of the PLP family. Nine PLPs were chosen and prepared by solid phase peptide synthesis. Their structural features were studied using solution NMR spectroscopy, and seven were found to possess regions of backbone order. Ordered regions consist of ß-turns, with some PLPs adopting two well-defined ß-turns within sequences as short as seven residues, which are largely the result of side chain interactions. Our data highlight that the sequence diversity within this family results in equally diverse structures. None of these nine PLPs showed antibacterial or antifungal activity.

Physical clustering of FLC alleles during Polycomb-mediated epigenetic silencing in vernalization.

Vernalization, the promotion of flowering by cold, involves Polycomb-mediated epigenetic silencing of FLOWERING LOCUS C (FLC). Cold progressively promotes cell-autonomous switching to a silenced state. Here, we used live-cell imaging of FLC-lacO to monitor changes in nuclear organization during vernalization. FLC-lacO alleles physically cluster during the cold and generally remain so after plants are returned to warm. Clustering is dependent on the Polycomb trans-factors necessary for establishment of the FLC silenced state but not on LIKE HETEROCHROMATIN PROTEIN 1, which functions to maintain silencing. These data support the view that physical clustering may be a common feature of Polycomb-mediated epigenetic switching mechanisms.

The macrocyclizing protease butelase 1 remains autocatalytic and reveals the structural basis for ligase activity.

Plant asparaginyl endopeptidases (AEPs) are expressed as inactive zymogens that perform maturation of seed storage protein upon cleavage-dependent autoactivation in the low-pH environment of storage vacuoles. The AEPs have attracted attention for their macrocyclization reactions, and have been classified as cleavage or ligation specialists. However, we have recently shown that the ability of AEPs to produce either cyclic or acyclic products can be altered by mutations to the active site region, and that several AEPs are capable of macrocyclization given favorable pH conditions. One AEP extracted from Clitoria ternatea seeds (butelase 1) is classified as a ligase rather than a protease, presenting an opportunity to test for loss of cleavage activity. Here, making recombinant butelase 1 and rescuing an Arabidopsis thaliana mutant lacking AEP, we show that butelase 1 retains cleavage functions in vitro and in vivo. The in vivo rescue was incomplete, consistent with some trade-off for butelase 1 specialization toward macrocyclization. Its crystal structure showed an active site with only subtle differences from cleaving AEPs, suggesting the many differences in its peptide-binding region are the source of its efficient macrocyclization. All considered, it seems that either butelase 1 has not fully specialized or a requirement for autocatalytic cleavage is an evolutionary constraint upon macrocyclizing AEPs.

Evidence for Ancient Origins of Bowman-Birk Inhibitors from Selaginella moellendorffii.

Bowman-Birk Inhibitors (BBIs) are a well-known family of plant protease inhibitors first described 70 years ago. BBIs are known only in the legume (Fabaceae) and cereal (Poaceae) families, but peptides that mimic their trypsin-inhibitory loops exist in sunflowers (Helianthus annuus) and frogs. The disparate biosynthetic origins and distant phylogenetic distribution implies these loops evolved independently, but their structural similarity suggests a common ancestor. Targeted bioinformatic searches for the BBI inhibitory loop discovered highly divergent BBI-like sequences in the seedless, vascular spikemoss Selaginella moellendorffii Using de novo transcriptomics, we confirmed expression of five transcripts in S. moellendorffii whose encoded proteins share homology with BBI inhibitory loops. The most highly expressed, BBI3, encodes a protein that inhibits trypsin. We needed to mutate two lysine residues to abolish trypsin inhibition, suggesting BBI3's mechanism of double-headed inhibition is shared with BBIs from angiosperms. As Selaginella belongs to the lycopod plant lineage, which diverged ∼200 to 230 million years before the common ancestor of angiosperms, its BBI-like proteins imply there was a common ancestor for legume and cereal BBIs. Indeed, we discovered BBI sequences in six angiosperm families outside the Fabaceae and Poaceae. These findings provide the evolutionary missing links between the well-known legume and cereal BBI gene families.

Cyclic Peptides in Seed of Annona muricata Are Ribosomally Synthesized.

Small, cyclic peptides are reported to have many bioactivities. In bacteria and fungi, they can be made by nonribosomal peptide synthetases, but in plants they are exclusively ribosomal. Cyclic peptides from the Annona genus possess cytotoxic and anti-inflammatory activities, but their biosynthesis is unknown. The medicinal soursop plant, Annona muricata, contains annomuricatins A (cyclo-PGFVSA) and B (cyclo-PNAWLGT). Here, using de novo transcriptomics and tandem mass spectrometry, we identify a suite of short transcripts for precursor proteins for 10 validated annomuricatins, 9 of which are novel. In their precursors, annomuricatins are preceded by an absolutely conserved Glu and each peptide sequence has a conserved proto-C-terminal Pro, revealing parallels with the segetalin orbitides from the seed of Vaccaria hispanica, which are processed through ligation by a prolyl oligopeptidase in a transpeptidation reaction.

Defining the Familial Fold of the Vicilin-Buried Peptide Family.

Plants and their seeds have been shown to be a rich source of cystine-stabilized peptides. Recently a new family of plant seed peptides whose sequences are buried within precursors for seed storage vicilins was identified. Members of this Vicilin-Buried Peptide (VBP) family are found in distantly related plant species including the monocot date palm, as well as dicotyledonous species like pumpkin and sesame. Genetic evidence for their widespread occurrence indicates that they are of ancient origin. Limited structural studies have been conducted on VBP family members, but two members have been shown to adopt a helical hairpin fold. We present an extensive characterization of VBPs using solution NMR spectroscopy, to better understand their structural features. Four peptides were produced by solid phase peptide synthesis and shown to favor a helix-loop-helix hairpin fold, as a result of the I-IV/II-III ladderlike connectivity of their disulfide bonds. Interhelical interactions, including hydrophobic contacts and salt bridges, are critical for the fold stability and control the angle at which the antiparallel α-helices interface. Activities reported for VBPs include trypsin inhibitory activity and inhibition of ribosomal function; however, their diverse structural features despite a common fold suggest that additional bioactivities yet to be revealed are likely.

Sequencing Orbitides by Acid-Mediated Ring Cleavage Followed by Tandem Mass Spectrometry.

Homodetic cyclic peptides have aroused interest because of their pharmacological potential. Sequencing cyclic peptides is difficult-Edman degradation is not possible as there is no N-terminus, NMR requires quantities that are hard to gather from native samples, and tandem mass spectrometry data are difficult to interpret due to the peptide ring opening at multiple points. Sequencing can be simplified by cleaving the peptide ring at a specific peptide bond. Partial acid hydrolysis is a possible solution, but to date sequencing by this method has only been demonstrated for linear peptides and cyclotides, which are larger cyclic peptides (∼30 amino acids) with three disulfide bonds. This study tests whether partial acid hydrolysis could be used to aid sequencing of Cys-less cyclic peptides with fewer than ten amino acid residues. We show that, with the right combination of temperature and acid, ring cleavage occurs and offers relatively simple MS/MS spectra amenable to sequencing. We describe how this method was used in our recent study in which we sequenced annomuricatin D (cyclo-GHSIFPPIP) from seeds of the soursop, Annona muricata. We found that orbitides can be linearized for MS/MS sequencing by incubation with 1.2 M HCl at 90 °C for 15-20 min. This fast, economical sequencing method will be useful to those studying small cyclic peptides lacking disulfide bonds, which are commonly found in many organisms, especially plants.

Targeting plant DIHYDROFOLATE REDUCTASE with antifolates and mechanisms for genetic resistance.

The folate biosynthetic pathway and its key enzyme dihydrofolate reductase (DHFR) is a popular target for drug development due to its essential role in the synthesis of DNA precursors and some amino acids. Despite its importance, little is known about plant DHFRs, which, like the enzymes from the malarial parasite Plasmodium, are bifunctional, possessing DHFR and thymidylate synthase (TS) domains. Here using genetic knockout lines we confirmed that either DHFR-TS1 or DHFR-TS2 (but not DHFR-TS3) was essential for seed development. Screening mutated Arabidopsis thaliana seeds for resistance to antimalarial DHFR-inhibitor drugs pyrimethamine and cycloguanil identified causal lesions in DHFR-TS1 and DHFR-TS2, respectively, near the predicted substrate-binding site. The different drug resistance profiles for the plants, enabled by the G137D mutation in DHFR-TS1 and the A71V mutation in DHFR-TS2, were consistent with biochemical studies using recombinant proteins and could be explained by structural models. These findings provide a great improvement in our understanding of plant DHFR-TS and suggest how plant-specific inhibitors might be developed, as DHFR is not currently targeted by commercial herbicides.

An Orbitide from Ratibida columnifera Seed Containing 16 Amino Acid Residues.

Cyclic peptides are abundant in plants and have attracted interest due to their bioactivity and potential as drug scaffolds. Orbitides are head-to-tail cyclic peptides that are ribosomally synthesized, post-translationally modified, and lack disulfide bonds. All known orbitides contain 5-12 amino acid residues. Here we describe PLP-53, a novel orbitide from the seed of Ratibida columnifera. PLP-53 consists of 16 amino acids, four residues larger than any known orbitide. NMR structural studies showed that, compared to previously characterized orbitides, PLP-53 is more flexible and, under the studied conditions, did not adopt a single ordered conformation based on analysis of NOEs and chemical shifts.

Two proteins for the price of one: Structural studies of the dual-destiny protein preproalbumin with sunflower trypsin inhibitor-1.

Seed storage proteins are both an important source of nutrition for humans and essential for seedling establishment. Interestingly, unusual napin-type 2S seed storage albumin precursors in sunflowers contain a sequence that is released as a macrocyclic peptide during post-translational processing. The mechanism by which such peptides emerge from linear precursor proteins has received increased attention; however, the structural characterization of intact precursor proteins has been limited. Here, we report the 3D NMR structure of the Helianthus annuus PawS1 (preproalbumin with sunflower trypsin inhibitor-1) and provide new insights into the processing of this remarkable dual-destiny protein. In seeds, PawS1 is matured by asparaginyl endopeptidases (AEPs) into the cyclic peptide SFTI-1 (sunflower trypsin inhibitor-1) and a heterodimeric 2S albumin. The structure of PawS1 revealed that SFTI-1 and the albumin are independently folded into well-defined domains separated by a flexible linker. PawS1 was cleaved in vitro with recombinant sunflower HaAEP1 and in situ using a sunflower seed extract in a way that resembled the expected in vivo cleavages. Recombinant HaAEP1 cleaved PawS1 at multiple positions, and in situ, its flexible linker was removed, yielding fully mature heterodimeric albumin. Liberation and cyclization of SFTI-1, however, was inefficient, suggesting that specific seed conditions or components may be required for in vivo biosynthesis of SFTI-1. In summary, this study has revealed the 3D structure of a macrocyclic precursor protein and provided important mechanistic insights into the maturation of sunflower proalbumins into an albumin and a macrocyclic peptide.

Stepwise Evolution of a Buried Inhibitor Peptide over 45 My.

The de novo evolution of genes and the novel proteins they encode has stimulated much interest in the contribution such innovations make to the diversity of life. Most research on this de novo evolution focuses on transcripts, so studies on the biochemical steps that can enable completely new proteins to evolve and the time required to do so have been lacking. Sunflower Preproalbumin with SFTI-1 (PawS1) is an unusual albumin precursor because in addition to producing albumin it also yields a potent, bicyclic protease-inhibitor called SunFlower Trypsin Inhibitor-1 (SFTI-1). Here, we show how this inhibitor peptide evolved stepwise over tens of millions of years. To trace the origin of the inhibitor peptide SFTI-1, we assembled seed transcriptomes for 110 sunflower relatives whose evolution could be resolved by a chronogram, which allowed dates to be estimated for the various stages of molecular evolution. A genetic insertion event in an albumin precursor gene ∼45 Ma introduced two additional cleavage sites for protein maturation and conferred duality upon PawS1-Like genes such that they also encode a small buried macrocycle. Expansion of this region, including two Cys residues, enlarged the peptide ∼34 Ma and made the buried peptides bicyclic. Functional specialization into a protease inhibitor occurred ∼23 Ma. These findings document the evolution of a novel peptide inside a benign region of a pre-existing protein. We illustrate how a novel peptide can evolve without de novo gene evolution and, critically, without affecting the function of what becomes the protein host.

Macrocyclization by asparaginyl endopeptidases.

Contents Summary 923 I. Introduction 923 II. Plant AEPs with macrocyclizing ability 924 III. Mechanism of macrocyclization by AEPs 925 IV. Conclusions 927 Acknowledgements 927 References 927 SUMMARY: Plant asparaginyl endopeptidases (AEPs) are important for the post-translational processing of seed storage proteins via cleavage of precursor proteins. Some AEPs also function as peptide bond-makers during the biosynthesis of several unrelated classes of cyclic peptides, namely the kalata-type cyclic peptides, PawS-Derived Peptides and cyclic knottins. These three families of gene-encoded peptides have different evolutionary origins, but all have recruited AEPs for their maturation. In the last few years, the field has advanced rapidly, with the biochemical characterization of three plant AEPs capable of peptide macrocyclization, and insights have been gained from the first AEP crystal structures, albeit mammalian ones. Although the biochemical studies have improved our understanding of the mechanism of action, the focus now is to understand what changes in AEP sequence and structure enable some plant AEPs to perform macrocyclization reactions.

Orientation and Location of the Cyclotide Kalata B1 in Lipid Bilayers Revealed by Solid-State NMR.

Cyclotides are ultra-stable cyclic disulfide-rich peptides from plants. Their biophysical effects and medically interesting activities are related to their membrane-binding properties, with particularly high affinity for phosphatidylethanolamine lipids. In this study we were interested in understanding the molecular details of cyclotide-membrane interactions, specifically with regard to the spatial orientation of the cyclotide kalata B1 from Oldenlandia affinis when embedded in a lipid bilayer. Our experimental approach was based on the use of solid-state 19F-NMR of oriented bilayers in conjunction with the conformationally restricted amino acid L-3-(trifluoromethyl)bicyclopent-[1.1.1]-1-ylglycine as an orientation-sensitive 19F-NMR probe. Its rigid connection to the kalata B1 backbone scaffold, together with the well-defined structure of the cyclotide, allowed us to calculate the protein alignment in the membrane directly from the orientation-sensitive 19F-NMR signal. The hydrophobic and polar residues on the surface of kalata B1 form well-separated patches, endowing this cyclotide with a pronounced amphipathicity. The peptide orientation, as determined by NMR, showed that this amphipathic structure matches the polar/apolar interface of the lipid bilayer very well. A location in the amphiphilic headgroup region of the bilayer was supported by 15N-NMR of uniformly labeled protein, and confirmed using solid-state 31P- and 2H-NMR. 31P-NMR relaxation data indicated a change in lipid headgroup dynamics induced by kalata B1. Changes in the 2H-NMR order parameter profile of the acyl chains suggest membrane thinning, as typically observed for amphiphilic peptides embedded near the polar/apolar bilayer interface. Furthermore, from the 19F-NMR analysis two important charged residues, E7 and R28, were found to be positioned equatorially. The observed location thus would be favorable for the postulated binding of E7 to phosphatidylethanolamine lipid headgroups. Furthermore, it may be speculated that this pair of side chains could promote oligomerization of kalata B1 through electrostatic intermolecular contacts via their complementary charges.

DNA Gyrase Is the Target for the Quinolone Drug Ciprofloxacin in Arabidopsis thaliana.

The Arabidopsis thaliana genome contains four genes that were originally annotated as potentially encoding DNA gyrase: ATGYRA, ATGYRB1, ATGYRB2, and ATGYRB3. Although we subsequently showed that ATGYRB3 does not encode a gyrase subunit, the other three genes potentially encode subunits of a plant gyrase. We also showed evidence for the existence of supercoiling activity in A. thaliana and that the plant is sensitive to quinolone and aminocoumarin antibiotics, compounds that target DNA gyrase in bacteria. However, it was not possible at that time to show whether the A. thaliana genes encoded an active gyrase enzyme, nor whether that enzyme is indeed the target for the quinolone and aminocoumarin antibiotics. Here we show that an A. thaliana mutant resistant to the quinolone drug ciprofloxacin has a point mutation in ATGYRA. Moreover we show that, as in bacteria, the quinolone-sensitive (wild-type) allele is dominant to the resistant gene. Further we have heterologously expressed ATGYRA and ATGYRB2 in a baculovirus expression system and shown supercoiling activity of the partially purified enzyme. Expression/purification of the quinolone-resistant A. thaliana gyrase yields active enzyme that is resistant to ciprofloxacin. Taken together these experiments now show unequivocally that A. thaliana encodes an organelle-targeted DNA gyrase that is the target of the quinolone drug ciprofloxacin; this has important consequences for plant physiology and the development of herbicides.

Evolutionary origins of a bioactive peptide buried within Preproalbumin.

The de novo evolution of proteins is now considered a frequented route for biological innovation, but the genetic and biochemical processes that lead to each newly created protein are often poorly documented. The common sunflower (Helianthus annuus) contains the unusual gene PawS1 (Preproalbumin with SFTI-1) that encodes a precursor for seed storage albumin; however, in a region usually discarded during albumin maturation, its sequence is matured into SFTI-1, a protease-inhibiting cyclic peptide with a motif homologous to unrelated inhibitors from legumes, cereals, and frogs. To understand how PawS1 acquired this additional peptide with novel biochemical functionality, we cloned PawS1 genes and showed that this dual destiny is over 18 million years old. This new family of mostly backbone-cyclic peptides is structurally diverse, but the protease-inhibitory motif was restricted to peptides from sunflower and close relatives from its subtribe. We describe a widely distributed, potential evolutionary intermediate PawS-Like1 (PawL1), which is matured into storage albumin, but makes no stable peptide despite possessing residues essential for processing and cyclization from within PawS1. Using sequences we cloned, we retrodict the likely stepwise creation of PawS1's additional destiny within a simple albumin precursor. We propose that relaxed selection enabled SFTI-1 to evolve its inhibitor function by converging upon a successful sequence and structure.

Natural structural diversity within a conserved cyclic peptide scaffold.

We recently isolated and described the evolutionary origin of a diverse class of small single-disulfide bonded peptides derived from Preproalbumin with SFTI-1 (PawS1) proteins in the seeds of flowering plants (Asteraceae). The founding member of the PawS derived peptide (PDP) family is the potent trypsin inhibitor SFTI-1 (sunflower trypsin inhibitor-1) from Helianthus annuus, the common sunflower. Here we provide additional structures and describe the structural diversity of this new class of small peptides, derived from solution NMR studies, in detail. We show that although most have a similar backbone framework with a single disulfide bond and in many cases a head-to-tail cyclized backbone, they all have their own characteristics in terms of projections of side-chains, flexibility and physiochemical properties, attributed to the variety of their sequences. Small cyclic and constrained peptides are popular as drug scaffolds in the pharmaceutical industry and our data highlight how amino acid side-chains can fine-tune conformations in these promising peptides.