Homology Modeling of Coagulase in Staphylococcus aureus

The close correlation between the ability of coagulase to clot blood plasma and their capacity to produce disease, and the corresponding absence of this property in nonpathogenic strains, have led to the assumption that the coagulase, plays important role in the pathogenesis of disease. Currently, crystal structure of coagulase in Staphylococcus aureus remains indefinable. Thus, the objectives of this research is to generate the three dimensional model of coagulase in S. aureus by using homology modeling approach. In this study, we used bioinformatics tools and databases such as BLAST (Basic Local Alignment Search Tool), GenBank, PDB (Protein Databank), and Discovery Studio to gain specific functional insights into coagulase. The model was validated using protein structure checking tools such as PROCHECK, Verify 3D and CE (Combinatorial Extension) for reliability. Therefore, structure prediction of coagulase in S. aureus can provide preliminary knowledge for understanding the function of the protein. The information from this finding will provide important information into the action and regulation mechanism of the coagulase protein in S. aureus.     

Structure and functional annotation of hypothetical proteins having putative Rubisco activase function from Vitis vinifera

Rubisco is a very large, complex and one of the most abundant proteins in the world and comprises up to 50% of all soluble protein in plants. The activity of Rubisco, the enzyme that catalyzes CO₂ assimilation in photosynthesis, is regulated by Rubisco activase (Rca). In the present study, we searched for hypothetical protein of Vitis vinifera which has putative Rubisco activase function. The Arabidopsis and tobacco Rubisco activase protein sequences were used as seed sequences to search against Vitis vinifera in UniprotKB database. The selected hypothetical proteins of Vitis vinifera were subjected to sequence, structural and functional annotation. Subcellular localization predictions suggested it to be cytoplasmic protein. Homology modelling was used to define the three-dimensional (3D) structure of selected hypothetical proteins of Vitis vinifera. Template search revealed that all the hypothetical proteins share more than 80% sequence identity with structure of green-type Rubisco activase from tobacco, indicating proteins are evolutionary conserved. The homology modelling was generated using SWISS-MODEL. Several quality assessment and validation parameters computed indicated that homology models are reliable. Further, functional annotation through PFAM, CATH, SUPERFAMILY, CDART suggested that selected hypothetical proteins of Vitis vinifera contain ATPase family associated with various cellular activities (AAA) and belong to the AAA+ super family of ring-shaped P-loop containing nucleoside triphosphate hydrolases. This study will lead to research in the optimization of the functionality of Rubisco which has large implication in the improvement of plant productivity and resource use efficiency.     

Homology modeling, comparative genomics and functional annotation of Mycoplasma genitalium hypothetical protein MG_237

Mycoplasma genitalium is a human pathogen associated with several sexually transmitted diseases. The complete genome of M. genitalium G37 has been sequenced and provides an opportunity to understand the pathogenesis and identification of therapeutic targets. However, complete understanding of bacterial function requires proper annotation of its proteins. The genome of M. genitalium consists of 475 proteins. Among these, 94 are without any known function and are described as 'hypothetical proteins'. We selected MG_237 for sequence and structural analysis using a bioinformatics approach. Primary and secondary structure analysis suggested that MG_237 is a hydrophilic protein containing a significant proportion of alpha helices, and subcellular localization predictions suggested it is a cytoplasmic protein. Homology modeling was used to define the three-dimensional (3D) structure of MG-237. A search for templates revealed that MG_237 shares 63% homology to a hypothetical protein of Mycoplasma pneumoniae, indicating this protein is evolutionary conserved. The refined 3D model was generated using (PS)(2)-v2 sever that incorporates MODELLER. Several quality assessment and validation parameters were computed and indicated that the homology model is reliable. Furthermore, comparative genomics analysis suggested MG_237 as non-homologous protein and involved in four different metabolic pathways. Experimental validation will provide more insight into the actual function of this protein in microbial pathways.     

Functional insights by comparison of modeled structures of 18kDa small heat shock protein and its mutant in Mycobacterium leprae

In this work we are proposing Homology modeled structures of Mycobacterium leprae 18kDa heat shock protein and its mutant. The more closely related structure of the small heat shock protein (sHSP) belonging to the eukaryotic species from wheat sHSP16.9 and 16.3kDa ACR1 protein from Mycobacterium tuberculosis were used as template structures. Each model contains an N-terminal domain, alpha-crystalline domain and a C-terminal tail. The models showed that a single point mutation from serine to proline at 52^nd position causes structural changes. The structural changes are observed in N-terminal region and alpha-crystalline domains. Serine in 52^nd position is observed in β4 strand and Proline in 52^nd position is observed in loop. The number of residues contributing α helix at N-terminal region varies in both models. In 18S more number of residues is present in α helix when compared to 18P. The loop regions between β3 and β4 strands of both models vary in number of residues present in it. Number of residues contributing β4 strand in both models vary. β6 strand is absent in both models. Major functional peptide region of alpha crystalline domains of both models varies. These differences observed in secondary structures support their distinct functional roles. It also emphasizes that a point mutation can cause structural variation.     

Structure modeling and antidiabetic activity of a seed protein of Momordica charantia in non-obese diabetic (NOD) mice

Momordica charantia is a well known medicinal plant used in the traditional medicinal system for the treatment of various diseases including diabetes mellitus. Recently, a novel protein termed as ADMc1 from the seed extract of M. charantia has been identified and isolated showing significant antihyperglycemic activity in type 1 diabetic rats in which diabetes was induced. However, the structure of this protein has not yet been analyzed. Homology modeling approach was used to generate a high quality protein 3D structure for the amino acid sequence of the ADMc1 protein in this study. The comparative assessment of secondary structures revealed ADMc1 as an all-alpha helix protein with random coils. Tertiary structure predicted on the template structure of Napin of B. Napus (PDB ID: 1SM7) with which the ADMc1 showed significant sequence similarity, was validated using protein structure validation tools like PROCHECK, WHAT_CHECK, VERIFY3D and ProSA. Arrangement of disulfide bridges formed by cysteine residues were predicted by the Dianna 1.1 server. The presence of multiple disulfide bond confers the stable nature of the ADMc1 protein. Further, the biological activity of the ADMc1 was assessed in non-obese diabetic (NOD) mice which are spontaneous model of type 1 diabetes. Significant reduction in the blood glucose levels of NOD mice was observed up to 8 h post administration of the rADMc1 protein. Overall, the structural characterizations with antihyperglycemic activity of this seed protein of Momordica charantia demonstrate its potential as an antidiabetic agent.     

Structure prediction and evolution of a halo-acid dehalogenase of Burkholderia mallei

Environmental pollutants containing halogenated organic compounds e.g. haloacid, can cause a plethora of health problems. The structural and functional analyses of the gene responsible of their degradation are an important aspect for environmental studies and are important to human well-being. It has been shown that some haloacids are toxic and mutagenic. Microorganisms capable of degrading these haloacids can be found in the natural environment. One of these, a soil-borne Burkholderia mallei posses the ability to grow on monobromoacetate (MBA). This bacterium produces a haloacid dehalogenase that allows the cell to grow on MBA, a highly toxic and mutagenic environmental pollutant. For the structural and functional analysis, a 346 amino acid encoding protein sequence of haloacid dehalogenase is retrieve from NCBI data base. Primary and secondary structure analysis suggested that the high percentage of helices in the structure makes the protein more flexible for folding, which might increase protein interactions. The consensus protein sub-cellular localization predictions suggest that dehalogenase protein is a periplasmic protein 3D2GO server, suggesting that it is mainly employed in metabolic process followed by hydrolase activity and catalytic activity. The tertiary structure of protein was predicted by homology modeling. The result suggests that the protein is an unstable protein which is also an important characteristic of active enzyme enabling them to bind various cofactors and substrate for proper functioning. Validation of 3D structure was done using Ramachandran plot ProsA-web and RMSD score. This predicted information will help in better understanding of mechanism underlying haloacid dehalogenase encoding protein and its evolutionary relationship.     

In silicio expression analysis of PKS genes isolated from Cannabis sativa L

Cannabinoids, flavonoids, and stilbenoids have been identified in the annual dioecious plant Cannabis sativa L. Of these, the cannabinoids are the best known group of this plant's natural products. Polyketide synthases (PKSs) are responsible for the biosynthesis of diverse secondary metabolites, including flavonoids and stilbenoids. Biosynthetically, the cannabinoids are polyketide substituted with terpenoid moiety. Using an RT-PCR homology search, PKS cDNAs were isolated from cannabis plants. The deduced amino acid sequences showed 51%-73% identity to other CHS/STS type sequences of the PKS family. Further, phylogenetic analysis revealed that these PKS cDNAs grouped with other non-chalcone-producing PKSs. Homology modeling analysis of these cannabis PKSs predicts a 3D overall fold, similar to alfalfa CHS2, with small steric differences on the residues that shape the active site of the cannabis PKSs.     

Annotation of hypothetical proteins orthologous in Pongo abelii and Sus scrofa

A hypothetical protein is predicted to be expressed from an open reading frame without known experimental evidence of translation. They constitute a substantial fraction of proteomes. Domain extraction from these hypothetical sequences helps to search for protein coding genes for protein structural and functional annotation. We describe the analysis of prediction data in a sequence dataset of hypothetical protein orthologs of Pongo abelii (orangutan) and Sus scrofa (pig). It should be noted that these orangutan-pig orthologs are also non-homologous to human proteins. These predicted data find application in the genome wide annotation of proteins in poorly understood genomes.

Abbreviations

PDB - Protein Data Bank, DEG - Database of Essential Genes, CDD - Conserved Domain Database, IUCN - International Union for Conservation of Nature.     

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.     

Tryparedoxin peroxidase of Leishmania braziliensis: homology modeling and inhibitory effects of flavonoids for anti-leishmanial activity

Inhibition of the Tryparedoxin peroxidase interaction has been becomes a new therapeutic strategy in leishmaniasis. Docking analysis was carried out to study the effects of quercetin and taxifolin on Tryparedoxin Peroxidase (TryP). Tryparedoxin peroxidase of Trypanosomatidae functions as antioxidants through their Peroxidase and peroxynitrite reductase activities. The 3D models of Tryparedoxin Peroxidase of Leishmania braziliensis (L. braziliensis TryP) was modeled using the template Tryparedoxin Peroxidase I from Leishmania Major (L. Major TryPI) (PDB ID: 3TUE). Further, we evaluated for TryP inhibitory activity of flavonoids such as quercetin and taxifolin using in silico docking studies. Docking results showed the binding energies of - 11.8601and -8.0851 for that quercetin and taxifolin respectively. Flavonoids contributed better L. braziliensis TryP inhibitory activity because of its structural parameters. Thus, from our in silico studies we identify that quercetin and taxifolin posses anti-leishmanial acitivities mediated through TryP inhibition mechanism.     

Homology modeling of Mycoplasma pneumoniae enolase and its molecular interaction with human plasminogen

Alpha (α)-enolase (e), a glycolytic enzyme, has an alternative role as a surface receptor of several bacteria mediating plasminogen (pg) binding. It is also recognized as a virulence factor of some pathogenic bacteria facilitating plasminogen activation and host cell invasion. A mycoplasmal α-enolase is also a plasminogen binding protein. Molecular interactions of enolase from Mycoplasma pneumoniae with host plasminogen would be useful for exploring the pathogen-host interaction. In an attempt to identify plasminogen binding sites of M. pneumoniae enolase, homology modeling and docking studies were conducted to obtain modeled structures of the M. pneumoniae enolase-plasminogen complex. The refined model was validated further by standard methods. Molecular docking revealed hydrogen bonding of eLys70-pgTyr50, eAsn165-pgThr66, eAla168-pgGlu21, eAsp17-pgLys70, and eAsn213-pgPro68/pgAsn69. Substantial decreases in accessible surface area (ASA) were observed and in concurrence with hydrogen bond pattern. These findings provide a detailed prediction of key residues that interact at the protein-protein interface. Our theoretical prediction is consistent with known biochemical data. The predicted interaction complex can be of great assistance in understanding structural insights, which is necessary to pathogen and host-component interaction. The ability of M. pneumoniae enolase to bind plasminogen may be indicative of an important role in invasion of this pathogen to host.     

Three dimensional (3D) structure prediction and substrate-protein interaction study of the chitin binding protein CBP24 from B. thuringiensis

Bacillus thuringiensis is an insecticidal bacterium whose chitinolytic system has been exploited to improve insect resistance in crops. In the present study, we studied the CBP24 from B. thuringiensis using homology modeling and molecular docking. The primary and secondary structure analyses showed CBP24 is a positively charged protein and contains single domain that belongs to family CBM33. The 3D model after refinement was used to explore the chitin binding characteristics of CBP24 using AUTODOCK. The docking analyses have shown that the surface exposed hydrophilic amino acid residues Thr-103, Lys-112 and Ser-162 interact with substrate through H-bonding. While, the amino acids resides Glu-39, Tyr-46, Ser-104 and Asn-109 were shown to have polar interactions with the substrate. The binding energy values evaluation of docking depicts a stable intermolecular conformation of the docked complex. The functional characterization of the CBP24 will elucidate the substrate-interaction pathway of the protein in specific and the carbohydrate binding proteins in general leading towards the exploration and exploitation of the prokaryotic substrate utilization pathways.     

Insights from the analysis of a predicted model of gp63 in Leishmania donovani

Leishmaniasis is a protozoal disease of human that occurs in most parts of the world. By considering the progress of bioinformatics in molecular modeling, major surface glycoprotein of Leishmania donovani (gp63) structure was modeled using homology modeling with high accuracy based on the X-ray crystal structure of the Leishmania major gp63 as a template, and then analyzed 3D structure of gp63 which can reveal exact facts about its structure and interaction. The objective of this study was to find folding and three dimensional structure of the gp63 as potent antigen for human. In this project, we applied the theory of evolution method, including comparative modeling and threading. This study presented a simple protocol for rapid and precise finding 3D structure of gp63 and investigation of its structural properties. The translated amino acid sequence showed that Leishmania donovani gp63 contains 590 amino acids precursor protein consisting of an NH₂-terminal signal peptide of 39 amino acids for membrane targeting, a pro region of 48 amino acids, the mature protein of 478 amino acids containing glycosylation and putative catalytic sites, and a COOH-terminal signal peptide of 25 amino acids for GPI attachment. Based on our model, the protein consists of three domains: the N-terminal, central and C-terminal domains. Additionally, these results could guide future structure-function analyses of gp63 protein.     

Modeling of human M1 aminopeptidases for in silico screening of potential Plasmodium falciparum alanine aminopeptidase (PfA-M1) specific inhibitors

Plasmodium falciparum alanine M1-aminopeptidase (PfA-M1) is a validated target for anti-malarial drug development. Presence ofsignificant similarity between PfA-M1 and human M1-aminopeptidases, particularly within regions of enzyme active site leads toproblem of non-specificity and off-target binding for known aminopeptidase inhibitors. Molecular docking based in silico screeningapproach for off-target binding has high potential but requires 3D-structure of all human M1-aminopeptidaes. Therefore, in thepresent study 3D structural models of seven human M1-aminopeptidases were developed. The robustness of docking parametersand quality of predicted human M1-aminopeptidases structural models was evaluated by stereochemical analysis and docking oftheir respective known inhibitors. The docking scores were in agreement with the inhibitory concentrations elucidated in enzymeassays of respective inhibitor enzyme combinations (r2≈0.70). Further docking analysis of fifteen potential PfA-M1 inhibitors(virtual screening identified) showed that three compounds had less docking affinity for human M1-aminopeptidases as comparedto PfA-M1. These three identified potential lead compounds can be validated with enzyme assays and used as a scaffold fordesigning of new compounds with increased specificity towards PfA-M1.