High temperature unfolding of Bacillus anthracis amidase-03 by molecular dynamics simulations

The stability of amidase-03 structure (a cell wall hydrolase protein) from Bacillus anthracis was studied using classical molecular dynamics (MD) simulation. This protein (GenBank accession number: {"type":"entrez-protein","attrs":{"text":"NP_844822","term_id":"30262445","term_text":"NP_844822"}}NP_844822) contains an amidase-03 domain which is known to exhibit the catalytic activity of N-acetylmuramoyl-L-alanine amidase (digesting MurNAc-Lalanine linkage of bacterial cell wall). The amidase-03 enzyme has stability at high temperature due to the core formed by the combination of several secondary structure elements made of β-sheets. We used root-mean-square-displacement (RMSD) of the simulated structure from its initial state to demonstrate the unfolding of the enzyme using its secondary structural elements. Results show that amidase-03 unfolds in transition state ensemble (TSE). The data suggests that α-helices unfold before β-sheets from the core during simulation.     

A molecular model for diacylglycerol acyltransferase from Mortierella ramanniana var. angulispora

Acyl CoA diacylglycerol acyltransferase (DGAT, EC 2.3.120) is recognized as a key player of cellular diacylglycerol metabolism. It catalyzes the terminal, yet the committed step in triacylglycerol synthesis using diacylglycerol and fatty acyl CoA as substrates. The protein sequence of diacylglycerol acyltransferse (DGAT) Type 2B in Moretierella ramanniana var. angulispora (Protein_ID = {"type":"entrez-protein","attrs":{"text":"AAK84180.1","term_id":"15099961"}}AAK84180.1) was retrieved from GenBank. However, a structure is not yet available for this sequence. The 3D structure of DGAT Type 2B was modeled using a template structure (PDB ID: 1K30) obtained from Protein databank (PDB) identified by searching with position specific iterative BLAST (PSI-BLAST). The template (PDB ID: 1K30) describes the structure of DGAT from Cucurbita moschata. Modeling was performed using Modeller 9v2 and protein model is hence generated. The DGAT type 2B protein model was subsequently docked with six inhibitors (sphingosine; trifluoroperazine; phosphatidic acid; lysophospatidylserine; KCl; 1, 2-diolein) using AutoDock (a molecular docking program). The binding of inhibitors to the protein model of DGAT type 2B is discussed.     

Insight from the structural molecular model of cytidylate kinase from Mycobacterium tuberculosis

Mycobacterium tuberculosis is a gram-positive bacterium causes tuberculosis in human. H37Rv strain is a pathogenic strain utilized for tuberculosis research. The cytidylate mono-phosphate (CMP) kinase of Mycobacterium tuberculosis belongs to the family nucleoside mono-phosphate kinase (NMK), this enzyme is required for the bacterial growth. Therefore, it is important to study the structural and functional features of this enzyme in the control of the disease. Hence, we developed the structural molecular model of the CMP kinase protein from Mycobacterium tuberculosis by homology modeling using the software MODELLER (9v10). Based on sequence similarity with protein of known structure (template) of Mycobacterium smegmatis (PDB ID: 3R20) was chosen from protein databank (PDB) by using BLASTp. The energy of constructed models was minimized and the qualities of the models were evaluated by PROCHECK and VERRIFY-3D. Resulted Ramachandran plot analysis showed that conformations for 100.00% of amino acids residues are within the most favored regions. A possible homologous deep cleft active site was identified in the Model using CASTp program. Amino acid composition and polarity of that protein was observed by CLC-Protein Workbench tool. Expasy''s Prot-param server and CYC_REC tool were used for physiochemical and functional characterization of the protein. Studied of secondary structure of that protein was carried out by computational program, ProFunc. The structure is finally submitted in Protein Model Database. The predicted model permits initial inferences about the unexplored 3D structure of the CMP kinase and may be promote in relational designing of molecules for structure-function studies.     

Surface Plasmon Resonance Biosensor for Detection of Bacillus anthracis, the Causative Agent of Anthrax from Soil Samples Targeting Protective Antigen

Bacillus anthracis, the causative agent of anthrax is one of the most important biological warfare agents. In this study, surface plasmon resonance (SPR) technology was used for indirect detection of B. anthracis by detecting protective antigen (PA), a common toxin produced by all live B. anthracis bacteria. For development of biosensor, a monoclonal antibody raised against B. anthracis PA was immobilized on carboxymethyldextran modified gold chip and its interaction with PA was characterized in situ by SPR and electrochemical impedance spectroscopy. By using kinetic evaluation software, K_D (equilibrium constant) and B_max (maximum binding capacity of analyte) were found to be 20 fM and 18.74, respectively. The change in Gibb’s free energy (∆G = −78.04 kJ/mol) confirmed the spontaneous interaction between antigen and antibody. The assay could detect 12 fM purified PA. When anthrax spores spiked soil samples were enriched, PA produced in the sample containing even a single spore of B. anthracis could be detected by SPR. PA being produced only by the vegetative cells of B. anthracis, confirms indirectly the presence of B. anthracis in the samples. The proposed method can be a very useful tool for screening and confirmation of anthrax suspected environmental samples during a bio-warfare like situation.     

A novel approach for large-scale polypeptide folding based on elastic networks using continuous optimization

We present a new computationally efficient method for large-scale polypeptide folding using coarse-grained elastic networks and gradient-based continuous optimization techniques. The folding is governed by minimization of energy based on Miyazawa-Jernigan contact potentials. Using this method we are able to substantially reduce the computation time on ordinary desktop computers for simulation of polypeptide folding starting from a fully unfolded state. We compare our results with available native state structures from Protein Data Bank (PDB) for a few de-novo proteins and two natural proteins, Ubiquitin and Lysozyme. Based on our simulations we are able to draw the energy landscape for a small de-novo protein, Chignolin. We also use two well known protein structure prediction software, MODELLER and GROMACS to compare our results. In the end, we show how a modification of normal elastic network model can lead to higher accuracy and lower time required for simulation.     

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.     

Identification of a Bacillus anthracis specific indel in the yeaC gene and development of a rapid pyrosequencing assay for distinguishing B. anthracis from the B. cereus group

Bacillus anthracis, the causative agent of anthrax, is a potential source of bioterrorism. The existing assays for its identification lack specificity due to the close genetic relationship it exhibits to other members of the B. cereus group. Our comparative analyses of protein sequences from Bacillus species have identified a 24 amino acid deletion in a conserved region of the YeaC protein that is uniquely present in B. anthracis. PCR primers based on conserved regions flanking this indel in the Bacillus cereus group of species (viz. Bacillus cereus, B. anthracis, B. thuringiensis, B. mycoides, B. weihenstephnensis and B. pseudomycoides) specifically amplified a 282 bp fragment from all six reference B. anthracis strains, whereas a 354 bp fragment was amplified from 15 other B. cereus group of species/strains. These fragments, due to large size difference, are readily distinguished by means of agarose gel electrophoresis. In contrast to the B. cereus group, no PCR amplification was observed with any of the non-B. cereus group of species/strains. This indel was also used for developing a rapid pyrosequencing assay for the identification of B. anthracis. Its performance was evaluated by examining the presence or absence of this indel in a panel of 81 B. cereus-like isolates from various sources that included 39 B. anthracis strains. Based upon the sequence data from the pyrograms, the yeaC indel was found to be a distinctive characteristic of various B. anthracis strains tested and not found in any other species/strains from these samples. Therefore, this B. anthracis specific indel provides a robust and highly-specific chromosomal marker for the identification of this high-risk pathogen from other members of the B. cereus group independent of a strain's virulence. The pyrosequencing platform also allows for the rapid and simultaneous screening of multiple samples for the presence of this B. anthracis-specific marker.     

Insight to pyrazinamide resistance inMycobacterium tuberculosisby molecular docking

Pyrazinamide (PZA) - an important drug in the anti-tuberculosis therapy, activated by an enzyme Pyrazinamidase (PZase). The basis of PZA resistance in Mycobacterium tuberculosis was owing to mutation in pncA gene coding for PZase. Homology modeling of PZase was performed using software Discovery Studio (DS) 2.0 based on the crystal structure of the PZase from Pyrococcus horikoshii (PDB code 1im5), in this study. The model comprises of one sheet with six parallel strands and seven helices with the amino acids Asp8, Asp49, Trp68, Lys96, Ala134, Thr135 and Cys138 at the active site. Five mutants were generated with Gly at position 8, Thr at position 96, Arg at position 104, Tyr and Ser at position 138. The Wild-type (WT) and five mutant models were docked with PZA. The results indicate that the mutants Lys96Thr, Ser104Arg Asp8Gly and Cys138Tyr may contribute to higher level drug resistance than Cys138Ser. These models provide the first in-silico evidence for the binding interaction of PZA with PZase and form the basis for rationalization of PZA resistance in naturally occurring pncA mutant strains of M. tuberculosis.     

Protein-protein docking on molecular models of Aspergillus niger RNase and human actin: novel target for anticancer therapeutics

The 3D models of human actin protein and A.niger RNase were designed using the templates ACTBIND (PDB ID: 3D3Z) and crystalline profilin-beta-actin (PDB ID: 2BTF), respectively in Modeller9v5. These models are testified using several validation methods including PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The stereo-chemical quality of the models was judged by Ramachandran plot with PROCHECK. The total quality G-factor −0.2, shows a good quality model. The ERRAT score for the human actin and A.niger RNase models are 86.104 and 84.615, respectively, fit well within the range of a high quality model. The ERRAT score for the templates 2BTF and 3D3Z are 91.111 and 97.391, respectively. The WHAT-IF evaluation justifies a reasonable homology model structure as none of the scores for each residue in the homology model is lower than −5.0. The energy-minimized model of human actin with PROSA reveals the Z-score value −10.52 between native conformations of the crystal structures. The VERIFY 3D average score is 0.36. All evidence suggests that the geometric quality of the backbone conformation, the residue interaction, the residue contact and the energy profile of the structures were well within the limits of reliable structures. The interaction energy of docking was calculated using the HEX server. The Etotal, lowest docked energy, and calculated RMSD values were −1.608 kcal mol⁻¹, -8.369 kcal mol⁻¹ and 0.617 Å, respectively. The study presented in the current project may be useful to design molecules that may have anticancer activity.

     

Peptides panned from a phage-displayed random peptide library are useful for the detection of Bacillus anthracis surrogates B. cereus 4342 and B. anthracis Sterne

Recent use of Bacillus anthracis as a bioweapon has highlighted the need for a sensitive monitoring system. Current bacterial detection tests use antibodies as bio-molecular recognition elements which have limitations with regard to time, specificity and sensitivity, creating the need for new and improved cost-effective high-affinity detection probes. In this study, we screened a commercially available bacteriophage-displayed random peptide library using Bacillus cereus 4342 cells as bait to identify peptides that could be used for detection of Bacillus. The method enabled us to identify two 12-amino acid consensus peptide sequences that specifically bind to B. cereus 4342 and B. anthracis Sterne, the nonpathogenic surrogates of B. anthracis strain. The two Bacillus-binding peptides (named BBP-1 and BBP-2) were synthesized with biotin tag to confirm their binding by four independent detection assays. Dot-blot analysis revealed that the peptides bind specifically to B. cereus 4342 and B. anthracis Sterne. Quantitative analysis of this interaction by ELISA and fluorometry demonstrated a detection sensitivity of 10² colony forming U/ml (CFU/ml) by both assays. When the peptides were used in combination with Qdots, the sensitivity was enhanced further by enabling detection of even a single bacterium by fluorescence microscopy. Immunoblot analysis and protein sequencing showed that BBP-1 and BBP-2 bound to the S-layer protein of B. anthracis Sterne. Overall, our findings validate the usefulness of synthetic versions of phage-derived peptides in combination with Qdot-liquid nanocrystals as high sensitivity bioprobes for various microbial detection platforms.     

Genomic comparison of 93 Bacillus phages reveals 12 clusters, 14 singletons and remarkable diversity

Background

The Bacillus genus of Firmicutes bacteria is ubiquitous in nature and includes one of the best characterized model organisms, B. subtilis, as well as medically significant human pathogens, the most notorious being B. anthracis and B. cereus. As the most abundant living entities on the planet, bacteriophages are known to heavily influence the ecology and evolution of their hosts, including providing virulence factors. Thus, the identification and analysis of Bacillus phages is critical to understanding the evolution of Bacillus species, including pathogenic strains.

Results

Whole genome nucleotide and proteome comparison of the 93 extant Bacillus phages revealed 12 distinct clusters, 28 subclusters and 14 singleton phages. Host analysis of these clusters supports host boundaries at the subcluster level and suggests phages as vectors for genetic transfer within the Bacillus cereus group, with B. anthracis as a distant member of the group. Analysis of the proteins conserved among these phages reveals enormous diversity and the uncharacterized nature of these phages, with a total of 4,922 protein families (phams) of which only 951 (19%) had a predicted function. In addition, 3,058 (62%) of phams were orphams (phams containing a gene product from a single phage). The most populated phams were those encoding proteins involved in DNA metabolism, virion structure and assembly, cell lysis, or host function. These included several genes that may contribute to the pathogenicity of Bacillus strains.

Conclusions

This analysis provides a basis for understanding and characterizing Bacillus phages and other related phages as well as their contributions to the evolution and pathogenicity of Bacillus cereus group bacteria. The presence of sparsely populated clusters, the high ratio of singletons to clusters, and the large number of uncharacterized, conserved proteins confirms the need for more Bacillus phage isolation in order to understand the full extent of their diversity as well as their impact on host evolution.     

Bacillusanthracis Surface-Layer Proteins Assemble by Binding to the Secondary Cell Wall Polysaccharide in a Manner that Requires csaB and tagO

Bacillus anthracis, the causative agent of anthrax, requires surface (S)-layer proteins for the pathogenesis of infection. Previous work characterized S-layer protein binding via the surface layer homology domain to a pyruvylated carbohydrate in the envelope of vegetative forms. The molecular identity of this carbohydrate and the mechanism of its display in the bacterial envelope are still unknown. Analyzing acid-solubilized, purified carbohydrates by mass spectrometry and NMR spectroscopy, we identify secondary cell wall polysaccharide (SCWP) as the ligand of S-layer proteins. In agreement with the model that surface layer homology domains bind to pyruvylated carbohydrate, SCWP was observed to be linked to pyruvate in a manner requiring csaB, the only structural gene known to be required for S-layer assembly. B. anthracis does not elaborate wall teichoic acids; however, its genome harbors tagO and tagA, genes responsible for the synthesis of the linkage unit that tethers teichoic acids to the peptidoglycan layer. The tagO gene appears essential for B. anthracis growth and complements the tagO mutant phenotypes of staphylococci. Tunicamycin-mediated inhibition of TagO resulted in deformed, S-layer-deficient bacilli. Together, these results suggest that tagO-mediated assembly of linkage units tethers pyruvylated SCWP to the B. anthracis envelope, thereby enabling S-layer assembly and providing for the pathogenesis of anthrax infections.     

Erratum to: genomic comparison of 93 Bacillus phages reveals 12 clusters, 14 singletons and remarkable diversity

Background

The Bacillus genus of Firmicutes bacteria is ubiquitous in nature and includes one of the best characterized model organisms, B. subtilis, as well as medically significant human pathogens, the most notorious being B. anthracis and B. cereus. As the most abundant living entities on the planet, bacteriophages are known to heavily influence the ecology and evolution of their hosts, including providing virulence factors. Thus, the identification and analysis of Bacillus phages is critical to understanding the evolution of Bacillus species, including pathogenic strains.

Results

Whole genome nucleotide and proteome comparison of the 83 extant, fully sequenced Bacillus phages revealed 10 distinct clusters, 24 subclusters and 15 singleton phages. Host analysis of these clusters supports host boundaries at the subcluster level and suggests phages as vectors for genetic transfer within the Bacillus cereus group, with B. anthracis as a distant member. Analysis of the proteins conserved among these phages reveals enormous diversity and the uncharacterized nature of these phages, with a total of 4,442 protein families (phams) of which only 894 (20%) had a predicted function. In addition, 2,583 (58%) of phams were orphams (phams containing a single member). The most populated phams were those encoding proteins involved in DNA metabolism, virion structure and assembly, cell lysis, or host function. These included several genes that may contribute to the pathogenicity of Bacillus strains.

Conclusions

This analysis provides a basis for understanding and characterizing Bacillus and other related phages as well as their contributions to the evolution and pathogenicity of Bacillus cereus group bacteria. The presence of sparsely populated clusters, the high ratio of singletons to clusters, and the large number of uncharacterized, conserved proteins confirms the need for more Bacillus phage isolation in order to understand the full extent of their diversity as well as their impact on host evolution.     

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.     

Molecular interactions between Bos taurus interferon-τ1c and human type I interferon receptor

Interferon (IFN)-τ secreted only by ruminant endometrium, helps in maternal recognition of pregnancy and exhibit antiviral and antiproliferative activity. Among different types of IFN-τ, IFN-τ1c and IFN- τ3a are the most highly expressed isoforms. In the present study structure of INF-τ1c was predicted using homology modelling. The best model was selected based on overall stereo-chemical quality. The generated 3D structure of the Interferon-τ1c protein of Bos taurus was predicted using the ovine interferon-τ (PDB ID: 1B5L_A) as template. The structure comprises of 5 α helices separated by loop regions, which is similar to the one predicted for other IFNs. Molecular interactions of bovine IFN-τ1c with human interferon Type 1 receptor (IFNAR1) was explored in an attempt to predict human IFNAR1 binding sites of IFN-τ1c.     

Affinity comparison of different THCA synthase to CBGA using modeling computational approaches

The Δ^9-Tetrahydrocannabinol (THCA) is the primary psychoactive compound of Cannabis Sativa. It is produced by Δ^1- Tetrahydrocannabinolic acid synthase (THCA) which catalyzes the oxidative cyclization of cannabigerolic acid (CBGA) the precursor of the THCA. In this study, we were interested by the three dimensional structure of THCA synthase protein. Generation of models were done by MODELLER v9.11 and homology modeling with Δ1-tetrahydrocannabinolic acid (THCA) synthase X ray structure (PDB code 3VTE) on the basis of sequences retrieved from GenBank. Procheck, Errat, and Verify 3D tools were used to verify the reliability of the six 3D models obtained, the overall quality factor and the Prosa Z-score were also used to check the quality of the six modeled proteins. The RMSDs for C-alpha atoms, main-chain atoms, side-chain atoms and all atoms between the modeled structures and the corresponding template ranged between 0.290 Å-1.252 Å, reflecting the good quality of the obtained models. Our study of the CBGA-THCA synthase docking demonstrated that the active site pocket was successfully recognized using computational approach. The interaction energy of CBGA computed in ‘fiber types’ proteins ranged between -4.1 95 kcal/mol and -5.95 kcal/mol whereas in the ‘drug type’ was about -7.02 kcal/mol to -7.16 kcal/mol, which maybe indicate the important role played by the interaction energy of CBGA in the determination of the THCA level in Cannabis Sativa L. varieties. Finally, we have proposed an experimental design in order to explore the binding energy source of ligand-enzyme in Cannabis Sativa and the production level of the THCA in the absence of any information regarding the correlation between the enzyme affinity and THCA level production. This report opens the doors to more studies predicting the binding site pocket with accuracy from the perspective of the protein affinity and THCA level produced in Cannabis Sativa.     

Filament Formation of the FtsZ/Tubulin-like Protein TubZ from the Bacillus cereus pXO1 Plasmid

Stable maintenance of low-copy-number plasmids requires partition (par) systems that consist of a nucleotide hydrolase, a DNA-binding protein, and a cis-acting DNA-binding site. The FtsZ/tubulin-like GTPase TubZ was identified as a partitioning factor of the virulence plasmids pBtoxis and pXO1 in Bacillus thuringiensis and Bacillus anthracis, respectively. TubZ exhibits high GTPase activity and assembles into polymers both in vivo and in vitro, and its “treadmilling” movement is required for plasmid stability in the cell. To investigate the molecular mechanism of pXO1 plasmid segregation by TubZ filaments, we determined the crystal structures of Bacillus cereus TubZ in apo-, GDP-, and guanosine 5′-3-O-(thio)triphosphate (GTPγS)-bound forms at resolutions of 2.1, 1.9, and 3.3 Å, respectively. Interestingly, the slowly hydrolyzable GTP analog GTPγS was hydrolyzed to GDP in the crystal. In the post-GTP hydrolysis state, GDP-bound B. cereus TubZ forms a dimer by the head-to-tail association of individual subunits in the asymmetric unit, which is similar to the protofilament formation of FtsZ and B. thuringiensis TubZ. However, the M loop interacts with the nucleotide-binding site of the adjacent subunit and stabilizes the filament structure in a different manner, which indicates that the molecular assembly of the TubZ-related par systems is not stringently conserved. Furthermore, we show that the C-terminal tail of TubZ is required for association with the DNA-binding protein TubR. Using a combination of crystallography, site-directed mutagenesis, and biochemical analysis, our results provide the structural basis of the TubZ polymer that may drive DNA segregation.