Expression and Characterization of Manduca sexta Stress Responsive Peptide-1; An Inducer of Antimicrobial Peptide Synthesis

In response to stress conditions such as wounding or infections in insects, several short peptides are processed to act as cytokines that induce AMP gene expression. To study their structure-activity relationship, immune inducibility, tissue specificity, stress responsiveness, and development relatedness, we chemically synthesized Manduca sexta stress response peptide-1, a 25-residue peptide with one disulfide bond (SRP1: FGVRVGTCPSGYVRRGTFCFPDDDY). Upon injection of the SRP1 into naïve larvae, several antimicrobial peptide genes were expressed at higher levels. The mRNA levels of SRP1 increased significantly in hemocytes and fat body after larvae were challenged with a mixture of bacteria and β-1,3-glucan. The expression patterns of SRP1 and its target genes are somewhat different from SRP2’s, suggesting overlapping yet distinct functions. We elucidated the 3D structure of SRP1 in solution by two-dimensional H-H NMR spectroscopy. The tertiary structure of SRP1 consists of two short β-strands at Y12−R15 and F18−F20, one type-II β-turn at R15−F18 in its well-defined core and is stabilized by a covalent disulfide bond between C8 and C19. The conformational ensemble of SRP1 from extensive atomistic simulation in explicit solvent (with 3.0 μs total effective sampling) shows high consistency with experimental intramolecular NOEs of the core region. The SRP1 core adopts a fold similar to the carboxyl-terminal subdomain of epidermal growth factor (EGF), suggesting that SRP1 may interact with EGF receptor-like molecules to trigger its biological function.


Introduction
Insects possess an effective innate immune system to fight against invading pathogens [1][2]. It is comprised of responses mediated by plasma factors made in fat body (pathogen recognition, melanization, antimicrobial peptides) and by blood cells known as hemocytes (phagocytosis, encapsulation). Extracellular serine protease pathways coordinate some of the responses, such as proteolytic activation of phenoloxidase, Spätzle, and plasmatocyte spreading peptide (PSP) precursors [3]. Phenoloxidase catalyzes the formation of reactive intermediates and melanin to kill and sequester pathogens, whereas Spätzle and PSP trigger the Toll pathway for antimicrobial peptide synthesis and cellular encapsulation, respectively [2,4]. In the holometabolous insects whose genomes are known, Spätzle and PSP homologs form two families of cytokines with as many as 7 and 11 members [5][6]. We are just beginning to understand the functions of these related peptides.
Because PSP and its homologs are mostly activated under stress conditions such as wounding, infection, ligation, and heat or cold shock, we chose stress response peptide (SRP) as a general term to describe these structurally related molecules [7]. These peptides stimulate hemocyte spreading and chemotaxis, block larval growth, induce antimicrobial peptide (AMP) synthesis, and regulate other physiological processes. In the tobacco hornworm Manduca sexta, two paralytic peptides [PP1 and PP2: ENFAGGCATG(Y/F)LRTADGRCKPTF] were isolated from hemolymph collected through a cut proleg [8]. Injection of chemically synthesized PP1 into larval hemocoel induced rapid paralysis, reduced bleeding, and stimulated plasmatocyte spreading and aggregation [9]. PP1 has an ordered structure comprised of a short antiparallel β-sheet at Y11−T14 and R18−P21, three β-turns at F3−G6, A8−Y11 and T14−G17, and a half turn at K20−F23 [10]. Genome analysis shows a single PP/PSP gene and ten homologs (uENF1, uENF2, SRP1−8) in M. sexta [6]. All of these genes encode predicted pre-pro-peptides, whose maturation is expected to require proteolytic cleavages. SRP2, which stimulates AMP expression, adopts a structure similar to PP1, and is composed of two short β-strands at R12−R15 and I18−V20, one type-I′ β-turn at R15−I18, and a half turn at C8−S10 [6]. The core structures of PP1 and SRP2 are each stabilized by a disulfide bond. To better understand structurefunction relationships of this family of cytokines, we investigated the immune inducibility, stress responsiveness, expression pattern of M. sexta SRP1 and carried out a twodimensional 1 H-1 H NMR to examine its solution structure. We observed that SRP1 has an ordered structure, including a short twisted anti-parallel β-sheet, a type-II β-turn, and is stabilized by one disulfide bond.

Rearing of Insects and Injection, Total RNA Preparation, and cDNA Synthesis
M. sexta eggs were purchased from Carolina Biological Supply (Burlington, NC) and larvae were reared on an artificial diet at 25°C [11]. Day 2, 5 th instar larvae were individually injected with a mixture of killed bacteria (2×10 7 Escherichia coli XL1-Blue and 20µg Micrococcus luteus, Sigma-Aldrich) and 20µg insoluble β-1,3-glucan from Alcaligenes faecalis (Sigma-Aldrich) suspended in water (30 µL per larva). Hemolymph was collected from cut prolegs of larvae (3 per group) after 24h and centrifuged at 5,000g for 5 min to harvest hemocytes. Fat body was then dissected from the bled larvae. Total RNA was prepared from pooled hemocytes and fat body using TRIZOL Reagent (Thermo Fisher Scientific) and was conducted in triplicate. Control hemocyte and fat body total RNA samples were prepared by the same method from day 3, 5 th instar naïve larvae (3 per group, 3 groups as biological replicates).
For tissue specific expression analysis, tissues were dissected from 5-10 day 3, 5 th instar naïve larvae (3-6 per group, 0 biological replicate) for RNA isolation. To study SRP1 expression after heat stress, day 2, 5th instar larvae (3 per group, 3 biological replicates) were incubated at 42°C for 1 h and then switched back to 25°C. Hemocytes, fat body and midgut tissues were collected 6 h later for total RNA preparation. The control insects (3 per group, 3 biological replicates) were day 2 naïve larvae continuously cultured at 25°C prior to hemocyte collection and tissue dissection. To examine developmental changes of SRP1 expression, hemocytes, fat body and midgut tissues were dissected from insects at stages ranging from day 2, 5 th instar larvae to day 1, pupae (3-6 insects per stage, 0 biological replicate) for total RNA preparation. Each total RNA sample (1 µg each) was incubated with 1× iScript Reverse Transcription Supermix (Bio-Rad) in a 10µL reaction at 42°C for 30min to synthesize cDNA. The reaction mixtures were then heated at 95°C for 3min to denature the enzyme.
To examine possible role of SRP1 in AMP expression regulation, 40µL of PBS (pH 7.4, 137mM NaCl, 2.7mM KCl, 10mM Na 2 HPO 4 , 1.8mM KH 2 PO4) or SRP1 (0.1µg/µL) in PBS were each injected into six day 2, 5 th instar larvae (3 control and 3 test). Fat body tissues were separately collected 6h later for RNA isolation and cDNA synthesis as described previously. This experiment was repeated twice (total of 3 biological replicates).

Synthesis of M. sexta Stress Response Peptide SRP1
M. sexta SRP1 was synthesized by automated protocols of stepwise solid-phase Fmoc chemistry on an automatic peptide synthesizer ABI model 431A (Applied Biosystems, Foster City, CA, USA). For peptide synthesis, CLEAR amide resin (0.3mmol/g; Peptides International, Louisville, KY) and N α -Fmoc amino acids (Anaspec Inc., San Jose, CA) were used. The linear peptides were cyclized by formation of a disulfide bond between C8 and C19 followed by purification on a high-pressure liquid chromatography (HPLC) system (Beckman Instruments, Inc., Fullerton, CA) using a Phenomenex reversed-phase C-18 column (Torrance, CA). The purified (>95%) peptide was eluted (1ml/min) using a linear gradient of 10-90% acetonitrile containing 0.1% trifluoroacetic acid. The characterization and mass Inducer of Antimicrobial Peptide Synthesis determination of the peptide was performed by matrixassisted-laser desorption time-of-flight (MALDI-TOF) mass spectrometric analysis (Bruker Ultraflex II spectrometer; Bruker Daltronics, Billerica, MA).

NMR Spectroscopy
The one-and two-dimension proton-proton NMR experiments described in these studies were carried out on 500 MHz Varian VNMR system (Agilent Technologies Inc., Palo Alto, CA) equipped with a 5mm cryogenic triple resonance inverse detection pulse field gradient probe. The NMR experimental data were acquired at 25°C. The 1.8 mM solution of the purified SRP1 in 650µL of 10% D 2 O and 90% triple distilled H 2 O was used for NMR data acquisition. A total of 4096 t2 and 512 t1 data points were acquired for the 2D 1 H-1 H Double Quantum Filtered Correlation Spectroscopy (DQF-COSY) [13] experiments with spectral width of 12ppm in each dimension, and 16 transient per increment. Whereas 2D 1 H-1 H Total Correlation Spectroscopy (TOCSY) [14] and Nuclear Overhauser Effect Spectroscopy (NOESY) [15] experiments were performed using 2048 t 2 and 256 t 1 data points. The MLEV-17 pulse sequence [14] with a spin-lock time of 80ms at B 1 field strength of 7kHz was used for acquisition of 2D 1 H-1 H TOCSY NMR data. Three different mixing times of 100, 300 and 500ms were used for acquisition 2D 1 H-1 H NOESY experimental data under similar conditions. All 2D NMR experiments were performed in hyper-complex phasesensitive mode. The 2D DQF-COSY NMR data sets were zero filled to 8K in the t2 dimension and 2K in the t 1 dimension respectively. For NMR data processing of all other 2D NMR experiments, the linear prediction and zero filling to 4K data points in t 2 dimension were used. When necessary, we used spectral Lorentzian-Gaussian apodization for enhancing the NMR spectral resolution. We used WATERGATE pulse scheme [16] during acquisition for suppression of the solvent peak (HOD). The residual water peak at 4.78ppm was used as reference for chemical shift assignments. NMR data was processed using VnmrJ3.2a (Agilent Technologies Inc., Palo Alto, CA) software and analysis was carried out using CCPNMR 2.4 [17] and Sparky 3.12 program [18]. The assignment of side chain proton resonances were achieved by comparing cross peaks in the 2D TOCSY spectra with those in the 2D NOESY spectra acquired under similar experimental conditions. The information for hydrogen bond was obtained through amide proton resonances either appearing after 12h or showing slow exchange in the 1D and 2D 1 H-1 H TOCSY spectra acquired in 100% D 2 O. The NOE cross peaks were classified as strong, medium, or weak intensity based on the observed number of contour lines for distance constraints calculation.

Structure Calculations
The inter-proton distance constraints were obtained using 2D NOESY spectra acquired with the different mixing times (100, 300 and 500ms). We used the 2D NOESY data acquired with mixing time of 300ms for distance constraint calculations and also confirmed the assigned cross peaks in the spectrum acquired with the 100ms mixing time. In order to rule out the peaks due to spin diffusion, we only considered cross peaks present in the NOESY spectra acquired with 100ms mixing time for distance constraint calculations. We assigned NOE-derived distance constraints into three distance ranges (1.8−2.5, 1.8−3.5, and 1.8−5.0Å) corresponding to strong, medium and weak cross peaks intensities. A correction to the upper bound distance constraint limits for methyl or non-stereo-specifically assigned methylene protons was done by adding 0.5Å appropriately for center averaging [19]. The CNS program version 1.3 [20] that uses a combination of simulated annealing in torsional space with restrained molecular dynamics (MD) in a Cartesian space was used for calculation of SRP1 structures with minimal restraint violation [21]. For structure calculations, the force constants were scaled throughout the procedure using default parameters of CNS. The molecular structure refinement protocol involved 100ps restrained MD simulation at 300K followed by energy minimization. For the NOE derived distance constraints and dihedral angles constraints the force constant value of 50kcal/mol/Å 2 and 100kcal/mol/Å 2 respectively were imposed during molecular dynamics calculations. The force constant of NOE derived distance constraints was increased to 100kcal/mol/Å 2 during final energy minimization. We calculated over 200 SRP1 structures using NOEs, dihedral angles and hydrogen bond derived distance constraints. An ensemble of three-dimensional NMR structures of SRP1 was obtained using 20 lowest energy structures with no restraint violation. The average structure was calculated to represent 3D solution NMR structure of M. sexta SRP1. The Visualized Molecular Dynamics (VMD) [22] and PROCHECK-NMR programs [23] were used to generate the figure and check the quality of the average structure respectively.

Explicit Solvent Molecular Dynamic (MD) Simulations
The 5 lowest energy NMR structures of SRP1 were used as the initial coordinates to initiate five independent MD simulations. Each of the 5 replicas was solvated in a cubic TIP3P explicit solvent water-box with an average volume of ~55 x 55 x 55Å 3 in CHARMM [24][25][26]. The system was then charge balanced with K + . Each of the solvated SRP1 structures was then energy minimized and equilibrated for 100 ps in the Charmm36 force field [27]. Five 600-ns production simulations at standard temperature and pressure were then performed using NAMD [28], together yielding a total of 3.0µs effective sampling of SRP1. Particle mesh Ewald was used for long-range electrostatic interactions and van der Waals interactions were smoothly switched off from 12 to 13Å [29]. Hydrogen mass was increased to 3 through bound heavy atom mass repartition [30]. The integration time-step was 4 fs, keeping all bonds to hydrogens atoms of fixed length with SHAKE [31]. Snapshots were saved at an interval of 10 ps for the duration of simulation time to be used for later analysis.

Clustering and NOE Analysis
All structural and clustering analysis were performed using the GROMACS program suite [32]. For clustering, snapshots at 1ns intervals were included, giving rise to 3004 unique snapshots together from all five replicas. GROMOS clustering was then employed using Cα RMSD with a cut-off of 3 Å [33]. Theoretical NOE distances were back calculated for each of the eight most populated clusters, which presumably represent key conformational states of SRP1. Back-calculated NOE-like distances were then compared to experimentally derived medium and long-range NOE upper bounds involving NH, CαH, CβH, and CγH protons used in NMR structure calculations (see Section 2.5).

Immune Inducibility, Tissue Specificity, Stress Responsiveness, and Developmental Profiles of SRP1 Expression
SRP1 is homologous to Hyphantria cunea Hdd23, a gene known to be upregulated by immune challenge, but with unknown function [34]. We examined the SRP1 mRNA levels in different tissues at various life stages represented by 52 previously obtained RNA-Seq datasets [35] ( Figure A1) and found that SRP1 is expressed at a lower level than SRP2 but substantially higher than SRP5, SRP7 and SRP8 (Cao and Jiang, unpublished data). Our previous transcriptome analysis showed that M. sexta SRP1 (CIFH contig 06581) was up-regulated 6.6-fold in fat body and 11.9-fold in hemocytes after immune challenge [36]. We further examined SRP1 expression profiles by qRT-PCR analysis. The mRNA levels of SRP1 were very low in hemocytes and fat body from naïve larvae and became significantly higher (p<0.05) in both tissues after larvae were injected with a mixture of bacteria and β-1,3-glucan ( Figure 1A). The SRP1 mRNA level in induced hemocytes was higher, relative to rpS3 mRNA, than that in induced fat body. In day 3, 5 th instar naïve larvae the SRP1 transcripts in integument and trachea were 3−4 fold more abundant than in nerve tissue and muscle ( Figure 1B). The mRNA levels in midgut and fat body were negligible. We tested if SRP1 expression is induced after larvae were heat shocked at 42°C for 1 h. At 6h after the treatment, we detected significant up-regulation of SRP1 expression in midgut and hemocytes ( Figure 1C). The highest mRNA level was detected in treated fat body but the p-value was >0.05 when compared with the control. We also measured the SRP1 mRNA levels in hemocytes, fat body and midgut during development from feeding stage 5 th instar larvae to pupae ( Figure 1D). The SRP1 transcripts in hemocytes increased three-fold in day one pupae, compared with expression in larval hemocytes. Two peaks of expression were detected in fat body, one in the early wandering stage (day 7) and the other at the early pupal stage (day 1). In similar stages, expression of AMP genes is also up-regulated in response to developmental signals [35]. 2,5,7,9,11 and pupa, day 1) were determined using the same method.

Induction of AMP Gene Expression by SRP1 Injection
In order to check whether SRP1 regulates AMP gene expression, we injected the SRP1 peptide solution (4.0 µg) or PBS as a control into day 2, 5 th instar larvae. We analyzed the test and control RNA samples by qRT-PCR and observed two of the seven AMP genes were significantly induced (p<0.05) at 6h after SRP1 injection (Figure 2). The diapausin-4 and attacins mRNA levels were 358.8 and 7.6-fold higher respectively than the control (p<0.05). These results suggest that SRP1 participates in regulation of immune responses in the tobacco hornworm and justifies a structural study of this cytokine [36].

Figure 2. Induction of AMP expression by injection of PBS (grey) or SRP1 (black). As described in Materials and methods, the test and control total RNA samples of fat body from the larvae injected with buffer or peptide were analyzed by qRT-PCR using the primers specific for the seven AMP genes.
After normalization against rpS3, the average SRP1 mRNA levels in the control samples were adjusted to one, and relative levels in the SRP1 injected group were calculated accordingly. The statistical significances were calculated using the Student's t-test and significant increases are shown as * (p<0.05) or ** (p<0.01)

NMR Resonance Assignments and Secondary Structure
Following the protocols published earlier [27][28][29], the analyses of 2D 1 H-1 H TOCSY and 2D 1 H-1 H NOESY NMR spectra provided intra-residue and inter-residue protonproton correlations on the amino acid spin system resonances and sequential proton chemical shift assignments respectively for SRP1. The TOCSY spectral data obtained at 25°C showed sufficient chemical shift dispersion for residues with limited spectral overlap. This enabled us to assign proton chemical shift of the spin systems resonances for all residues and indicated the presence of a well-defined secondary and tertiary structure for this peptide. In these assignments, CαH (i) -CδH (i+1: P) (dαδ) or CαH (i) -CαH (i+1: P) (dαα) NOEs instead of dαN(i, i+1) were used to assign trans or cis conformations for P9 and P21. We observed that both prolines (P9 and P21) in SRP1 showing strong dαδ NOEs suggesting trans configuration. The absence of dαα NOE cross peaks in the 2D 1 H-1 H NOESY spectra indicated a lack of any Pro in cis configuration. The proton chemical shifts assignments data for M. sexta SRP1 have been deposited in BioMagResBank (http://www.bmrb.wisc.edu) under the accession number 30507.  Figure 3 displays NH-CαH, NH-CβH, NH-CγH, and NH-CδH cross peak assignments of individual amino acid spin system in the fingerprint region of the TOCSY spectrum used for residue identification. The 2D 1 H-1 H NOESY spectral analysis allowed the identification of dαN(i, i+1) cross peaks for sequential assignments and the disambiguation of repeated amino acid residues. Further analysis of the NOESY spectra acquired with different mixing times characterized short, medium and long-range connectivity including dαN, dβN and dNN that provided structural information. An example of 2D 1 H-1 H NOESY spectrum fingerprint region showing dNN connectivity is provided in Supporting Information ( Figure A2). Figure 4 illustrates all the sequential and medium range distance NOE connectivity as well as deuterium exchange resistant data for residues of SRP1. A number of non-sequential [weak dNN(i, i+3) NOE between R15 and F18 and weak dNN(i, i+7) NOE between V13 and F20] NOEs suggested an antiparallel β-sheet arrangement. Residues Y12, V13, R14, R15, F18, C19 and F20 were likely to form two short β-strands. This secondary structural observation was further supported by the down field shift of Cα proton resonances in chemical shift index (∆CαH) values ( Figure A3) and larger (9.5 -15Hz) coupling constant ( 3 JNH-CαH) values for these residues. The random coil values for the chemical shift of CαH protons of residues were taken from the published data [30]. The SRP1 structure consisted of two short β-strands at region Y12-R15 [weak dαN(i, i+1) NOEs from Y12 to R14] and F18-F20 [medium dαN(i, i+1) NOE between C19 and F20]. These β-strands flanked a type-II β-turn between residues R15 -F18 [weak dNN(i, i+3) NOE from R15 to F18; weak dαN(i, i+1) NOE from residues R15 to F18; weak dNN(i, i+1) NOE between G16 and T17; medium dNN(i, i+1) NOE between T17 and F18]. The short stretch of twisted anti-parallel β-sheet was confirmed by observation of hydrogen bonds between V13 -F20. Low number of NOE cross peaks for N-terminal residues F1 -T7 suggested an extended conformation for this region. However, the C-terminal residues P21 -Y25 displayed turn-like characteristics [weak dαN(i, i+3) NOE between P21 and D24 and dNN(i, i+2) NOE between D23 and Y25].

Tertiary Structure of M. sexta SRP1
We determined the tertiary structure of SRP1 by hybrid distance geometry-simulated annealing molecular dynamics calculations using 284 constraints including 139 intraresidue, 63 sequential, 29 medium range, 21 long range NOE and 6 hydrogen bond distance constraints as well as 22 dihedral angle constraints. Two additional distance constraints were used between Sγ and Cβ atoms of C8 and C19 to account for normal disulfide bond geometry (2.0 Å < dS-S < 2.5 Å, 2.2 Å < dCβ-S < 3.5 Å). Table 1 summarized structural statistics of the 20 lowest energy structures for SRP1 in water. The calculated structures converged to a single family of conformers regardless of the starting structure. For these structure calculations, randomly chosen structures with extended conformation were selected as starting structures. From a total of 200 conformers calculated, 20 structures with lowest energy were selected for further analysis. The torsion angle dynamics of all the 20 selected conformers were consistent with a pair of hydrogen bonds involving the backbone carbonyl of V13 and amide proton of F20, and vice versa. Conformational analysis of these 20 structures suggested two additional hydrogen bonds (F18 backbone carbonyl to backbone amide proton of R15, and R15 backbone carbonyl to backbone amide proton of F18). Figure 5A displays superimposition of the 20 lowest energy structure of SRP1. Alignments along residues C8-P21, forming the core region for SRP1, of the 20 lowest energy structures showed backbone root mean square deviation (RMSD) value of 0.28Å.
These structures showed good agreement with experimental data with no distance violation larger than 0.3Å and no dihedral angle violation larger than 5°. A Ramachandran plot analysis using Procheck-NMR [23] showed that 79.5% of residues were in most favored stereochemically allowed region and 18.4% were in additional allowed regions ( Table 1). The schematic cartoon diagram of the energy minimized average tertiary structure of M. sexta SRP1's core region (residues C8-P21) is shown in Figure 5B & C. Two hydrogen bonds that stabilized the tertiary structure of SRP1 were confirmed by the slowly exchanged amide protons from V13 and F20.  As a means of structural refinement of the calculated NMR ensemble, 600-ns MD simulations seeded with the five lowest energy structures from NMR were conducted in explicit solvent. A total of 3004 unique snapshots taken at 1ns intervals were collected from these simulations. The average structures and populations of the eight largest clusters of the simulated ensemble are shown in Figure 6, which shows a strong degree of agreement in positions and relative integrity of the β-strands throughout simulation with respect to the NMR structure. It also elucidates some conformational heterogeneity present for both N-and Ctermini. We have further performed NOE violation analysis on the sub-ensembles of structures that belong to the eight largest clusters. The NOE violation analysis results are summarized in Figure 7. These results suggest that constraints pertaining to the core region residues (C8-P21) are in strong agreement with NOE observations from NMR. However, there are NOEs that only satisfy in smaller subpopulations of the ensemble pertaining to residues D23 and D24. This likely reflects SRP1's conformational hypervariability in the C-terminus.  The analysis of the tertiary structure of SRP1 suggested that in aqueous solution this peptide possessed a very welldefined core region with a short twisted anti-parallel β-sheet that was stabilized by a type-II β-turn and a disulfide bond. The C-terminus displayed turn-like characteristics, possibly due to electrostatic repulsion of the aspartate residues, while the N-terminus was largely unstructured with substantial flexibility. Overall tertiary structure was stabilized by a combination of the hydrogen bonds within the β-sheet structure and the covalent disulfide linkage. Within the core region of the molecule, one side of the structure was mainly composed of charged residues R14 and R15. The opposite side of the structure was composed of hydrophobic residues, including Y12, F18 and F20. The side chain of Y12 showed NOE contacts to the side chains of C8, P9, and C19 (data not shown). This suggested residue Y12 is embedded, as observed for residues at a similar position in MsPP1 and MsSRP2 [6,10]. These structural characteristics, similar to other ENF peptides [10], indicated that the SRP1 core region adopted a fold very similar to the carboxyl-terminal region of the vertebrate EGF, though they had very low sequence similarities. Interestingly, plasmatocyte-spreading peptide 1 (PSP1) from insect Pseudoplusia includens showing very similar structural characteristics of its core region also adopted a fold with close similarity to the C-terminal subdomain of EGF-like module 5 of human thrombomodulin (hTM5) [37]. The sequence identity between PSP1 and SRP1 was limited to mainly cysteine, glycine, tyrosine residues and structures of human epidermal growth factor (hEGF) and the hTM5 shared the general structural features of the family of EGF-like domains [38,39]. These findings suggested that SRP1 interactions with its receptor may be analogous to EGF-receptor interactions to trigger the biological function of this stress response peptide.

Conclusion
In the M. sexta genome there are 11 encoded small prepropeptides in the SRP family that function as cytokines. Appropriate maturation of the peptide involves a combination of intracellular and extracellular serine proteases. The tissue expression profile of the SRP1 gene during development shows discrete tissue specific and temporal regulation of the mRNA transcripts. Induction of SRP1 gene expression by biotic stressors was observed in hemocytes and fat body, while heat shock induced gene expression in hemocytes and midgut. Injection of mature SRP1 demonstrates that it stimulates expression of specific AMP genes. These targeted genes are unique and somewhat different from AMP genes induced by SRP2, indicating a distinct function of SRP1. Our results suggest the existence of a proSRP1 activating protease(s), an SRP1 receptor, and a unique role of SRP1 within the immune signal transduction pathway. The structural information of SRP1 will provide insight into our understanding of the interactions with its putative receptor.  [40]. Bars represent difference between published values and values observed during NMR acquisition for SRP1.