Modeling and phylogenetic analysis of cytosolic ascorbate peroxidase (OsAPX1) from rice reveal signature motifs that may play a role in stress tolerance

Ascorbate peroxidase (APX) is a crucial, haeme-containing enzyme of the ascorbate glutathione cycle that detoxifies reactive oxygen species in plants by catalyzing the conversion of hydrogen peroxide to water using ascorbate as a specific electron donor. Different APX isoforms are present in discrete subcellular compartments in rice and their expression is stress regulated. We revealed the homology model of OsAPX1 protein using the crystal structure of soybean GmAPX1 (PDB ID: 2XIF) as template by Modeller 9.12. The resultant OsAPX1 model structure was refined by PROCHECK, ProSA, Verify3D and RMSD that indicated the model structure is reliable with 83 % amino acid sequence identity with template, RMSD (1.4 Å), Verify3D (86.06 %), Zscores (-8.44) and Ramachandran plot analysis showed that conformations for 94.6% of amino acid residues are within the most favoured regions. Investigation revealed two conserved signatures for haeme ligand binding and peroxidase activity in the alpha helical region that may play a significant role during stress.     

In silico exploration of inhibition potency of maslinic acid,a plant based triterpene,against Helicoverpa armigera serine protease

Gut proteases are accountable for survival of Helicoverpa armigera on protein rich parts of plant devastating many important agricultural crops. The aim of present study was to identify potential natural compounds having inhibitory potency against Helicoverpa armigera gut proteases. We have modeled structure of H. armigera serine protease (UniProt ID: O18447) and analyzed its interactions with maslinic acid (Zinc ID: ZINC38140521). A 3D model was generated using bovine trypsin in complex with analogues of sunflower inhibitor 1 as template with the help of Chimera Modeler 1.11. The PROCHECK and Modfold analysis have revealed 81.8% of residue in favored region. The POOL and COACH analysis have revealed 18 amino acids in the active site. In the 10 ns MD simulations of modeled structure, the RMSD of the protein backbone increased slightly and later stabilized from 7 ns to 10 ns. The modeled structure was stabilized at gyration distance of about 1.65 nm at 7 ns. Potential hit compounds from the ZINC database identified in this study showed good inhibitory bindings with modeled structure. Among these compounds maslinic acid, a plant based pentacyclic triterpenes was found to be potent lead compound with good binding affinity (−9.5 kcal/mol). RMSD profile was <0.45 nm for complex with stabilization at about 18,000 ps (18 nm) suggesting stable interaction. This work demonstrates reasonable in silico inhibitory action of maslinic acid against H. armigera serine protease and depicts utility of in silico methodologies for designing competent strategies against dreaded insect pests like H. armigera.     

Novel oxidatively stable subtilisin-like serine proteases from alkaliphilic Bacillus spp.: enzymatic properties, sequences, and evolutionary relationships

The genes for five subtilisin-like serine proteases from alkaliphilic strains of Bacillus exhibiting resistance to oxidative inactivation were cloned and sequenced. The deduced amino acid sequences of the enzymes were highly homologous (greater than 88% identity). They were composed of 638 or 639 amino acids, including a possible approximately 200-amino acid prepro-peptide, and unique stretches of approximately 160 amino acids were found in the C-terminal regions. The molecular masses of mature enzymes (433 or 434 amino acids) were approximately 45 kDa for all. Amino acid sequence comparison and phylogenetic analysis indicated that these enzymes are far removed from other known subtilisins in the line of molecular evolution. We propose that these novel proteases be categorized as a new class of subtilisins, named oxidatively stable, alkaline protease.     

Molecular cloning, sequence and structural analysis of dehairing Mn(2+) dependent alkaline serine protease (MASPT) of Bacillus pumilus TMS55

Leather industries release a large amount of pollution-causing chemicals which creates one of the major industrial pollutions. The development of enzyme based processes as a potent alternative to pollution-causing chemicals is useful to overcome this issue. Proteases are enzymes which have extensive applications in leather processing and in several bioremediation processes due to their high alkaline protease activity and dehairing efficacy. In the present study, we report cloning, characterization of a Mn2+ dependent alkaline serine protease gene (MASPT) of Bacillus pumilus TMS55. The gene encoding the protease from B. pumilus TMS55 was cloned and its nucleotide sequence was determined. This gene has an open reading frame (ORF) of 1,149 bp that encodes a polypeptide of 383 amino acid residues. Our analysis showed that this polypeptide is composed of 29 residues N-terminal signal peptide, a propeptide of 79 residues and a mature protein of 275 amino acids. We performed bioinformatics analysis to compare MASPT enzyme with other proteases. Homology modeling was employed to model three dimensional structure for MASPT. Structural analysis showed that MASPT structure is composed of nine α-helices and nine β-strands. It has 3 catalytic residues and 14 metal binding residues. Docking analysis showed that residues S223, A260, N263, T328 and S329 interact with Mn2+. This study allows initial inferences about the structure of the protease and will allow the rational design of its derivatives for structure-function studies and also for further improvement of the enzyme.     

地衣芽孢杆菌2709碱性蛋白酶编码序列的PCR扩增、克隆及序列测定

在地衣芽孢杆菌NCIB 6816菌株碱性蛋白酶基因已知序列的基础上,通过设计合适的引物,利用PCR(Polymerase Chain Reaction)技术从地衣芽孢杆菌2709菌株的柒色体DNA中扩增了2709碱性蛋白酶的编码序列。对两个克隆的PCR片段的全序列分析结果显示,2709碱性蛋白酶的编码序列同相应的NCIB 6816序列相比有3％左右的碱基组成差异。由此推定的2709碱性蛋白酶的氨基酸序列肯定了2709碱性蛋白酶属典型的subtilisin Carlsberg类,同时还表明来源于不同地衣芽孢杆菌菌株的subtilisin Carlsberg存在着若干氨基酸组成上的差异。     

Molecular cloning, sequencing, and expression of a fibrinolytic serine-protease gene from the earthworm Lumbricus rubellus

The full-length cDNA of the lumbrokinase fraction 6 (F6) protease gene of Lumbricus rubellus was amplified using an mRNA template, sequenced and expressed in E. coli cells. The F6 protease gene consisted of pro- and mature sequences by gene sequence analysis, and the protease was translated and modified into active mature polypeptide by N-terminal amino acid sequence analysis of the F6 protease. The pro-region of F6 protease consisted of the 44 residues from methionine-1 to lysine-44, and the mature polypeptide sequence (239 amino acid residues and one stop codon; 720 bp) started from isoleucine-45 and continued to the terminal residue. F6 protease gene clones having pro-mature sequence and mature sequence produced inclusion bodies in E. coli cells. When inclusion bodies were orally administrated rats, generated thrombus weight in the rat's venous was reduced by approximately 60 % versus controls. When the inclusion bodies were solubilized in pepsin and/or trypsin solutions, the solubilized enzymes showed hemolytic activity in vitro. It was concluded the F6 protease has hemolytic activity, and that it is composed of pro- and mature regions.     

The primary structure of staphylococcal protease.

The amino acid sequence of staphylococcal protease has been determined by analysis of tryptic peptides obtained from cyanogen bromide fragments. Selected peptides obtained from digests with staphylococcal protease, thermolysin, and chymotrypsin provided the information necessary to align the tryptic peptides and the cyanogen bromide fragments. The protease is a single polypeptide chain of some 250 amino acids and is devoid of sulfhydryl groups. The COOH-terminal tryptic peptide of of the protease molecule contains some 43 residues, most of which are aspartic acids, asparagines, and prolines. The amino acid sequence of this peptide was not determined. The primary structure near the active serine residue indicates that staphylococcal protease is related to the pancreatic serine proteases. However, it has little or no additional sequence homologies with these enzymes except for the regions near histidine-50 and aspartic acid - 91. These regions have striking similarities with the corresponding regions of protease B and the trypsin-like enzyme of Streptomyces griseus.     

Consensus alignment for reliable framework prediction in homology modeling

MOTIVATION: Even the best sequence alignment methods frequently fail to correctly identify the framework regions for which backbones can be copied from the template into the target structure. Since the underprediction and, more significantly, the overprediction of these regions reduces the quality of the final model, it is of prime importance to attain as much as possible of the true structural alignment between target and template. RESULTS: We have developed an algorithm called Consensus that consistently provides a high quality alignment for comparative modeling. The method follows from a benchmark analysis of the 3D models generated by ten alignment techniques for a set of 79 homologous protein structure pairs. For 20-to-40% of the targets, these methods yield models with at least 6 A root mean square deviation (RMSD) from the native structure. We have selected the top five performing methods, and developed a consensus algorithm to generate an improved alignment. By building on the individual strength of each method, a set of criteria was implemented to remove the alignment segments that are likely to correspond to structurally dissimilar regions. The automated algorithm was validated on a different set of 48 protein pairs, resulting in 2.2 A average RMSD for the predicted models, and only four cases in which the RMSD exceeded 3 A. The average length of the alignments was about 75% of that found by standard structural superposition methods. The performance of Consensus was consistent from 2 to 32% target-template sequence identity, and hence it can be used for accurate prediction of framework regions in homology modeling.     

A cumulative specificity model for proteases from human immunodeficiency virus types 1 and 2, inferred from statistical analysis of an extended substrate data base.

Statistical analysis of an expanded data base of regions in viral polyproteins and in non-viral proteins that are sensitive to hydrolysis by the protease from human immunodeficiency virus (HIV) type 1 has generated a model which characterizes the substrate specificity of this retroviral enzyme. The model leads to an algorithm for predicting protease-susceptible sites from primary structure. Amino acids in each of the sites from P4 to P4' are tabulated for 40 protein substrates, and the frequency of occurrence for each residue is compared to the natural abundance of that amino acid in a selected data set of globular proteins. The results suggest that the highest stringency for particular amino acid residues is at the P2, P1, and P2' positions of the substrate. The broad specificity of the HIV-1 protease appears to be a consequence of its being able to bind productively substrates in which interactions with only a few Pi or Pi' side-chains need be optimized. The analysis, extended to 22 protein segments cleaved by the HIV-2 protease, delineates marked differences in specificity from that of the HIV-1 enzyme.     

地衣芽孢杆菌YP1A耐有机溶剂蛋白酶基因的克隆与功能表达

根据B．licheniformis YP1A来源的碱性蛋白酶具有的高强度耐有机溶剂性能及相关数据库分析,采用PCR克隆B．licheniformis YP1A耐有机溶剂碱性蛋白酶基因,序列分析显示该基因（1264bp）包含启动子与编码380个氨基酸的开放阅读框（ORF）,ORF包括信号肽、前肽及编码254个氨基酸的成熟肽序列。相关基因分析表明,YP1A耐有机溶剂碱性蛋白酶基因与地衣芽孢杆菌ATCC14580的碱性蛋白酶基因仅有6个氨基酸残基差异：构建2种含YP1A碱性蛋白酶CDS的组成型穿梭表达载体pHY／aprYP与pHY／aprP43,前者采用YP1A蛋白酶自带的启动子,后者则采用来自于质粒pP43NMK的P43强启动子。利用这2种表达载体在枯草芽孢杆菌WB800中成功进行蛋白酶的功能表达．其中P43强启动子的表达能力明显优于碱性蛋白酶自带的启动子,表达的蛋白酶比酶活为395U／ml。重组菌表达的碱性蛋白酶在体积分数50％的亲水及疏水有机溶剂中表现出了很好的耐受性,验证了克隆基因为地衣芽孢杆菌YP1A的高强度耐有机溶剂碱性蛋白酶基因.     

Three-dimensional structure of heat shock protein 90 from <Emphasis Type="Italic">Plasmodium falciparum</Emphasis>: molecular modelling approach to rational drug design against malaria

We have recently implicated heat shock protein 90 from Plasmodium falciparum (PfHsp90) as a potential drug target against malaria. Using inhibitors specific to the nucleotide binding domain of Hsp90, we have shown potent growth inhibitory effects on development of malarial parasite in human erythrocytes. To gain better understanding of the vital role played by PfHsp90 in parasite growth,we have modeled its three dimensional structure using recently described full length structure of yeast Hsp90.S equence similarity found between PfHsp90 and yeast Hsp90 allowed us to model the core structure with high confidence. The superimposition of the predicted structure with that of the template yeast Hsp90 structure reveals an RMSD of 3.31 Angstrom. The N-terminal and middle domains showed the least RMSD (1.76 Angstrom) while the more divergent C-terminus showed a greater RMSD (2.84 Angstrom) with respect to the template. The structure shows overall conservation of domains involved in nucleotide binding, ATPase activity, co-chaperone binding as well as inter-subunit interactions. Important co-chaperones known to modulate Hsp90 function in other eukaryotes are conserved in malarial parasite as well. An acidic stretch of amino acids found in the linker region, which is uniquely extended in PfHsp90 could not be modeled in this structure suggesting a flexible conformation. Our results provide a basis to compare the overall structure and functional pathways dependent on PfHsp90 in malarial parasite. Further analysis of differences found between human and parasite Hsp90 may make it possible to design inhibitors targeted specifically against malaria.     

Molecular cloning,nucleotide sequence and expression of the structural gene for a thermostable alkaline protease from Bacillus sp. no. AH-101

Summary Alkaliphilic Bacillus sp. no. AH-101 produces an extremely thermostable alkaline serine protease that has a high optimum pH (pH 12–13) and shows keratinolytic activity. The gene encoding this protease was cloned in Escherichia coli and expressed in B. subtilis. The cloned protease was identical to the AH-101 protease in its optimum pH and thermostability at high alkaline pH. An open reading frame of 1083 bases, identified as the protease gene, was preceded by a putative Shine-Dalgarno sequence (AAAGGAGG) with a spacing of 11 bases. The deduced amino acid sequence revealed a pre-pro-peptide of 93 residues followed by the mature protease comprising 268 residues. AH-101 protease showed slightly higher homology to alkaline proteases from alkaliphilic bacilli (61.2% and 65.3%) than to those from neutrophilic bacilli (54.9–56.7%). Also AH-101 protease and other proteases from alkaliphilic bacilli shared common amino acid changes and a four amino acid deletion when compared to the proteases from neutrophilic bacilli. AH-101 protease, however, was distinct among the proteases from alkaliphilic bacilli in showing the lowest homology to the others.Correspondence to: H. Takami     

Statistical coupling analysis of aspartic proteinases based on crystal structures of the Trichoderma reesei enzyme and its complex with pepstatin A

The crystal structures of an aspartic proteinase from Trichoderma reesei (TrAsP) and of its complex with a competitive inhibitor, pepstatin A, were solved and refined to crystallographic R-factors of 17.9% (R_free = 21.2%) at 1.70 Å resolution and 15.8% (R_free = 19.2%) at 1.85 Å resolution, respectively. The three-dimensional structure of TrAsP is similar to structures of other members of the pepsin-like family of aspartic proteinases. Each molecule is folded in a predominantly β-sheet bilobal structure with the N-terminal and C-terminal domains of about the same size. Structural comparison of the native structure and the TrAsP-pepstatin complex reveals that the enzyme undergoes an induced-fit, rigid-body movement upon inhibitor binding, with the N-terminal and C-terminal lobes tightly enclosing the inhibitor. Upon recognition and binding of pepstatin A, amino acid residues of the enzyme active site form a number of short hydrogen bonds to the inhibitor that may play an important role in the mechanism of catalysis and inhibition. The structures of TrAsP were used as a template for performing statistical coupling analysis of the aspartic protease family. This approach permitted, for the first time, the identification of a network of structurally linked residues putatively mediating conformational changes relevant to the function of this family of enzymes. Statistical coupling analysis reveals coevolved continuous clusters of amino acid residues that extend from the active site into the hydrophobic cores of each of the two domains and include amino acid residues from the flap regions, highlighting the importance of these parts of the protein for its enzymatic activity.     

Molecular cloning and homology modelling of a subtilisin-like serine protease from the marine fungus, <Emphasis Type="Italic">Engyodontium album</Emphasis> BTMFS10

An alkaline protease gene (Eap) was isolated for the first time from a marine fungus, Engyodontium album. Eap consists of an open reading frame of 1,161 bp encoding a prepropeptide consisting of 387 amino acids with a calculated molecular mass of 40.923 kDa. Homology comparison of the deduced amino acid sequence of Eap with other known proteins indicated that Eap encode an extracellular protease that belongs to the subtilase family of serine protease (Family S8). A comparative homology model of the Engyodontium album protease (EAP) was developed using the crystal structure of proteinase K. The model revealed that EAP has broad substrate specificity similar to Proteinase K with preference for bulky hydrophobic residues at P1 and P4. Also, EAP is suggested to have two disulfide bonds and more than two Ca²⁺ binding sites in its 3D structure; both of which are assumed to contribute to the thermostable nature of the protein.     

Rat submaxillary gland serine protease, tonin. Structure solution and refinement at 1.8 A resolution

Tonin is a mammalian serine protease that is capable of generating the vasoconstrictive agent, angiotensin II, directly from its precursor protein, angiotensinogen, a process that normally requires two enzymes, renin and angiotensin-converting enzyme. The X-ray crystallographic structure determination and refinement of tonin at 1.8 A resolution and the analysis of the resulting model are reported. The initial phases were obtained by the method of molecular replacement using as the search model the structure of bovine trypsin. The refined model of tonin consists of 227 amino acid residues out of the 235 in the complete molecule, 149 water molecules, and one zinc ion. The R-factor (R = sigma Fo - Fc/sigma Fo) is 0.196 for the 14,997 measured data between 8 and 1.8 A resolution with I greater than or equal to sigma (I). It is estimated that the overall root-mean-square error in the coordinates is about 0.3 A. The structure of tonin that has been determined is not in its active conformation, but one that has been perturbed by the binding of Zn2+ in the active site. Zn2+ was included in the buffer to aid the crystallization. Nevertheless, the structure of tonin that is described is for the most part similar to its native form as indicated by the close tertiary structural homology with kallikrein. The differences in the structures of the two enzymes are concentrated in several loop regions; these structural differences are probably responsible for the differences in their reactivities and specificities.     

Critical assessment of important regions in the subunit association and catalytic action of the severe acute respiratory syndrome coronavirus main protease

The severe acute respiratory syndrome (SARS) coronavirus (CoV) main protease represents an attractive target for the development of novel anti-SARS agents. The tertiary structure of the protease consists of two distinct folds. One is the N-terminal chymotrypsin-like fold that consists of two structural domains and constitutes the catalytic machinery; the other is the C-terminal helical domain, which has an unclear function and is not found in other RNA virus main proteases. To understand the functional roles of the two structural parts of the SARS-CoV main protease, we generated the full-length of this enzyme as well as several terminally truncated forms, different from each other only by the number of amino acid residues at the C- or N-terminal regions. The quaternary structure and K(d) value of the protease were analyzed by analytical ultracentrifugation. The results showed that the N-terminal 1-3 amino acid-truncated protease maintains 76% of enzyme activity and that the major form is a dimer, as in the wild type. However, the amino acids 1-4-truncated protease showed the major form to be a monomer and had little enzyme activity. As a result, the fourth amino acid seemed to have a powerful effect on the quaternary structure and activity of this protease. The last C-terminal helically truncated protease also exhibited a greater tendency to form monomer and showed little activity. We concluded that both the C- and the N-terminal regions influence the dimerization and enzyme activity of the SARS-CoV main protease.     

Comparative molecular model building of two serine proteinases from cytotoxic T lymphocytes

Two genes that are expressed when precursor cytotoxic T lymphocytes are transformed to T killer cells have been cloned and sequenced. The derived amino acid sequences, coding for cytotoxic cell protease 1 (CCP1) and Hannuka factor (HF) are highly homologous to members of the serine proteinase family. Comparative molecular model building using the known three-dimensional structures and the derived amino acid sequences of the lymphocyte enzymes has provided useful structural information, especially in predicting the conformations of the substrate binding sites. In applying this modelling procedure, we used the X-ray structures of four serine proteinases to provide a structurally based sequence alignment: alpha-chymotrypsin (CHT), bovine trypsin (BT), Streptomyces griseus trypsin (SGT), and rat mast cell protease 2 (RMCP2). The root mean square differences in alpha-carbon atom positions among these four structures when compared in a pairwise fashion range from 0.79 to 0.97 A for structurally equivalent residues. The sequences of the two lymphocyte enzymes were then aligned to these proteinases using chemical criteria and the superimposed X-ray structures as guides. The alignment showed that the sequence of CCP1 was most similar to RMCP2, whereas HF has regions of homology with both RMCP2 and BT. With RMCP2 as a template for CCP1 and the two enzymes RMCP2 and BT as templates for HF, the molecular models were constructed. Intramolecular steric clashes that resulted from the replacement of amino acid side chains of the templates by the aligned residues of CCP1 and HF were relieved by adjustment of the side chain conformational angles in an interactive computer graphics device. This process was followed by energy minimization of the enzyme model to optimize the stereochemical geometry and to relieve any remaining unacceptably close nonbonded contacts. The resulting model of CCP1 has an arginine residue at position 226 in the specificity pocket, thereby predicting a substrate preference for P1 aspartate or glutamate residues. The model also predicts favorable binding for a small hydrophobic residue at the P2 position of the substrate. The primary specificity pocket of HF resembles that of BT and therefore predicts a lysine or arginine preference for the P1 residue. The arginine at position 99 in the model of HF suggests a preference for aspartate or glutamate side chains in the P2 position of the substrate. Both CCP1 and HF have a free cysteine in the segment of polypeptide 88 to 93.(ABSTRACT TRUNCATED AT 400 WORDS)     

海洋细菌QD80低温碱性蛋白酶的基因克隆和性质

对海洋细菌QD80所产低温碱性蛋白酶进行了基因克隆和序列分析,对此酶的性质进行了初步研究.此酶基因开放阅读框架为1377bp,分子量为49.9kD.此序列上游-8bp处为该基因的SD序列,-10区和-35区分别有5′TAGAAT3′和5′TTGACC3′的保守序列.该酶最适pH为9.5,最适反应温度为30℃,在10℃酶活力仍能保持30%以上.该酶对氧化剂H2O2的抗氧化作用明显,浓度达到4gL时酶活仍保留85%.该蛋白酶的低温适应性和抗氧化特性将对其在低温洗涤领域的应用提供广泛的潜在应用价值.     

Isolation of subtilisin pro-sequence mutations that affect formation of active protease by localized random polymerase chain reaction mutagenesis

In order to analyze the role of the pro-sequence in folding of the alkaline serine protease subtilisin, localized random mutagenesis using the polymerase chain reaction with Taq DNA polymerase was employed to obtain mutations in the pro-sequence which prevent production of active protease. The unique aspect of this procedure is that random mutations can be easily generated in vitro over large but defined regions of a specific gene. The method was applied to a 458-base pair fragment encompassing the coding region of the pro-sequence of subtilisin, a region of the protein which has been shown to be required for proper folding. Protease-deficient mutants containing a variety of amino acid substitutions were isolated with a frequency of 4.3%. From analysis of these mutants, four independent amino acid substitution mutations in the pro-sequence were identified. The present results demonstrate that polymerase chain reaction is an efficient and simple method for obtaining random mutations within a localized region of a given gene.     

Sequence analysis of NS3 protease gene in clinical strains of hepatitis C virus

The amino terminal region of the non structural gene 3 (NS3) of hepatitis C virus (HCV) is a chymotripsinlike serine-protease responsible for cleavage of the non structural proteins of Hepatitis C virus (HCV). In order to investigate the genetic variation of this region, we developed a nested PCR to obtain NS3 protease sequences from 54 patients chronically infected with HCV genotypes 1a, 1b and 3, respectively. Comparison of nucleotide and amino acids sequences of NS3 protease domain with consensus sequence obtained within the same genotype, showed 3.73% nucleotide divergence and 1.64% amino acid divergence in isolates of genotype 3a, whereas isolates 1a exhibited 4.45% nucleotide and 4% amino acid change, respectively. Finally, NS3 sequence from 1b isolates revealed 6.47% nucleotide and 3.5 % aa changes. Comparison of consensus amino acid sequences derived from isolates 1a, 1b and 3, with the HCV prototypes showed a low amino acid sequence diversity. However, the consensus sequence of HCV genotype 3 isolates showed an amino acid changed from the prototype, that was located within a region important for enzyme structure and activity. These results indicated that the NS3 protease gene is highly conserved within the same HCV genotype. The domains involved in enzyme function were highly conserved in 1a and 1b strains, whereas consensus sequence of isolates 3a showed that the majority of these strains were not perfectly conserved in one of such regions. These findings altogether suggested that the NS3 protease enzyme of HCV may constitute an important target for antiviral therapy, but the NS3 protease variability of isolates 3 within a region that is a potential target for antiviral therapy could pose a problem for structure based drug development.