Structure and multiple conformations of the kunitz-type domain from human type VI collagen alpha3(VI) chain in solution

A novel family of defensin-like peptides from Hermetia illucens with antibacterial properties

BACKGROUND

The world faces a major infectious disease challenge. Interest in the discovery, design, or development of antimicrobial peptides (AMPs) as an alternative approach for the treatment of bacterial infections has increased. Insects are a good source of AMPs which are the main effector molecules of their innate immune system. Black Soldier Fly Larvae (BSFL) are being developed for large-scale rearing for food sustainability, waste reduction and as sustainable animal and fish feed. Bioinformatic studies have suggested that BSFL have the largest number of AMPs identified in insects. However, most AMPs identified in BSF have not yet undergone antimicrobial evaluation but are promising leads to treat critical infections.

RESULTS

Jg7197.t1, Jg7902.t1 and Jg7904.t1 were expressed into the haemolymph of larvae following infection with Salmonella enterica serovar Typhimurium and were predicted to be AMPs using the computational tool ampir. The genes encoding these proteins were within 2 distinct clusters in chromosome 1 of the BSF genome. Following removal of signal peptides, predicted structures of the mature proteins were superimposed, highlighting a high degree of structural conservation. The 3 AMPs share primary sequences with proteins that contain a Kunitz-binding domain; characterised for inhibitory action against proteases, and antimicrobial activities. An in vitro antimicrobial screen indicated that heterologously expressed SUMO-Jg7197.t1 and SUMO-Jg7902.t1 did not show activity against 12 bacterial strains. While recombinant SUMO-Jg7904.t1 had antimicrobial activity against a range of Gram-negative and Gram-positive bacteria, including the serious pathogen Pseudomonas aeruginosa.

CONCLUSIONS

We have cloned and purified putative AMPs from BSFL and performed initial in vitro experiments to evaluate their antimicrobial activity. In doing so, we have identified a putative novel defensin-like AMP, Jg7904.t1, encoded in a paralogous gene cluster, with antimicrobial activity against P. aeruginosa.

Lack of evidence for a role of anthrax toxin receptors as surface receptors for collagen VI and for its cleaved-off C5 domain/endotrophin

The microfibril-forming collagen VI is proteolytically cleaved and it was proposed that the released C-terminal Kunitz domain (C5) of the α3 chain is an adipokine important for tumor progression and fibrosis. Designated “endotrophin,” C5 is a potent biomarker for fibroinflammatory diseases. However, the biochemical mechanisms behind endotrophin activity were not investigated. Earlier, anthrax toxin receptor 1 was found to bind C5, but this potential interaction was not further studied. Given the proposed physiological role of endotrophin, we aimed to determine how the signal is transmitted. Surprisingly, we could not detect any interaction between endotrophin and anthrax toxin receptor 1 or its close relative, anthrax toxin receptor 2. Moreover, we detect no binding of fully assembled collagen VI to either receptor. We also studied the collagen VI receptor NG2 (CSPG4) and confirmed that NG2 binds assembled collagen VI, but not cleaved C5/endotrophin. A cellular receptor for C5/endotrophin, therefore, still remains elusive.

•

ANTXR1 does not support collagen VI or C5/endotrophin binding to the cell surface

•

ANTXR2 does not support collagen VI or C5/endotrophin binding to the cell surface

•

NG2/CSPG4 supports collagen VI, but not C5/endotrophin binding to the cell surface

Analysis of High Molecular Mass Compounds from the Spider Pamphobeteus verdolaga Venom Gland. A Transcriptomic and MS ID Approach

Nowadays, spider venom research focuses on the neurotoxic activity of small peptides. In this study, we investigated high-molecular-mass compounds that have either enzymatic activity or housekeeping functions present in either the venom gland or venom of Pamphobeteus verdolaga. We used proteomic and transcriptomic-assisted approaches to recognize the proteins sequences related to high-molecular-mass compounds present in either venom gland or venom. We report the amino acid sequences (partial or complete) of 45 high-molecular-mass compounds detected by transcriptomics showing similarity to other proteins with either enzymatic activity (i.e., phospholipases A₂, kunitz-type, hyaluronidases, and sphingomyelinase D) or housekeeping functions involved in the signaling process, glucanotransferase function, and beta-N-acetylglucosaminidase activity. MS/MS analysis showed fragments exhibiting a resemblance similarity with different sequences detected by transcriptomics corresponding to sphingomyelinase D, hyaluronidase, lycotoxins, cysteine-rich secretory proteins, and kunitz-type serine protease inhibitors, among others. Additionally, we report a probably new protein sequence corresponding to the lycotoxin family detected by transcriptomics. The phylogeny analysis suggested that P. verdolaga includes a basal protein that underwent a duplication event that gave origin to the lycotoxin proteins reported for Lycosa sp. This approach allows proposing an evolutionary relationship of high-molecular-mass proteins among P. verdolaga and other spider species.

Neuropeptide signalling systems - An underexplored target for venom drug discovery

Neuropeptides are signalling molecules mainly secreted from neurons that act as neurotransmitters or peptide hormones to affect physiological processes and modulate behaviours. In humans, neuropeptides are implicated in numerous diseases and understanding their role in physiological processes and pathologies is important for therapeutic development. Teasing apart the (patho)physiology of neuropeptides remains difficult due to ligand and receptor promiscuity and the complexity of the signalling pathways. The current approach relies on a pharmacological toolbox of agonists and antagonists displaying high selectivity for independent receptor subtypes, with the caveat that only few selective ligands have been discovered or developed. Animal venoms represent an underexplored source for novel receptor subtype-selective ligands that could aid in dissecting human neuropeptide signalling systems. Multiple endogenous-like neuropeptides as well as peptides acting on neuropeptide receptors are present in venoms. In this review, we summarise current knowledge on neuropeptides and discuss venoms as a source for ligands targeting neuropeptide signalling systems.

Evolutionary Aspects of the Structural Convergence and Functional Diversification of Kunitz-Domain Inhibitors

Kunitz-type domains are ubiquitously found in natural systems as serine protease inhibitors or animal toxins in venomous animals. Kunitz motif is a cysteine-rich peptide chain of ~ 60 amino acid residues with alpha and beta fold, stabilized by three conserved disulfide bridges. An extensive dataset of amino acid variations is found on sequence analysis of various Kunitz peptides. Kunitz peptides show diverse biological activities like inhibition of proteases of other classes and/or adopting a new function of blocking or modulating the ion channels. Based on the amino acid residues at the functional site of various Kunitz-type inhibitors, it is inferred that this 'flexibility within the structural rigidity' is responsible for multiple biological activities. Accelerated evolution of functional sites in response to the co-evolving molecular targets of the hosts of venomous animals or parasites, gene sharing, and gene duplication have been discussed as the most likely mechanisms responsible for the functional heterogeneity of Kunitz-domain inhibitors.

Enhanced dynamics of conformationally heterogeneous T7 bacteriophage lysozyme native state attenuates its stability and activity

Proteins are dynamic in nature and exist in a set of equilibrium conformations on various timescale motions. The flexibility of proteins governs various biological functions, and therefore elucidation of such functional dynamics is essential. In this context, we have studied the structure–dynamics–stability–activity relationship of bacteriophage T7 lysozyme/endolysin (T7L) native-state ensemble in the pH range of 6–8. Our studies established that T7L native state is conformationally heterogeneous, as several residues of its C-terminal half are present in two conformations (major and minor) in the slow exchange time scale of nuclear magnetic resonance (NMR). Structural and dynamic studies suggested that the residues belonging to minor conformations do exhibit native-like structural and dynamic features. Furthermore, the NMR relaxation experiments unraveled that the native state is highly dynamic and the dynamic behavior is regulated by the pH, as the pH 6 conformation exhibited enhanced dynamics compared with pH 7 and 8. The stability measurements and cell-based activity studies on T7L indicated that the native protein at pH 6 is ∼2 kcal less stable and is ∼50% less active than those of pH 7 and 8. A comprehensive analysis of the T7L active site, unfolding initiation sites and the residues with altered dynamics outlined that the attenuation of stability and activity is a resultant of its enhanced dynamic properties, which, in turn, can be attributed to the protonation/deprotonation of its partially buried His residues. Our study on T7L structure–dynamics–activity paradigm could assist in engineering novel amidase-based endolysins with enhanced activity and stability over a broad pH range.

The Anemonia viridis Venom: Coupling Biochemical Purification and RNA-Seq for Translational Research

Blue biotechnologies implement marine bio-resources for addressing practical concerns. The isolation of biologically active molecules from marine animals is one of the main ways this field develops. Strikingly, cnidaria are considered as sustainable resources for this purpose, as they possess unique cells for attack and protection, producing an articulated cocktail of bioactive substances. The Mediterranean sea anemone Anemonia viridis has been studied extensively for years. In this short review, we summarize advances in bioprospecting of the A. viridis toxin arsenal. A. viridis RNA datasets and toxin data mining approaches are briefly described. Analysis reveals the major pool of neurotoxins of A. viridis, which are particularly active on sodium and potassium channels. This review therefore integrates progress in both RNA-Seq based and biochemical-based bioprospecting of A. viridis toxins for biotechnological exploitation.

Defining the hierarchical organisation of collagen VI microfibrils at nanometre to micrometre length scales

Extracellular matrix microfibrils are critical components of connective tissues with a wide range of mechanical and cellular signalling functions. Collagen VI is a heteromeric network-forming collagen which is expressed in tissues such as skin, lung, blood vessels and articular cartilage where it anchors cells into the matrix allowing for transduction of biochemical and mechanical signals. It is not understood how collagen VI is arranged into microfibrils or how these microfibrils are arranged into tissues. Therefore we have characterised the hierarchical organisation of collagen VI across multiple length scales. The frozen hydrated nanostructure of purified collagen VI microfibrils was reconstructed using cryo-TEM. The bead region has a compact hollow head and flexible tail regions linked by the collagenous interbead region. Serial block face SEM imaging coupled with electron tomography of the pericellular matrix (PCM) of murine articular cartilage revealed that the PCM has a meshwork-like organisation formed from globular densities ∼30nm in diameter. These approaches can characterise structures spanning nanometer to millimeter length scales to define the nanostructure of individual collagen VI microfibrils and the micro-structural organisation of these fibrils within tissues to help in the future design of better mimetics for tissue engineering.

STATEMENT OF SIGNIFICANCE

Cartilage is a connective tissue rich in extracellular matrix molecules and is tough and compressive to cushion the bones of joints. However, in adults cartilage is poorly repaired after injury and so this is an important target for tissue engineering. Many connective tissues contain collagen VI, which forms microfibrils and networks but we understand very little about these assemblies or the tissue structures they form. Therefore, we have use complementary imaging techniques to image collagen VI microfibrils from the nano-scale to the micro-scale in order to understand the structure and the assemblies it forms. These findings will help to inform the future design of scaffolds to mimic connective tissues in regenerative medicine applications.

Purification, characterization and molecular cloning of chymotrypsin inhibitor peptides from the venom of Burmese Daboia russelii siamensis

One novel Kunitz BPTI-like peptide designated as BBPTI-1, with chymotrypsin inhibitory activity was identified from the venom of Burmese Daboia russelii siamensis. It was purified by three steps of chromatography including gel filtration, cation exchange and reversed phase. A partial N-terminal sequence of BBPTI-1, HDRPKFCYLPADPGECLAHMRSF was obtained by automated Edman degradation and a Ki value of 4.77nM determined. Cloning of BBPTI-1 including the open reading frame and 3' untranslated region was achieved from cDNA libraries derived from lyophilized venom using a 3' RACE strategy. In addition a cDNA sequence, designated as BBPTI-5, was also obtained. Alignment of cDNA sequences showed that BBPTI-5 exhibited an identical sequence to BBPTI-1 cDNA except for an eight nucleotide deletion in the open reading frame. Gene variations that represented deletions in the BBPTI-5 cDNA resulted in a novel protease inhibitor analog. Amino acid sequence alignment revealed that deduced peptides derived from cloning of their respective precursor cDNAs from libraries showed high similarity and homology with other Kunitz BPTI proteinase inhibitors. BBPTI-1 and BBPTI-5 consist of 60 and 66 amino acid residues respectively, including six conserved cysteine residues. As these peptides have been reported to have influence on the processes of coagulation, fibrinolysis and inflammation, their potential application in biomedical contexts warrants further investigation.

Papilin, a novel component of basement membranes, in relation to ADAMTS metalloproteases and ECM development

Papilins are homologous, secreted extracellular matrix proteins which share a common order of protein domains. They occur widely, from nematodes to man, and can differ in the number of repeats of a given type of domain. Within one species the number of repeats can vary by differential RNA splicing. A distinctly conserved cassette of domains at the amino-end of papilins is homologous with a cassette of protein domains at the carboxyl-end of the ADAMTS subgroup of secreted, matrix-associated metalloproteases. Papilins primarily occur in basement membranes. Papilins interact with several extracellular matrix components and ADAMTS enzymes. Papilins are essential for embryonic development of Drosophila melanogaster and Caenorhabditis elegans.

Solution structure of a Kunitz-type chymotrypsin inhibitor isolated from the elapid snake Bungarus fasciatus

Bungarus fasciatus fraction IX (BF9), a chymotrypsin inhibitor, consists of 65 amino acid residues with three disulfide bridges. It was isolated from the snake venom of B. fasciatus by ion-exchange chromatography and belongs to the bovine pancreatic trypsin inhibitor (BPTI)-like superfamily. It showed a dissociation constant of 5.8 x 10(-8) m with alpha-chymotrypsin as measured by a BIAcore binding assay system. The isothermal titration calorimetry revealed a 1:1 binding stoichiometry between this inhibitor and chymotrypsin and apparently no binding with trypsin. We further used CD and NMR to determine the solution structure of this venom-derived chymotrypsin inhibitor. The three-dimensional NMR solution structures of BF9 were determined on the basis of 582 restraints by simulated annealing and energy minimization calculations. The final set of 10 NMR structures was well defined, with average root mean square deviations of 0.47 A for the backbone atoms in the secondary structure regions and 0.86 A for residues The side chains of Phe(23), Tyr(24), Tyr(25), Phe(35), and Phe(47) exhibited many long-range nuclear Overhauser effects and were the principal components of the hydrophobic core in BF9. To gain insight into the structure-function relationships among proteins in the BPTI-like superfamily, we compared the three-dimensional structure of BF9 with three BPTI-like proteins that possess distinct biological functions. These proteins possessed similar secondary structure elements, but the loop regions and beta-turn were different from one another. Based on residues at the functional site of each protein, we suggest that the flexibility, rigidity, and variations of the amino acid residues in both the loop and beta-turn regions are related to their biological functions.

High-resolution structure of bovine pancreatic trypsin inhibitor with altered binding loop sequence

A mutant of bovine pancreatic trypsin inhibitor (BPTI) has been constructed and expressed in Escherichia coli in order to probe the kinetic and structural consequences of truncating the binding loop residues to alanine. In addition to two such mutations (Thr11Ala and Pro13Ala), it has a conservative Lys15Arg substitution at position P(1) and an unrelated Met52Leu change. In spite of the binding loop modification, the affinity for trypsin is only 30 times lower than that of the wild-type protein. At pH 7.5 the protein can be crystallized on the time-scale of hours, yielding very stable crystals of a new (tetragonal) form of BPTI. Conventional source X-ray data collected to 1.4 A at room temperature allowed anisotropic structure refinement characterized by R=0.1048. The structure reveals all 58 residues, including the complete C terminus, which is in a salt-bridge contact with the N terminus. The Cys14-Cys38 disulfide bridge is observed in two distinct chiralities. This bridge, together with an internal water molecule, contributes to the stabilization of the binding loop. The Ala mutations have only an insignificant and localized effect on the binding loop, which retains its wild-type conformation (maximum deviation of loop C(alpha) atoms of 0.7 A at Ala13). Four (instead of the typical three) additional water molecules are buried in an internal cleft and connected to the surface via a sulfate anion. Three more SO(4)(2-) anions are seen in the electron density, one of them located on a 2-fold axis. It participates in the formation of a dimeric structure between symmetry-related BPTI molecules, in which electrostatic and hydrogen bonding interactions resulting from the mutated Lys15Arg substitution are of central importance. This dimeric interaction involves direct recognition loop-recognition loop contacts, part of which are hydrophobic interactions of the patches created by the alanine mutations. Another 2-fold symmetric interaction between the BPTI molecules involves the formation of an antiparallel intermolecular beta-sheet that, together with the adjacent intramolecular beta-hairpin loops, creates a four-stranded structure.

Interscaffolding additivity: binding of P1 variants of bovine pancreatic trypsin inhibitor to four serine proteases

Different families of protein inhibitors of serine proteases share similar conformation of the enzyme-binding loop, while their scaffoldings are completely different. In the enzyme-inhibitor complex, the P1position of the loop makes numerous contacts within the S1pocket and significantly influences the energy of the interaction. Here, we determine the association energies (DeltaGavalues) for the interaction of coded P1variants of bovine pancreatic trypsin inhibitor (BPTI) with bovine beta-trypsin (BT), anionic salmon trypsin (AST), bovine alpha-chymotrypsin (BCHYM), and human neutrophil elastase (HNE). The respective DeltaGaranges are 15, 13, 9, and 8 kcal mol-1(1 cal=4.18 J). Next, through interscaffolding additivity cycles, we compare our set of DeltaGavalues determined for BCHYM and HNE with similar data sets available in the literature for three other inhibitor families. The analysis of the cycles shows that 27 to 83 % of cycles fulfil the criteria of additvity. In one particular case (comparison of associations of P1variants of BPTI and OMTKY3 with BCHYM) there is a structural basis for strongly non-additive behaviour. We argue that the interscaffolding additvity depends on sequential and conformational similarities of sites where the mutation(s) are introduced and on the particular substitution. In the second interscaffolding analysis, we compare binding of the same P1mutants to BT and AST. The high correlation coefficient shows that both trypsins recognize with comparable strength the non-cognate side-chains. However, the cognate Arg and Lys side-chains are recognized significantly more strongly by AST.

Assignment of side-chain conformation using adiabatic energy mapping, free energy perturbation, and molecular dynamic simulations

NMR spectroscopic analysis of the C-terminal Kunitz domain fragment (alpha3(VI)) from the human alpha3-chain of type VI collagen has revealed that the side chain of Trp21 exists in two unequally populated conformations. The major conformation (M) is identical to the conformation observed in the X-ray crystallographic structure, while the minor conformation (m) cannot structurally be resolved in detail by NMR due to insufficient NOE data. In the present study, we have applied: (1) rigid and adiabatic mapping, (2) free energy simulations, and (3) molecular dynamic simulations to elucidate the structure of the m conformer and to provide a possible pathway of the Trp21 side chain between the two conformers. Adiabatic energy mapping of conformations of the Trp21 side chain obtained by energy minimization identified two energy minima: One corresponding to the conformation of Trp21 observed in the X-ray crystallographic structure and solution structure of alpha3(VI) (the M conformation) and the second corresponding to the m conformation predicted by NMR spectroscopy. A transition pathway between the M and m conformation is suggested. The free-energy difference between the two conformers obtained by the thermodynamic integration method is calculated to 1.77+/-0.7 kcal/mol in favor of the M form, which is in good agreement with NMR results. Structural and dynamic properties of the major and minor conformers of the alpha3(VI) molecule were investigated by molecular dynamic. Essential dynamics analysis of the two resulting 800 ps trajectories reveals that when going from the M to the m conformation only small, localized changes in the protein structure are induced. However, notable differences are observed in the mobility of the binding loop (residues Thr13-Ile18), which is more flexible in the m conformation than in the M conformation. This suggests that the reorientation of Trp2 might influence the inhibitory activity against trypsin, despite the relative large distance between the binding loop and Trp21.

Structure of single-disulfide variants of bovine pancreatic trypsin inhibitor (BPTI) as probed by their binding to bovine beta-trypsin

Native bovine pancreatic trypsin inhibitor (BPTI) contains three disulfide bonds: Cys5-Cys55, Cys14-Cys38 and Cys30-Cys51. Correct cysteine pairing, native structure, and full anti-proteinase activity can be restored in the process of oxidative refolding of reduced BPTI. Oxidative refolding starts with the formation of single disulfide intermediates. All 15 single-disulfide variants of BPTI (three native and 12 non-native combinations) have been expressed in Escherichia coli. In each variant the remaining four cysteine residues were replaced by alanine. Four of these variants are shown here to inhibit bovine beta-trypsin: three of them contain native and one non-native (Cys5-Cys51) disulfide. All but one (Cys5-Cys55) variant are slowly digested by the enzyme, therefore measurements were performed at pH 4.0, at which trypsin activity is low. Binding constants of these four single disulfide variants were at least two orders of magnitude lower than for native BPTI. Remarkably, in some of the variants the binding constants were found to be higher for the reduced rather than for the oxidized form of the variant. Also for the fully reduced native BPTI, determined here, the binding constant is of considerable value. Two sets of control experiments demonstrated that the binding of reduced native BPTI to trypsin is specific. In the first, mutation of Lys15 (P1 position) in the binding loop abolished binding of the reduced forms to trypsin. In the second, the binding of reduced native BPTI to anhydrotrypsin yielded the expected UV difference spectra. In general, the results obtained indicate that the inhibitor activity can be induced even in the reduced protein. This activity is not a local effect, such as the nature of residues surrounding the binding loop, but rather is induced by residual structure in the unfolded protein. This structure has been shown to consist of a set of hydrophobic residues and the data presented here indicate that reduced cysteine residues provide further stabilization of such a hydrophobic cluster. On the other hand, improper pairing of the cysteine residues in non-native single disulfide variants destabilizes the enzyme-inhibitor complex by inducing deformations of the binding loop region.

Elucidation of the origin of multiple conformations of the human alpha 3-chain type VI collagen C-terminal Kunitz domain: the reorientation of the Trp21 ring

The human alpha 3-chain type VI collagen C-terminal Kunitz domain fragment (alpha 3(VI)) has been studied by two dimensional 1H-1H and 1H-13C NMR spectroscopy at 303 K. It is shown that the secondary structure of the protein is strikingly similar to that of BPTI, and a number of unusual H alpha chemical shifts, which are highly conserved in Kunitz-domain proteins, are also observed for a alpha 3(VI). Furthermore, a series of exchange cross peaks observed in 1H-1H spectra shows that a large number of protons in the central beta-sheet exist in two different chemical environments, corresponding to two unequally populated conformations that are slowly exchanging on the NMR time scale. Several protons, including Ser47(53) H alpha, Arg32(28) H(gamma 1) and H(gamma 2), and GLN48(54) H(beta 2), all located in the vicinity of the Trp21(27) ring in the crystal structure of alpha 3(VI) [Arnoux, B. et al. (1995) J. Mol. Biol., 246, 609-617], have very different chemical shifts in the two conformations, the most affected being Gln48(54) H(beta 2) (delta sigma = 3 ppm), which is placed directly above the Trp21(27) ring in the crystal structure of alpha 3(VI). It should be concluded that the origin of the multiple conformations of the central beta-sheet is a reorientation of the Trp21(27) ring. From the intensities of corresponding signals in the two conformations, the populations, the population of the minor conformation was found to be 6.4 +/- 0.2% of that of the major conformation, while a rate constant kM = 1.01 +/- 0.05 s-1 for the major to minor interconversion was obtained from a series of NOESY spectra with different mixing times. In addition, it is shown that Cys14(20)-Cys38(44) disulfide bond isomerization, previously observed in BPTI [Otting, G. et al. (1993) Biochemistry, 32, 3571-3582], is also likely to occur in alpha 3(VI).

Conversion of the Kunitz-type module of collagen VI into a highly active trypsin inhibitor by site-directed mutagenesis

The recombinant Kunitz protease inhibitor module (domain C5) of human collagen alpha 3(VI) chain was previously shown to lack inhibitory activity for proteases with trypsin-like specificity and some other proteases. We have now prepared mutants in the binding loop region including the P1' site (D2889-->A), the P2' site (F2890-->R) and the P3 site (T2886-->P) and in a more remote region (W2907-->V) either as individual substitutions or combinations of them. These mutants were analyzed for their kinetics of binding to trypsin by surface plasmon resonance and for their capacity to inhibit various proteases. Single substitutions (D-->A, T-->P, W-->V) showed an effect only for D->A which bound to trypsin with Kd = 0.25 microM. A 25-100-fold increase in affinity was observed for the double mutants T-->P/D-->A and F-->R/D-->A and approached the affinity of aprotinin (Kd approximately 0.01 nM) in two different triple mutants. These affinities correlated well with the inhibitory capacities of the mutants for trypsin in the cleavage of a large protein and a small peptide substrate. A similar but not completely identical improvement in inhibitory capacity was also observed for leucocyte elastase but not for thrombin. These data could be interpreted in terms of steric interferences or lack of hydrogen bonding of a few critical residues based on three-dimensional structures available for the C5 domain.