Structural homology guided alignment of cysteine rich proteins

Structural and functional basis of VLDLR usage by Eastern equine encephalitis virus

The very-low-density lipoprotein receptor (VLDLR) comprises eight LDLR type A (LA) domains and supports entry of distantly related alphaviruses, including Eastern equine encephalitis virus (EEEV) and Semliki Forest virus (SFV). Here, by resolving multiple cryo-electron microscopy structures of EEEV-VLDLR complexes and performing mutagenesis and functional studies, we show that EEEV uses multiple sites (E1/E2 cleft and E2 A domain) to engage more than one LA domain simultaneously. However, no single LA domain is necessary or sufficient to support efficient EEEV infection. Whereas all EEEV strains show conservation of two VLDLR-binding sites, the EEEV PE-6 strain and a few other EEE complex members feature a single amino acid substitution that enables binding of LA domains to an additional site on the E2 B domain. These structural and functional analyses informed the design of a minimal VLDLR decoy receptor that neutralizes EEEV infection and protects mice from lethal challenge.

Hyperpolarisation of Mitochondrial Membranes Is a Critical Component of the Antifungal Mechanism of the Plant Defensin, Ppdef1

Plant defensins are a large family of small cationic proteins with diverse functions and mechanisms of action, most of which assert antifungal activity against a broad spectrum of fungi. The partial mechanism of action has been resolved for a small number of members of plant defensins, and studies have revealed that many act by more than one mechanism. The plant defensin Ppdef1 has a unique sequence and long loop 5 with fungicidal activity against a range of human fungal pathogens, but little is known about its mechanism of action. We screened the S. cerevisiae non-essential gene deletion library and identified the involvement of the mitochondria in the mechanism of action of Ppdef1. Further analysis revealed that the hyperpolarisation of the mitochondrial membrane potential (MMP) activates ROS production, vacuolar fusion and cell death and is an important step in the mechanism of action of Ppdef1, and it is likely that a similar mechanism acts in Trichophyton rubrum.

Structural and functional basis of VLDLR receptor usage by Eastern equine encephalitis virus

The very low-density lipoprotein receptor (VLDLR) is comprised of eight LDLR type A (LA) domains and supports entry of distantly related Eastern equine encephalitis (EEEV) and Semliki Forest (SFV) alphaviruses. Here, by resolving multiple cryo-electron microscopy structures of EEEV-VLDLR complexes and performing mutagenesis and functional studies, we show that EEEV uses multiple sites (E1/E2 cleft and E2 A domain) to engage different LA domains simultaneously. However, no single LA domain is necessary or sufficient to support efficient EEEV infection, highlighting complexity in domain usage. Whereas all EEEV strains show conservation of two VLDLR binding sites, the EEEV PE-6 strain and other EEE complex members feature a single amino acid substitution that mediates binding of LA domains to an additional site on the E2 B domain. These structural and functional analyses informed the design of a minimal VLDLR decoy receptor that neutralizes EEEV infection and protects mice from lethal challenge.

Comprehensive computational analysis of the SRK-SP11 molecular interaction underlying self-incompatibility in Brassicaceae using improved structure prediction for cysteine-rich proteins

Plants employ self-incompatibility (SI) to promote cross-fertilization. In Brassicaceae, this process is regulated by the formation of a complex between the pistil determinant S receptor kinase (SRK) and the pollen determinant S-locus protein 11 (SP11, also known as S-locus cysteine-rich protein, SCR). In our previous study, we used the crystal structures of two eSRK-SP11 complexes in Brassica rapa S8 and S9 haplotypes and nine computationally predicted complex models to demonstrate that only the SRK ectodomain (eSRK) and SP11 pairs derived from the same S haplotype exhibit high binding free energy. However, predicting the eSRK-SP11 complex structures for the other 100 + S haplotypes and genera remains difficult because of SP11 polymorphism in sequence and structure. Although protein structure prediction using AlphaFold2 exhibits considerably high accuracy for most protein monomers and complexes, 46% of the predicted SP11 structures that we tested showed < 75 mean per-residue confidence score (pLDDT). Here, we demonstrate that the use of curated multiple sequence alignment (MSA) for cysteine-rich proteins significantly improved model accuracy for SP11 and eSRK-SP11 complexes. Additionally, we calculated the binding free energies of the predicted eSRK-SP11 complexes using molecular dynamics (MD) simulations and observed that some Arabidopsis haplotypes formed a binding mode that was critically different from that of B. rapa S8 and S9. Thus, our computational results provide insights into the haplotype-specific eSRK-SP11 binding modes in Brassicaceae at the residue level. The predicted models are freely available at Zenodo, https://doi.org/10.5281/zenodo.8047768.

Reduction in PLANT DEFENSIN 1 expression in Arabidopsis thaliana results in increased resistance to pathogens and zinc toxicity

Ectopic expression of defensins in plants correlates with their increased capacity to withstand abiotic and biotic stresses. This applies to Arabidopsis thaliana, where some of the seven members of the PLANT DEFENSIN 1 family (AtPDF1) are recognised to improve plant responses to necrotrophic pathogens and increase seedling tolerance to excess zinc (Zn). However, few studies have explored the effects of decreased endogenous defensin expression on these stress responses. Here, we carried out an extensive physiological and biochemical comparative characterization of (i) novel artificial microRNA (amiRNA) lines silenced for the five most similar AtPDF1s, and (ii) a double null mutant for the two most distant AtPDF1s. Silencing of five AtPDF1 genes was specifically associated with increased aboveground dry mass production in mature plants under excess Zn conditions, and with increased plant tolerance to different pathogens - a fungus, an oomycete and a bacterium, while the double mutant behaved similarly to the wild type. These unexpected results challenge the current paradigm describing the role of PDFs in plant stress responses. Additional roles of endogenous plant defensins are discussed, opening new perspectives for their functions.

Killer Knots: Molecular Evolution of Inhibitor Cystine Knot Toxins in Wandering Spiders (Araneae: Ctenidae)

Venom expressed by the nearly 50,000 species of spiders on Earth largely remains an untapped reservoir of a diverse array of biomolecules with potential for pharmacological and agricultural applications. A large fraction of the noxious components of spider venoms are a functionally diverse family of structurally related polypeptides with an inhibitor cystine knot (ICK) motif. The cysteine-rich nature of these toxins makes structural elucidation difficult, and most studies have focused on venom components from the small handful of medically relevant spider species such as the highly aggressive Brazilian wandering spider Phoneutria nigriventer. To alleviate difficulties associated with the study of ICK toxins in spiders, we devised a comprehensive approach to explore the evolutionary patterns that have shaped ICK functional diversification using venom gland transcriptomes and proteomes from phylogenetically distinct lineages of wandering spiders and their close relatives. We identified 626 unique ICK toxins belonging to seven topological elaborations. Phylogenetic tests of episodic diversification revealed distinct regions between cysteine residues that demonstrated differential evidence of positive or negative selection, which may have structural implications towards the specificity and efficacy of these toxins. Increased taxon sampling and whole genome sequencing will provide invaluable insights to further understand the evolutionary processes that have given rise to this diverse class of toxins.

Orientational Preferences of GPI-Anchored Ly6/uPAR Proteins

Ly6/uPAR proteins regulate many essential functions in the nervous and immune systems and epithelium. Most of these proteins contain single β-structural LU domains with three protruding loops and are glycosylphosphatidylinositol (GPI)-anchored to a membrane. The GPI-anchor role is currently poorly studied. Here, we investigated the positional and orientational preferences of six GPI-anchored proteins in the receptor-unbound state by molecular dynamics simulations. Regardless of the linker length between the LU domain and GPI-anchor, the proteins interacted with the membrane by polypeptide parts and N-/O-glycans. Lynx1, Lynx2, Lypd6B, and Ly6H contacted the membrane by the loop regions responsible for interactions with nicotinic acetylcholine receptors, while Lypd6 and CD59 demonstrated unique orientations with accessible receptor-binding sites. Thus, GPI-anchoring does not guarantee an optimal 'pre-orientation' of the LU domain for the receptor interaction.

Histidine-Rich Defensins from the Solanaceae and Brasicaceae Are Antifungal and Metal Binding Proteins

Plant defensins are best known for their antifungal activity and contribution to the plant immune system. The defining feature of plant defensins is their three-dimensional structure known as the cysteine stabilized alpha-beta motif. This protein fold is remarkably tolerant to sequence variation with only the eight cysteines that contribute to the stabilizing disulfide bonds absolutely conserved across the family. Mature defensins are typically 46–50 amino acids in length and are enriched in lysine and/or arginine residues. Examination of a database of approximately 1200 defensin sequences revealed a subset of defensin sequences that were extended in length and were enriched in histidine residues leading to their classification as histidine-rich defensins (HRDs). Using these initial HRD sequences as a query, a search of the available sequence databases identified over 750 HRDs in solanaceous plants and 20 in brassicas. Histidine residues are known to contribute to metal binding functions in proteins leading to the hypothesis that HRDs would have metal binding properties. A selection of the HRD sequences were recombinantly expressed and purified and their antifungal and metal binding activity was characterized. Of the four HRDs that were successfully expressed all displayed some level of metal binding and two of four had antifungal activity. Structural characterization of the other HRDs identified a novel pattern of disulfide linkages in one of the HRDs that is predicted to also occur in HRDs with similar cysteine spacing. Metal binding by HRDs represents a specialization of the plant defensin fold outside of antifungal activity.

Orientation-dependent toxic effect of human papillomavirus type 33 long control region DNA in Escherichia coli cells

The functional analysis of human papillomavirus (HPV) sequence variation requires the molecular cloning of different genomic regions of virus variants. In this study, we report an unexpected difficulty experienced when trying to clone HPV33 long control region (LCR) variants in Escherichia coli. Standard cloning strategies proved to be inappropriate to clone HPV33 LCR variants in the forward orientation into a eukaryotic reporter vector (pGL2-Basic). However, by slight modification of culture conditions (incubation at 25 °C instead of 37 °C), constructs containing the HPV33 LCR variants in the forward orientation were obtained. Transformation experiments performed with different HPV33 LCR constructs indicated that there is a sequence element in the 5′ LCR of HPV33 causing temperature-dependent toxic effect in E. coli. Sequence analysis revealed the presence of an open reading frame (ORF) in the 5′ part of HPV33 LCR potentially encoding a 116-amino acid polypeptide. Protein structure prediction suggested that this putative protein might have a structural similarity to transmembrane proteins. Even a low-level expression of this protein may cause significant toxicity in the host bacteria. In silico analysis of the LCR of HPV33 and some other HPV types belonging to the species Alphapapillomavirus 9 (HPV31, 35 and 58) seemed to support the assumption that the ORFs found in the 5′ LCR of these HPVs are protein-coding sequences. Further studies should be performed to prove that these putative proteins are really expressed in the infected host cells and to identify their function.

Mapping the chemical and sequence space of the ShKT superfamily

The ShKT superfamily is widely distributed throughout nature and encompasses a wide range of documented functions and processes, from modulation of potassium channels to involvement in morphogenesis pathways. Cysteine-rich secretory proteins (CRISPs) contain a cysteine-rich domain (CRD) at the C-terminus that is similar in structure to the ShK fold. Despite the structural similarity of the CRD and ShK-like domains, we know little of the sequence-function relationships in these families. Here, for the first time, we examine the evolution of the biophysical properties of sequences within the ShKT superfamily in relation to function, with a focus on the ShK-like superfamily. ShKT data were sourced from published sequences in the protein family database, in addition to new ShK-like sequences from the Australian speckled anemone (Oulactis sp.). Our analysis clearly delineates the ShK-like family from the CRDs of CRISP proteins. The four CRISP subclusters separate out into the main phyla of Mammalia, Insecta and Reptilia. The ShK-like family is in turn composed of seven subclusters, the largest of which contains members from across the eukaryotes, with a continuum of intermediate properties. Smaller sub-clusters contain specialised members such as nematode ShK-like sequences. Several of these ShKT sub-clusters contain no functionally characterised sequences. This chemical space analysis should be useful as a guide to select sequences for functional studies and to gain insight into the evolution of these highly divergent sequences with an ancient conserved fold.

The evolution, function and mechanisms of action for plant defensins

Plant defensins are an extensive family of small cysteine rich proteins characterised by a conserved cysteine stabilised alpha beta protein fold which resembles the structure of insect and vertebrate defensins. However, secondary structure and disulphide topology indicates two independent superfamilies of defensins with similar structures that have arisen via an extreme case of convergent evolution. Defensins from plants and insects belong to the cis-defensin superfamily whereas mammalian defensins belong to the trans-defensin superfamily. Plant defensins are produced by all species of plants and although the structure is highly conserved, the amino acid sequences are highly variable with the exception of the cysteine residues that form the stabilising disulphide bonds and a few other conserved residues. The majority of plant defensins are components of the plant innate immune system but others have evolved additional functions ranging from roles in sexual reproduction and development to metal tolerance. This review focuses on the antifungal mechanisms of plant defensins. The activity of plant defensins is not limited to plant pathogens and many of the described mechanisms have been elucidated using yeast models. These mechanisms are more complex than simple membrane permeabilisation induced by many small antimicrobial peptides. Common themes that run through the characterised mechanisms are interactions with specific lipids, production of reactive oxygen species and induction of cell wall stress. Links between sequence motifs and functions are highlighted where appropriate. The complexity of the interactions between plant defensins and fungi helps explain why this protein superfamily is ubiquitous in plant innate immunity.

Convergent evolution of defensin sequence, structure and function

Defensins are a well-characterised group of small, disulphide-rich, cationic peptides that are produced by essentially all eukaryotes and are highly diverse in their sequences and structures. Most display broad range antimicrobial activity at low micromolar concentrations, whereas others have other diverse roles, including cell signalling (e.g. immune cell recruitment, self/non-self-recognition), ion channel perturbation, toxic functions, and enzyme inhibition. The defensins consist of two superfamilies, each derived from an independent evolutionary origin, which have subsequently undergone extensive divergent evolution in their sequence, structure and function. Referred to as the cis- and trans-defensin superfamilies, they are classified based on their secondary structure orientation, cysteine motifs and disulphide bond connectivities, tertiary structure similarities and precursor gene sequence. The utility of displaying loops on a stable, compact, disulphide-rich core has been exploited by evolution on multiple occasions. The defensin superfamilies represent a case where the ensuing convergent evolution of sequence, structure and function has been particularly extreme. Here, we discuss the extent, causes and significance of these convergent features, drawing examples from across the eukaryotes.

AlignStat: a web-tool and R package for statistical comparison of alternative multiple sequence alignments

BACKGROUND

Alternative sequence alignment algorithms yield different results. It is therefore useful to quantify the similarities and differences between alternative alignments of the same sequences. These measurements can identify regions of consensus that are likely to be most informative in downstream analysis. They can also highlight systematic differences between alignments that relate to differences in the alignment algorithms themselves.

RESULTS

Here we present a simple method for aligning two alternative multiple sequence alignments to one another and assessing their similarity. Differences are categorised into merges, splits or shifts in one alignment relative to the other. A set of graphical visualisations allow for intuitive interpretation of the data.

CONCLUSIONS

AlignStat enables the easy one-off online use of MSA similarity comparisons or into R pipelines. The web-tool is available at AlignStat.Science.LaTrobe.edu.au. The R package, readme and example data are available on CRAN and GitHub.com/TS404/AlignStat.

The plant defensin NaD1 introduces membrane disorder through a specific interaction with the lipid, phosphatidylinositol 4,5 bisphosphate

Plant defensins interact with phospholipids in bilayers as part of their cytotoxic activity. Solanaceous class II defensins with the loop 5 sequence pattern "S-[KR]-[ILVQ]-[ILVQ]-[KR]-[KR]" interact with PI(4,5)P2. Here, the prototypical defensin of this class, NaD1, is used to characterise the biophysical interactions between these defensins and phospholipid bilayers. Binding of NaD1 to bilayers containing PI(4,5)P2 occurs rapidly and the interaction is very strong. Dual polarisation interferometry revealed that NaD1 does not dissociate from bilayers containing PI(4,5)P2. Binding of NaD1 to bilayers with or without PI(4,5)P2 induced disorder in the bilayer. However, permeabilisation assays revealed that NaD1 only permeabilised liposomes with PI(4,5)P2 in the bilayer, suggesting a role for this protein-lipid interaction in the plasma membrane permeabilising activity of this defensin. No defensins in the available databases have the PI(4,5)P2 binding sequence outside the solanaceous class II defensins, leading to the hypothesis that PI(4,5)P2 binding co-evolved with the C-terminal propeptide to protect the host cell against the effects of the tight binding of these defensins to their cognate lipid as they travel along the secretory pathway. This data has allowed us to develop a new model to explain how this class of defensins permeabilises plasma membranes to kill target cells.