Leishmania species are the causative agents for Leishmaniasis which is one of the neglected tropical diseases causing 70,000 deaths worldwide each year. Squalene synthase enzyme plays a vital role in sterol metabolism which is essential for Leishmania parasite viability. Therefore squalene synthase of Leishmania donovani is a therapeutic target to inhibit growth of parasite. The 3D model of Leishmania donovani Squalene Synthase (LdSQS) was generated by homology modeling and validated through PROCHECK, ERRAT, VERIFY3D and PROSA tools. Virtual screening of the protein was performed by AutoDock with reported inhibitor, E5700 and two natural alkaloids. Molecular interactions were explored to understand the nature of intermolecular bonds between active ligand and the protein binding site residues using UCSF Chimera and PLIP server. The reported inhibitor showed the best binding affinity (-9.75 kcal/mol) closely followed by Ancistrotanzanine B (-9.55 kcal/mol) and Holamine (-8.79 kcal/mol). Ancistrotanzanine B showed low binding energy and permissible ADMET properties. Based on the present study, homology model of LdSQS and Ancistrotanzanine B can be used to design inhibitors with antileishmanial activity.

Chitinase is responsible for insect chitin hydrolyzation, which is a key process in insect molting and pupation. However, little is known about the chitinase of Spodoptera exigua (SeChi). In this study, based on the SeChi gene (ADI24346) identified in our laboratory, we constructed the recombinant baculovirus P-Chi for the expression of recombinant SeChi (rSeChi) in Hi5 cells. The rSeChi was purified by chelate affinity chromatography, and the purified protein showed activity comparable with that of a commercial SgChi, suggesting that we harvested active SeChi for the first time. The purified protein was subsequently tested for enzymatic properties and revealed to exhibit its highest activity at pH 8 and 40 C. Using homology modeling and molecular docking techniques, the three-dimensional model of SeChi was constructed and screened for inhibitors. In two rounds of screening, twenty compounds were selected. With the purified rSeChi, we tested each of the twenty compounds for inhibitor activity against rSeChi, and seven compounds showed obvious activity. This study provided new information for the chitinase of beet armyworm and for chitinase inhibitor development.

Unavailability of vaccines and effective drugs are primarily responsible for the growing menace of cryptosporidiosis. This study has incorporated a bioinformatics-based screening approach to explore potential vaccine candidates and novel drug targets in Cryptosporidium parvum proteome. A systematic strategy was defined for comparative genomics, orthology with related Cryptosporidium species, prioritization parameters and MHC class I and II binding promiscuity. The approach reported cytoplasmic protein cgd7_1830, a signal peptide protein, as a novel drug target. SWISS-MODEL online server was used to generate the 3D model of the protein and was validated by PROCHECK. The model has been subjected to in silico docking study with screened potent lead compounds from the ZINC database, PubChem and ChEMBL database using Flare software package of Cresset®. Furthermore, the approach reported protein cgd3_1400, as a vaccine candidate. The predicted B- and T-cell epitopes on the proposed vaccine candidate with highest scores were also subjected to docking study with MHC class I and II alleles using ClusPro web server. Results from this study could facilitate selection of proteins which could serve as drug targets and vaccine candidates to efficiently tackle the growing threat of cryptosporidiosis.

The trematode Fasciola hepatica secretes a number of cathepsin L-like proteases that are proposed to be involved in feeding, migration, and immune evasion by the parasite. To date, six full cDNA sequences encoding cathepsin L preproproteins have been identified. Previous studies have demonstrated that one of these cathepsins (L2) is unusual in that it is able to cleave substrates with a proline in the P2 position, translating into an unusual ability (for a cysteine proteinase) to clot fibrinogen. In this study, we report the sequence of a novel cathepsin (L5) and compare the substrate specificity of a recombinant enzyme with that of recombinant cathepsin L2. Despite sharing 80% sequence identity with cathepsin L2, cathepsin L5 does not exhibit substantial catalytic activity against substrates containing proline in the P2 position. Molecular modeling studies suggested that a single amino acid change (L69Y) in the mature proteinases may account for the difference in specificity at the S2 subsite. Recombinant cathepsin L5/L69Y was expressed in yeast and a substantial increase in the ability of this variant to accommodate substrates with a proline residue in the P2 position was observed. Thus, we have identified a single amino acid substitution that can substantially influence the architecture of the S2 subsite of F. hepatica cathepsin L proteases.

Metabotropic glutamate receptors (mGluRs) belong to the family 3 of G-protein-coupled receptors. On these proteins, agonist binding on the extracellular domain leads to conformational changes in the 7-transmembrane domains required for G-protein activation. To elucidate the structural features that might be responsible for such an activation mechanism, we have generated models of the amino terminal domain (ATD) of type 4 mGluR (mGlu4R). The fold recognition search allowed the identification of three hits with a low sequence identity, but with high secondary structure conservation: leucine isoleucine valine-binding protein (LIVBP) and leucine-binding protein (LBP) as already known, and acetamide-binding protein (AmiC). These proteins are characterized by a bilobate structure in an open state for LIVBP/LBP and a closed state for AmiC, with ligand binding in the cleft. Models for both open and closed forms of mGlu4R ATD have been generated. ACPT-I (1-aminocyclopentane 1,3,4-tricarboxylic acid), a selective agonist, has been docked in the two models. In the open form, ACPT-I is only bound to lobe I through interactions with Lys74, Arg78, Ser159, and Thr182. In the closed form, ACPT-I is trapped between both lobes with additional binding to Tyr230, Asp312, Ser313, and Lys317 from lobe II. These results support the hypothesis that mGluR agonists bind a closed form of the ATDs, suggesting that such a conformation of the binding domain corresponds to the active conformation.

We have modeled the structure of human lymphotactin (hLpnt), by homology modeling and molecular dynamics simulations. This chemokine is unique in having a single disulfide bond and a long C-terminal tail. Because other structural classes of chemokines have two pairs of Cys residues, compared to one in Lpnt, and because it has been shown that both disulfide bonds are required for stability and function, the question arises how the Lpnt maintains its structural integrity. The initial structure of hLpnt was constructed by homology modeling. The first 63 residues in the monomer of hLpnt were modeled using the structure of the human CC chemokine, RANTES, whose sequence appeared most similar. The structure of the long C-terminal tail, missing in RANTES, was taken from the human muscle fatty-acid binding protein. In a Protein Data Bank search, this protein was found to contain a sequence that was most homologous to the long tail. Consequently, the modeled hLpnt C-terminal tail consisted of both α-helical and β-motifs. The complete model of the hLpnt monomer consisted of two α-helices located above the five-stranded β-sheet. Molecular dynamics simulations of the solvated initial model have indicated that the stability of the predicted fold is related to the geometry of Pro78. The five-stranded β-sheet appeared to be preserved only when Pro78 was modeled in the cis conformation. Simulations were also performed both for the C-terminal truncated forms of the hLpnt that contained one or two (CC chemokine-like) disulfide bonds, and for the chicken Lpnt (cLpnt). Our MD simulations indicated that the turn region (T30–G34) in hLpnt is important for the interactions with the receptor, and that the long C-terminal region stabilizes both the turn (T30–G34) and the five-stranded β-sheet. The major conclusion from our theoretical studies is that the lack of one disulfide bond and the extension of the C-terminus in hLptn are mutually complementary. It is very likely that removal of two Cys residues sufficiently destabilizes the structure of a chemokine molecule, particularly the core β-sheet, to abolish its biological function. However, this situation is rectified by the long C-terminal segment. The role of this long region is most likely to stabilize the first β-turn region and α-helix H1, explaining how this chemokine can function with a single disulfide bond.

In Saccharomyces cerevisiae, Prp17p is required for the efficient completion of the second step of pre-mRNA splicing. The function and interacting factors for this protein have not been elucidated. We have performed a mutational analysis of yPrp17p to identify protein domains critical for function. A series of deletions were made throughout the region spanning the N-terminal 158 amino acids of the protein, which do not contain any identified structural motifs. The C-terminal portion (amino acids 160–455) contains a WD domain containing seven WD repeats. We determined that a minimal functional Prp17p consists of the WD domain and 40 amino acids N-terminal to it. We generated a three-dimensional model of the WD repeats in Prp17p based on the crystal structure of the β-transducin WD domain. This model was used to identify potentially important amino acids for in vivo functional characterization. Through analysis of mutations in four different loops of Prp17p that lie between β strands in the WD repeats, we have identified four amino acids, 235TETG238, that are critical for function. These amino acids are predicted to be surface exposed and may be involved in interactions that are important for splicing. Temperature-sensitive prp17 alleles with mutations of these four amino acids are defective for the second step of splicing and are synthetically lethal with a U5 snRNA loop I mutation, which is also required for the second step of splicing. These data reinforce the functional significance of this region within the WD domain of Prp17p in the second step of splicing.

Several recent publications illustrated advantages of using sequence profiles in recognizing distant homologies between proteins. At the same time, the practical usefulness of distant homology recognition depends not only on the sensitivity of the algorithm, but also on the quality of the alignment between a prediction target and the template from the database of known proteins. Here, we study this question for several supersensitive protein algorithms that were previously compared in their recognition sensitivity (Rychlewski et al., 2000). A database of protein pairs with similar structures, but low sequence similarity is used to rate the alignments obtained with several different methods, which included sequence–sequence, sequence–profile, and profile–profile alignment methods. We show that incorporation of evolutionary information encoded in sequence profiles into alignment calculation methods significantly increases the alignment accuracy, bringing them closer to the alignments obtained from structure comparison.

In general, alignment quality is correlated with recognition and alignment score significance. For every alignment method, alignments with statistically significant scores correlate with both correct structural templates and good quality alignments. At the same time, average alignment lengths differ in various methods, making the comparison between them difficult. For instance, the alignments obtained by FFAS, the profile–profile alignment algorithm developed in our group are always longer that the alignments obtained with the PSI-BLAST algorithms. To address this problem, we develop methods to truncate or extend alignments to cover a specified percentage of protein lengths. In most cases, the elongation of the alignment by profile–profile methods is reasonable, adding fragments of similar structure. The examples of erroneous alignment are examined and it is shown that they can be identified based on the model quality.

Previously, we determined the DNA and amino acid sequences as well as biochemical and biophysical properties of a series of fungal phytases. The amino acid sequences displayed 49–68% identity between species, and the catalytic properties differed widely in terms of specific activity, substrate specificity, and pH optima. With the ultimate goal to combine the most favorable properties of all phytases in a single protein, we attempted, in the present investigation, to increase the specific activity of Aspergillus fumigatus phytase. The crystal structure of Aspergillus niger NRRL 3135 phytase known at 2.5 Å resolution served to specify all active site residues. A multiple amino acid sequence alignment was then used to identify nonconserved active site residues that might correlate with a given favorable property of interest. Using this approach, Gln27 of A. fumigatus phytase (amino acid numbering according to A. niger phytase) was identified as likely to be involved in substrate binding and/or release and, possibly, to be responsible for the considerably lower specific activity (26.5 vs. 196 U·[mg protein]−1 at pH 5.0) of A. fumigatus phytase when compared to Aspergillus terreus phytase, which has a Leu at the equivalent position. Site-directed mutagenesis of Gln27 of A. fumigatus phytase to Leu in fact increased the specific activity to 92.1 U·(mg protein)−1, and this and other mutations at position 27 yielded an interesting array of pH activity profiles and substrate specificities. Analysis of computer models of enzyme–substrate complexes suggested that Gln27 of wild-type A. fumigatus phytase forms a hydrogen bond with the 6-phosphate group of myo-inositol hexakisphosphate, which is weakened or lost with the amino acid substitutions tested. If this hydrogen bond were indeed responsible for the differences in specific activity, this would suggest product release as the rate-limiting step of the A. fumigatus wild-type phytase reaction.

Using a combination of theoretical sequence structure recognition predictions and experimental disulfide bond assignments, a three-dimensional (3D) model of human interleukin-7 (hIL-7) was constructed that predicts atypical surface chemistry in helix D that is important for receptor activation. A 3D model of hIL-7 was built using the X-ray crystal structure of interleukin-4 (IL-4) as a template (Walter MR et al., 1992, J Mol Biol. 224:1075–1085; Walter MR et al., 1992, J Biol Chem 267:20371–20376). Core secondary structures were constructed from sequences of hIL-7 predicted to form helices. The model was constructed by superimposing IL-7 helices onto the IL-4 template and connecting them together in an up–up down–down topology. The model was finished by incorporating the disulfide bond assignments (Cys3, Cys142), (Cys35, Cys130), and (Cys48, Cys93), which were determined by MALDI mass spectroscopy and site-directed mutagenesis (Cosenza L, Sweeney E, Murphy JR, 1997, J Biol Chem 272:32995–33000). Quality analysis of the hIL-7 model identified poor structural features in the carboxyl terminus that, when further studied using hydrophobic moment analysis, detected an atypical structural property in helix D, which contains Cys130 and Cys142. This analysis demonstrated that helix D had a hydrophobic surface exposed to bulk solvent that accounted for the poor quality of the model, but was suggestive of a region in IL-7 that maybe important for protein interactions. Alanine (Ala) substitution scanning mutagenesis was performed to test if the predicted atypical surface chemistry of helix D in the hIL-7 model is important for receptor activation. This analysis resulted in the construction, purification, and characterization of four hIL-7 variants, hIL-7(K121A), hIL-7(L136A), hIL-7(K140A), and hIL-7(W143A), that displayed reduced or abrogated ability to stimulate a murine IL-7 dependent pre-B cell proliferation. The mutant hIL-7(W143A), which is biologically inactive and displaces [125I]-hIL-7, is the first reported IL-7R system antagonist.

Protein phosphatase 1 (PP1) is widely distributed among tissues and species and acts as a regulator of many important cellular processes. By targeting the catalytic part of PP1 (PP1C) toward particular loci and substrates, regulatory subunits constitute key elements conferring specificity to the holoenzyme. Here, we report the identification of an (α/β)8-barrel-like structure within the N-ter stretch of the human PP1 regulatory subunit hGM, which is part of the family of diverse proteins associated with glycogen metabolism. Protein homology modeling gave rise to a three-dimensional (3D) model for the 381 N-ter residue stretch of hGM, based on sequence similarity with Streptomyces olivochromogenes xylose isomerase, identified by using FASTA. The alignment was subsequently extended by using hydrophobic cluster analysis. The homology-derived model includes the putative glycogen binding area located within the 142–230 domain of hGM as well as a structural characterization of the PP1C interacting domain (segment 51–67). Refinement of the latter by molecular dynamics afforded a topology that is in agreement with previous X-ray studies (Egloff et al., 1997). Finite difference Poisson–Boltzmann calculations performed on the interacting domains of PP1C and hGM confirm the complementarity of the local electrostatic potentials of the two partners. This work highlights the presence of a conserved fold among distant species (mammalian, Caenorhabditis elegans, yeast) and, thus, emphasizes the involvement of PP1 in crucial basic cellular functions.

The moderate thermophilic eubacterium Alicyclobacillus (formerly Bacillus) acidocaldarius expresses a thermostable carboxylesterase (esterase 2) belonging to the hormone-sensitive lipase (HSL)-like group of the esterase/lipase family. Based on secondary structures predictions and a secondary structure-driven multiple sequence alignment with remote homologous protein of known three-dimensional (3D) structure, we previously hypothesized for this enzyme the α/β-hydrolase fold typical of several lipases and esterases and identified Ser155, Asp252, and His282 as the putative members of the catalytic triad. In this paper we report the construction of a 3D model for this enzyme based on the structure of mouse acetylcholinesterase complexed with fasciculin. The model reveals the topological organization of the fold corroborating our predictions. As regarding the active-site residues, Ser155, Asp252, and His282 are located close to each other at hydrogen bond distances. Their catalytic role was here probed by biochemical and mutagenic studies. Moreover, on the basis of the secondary structure-driven multiple sequence alignment and the 3D structural model, a residue supposed important for catalysis, Gly84, was mutated to Ser. The activity of the mutated enzyme was drastically reduced. We propose that Gly84 is part of a putative “oxyanion hole” involved in the stabilization of the transition state similar to the C group of the esterase/lipase family.

The thioredoxin action upon the 2-oxoacid dehydrogenase complexes is investigated by using different thioredoxins, both wild-type and mutated. The attacking cysteine residue of thioredoxin is established to be essential for the thioredoxin-dependent activation of the complexes. Mutation of the buried cysteine residue to serine is not crucial for the activation, but prevents inhibition of the complexes, exhibited by the Clamydomonas reinhardtii thioredoxin m disulfide. Site-directed mutagenesis of D26, W31, F/W12, and Y/A70 (the Escherichia coli thioredoxin numbering is employed for all the thioredoxins studied) indicates that both the active site and remote residues of thioredoxin are involved in its interplay with the 2-oxoacid dehydrogenase complexes. Sequences of 11 thioredoxin species tested biochemically are aligned. The thioredoxin residues at the contact between the α3/310 and α1 helices, the length of the α1 helix and the charges in the α2-β3 and β4-β5 linkers are found to correlate with the protein influence on the 2-oxoacid dehydrogenase complexes (the secondary structural elements of thioredoxin are defined according to Eklund H et al., 1991, Proteins 11:13–28). The distribution of the charges on the surface of the thioredoxin molecules is analyzed. The analysis reveals the species specific polarization of the thioredoxin active site surroundings, which corresponds to the efficiency of the thioredoxin interplay with the 2-oxoacid dehydrogenase systems. The most effective mitochondrial thioredoxin is characterized by the strongest polarization of this area and the highest value of the electrostatic dipole vector of the molecule. Not only the magnitude, but also the orientation of the dipole vector show correlation with the thioredoxin action. The dipole direction is found to be significantly influenced by the charges of the residues 13/14, 51, and 83/85, which distinguish the activating and inhibiting thioredoxin disulfides.