Pepsin-like Aspartic Proteases Research Articles

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin.

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

Toxoplasma gondii aspartic protease 5: structural basis of substrate binding and inhibition mechanism

Toxoplasma gondii, a worldwide prevalent parasite is responsible for causing toxoplasmosis in almost all warm-blooded animals, including humans. Golgi-resident T. gondii aspartic protease 5 (TgASP5) plays an essential role in the maturation and export of the effector proteins those modulate the host immune system to establish a successful infection. Hence, inhibiting this enzyme can be a possible therapeutic strategy against toxoplasmosis. This is the first report of the detailed structural investigations of the TgASP5 mature enzyme using molecular modeling and an all-atom simulation approach which provide in-depth knowledge of the active site architecture of TgASP5. The analysis of the binding mode of the TEXEL (Toxoplasma EXport Element) substrate to TgASP5 highlighted the importance of the active site residues. Ser505, Ala776 and Tyr689 in the S2 binding pocket are responsible for the specificity towards Arg at the P2 position of TEXEL substrate. The molecular basis of inhibition by the only known inhibitor RRLStatine has been identified, and our results show that it blocks the active site by forming a hydrogen bond with a catalytic aspartate. Besides that, known aspartic protease inhibitors were screened against TgASP5 using docking, MD simulations and MM-PBSA binding energy calculations. The top-ranked inhibitors (SC6, ZY1, QBH) showed higher binding energy than RRLStatine. Understanding the structural basis of substrate recognition and the binding mode of these inhibitors will help to develop potent mechanistic inhibitors against TgASP5. This study will also provide insights into the structural features of pepsin-like aspartic proteases from other apicomplexan parasites for developing antiparasitic agents. Communicated by Ramaswamy H. Sarma

Computer-Aided Prediction of the Interactions of Viral Proteases with Antiviral Drugs: Antiviral Potential of Broad-Spectrum Drugs.

Human society is facing the threat of various viruses. Proteases are promising targets for the treatment of viral infections. In this study, we collected and profiled 170 protease sequences from 125 viruses that infect humans. Approximately 73 of them are viral 3-chymotrypsin-like proteases (3CLpro), and 11 are pepsin-like aspartic proteases (PAPs). Their sequences, structures, and substrate characteristics were carefully analyzed to identify their conserved nature for proposing a pan-3CLpro or pan-PAPs inhibitor design strategy. To achieve this, we used computational prediction and modeling methods to predict the binding complex structures for those 73 3CLpro with 4 protease inhibitors of SARS-CoV-2 and 11 protease inhibitors of HCV. Similarly, the complex structures for the 11 viral PAPs with 9 protease inhibitors of HIV were also obtained. The binding affinities between these compounds and proteins were also evaluated to assess their pan-protease inhibition via MM-GBSA. Based on the drugs targeting viral 3CLpro and PAPs, repositioning of the active compounds identified several potential uses for these drug molecules. As a result, Compounds 1-2, modified based on the structures of Ray1216 and Asunaprevir, indicate potential inhibition of DENV protease according to our computational simulation results. These studies offer ideas and insights for future research in the design of broad-spectrum antiviral drugs.

Artificial intelligence based de-novo design for novel Plasmodium falciparum plasmepsin (PM) X inhibitors

Plasmodium falciparum is the leading cause of malaria with 627,000 deaths annually. Invasion and egress are critical stages for successful infection of the host yet depend on proteins that are extensively pre-processed by various maturases. Plasmepsins (Plasmodium pepsins, abbreviated PM, I-X) are pepsin-like aspartic proteases that are involved in almost all stages of the life cycle. The goal of this study was to use de-novo generative modeling techniques to create novel potential PfPMX inhibitors. A total of 4325 compounds were virtually screened by structural-based docking methods. The obtained hits were utilized to refine a structure-based Ligand Neural Network (L-Net) generative model to generate related compounds. The obtained optimal L-Net Compounds with smina scores ≤ −5.00KCalmol-1 and QED ≥ 0.35 were further taken for amplification utilizing Ligand Based Transformer modeling using Deep generative learning (Drug Explorer/DrugEx). The resulting hits were then subjected to XP Glide conventional Molecular docking and QikProp ADMET screening; molecules with XP Docking score ≤ −7.00KCalmol-1 were retained. Based on their Glide ligand efficiency, originality, and uniqueness, 30 compounds were chosen for binding affinity and MM_GBSA energy determination. Following Induced Fit docking (IFD), 7 compounds were taken for 50 ns MD simulations and FEP/MD calculations. This study reported novel potential PfPMX inhibitors with acceptable ADMET profiles and reasonable synthetic accessibility scores, as well as sufficient docking scores against other PMs were generated. The PfPMX inhibitors reported in this article are promising antimalarials for the next stages of drug development, and the first of their kind to be investigated thoroughly. Communicated by Ramaswamy H. Sarma

Histo-aspartic protease (HAP): a unique pepsin-like aspartic protease with catalytic histidine and aspartate

Histo-aspartic protease (HAP): a unique pepsin-like aspartic protease with catalytic histidine and aspartate

Bioprospection of Selected Plant Secondary Metabolites as Modulators of the Proteolytic Activity of Plasmodium falciparum Plasmepsin V.

Malaria is a devastating disease, and its management is only achieved through chemotherapy. However, resistance to available medication is still a challenge; therefore, there is an urgent need for the discovery and development of therapeutics with a novel mechanism of action to counter the resistance scourge consistent with the currently available antimalarials. Recently, plasmepsin V was validated as a therapeutic target for the treatment of malaria. The pepsin-like aspartic protease anchored in the endoplasmic reticulum is responsible for the trafficking of parasite-derived proteins to the erythrocytic surface of the host cells. In this study, a small library of compounds was preliminarily screened in vitro to identify novel modulators of Plasmodium falciparum plasmepsin V (PfPMV). The results obtained revealed kaempferol, quercetin, and shikonin as possible PfPMV inhibitors, and these compounds were subsequently probed for their inhibitory potentials using in vitro and in silico methods. Kaempferol and shikonin noncompetitively and competitively inhibited the specific activity of PfPMV in vitro with IC50 values of 22.4 and 43.34 μM, respectively, relative to 62.6 μM obtained for pepstatin, a known aspartic protease inhibitor. Further insight into the structure-activity relationship of the compounds through a 100 ns molecular dynamic (MD) simulation showed that all the test compounds had a significant affinity for PfPMV, with quercetin (-36.56 kcal/mol) being the most prominent metabolite displaying comparable activity to pepstatin (-35.72 kcal/mol). This observation was further supported by the compactness and flexibility of the resulting complexes where the compounds do not compromise the structural integrity of PfPMV but rather stabilized and interacted with the active site amino acid residues critical to PfPMV modulation. Considering the findings in this study, quercetin, kaempferol, and shikonin could be proposed as novel aspartic protease inhibitors worthy of further investigation in the treatment of malaria.

The nepenthesin insert in the Plasmodium falciparum aspartic protease plasmepsin V is necessary for enzyme function

Plasmepsin V (PM V) is a pepsin-like aspartic protease essential for growth of the malarial parasite Plasmodium falciparum. Previous work has shown PM V to be an endoplasmic reticulum–resident protease that processes parasite proteins destined for export into the host cell. Depletion or inhibition of the enzyme is lethal during asexual replication within red blood cells as well as during the formation of sexual stage gametocytes. The structure of the Plasmodium vivax PM V has been characterized by X-ray crystallography, revealing a canonical pepsin fold punctuated by structural features uncommon to secretory aspartic proteases; however, the function of this unique structure is unclear. Here, we used parasite genetics to probe these structural features by attempting to rescue lethal PM V depletion with various mutant enzymes. We found an unusual nepenthesin 1-type insert in the PM V gene to be essential for parasite growth and PM V activity. Mutagenesis of the nepenthesin insert suggests that both its amino acid sequence and one of the two disulfide bonds that undergird its structure are required for the insert’s role in PM V function. Furthermore, molecular dynamics simulations paired with Markov state modeling suggest that mutations to the nepenthesin insert may allosterically affect PM V catalysis through multiple mechanisms. Taken together, these data provide further insights into the structure of the P. falciparum PM V protease.

Flap Dynamics in Pepsin-Like Aspartic Proteases: A Computational Perspective Using Plasmepsin-II and BACE-1 as Model Systems.

The flexibility of β hairpin structure known as the flap plays a key role in catalytic activity and substrate intake in pepsin-like aspartic proteases. Most of these enzymes share structural and sequential similarity. In this study, we have used apo Plm-II and BACE-1 as model systems. In the apo form of the proteases, a conserved tyrosine residue in the flap region remains in a dynamic equilibrium between the normal and flipped states through rotation of the χ1 and χ2 angles. Independent MD simulations of Plm-II and BACE-1 remained stuck either in the normal or flipped state. Metadynamics simulations using side-chain torsion angles (χ1 and χ2 of tyrosine) as collective variables sampled the transition between the normal and flipped states. Qualitatively, the two states were predicted to be equally populated. The normal and flipped states were stabilized by H-bond interactions to a tryptophan residue and to the catalytic aspartate, respectively. Further, mutation of tyrosine to an amino-acid with smaller side-chain, such as alanine, reduced the flexibility of the flap and resulted in a flap collapse (flap loses flexibility and remains stuck in a particular state). This is in accordance with previous experimental studies, which showed that mutation to alanine resulted in loss of activity in pepsin-like aspartic proteases. Our results suggest that the ring flipping associated with the tyrosine side-chain is the key order parameter that governs flap dynamics and opening of the binding pocket in most pepsin-like aspartic proteases.

Structures of plasmepsin X from Plasmodium falciparum reveal a novel inactivation mechanism of the zymogen and molecular basis for binding of inhibitors in mature enzyme.

Plasmodium falciparum plasmepsin X (PfPMX), involved in the invasion and egress of this deadliest malarial parasite, is essential for its survival and hence considered as an important drug target. We report the first crystal structure of PfPMX zymogen containing a novel fold of its prosegment. A unique twisted loop from the prosegment and arginine 244 from the mature enzyme is involved in zymogen inactivation; such mechanism, not previously reported, might be common for apicomplexan proteases similar to PfPMX. The maturation of PfPMX zymogen occurs through cleavage of its prosegment at multiple sites. Our data provide thorough insights into the mode of binding of a substrate and a potent inhibitor 49c to PfPMX. We present molecular details of inactivation, maturation, and inhibition of PfPMX that should aid in the development of potent inhibitors against pepsin-like aspartic proteases from apicomplexan parasites.

Improving the alkaline stability of pepsin through rational protein design using renin, an alkaline-stable aspartic protease, as a structural and functional reference

The present study sought to identify the structural determinants of aspartic protease structural stability and activity at elevated pH. Various hypotheses have been published regarding the features responsible for the unusual alkaline structural stability of renin, however, few structure-function studies have verified these claims. Using pepsin as a model system, and renin as a template for functional and structural alkaline stability, a rational re-design of pepsin was undertaken to identify residues contributing to the alkaline instability of pepsin-like aspartic proteases in regards to both structure and function. We constructed 13 mutants based on this strategy. Among them, mutants D159 L and D60A led to an increase in activity at elevated pH levels (p ≤ 0.05) and E4V and H53F were shown to retain native-like structure at elevated pH (p ≤ 0.05). Previously suggested carboxyl groups Asp11, Asp118, and Glu13 were individually shown not to be responsible for the structural instability or lack of activity at neutral pH in pepsin. The importance of the β-barrel to structural stability was highlighted as the majority of the stabilizing residues identified, and 39% of the weakly conserved residues in the N-terminal lobe, were located in β-sheet strands of the barrel. The results of the present study indicate that alkaline stabilization of pepsin will require reduction of electrostatic repulsions and an improved understanding of the role of the hydrogen bonding network of the characteristic β-barrel.

Pepsin-like aspartic proteases (PAPs) as model systems for combining biomolecular simulation with biophysical experiments.

Pepsin-like aspartic proteases (PAPs) are a class of aspartic proteases which shares tremendous structural similarity with human pepsin. One of the key structural features of PAPs is the presence of a β-hairpin motif otherwise known as flap. The biological function of the PAPs is highly dependent on the conformational dynamics of the flap region. In apo PAPs, the conformational dynamics of the flap is dominated by the rotational degrees of freedom associated with χ1 and χ2 angles of conserved Tyr (or Phe in some cases). However it is plausible that dihedral order parameters associated with several other residues might play crucial roles in the conformational dynamics of apo PAPs. Due to their size, complexities associated with conformational dynamics and clinical significance (drug targets for malaria, Alzheimer's disease etc.), PAPs provide a challenging testing ground for computational and experimental methods focusing on understanding conformational dynamics and molecular recognition in biomolecules. The opening of the flap region is necessary to accommodate substrate/ligand in the active site of the PAPs. The BIG challenge is to gain atomistic details into how reversible ligand binding/unbinding (molecular recognition) affects the conformational dynamics. Recent reports of kinetics (Ki, Kd) and thermodynamic parameters (ΔH, TΔS, and ΔG) associated with macro-cyclic ligands bound to BACE1 (belongs to PAP family) provide a perfect challenge (how to deal with big ligands with multiple torsional angles and select optimum order parameters to study reversible ligand binding/unbinding) for computational methods to predict binding free energies and kinetics beyond typical test systems e.g. benzamide–trypsin. In this work, i reviewed several order parameters which were proposed to capture the conformational dynamics and molecular recognition in PAPs. I further highlighted how machine learning methods can be used as order parameters in the context of PAPs. I then proposed some open ideas and challenges in the context of molecular simulation and put forward my case on how biophysical experiments e.g. NMR, time-resolved FRET etc. can be used in conjunction with biomolecular simulation to gain complete atomistic insights into the conformational dynamics of PAPs.

Comparative bioinformatic and structural analyses of pepsin and renin

Pepsin, the archetypal pepsin-like aspartic protease, is irreversibly denatured when exposed to neutral pH conditions whereas renin, a structural homologue of pepsin, is fully stable and optimally active in the same conditions despite sharing highly similar enzyme architecture. To gain insight into the structural determinants of differential aspartic protease pH stability, the present study used comparative bioinformatic and structural analyses. In pepsin, an abundance of polar and aspartic acid residues were identified, a common trait with other acid-stable enzymes. Conversely, renin was shown to have increased levels of basic amino acids. In both pepsin and renin, the solvent exposure of these charged groups was high. Having similar overall acidic residue content, the solvent-exposed basic residues may allow for extensive salt bridge formation in renin, whereas in pepsin, these residues are protonated and serve to form stabilizing hydrogen bonds at low pH. Relative differences in structure and sequence in the turn and joint regions of the β-barrel and ψ-loop in both the N- and C-terminal lobes were identified as regions of interest in defining divergent pH stability. Compared to the structural rigidity of renin, pepsin has more instability associated with the N-terminus, specifically the B/C connector. By contrast, renin exhibits greater C-terminal instability in turn and connector regions. Overall, flexibility differences in connector regions, and amino acid composition, particularly in turn and joint regions of the β-barrel and ψ-loops, likely play defining roles in determining pH stability for renin and pepsin.

Activation mechanism of plasmepsins, pepsin-like aspartic proteases from Plasmodium, follows a unique trans-activation pathway.

Plasmodium parasites that cause malaria produce plasmepsins (PMs), pepsin-like aspartic proteases that are important antimalarial drug targets due to their role in host hemoglobin degradation. The enzymes are synthesized as inactive zymogens (pro-PMs), and the mechanism of their conversion to the active, mature forms has not been clearly elucidated. Our structural investigations of vacuolar pro-PMs with truncated prosegment (pro-tPMs) reveal that the formation of the S-shaped dimer is their innate property. Further structural studies, biochemical analysis, and molecular dynamics simulations indicate that disruption of the Tyr-Asp loop (121p-4), coordinated with the movement of the loop L1 (237-247) and helix H2 (101p-113p), is responsible for the extension of the pro-mature region (harboring the cleavage site). Consequently, under acidic pH conditions, these structural changes result in the dissociation of the dimers to monomers and the protonation of the residues in the prosegment prompts its unfolding. Subsequently, we demonstrated that the active site of the monomeric pro-tPMs with the unfolded prosegment is accessible for peptide substrate binding; in contrast, the active site is blocked in folded prosegment form of pro-tPMs. Thus, we propose a novel mechanism of auto-activation of vacuolar pro-tPMs that under acidic conditions can form a catalytically competent active site. One monomer cleaves the prosegment of the other one through a trans-activation process, resulting in formation of mature enzyme. As a result, once a mature enzyme is generated, it leads to the complete conversion of all the inactive pro-tPMs to their mature form. DATABASE: Atomic coordinates and structure factors have been submitted in the Protein Data Bank (PDB) under the PDB IDs 6KUB, 6KUC, and 6KUD.

Understanding the structural basis of substrate recognition by Plasmodium falciparum plasmepsin V to aid in the design of potent inhibitors.

Plasmodium falciparum plasmepsin V (PfPMV) is an essential aspartic protease required for parasite survival, thus, considered as a potential drug target. This study reports the first detailed structural analysis and molecular dynamics simulation of PfPMV as an apoenzyme and its complexes with the substrate PEXEL as well as with the inhibitor saquinavir. The presence of pro-peptide in PfPMV may not structurally hinder the formation of a functionally competent catalytic active site. The structure of PfPMV-PEXEL complex shows that the unique positions of Glu179 and Gln222 are responsible for providing the specificity of PEXEL substrate with arginine at P3 position. The structural analysis also reveals that the S4 binding pocket in PfPMV is occupied by Ile94, Ala98, Phe370 and Tyr472, and therefore, does not allow binding of pepstatin, a potent inhibitor of most pepsin-like aspartic proteases. Among the screened inhibitors, the HIV-1 protease inhibitors and KNI compounds have higher binding affinities for PfPMV with saquinavir having the highest value. The presence of a flexible group at P2 and a bulky hydrophobic group at P3 position of the inhibitor is preferred in the PfPMV substrate binding pocket. Results from the present study will aid in the design of potent inhibitors of PMV.

The prosegment catalyzes native folding of Plasmodium falciparum plasmepsin II

Plasmepsin II is a malarial pepsin-like aspartic protease produced as a zymogen containing an N-terminal prosegment domain that is removed during activation. Despite structural similarities between active plasmepsin II and pepsin, their prosegments adopt different conformations in the respective zymogens. In contrast to pepsinogen, the proplasmepsin II prosegment is 80 residues longer, contains a transmembrane region and is non-essential for recombinant expression in an active form, thus calling into question the prosegment's precise function. The present study examines the role of the prosegment in the folding mechanism of plasmepsin II. Both a shorter (residues 77–124) and a longer (residues 65–124) prosegment catalyze plasmepsin II folding at rates more than four orders of magnitude faster compared to folding without prosegment. Native plasmepsin II is kinetically trapped and requires the prosegment both to catalyze folding and to shift the folding equilibrium towards the native conformation. Thus, despite low sequence identity and distinct zymogen conformations, the folding landscapes of plasmepsin II and pepsin, both with and without prosegment, are qualitatively identical. These results imply a conserved and unusual feature of the pepsin-like protease topology that necessitates prosegment-assisted folding.

Structure of RC1339/APRc from Rickettsia conorii, a retropepsin-like aspartic protease.

The crystal structures of two constructs of RC1339/APRc from Rickettsia conorii, consisting of either residues 105-231 or 110-231 followed by a His tag, have been determined in three different crystal forms. As predicted, the fold of a monomer of APRc resembles one-half of the mandatory homodimer of retroviral pepsin-like aspartic proteases (retropepsins), but the quaternary structure of the dimer of APRc differs from that of the canonical retropepsins. The observed dimer is most likely an artifact of the expression and/or crystallization conditions since it cannot support the previously reported enzymatic activity of this bacterial aspartic protease. However, the fold of the core of each monomer is very closely related to the fold of retropepsins from a variety of retroviruses and to a single domain of pepsin-like eukaryotic enzymes, and may represent a putative common ancestor of monomeric and dimeric aspartic proteases.

Foldase and inhibitor functionalities of the pepsinogen prosegment are encoded within discrete segments of the 44 residue domain

Pepsin is initially produced as the zymogen pepsinogen, containing a 44 residue prosegment (PS) domain. When folded without the PS, pepsin forms a thermodynamically stable denatured state (refolded pepsin, Rp). To guide native folding, the PS binds to Rp, stabilizes the folding transition state, and binds tightly to native pepsin (Np), thereby driving the folding equilibrium to favor the native state. It is unknown whether these functionalities of the PS are encoded within the entire sequence or within discrete segments. PS residues 1p–29p correspond to a highly conserved region in pepsin-like aspartic proteases and we hypothesized that this segment is critical to PS-catalyzed folding. This notion was tested in the present study by characterizing the ability of various truncated PS peptides to bind Rp, catalyze folding from Rp to Np, and to inhibit Np. Four PS truncations were examined, corresponding to PS residues 1p–16p (PS1–16), 1p–29p (PS1–29), 17p–44p (PS17–44) and 30p–44p (PS30–44). The three PS functionalities could be ascribed primarily to discrete regions within the highly conserved motif: 1p–16p dictated Rp binding, 17p–29p dictated Np binding/inhibition, while the entire 1p–29p dictated transition state binding/catalyzing folding. Conversely, PS30–44 played no obvious role in PS-catalyzed folding; it is hypothesized that this more variable region may serve as a linker between PS1–29 and the mature domain. The high sequence conservation of PS1–29 and its role in catalyzing pepsin folding strongly suggest that there is a conserved PS-catalyzed folding mechanism shared by pepsin-like aspartic proteases with this motif.

Structure of γ-conglutin: insight into the quaternary structure of 7S basic globulins from legumes.

γ-Conglutin from lupin seeds is an unusual 7S basic globulin protein. It is capable of reducing glycaemia in mammals, but the structural basis of this activity is not known. γ-Conglutin shares a high level of structural homology with glycoside hydrolase inhibitor proteins, although it lacks any kind of inhibitory activity against plant cell-wall degradation enzymes. In addition, γ-conglutin displays a less pronounced structural similarity to pepsin-like aspartic proteases, but it is proteolytically dysfunctional. Only one structural study of a legume 7S basic globulin, that isolated from soybean, has been reported to date. The quaternary assembly of soybean 7S basic globulin (Bg7S) is arranged as a cruciform-shaped tetramer comprised of two superposed dimers. Here, the crystal structure of γ-conglutin isolated from Lupinus angustifolius seeds (LangC) is presented. The polypeptide chain of LangC is post-translationally cleaved into α and β subunits but retains its covalent integrity owing to a disulfide bridge. The protomers of LangC undergo an intricate quaternary assembly, resulting in a ring-like hexamer with noncrystallographic D3 symmetry. The twofold-related dimers are similar to those in Bg7S but their assembly is different as a consequence of mutations in a β-strand that is involved in intermolecular β-sheet formation in γ-conglutin. Structural elucidation of γ-conglutin will help to explain its physiological role, especially in the evolutionary context, and will guide further research into the hypoglycaemic activity of this protein in humans, with potential consequences for novel antidiabetic therapies.

Allosteric inhibition of BACE1 by an exosite-binding antibody

β-Secretase (BACE1) is a membrane-anchored pepsin-like aspartic protease and is the rate-limiting enzyme in the β-amyloidogenic pathway. Thus, inhibitors of BACE1 activity have therapeutic potential for Alzheimer's disease. While much effort has focused on small molecule active site inhibitors, recent exploration of BACE1 inhibition by peptides and antibodies has revealed exosites that can regulate enzymatic activity. This type of allosteric regulation by proteinaceous factors, while frequently found in serine and cysteine proteases, is rarely seen in aspartic proteases. A crystal structure of the anti-BACE1/enzyme complex shows altered structural features and dynamic characteristics near the substrate-binding cleft. This binding mode, along with the enzymatic inhibition pattern, suggests that anti-BACE1 functions through an allosteric inhibition mechanism.

Pepsin-like aspartic protease (Sc-ASP155) cloning, molecular characterization and gene expression analysis in developmental stages of nematode Steinernema carpocapsae

Steinernema carpocapsae is an insect parasitic nematode associated with the bacterium Xenorhabdus nematophila. These symbiotic complexes are virulent against the insect host. Many protease genes were shown previously to be induced during parasitism, including one predicted to encode an aspartic protease, which was cloned and analyzed in this study. A cDNA encoding Sc-ASP155 was cloned based on the EST fragment. The full-length cDNA of Sc-ASP155 consists of 955 nucleotides with multiple domains, including a signal peptide (aa1–15), a pro-peptide region (aa16–45), and a typical catalytic aspartic domain (aa71–230). The putative 230 amino acid residues have a calculated molecular mass of 23,812Da and a theoretical pI of 5.01. Sc-ASP155 blastp analysis showed 40–62% amino acid sequence identity to aspartic proteases from parasitic and free-living nematodes. Expression analysis showed that the sc-asp155 gene was up-regulated during the initial parasitic stage, especially in L3 gut and 6h induced nematodes. Sequence comparison revealed that Sc-ASP155 was a member of an aspartic protease family and phylogenetic analysis indicated that Sc-ASP155 was clustered with Sc-ASP113. In situ hybridization showed that sc-asp155 was expressed in subventral cells. Additionally, we determined that sc-asp155 is a single-copy gene in S. carpocapsae. Homology modeling showed that Sc-ASP155 adopts a typical aspartic protease structure. The up-regulated Sc-ASP155 expression revealed that this protease could play a role in the parasitic process. In this study, we have cloned the gene and determined the expression of the pepsin-like aspartic protease Sc-ASP155 in S. carpocapsae.