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.     

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 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.     

Homology modeling of Ferredoxin-nitrite reductase from Arabidopsis thaliana

Different subtypes of Influenza A virus are associated with species specific, zoonotic or pandemic Influenza. The cause of its severity underlies in complicated evolution of its segmented RNA genome. Although genetic shift and genetic drift are well known in the evolution of this virus, we reported the significant role of unique RNA palindromes in its evolution. Our computational approach identified the existence of unique palindromes in each subtype of Influenza A virus with its absence in Influenza B relating the fact of virulence and vigorous genetic hitchhiking in Influenza A. The current study focused on the re-assortment event responsible for the emergence of pandemic-2009 H1N1 virus, which is associated with outgrow of new palindrome and in turn, changing its RNA structure. We hypothesize that the change in RNA structure due to the presence of palindrome facilitates the event of re-assortment in Influenza A. Thus the evolutionary process of Influenza A is much more complicated as previously known, and that has been demonstrated in this study.     

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.     

Homology modeling of the Abl-SH3 domain

A tertiary structure model of the Abl-SH3 domain is predicted by using homology modeling techniques coupled to molecular dynamics simulations. Two template proteins were used, Fyn-SH3 and Spc-SH3. The refined model was extensively checked for errors using criteria based on stereochemistry, packing, solvation free-energy, accessible surface areas, and contact analyses. The different checking methods do not totally agree, as each one evaluates a different characteristic of protein structures. Several zones of the protein are more susceptible to incorporating errors. These include residues 13, 15, 35, 39, 45, 46, 50, and 60. An interesting finding is that the measurement of the Cα chirality correlated well with the rest of the criteria, suggesting that this parameter might be a good indicator of correct local conformation. Deviations of more than 4 degrees may be indicative of poor local structure. © 1994 Wiley-Liss, Inc.     

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.     

Homology modeling and consensus protein disorder prediction of human filamin

Filamins are dimeric actin-binding proteins participating in the organization of the actin-based cytoskeleton. Their modular domain organization is made up of an N-terminal actin-binding domain composed of two CH domains followed by flexible rod regions that consist of 24 Ig-like domains. Homology modeling was used to model human filamin using Modeller 9v5. The resulting model assessed by Verify 3D and PROCHECK showed that the final model is reliable. The conformational disorder prediction of human filamin residues were also mapped on the validated structure of human filamin. Prediction of protein disorder in filamin structures will help structural biologists to find suitable targets to be analyzed and for understanding protein function.     

Homology modeling, docking studies and functional analysis of various azoreductase accessory interacting proteins of Nostoc sp.PCC7120

Azo dyes have become a threat to public health because of its toxicity and carcinogenicity. Azoreductase enzyme plays a pivotal role in the degradation of azodyes released by industrial effluents and other resources. The degradation pathway has to be studied in detail for increasing the activity of azoreductase and for better degradation of azo dyes. But the data available on cyanobacterial azoreductase enzyme and its degradation pathway are still very less. Therefore the present work explored the azoreductase pathway of the cyanobacterium Nostoc sp. PCC7120 for better understanding of the degradation pathway and the other accessory interacting proteins involved. The accessory interacting proteins of azoreductase from cyanobacterium Nostoc sp. PCC7120 were obtained from STRING database. The proteins do not have a comprehensive three dimensional structure and are hypothetical. The secondary structure and functional analysis indicated that the proteins are all soluble proteins, without disulphide bonds and have alpha helices only. The structural prediction and docking study showed that alr2106, alr1063 and alr2326 have best docking result which tally with the STRING database confidence score and thus these proteins could possibly enhance the azoreductase activity and better dye degradation. These results will pave way for further increase in azoreductase activity and for better understanding of the dye degradation pathway.     

Solution structure of HI1506, a novel two-domain protein from Haemophilus influenzae

HI1506 is a 128-residue hypothetical protein of unknown function from Haemophilus influenzae. It was originally annotated as a shorter 85-residue protein, but a more detailed sequence analysis conducted in our laboratory revealed that the full-length protein has an additional 43 residues on the C terminus, corresponding with a region initially ascribed to HI1507. As part of a larger effort to understand the functions of hypothetical proteins from Gram-negative bacteria, and H. influenzae in particular, we report here the three-dimensional solution NMR structure for the corrected full-length HI1506 protein. The structure consists of two well-defined domains, an alpha/beta 50-residue N-domain and a 3-alpha 32-residue C-domain, separated by an unstructured 30-residue linker. Both domains have positively charged surface patches and weak structural homology with folds that are associated with RNA binding, suggesting a possible functional role in binding distal nucleic acid sites.     

Homology modeling, molecular docking and electrostatic potential analysis of MurF ligase from Klebsiella pneumonia

In spite of availability of moderately protective vaccine and antibiotics, new antibacterial agents are urgently needed to decrease the global incidence of Klebsiella pneumonia infections. MurF ligase, a key enzyme, which participates in the bacterial cell wall assembly, is indispensable to existence of K. pneumonia. MurF ligase lack mammalian vis-à-vis and have high specificity, uniqueness, and occurrence only in eubacteria, epitomizing them as promising therapeutic targets for intervention. In this study, we present a unified approach involving homology modeling and molecular docking studies on MurF ligase enzyme. As part of this study, a homology model of K. pneumonia (MurF ligase) enzyme was predicted for the first time in order to carry out structurebased drug design. The accuracy of the model was further validated using different computational approaches. The comparative molecular docking study on this enzyme was undertaken using different phyto-ligands from Desmodium sp. and a known antibiotic Ciprofloxacin. The docking analysis indicated the importance of hotspots (HIS 281 and ASN 282) within the MurF binding pocket. The Lipinski's rule of five was analyzed for all ligands considered for this study by calculating the ADME/Tox, drug likeliness using Qikprop simulation. Only ten ligands were found to comply with the Lipinski rule of five. Based on the molecular docking results and Lipinki values 6-Methyltetrapterol A was confirmed as a promising lead compound. The present study should therefore play a guiding role in the experimental design and development of 6-Methyltetrapterol A as a bactericidal agent.     

Homology modeling of an antifungal metabolite plipastatin synthase from the Bacillus subtilis 168

Lipopeptides have a widespread role in different pathways of Bacillus subtilis; they can act as antagonists, spreader and immunostimulators. Plipastatin, an antifungal antibiotic, is one of the most important lipopeptide nonribosomly produced by Bacillus subtilis. Plipastatin has strong fungitoxic activity and involve in inhibition of phospholipase A2 and biofilm formation. For better understanding of the molecule and pathway by which lipopeptide plipastatin is synthesized, we present a computationally predicted structure of plipastatin using homology modeling. Primary and secondary structure analysis suggested that ppsD is a hydrophilic protein containing a significant proportion of alpha helices, and subcellular localization predictions suggested it is a cytoplasmic protein. The tertiary structure of protein (plipastatin synthase subunit D) was predicted by homology modeling. The results suggest a flexible structure which is also an important characteristic of active enzymes enabling them to bind various cofactors and substrates for proper functioning. Validation of 3D structure was done using Ramachandran plot ProsA-web and QMEAN score.This predicted information will help in better understanding of mechanisms underlying plipastatin synthase subunit D synthesis. Plipastatin can be used as an inhibitor of various fungal diseases in plants.     

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on Cα atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 Å, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low. © 1995 Wiley-Liss, Inc.     

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.     

Homology modeling of an immunoglobulin-like domain in the Saccharomyces cerevisiae adhesion protein alpha-agglutinin.

The Saccharomyces cerevisiae adhesion protein alpha-agglutinin is expressed by cells of alpha mating type. On the basis of sequence similarities, alpha-agglutinin has been proposed to contain variable-type immunoglobulin-like (IgV) domains. The low level of sequence similarity to IgV domains of known structure made homology modeling using standard sequence-based alignment algorithms impossible. We have therefore developed a secondary structure-based method that allowed homology modeling of alpha-aggulutinin domain III, the domain most similar to IgV domains. The model was assessed and where necessary refined to accommodate information obtained by biochemical and molecular genetic approaches, including the positions of a disulfide bond, glycosylation sites, and proteolytic sites. The model successfully predicted surface exposure of glycosylation and proteolytic sites, as well as identifying residues essential for binding activity. One side of the domain was predicted to be covered by carbohydrate residues. Surface accessibility and volume packing analyses showed that the regions of the model that have greatest sequence dissimilarity from the IgV consensus sequence are poorly structured in the biophysical sense. Nonetheless, the utility of the model suggests that these alignment and testing techniques should be of general use for building and testing of models of proteins that share limited sequence similarity with known structures.     

Refinement of protein structures by iterative comparative modeling and CryoEM density fitting

We developed a method for structure characterization of assembly components by iterative comparative protein structure modeling and fitting into cryo-electron microscopy (cryoEM) density maps. Specifically, we calculate a comparative model of a given component by considering many alternative alignments between the target sequence and a related template structure while optimizing the fit of a model into the corresponding density map. The method relies on the previously developed Moulder protocol that iterates over alignment, model building, and model assessment. The protocol was benchmarked using 20 varied target-template pairs of known structures with less than 30% sequence identity and corresponding simulated density maps at resolutions from 5A to 25A. Relative to the models based on the best existing sequence profile alignment methods, the percentage of C(alpha) atoms that are within 5A of the corresponding C(alpha) atoms in the superposed native structure increases on average from 52% to 66%, which is half-way between the starting models and the models from the best possible alignments (82%). The test also reveals that despite the improvements in the accuracy of the fitness function, this function is still the bottleneck in reducing the remaining errors. To demonstrate the usefulness of the protocol, we applied it to the upper domain of the P8 capsid protein of rice dwarf virus that has been studied by cryoEM at 6.8A. The C(alpha) root-mean-square deviation of the model based on the remotely related template, bluetongue virus VP7, improved from 8.7A to 6.0A, while the best possible model has a C(alpha) RMSD value of 5.3A. Moreover, the resulting model fits better into the cryoEM density map than the initial template structure. The method is being implemented in our program MODELLER for protein structure modeling by satisfaction of spatial restraints and will be applicable to the rapidly increasing number of cryoEM density maps of macromolecular assemblies.     

Homology modeling of histidine-containing phosphocarrier protein and eosinophil-derived neurotoxin: Construction of models and comparison with experiment

Homology modeling methods have been used to construct models of two proteins—the histidine-containing phosphocarrier protein (HPr) from Mycoplasma capricolum and human eosinophil-derived neurotoxin (EDN). Comparison of the models with the subsequently determined X-ray crystal structures indicates that the core regions of both proteins are reasonably well reproduced, although the template structures are closer to the X-ray structures in these regions—possible enhancements are discussed. The conformations of most of the side chains in the core of HPr are well reproduced in the modeled structure. As expected, the conformations of surface side chains in this protein differ significantly from the X-ray structure. The loop regions of EDN were incorrectly modeled—reasons for this and possible enhancements are discussed. © 1995 Wiley-Liss, Inc.     

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.     

Sampling of near-native protein conformations during protein structure refinement using a coarse-grained model, normal modes, and molecular dynamics simulations

Protein structure refinement from comparative models with the goal of predicting structures at near-experimental accuracy remains an unsolved problem. Structure refinement might be achieved with an iterative protocol where the most native-like structure from a set of decoys generated from an initial model in one cycle is used as the starting structure for the next cycle. Conformational sampling based on the coarse-grained SICHO model, atomic level of detail molecular dynamics simulations, and normal-mode analysis is compared in the context of such a protocol. All of the sampling methods can achieve significant refinement close to experimental structures, although the distribution of structures and the ability to reach native-like structures differs greatly. Implications for the practical application of such sampling methods and the requirements for scoring functions in an iterative refinement protocol are analyzed in the context of theoretical predictions for the distribution of protein-like conformations with a random sampling protocol.     

Structure and action of the N-oxygenase AurF from Streptomyces thioluteus

Nitro groups are found in a number of bioactive compounds. Most of them arise by a stepwise mono-oxygenation of amino groups. One of the involved enzymes is AurF participating in the biosynthesis of aureothin. Its structure was established at 2.1 A resolution showing a homodimer with a binuclear manganese cluster. The enzyme preparation, which yielded the analyzed crystals, showed activity using in vitro and in vivo assays. Chain fold and cluster are homologous with ribonucleotide reductase subunit R2 and related enzymes. The two manganese ions and an iron content of about 15% were established by anomalous X-ray diffraction. A comparison of the cluster with more common di-iron clusters suggested an additional histidine in the coordination sphere to cause the preference for manganese over iron. There is no oxo-bridge. The substrate p-amino-benzoate was modeled into the active center. The model is supported by mutant activity measurements. It shows the geometry of the reaction and explains the established substrate spectrum.