In silico characterization and structural modeling of bacterial metalloprotease of family M4

Lactiplantibacillus plantarum LM1001 Improves Digestibility of Branched-Chain Amino Acids in Whey Proteins and Promotes Myogenesis in C2C12 Myotubes

Lactiplantibacillus plantarum is a valuable potential probiotic species with various proven health-beneficial effects. L. plantarum LM1001 strain was selected among ten strains of L. plantarum based on proteolytic activity on whey proteins. L. plantarum LM1001 produced higher concentrations of total free amino acids and branched-chain amino acids (Ile, Leu, and Val) than other L. plantarum strains. Treatment of C2C12 myotubes with whey protein culture supernatant (1%, 2% and 3%, v/v) using L. plantarum LM1001 significantly increased the expression of myogenic regulatory factors, such as Myf-5, MyoD, and myogenin, reflecting the promotion of myotubes formation (p<0.05). L. plantarum LM1001 displayed β-galactosidase activity but did not produce harmful β-glucuronidase. Thus, the intake of whey protein together with L. plantarum LM1001 has the potential to aid protein digestion and utilization.

Copper Forms a PPII Helix-Like Structure with the Catalytic Domains of Bacterial Zinc Metalloproteases

The rapid spread of antibiotic-resistant bacteria continuously raises concerns about the future ineffectiveness of current antimicrobial treatments against infectious diseases. To address this problem, new therapeutic strategies and antimicrobial drugs with unique modes of action are urgently needed. Inhibition of metalloproteases, bacterial virulence factors, is a promising target for the development of antibacterial treatments. In this study, the interaction among Zn(II), Cu(II), and the metal-binding domains of two metalloproteases, AprA (Pseudomonas aureginosa) and CpaA (Acinetobacter baumanii), was investigated. The objective was to determine the coordination sphere of Zn(II) with a peptide model of two zinc-dependent metalloproteases. Additionally, the study explored the formation of Cu(II) complexes with the domains, as Cu(II) has been shown to inhibit metalloproteases. The third aim was to understand the role of nonbinding amino acids in stabilizing the metal complexes formed by these proteases. This work identified specific coordination patterns (HExxHxxxxxH) for both Zn(II) and Cu(II) complexes, with AprA and CpaA exhibiting a higher affinity for Cu(II) compared to Zn(II). The study also found that the CpaA domain has greater stability for both Zn(II) and Cu(II) complexes compared to AprA. The nonbinding amino acids of CpaA surrounding the metal ion contribute to the increased thermodynamic stability of the metal-peptide complex through various intramolecular interactions. These interactions can also influence the secondary structures of the peptides. The presence of certain amino acids, such as tyrosine, arginine, and glutamic acid, and their interactions contribute to the stability and, only in the case of Cu(II) complexes, the formation of a rare protein structure called a left-handed polyproline II helix (PPII), which is known to play a role in the stability and function of various proteins. These findings provide valuable insights into the coordination chemistry of bacterial metalloproteases and expand our understanding of potential mechanisms for inhibiting these enzymes.

Genetic insights into obesity: in silico identification of pathogenic SNPs in MBOAT4 gene and their structural molecular dynamics consequences

Membrane Bound O-Acyltransferase Domain-Containing 4 (MBOAT4) protein catalyzes ghrelin acylation, leading to prominent ghrelin activity, hence characterizing its role as an anti-obesity target. We extracted 625 exonic SNPs from the ENSEMBL database and one phenotype-based missense mutation associated with obesity (A46T) from the HGMD (Human Gene Mutation Database). These were differentiated on deleterious missense SNPs of the MBOAT4 gene through MAF (minor allele frequency: <0.01) cut-off criteria in relation to some bioinformatics-based supervised machine learning tools. We found 8 rare-coding and harmful missense SNPs. The consensus classifier (PredictSNP) tool predicted that the SNP (G57S, C: rs561065025) was the most pathogenic. Several trained in silico algorithms have predicted decreased protein stability [ΔΔG (kcal/mol)] function in the presence of these rare-coding pathogenic mutations in the MBOAT4 gene. Then, a stereochemical quality check (i.e. validation and assessment) of the 3D model was performed, followed by a blind cavity docking approach, used to search for druggable cavities and molecular interactions with citrus flavonoids of the Rutaceae family, ranked with energetic estimations. Significant interactions with Phloretin 3',5'-Di-C-Glucoside were also observed at R304, W306, N307, A311, L314 and H338 with (iGEMDOCK: -95.82 kcal/mol and AutoDock: -7.80 kcal/mol). The RMSD values and other variables of MD simulation analyses on this protein further validated its significant interactions with the above flavonoids. The MBOAT4 gene and its molecular interactions could serve as an interventional future anti-obesity target. The current study's findings will benefit future prospects for large population-based studies and drug development, particularly for generating personalized medicine.Communicated by Ramaswamy H. Sarma.

In Silico Evaluation, Phylogenetic Analysis, and Structural Modeling of the Class II Hydrophobin Family from Different Fungal Phytopathogens

The class II hydrophobin group (HFBII) is an extracellular group of proteins that contain the HFBII domain and eight conserved cysteine residues. These proteins are exclusively secreted by fungi and have multiple functions with a probable role as effectors. In the present study, a total of 45 amino acid sequences of hydrophobin class II proteins from different phytopathogenic fungi were retrieved from the NCBI database. We used the integration of well-designed bioinformatic tools to characterize and predict their physicochemical parameters, novel motifs, 3D structures, multiple sequence alignment (MSA), evolution, and functions as effector proteins through molecular docking. The results revealed new features for these protein members. The ProtParam tool detected the hydrophobicity properties of all proteins except for one hydrophilic protein (KAI3335996.1). Out of 45 proteins, six of them were detected as GPI-anchored proteins by the PredGPI server. Different 3D structure templates with high pTM scores were designed by Multifold v1, AlphaFold2, and trRosetta. Most of the studied proteins were anticipated as apoplastic effectors and matched with the ghyd5 gene of Fusarium graminearum as virulence factors. A protein-protein interaction (PPI) analysis unraveled the molecular function of this group as GTP-binding proteins, while a molecular docking analysis detected a chitin-binding effector role. From the MSA analysis, it was observed that the HFBII sequences shared conserved 2 Pro (P) and 2 Gly (G) amino acids besides the known eight conserved cysteine residues. The evolutionary analysis and phylogenetic tree provided evidence of episodic diversifying selection at the branch level using the aBSREL tool. A detailed in silico analysis of this family and the present findings will provide a better understanding of the HFBII characters and evolutionary relationships, which could be very useful in future studies.

In Silico Functional Characterization of a Hypothetical Protein From Pasteurella Multocida Reveals a Novel S-Adenosylmethionine-Dependent Methyltransferase Activity

Genomes may now be sequenced in a matter of weeks, leading to an influx of "hypothetical" proteins (HP) whose activities remain a mystery in GenBank. The information included inside these genes has quickly grown in prominence. Thus, we selected to look closely at the structure and function of an HP (AFF25514.1; 246 residues) from Pasteurella multocida (PM) subsp. multocida str. HN06. Possible insights into bacterial adaptation to new environments and metabolic changes might be gained by studying the functions of this protein. The PM HN06 2293 gene encodes an alkaline cytoplasmic protein with a molecular weight of 28352.60 Da, an isoelectric point (pI) of 9.18, and an overall average hydropathicity of around -0.565. One of its functional domains, tRNA (adenine (37)-N6)-methyltransferase TrmO, is a S-adenosylmethionine (SAM)-dependent methyltransferase (MTase), suggesting that it belongs to the Class VIII SAM-dependent MTase family. The tertiary structures represented by HHpred and I-TASSER models were found to be flawless. We predicted the model's active site using the Computed Atlas of Surface Topography of Proteins (CASTp) and FTSite servers, and then displayed it in 3 dimensional (3D) using PyMOL and BIOVIA Discovery Studio. Based on molecular docking (MD) results, we know that HP interacts with SAM and S-adenosylhomocysteine (SAH), 2 crucial metabolites in the tRNA methylation process, with binding affinities of 7.4 and 7.5 kcal/mol, respectively. Molecular dynamic simulations (MDS) of the docked complex, which included only modest structural adjustments, corroborated the strong binding affinity of SAM and SAH to the HP. Evidence for HP's possible role as an SAM-dependent MTase was therefore given by the findings of Multiple sequence alignment (MSA), MD, and molecular dynamic modeling. These in silico data suggest that the investigated HP might be used as a useful adjunct in the investigation of Pasteurella infections and the development of drugs to treat zoonotic pasteurellosis.

Bacterial Zinc Metalloenzyme Inhibitors: Recent Advances and Future Perspectives

Human deaths caused by Gram-negative bacteria keep rising due to the multidrug resistance (MDR) phenomenon. Therefore, it is a priority to develop novel antibiotics with different mechanisms of action. Several bacterial zinc metalloenzymes are becoming attractive targets since they do not show any similarities with the human endogenous zinc-metalloproteinases. In the last decades, there has been an increasing interest from both industry and academia in developing new inhibitors against those enzymes involved in lipid A biosynthesis, and bacteria nutrition and sporulation, e.g., UDP-[3-O-(R)-3-hydroxymyristoyl]-N-acetylglucosamine deacetylase (LpxC), thermolysin (TLN), and pseudolysin (PLN). Nevertheless, targeting these bacterial enzymes is harder than expected and the lack of good clinical candidates suggests that more effort is needed. This review gives an overview of bacterial zinc metalloenzyme inhibitors that have been synthesized so far, highlighting the structural features essential for inhibitory activity and the structure-activity relationships. Our discussion may stimulate and help further studies on bacterial zinc metalloenzyme inhibitors as possible novel antibacterial drugs.

Comparative proteomics study of exosomes in Vibrio harveyi and Vibrio anguillarum

Exosomes are a class of extracellular vesicles released by bacteria and contain diverse biomolecules. In this study, we isolated exosomes from Vibrio harveyi and Vibrio anguillarum, which are both serious pathogens in mariculture, using a supercentrifugation method, and the proteins in the exosomes of these two vibrios were analyzed by LC-MS/MS proteomics. Exosome proteins released by V. harveyi and V. anguillarum were different; they not only contained virulence factors (such as lipase and phospholipase in V. harveyi, metalloprotease and hemolysin in V. anguillarum), but also participated in the important life activities of bacteria (such as fatty acid biosynthesis, biosynthesis of antibiotics, carbon metabolism). Subsequently, to verify whether the exosomes participated in bacterial toxicity, after Ruditapes philippinarum was challenged with V. harveyi and V. anguillarum, the corresponding genes of virulence factors from exosomes screened by proteomics were tested by quantitative real-time PCR. All the genes detected were upregulated which suggested that exosomes were involved in vibrio toxicity. The results could provide an effective proteome database for decoding the pathogenic mechanism of vibrios from the exosome perspective.

Identification and Functional Characterization of General Odorant Binding Proteins in Orthaga achatina

The olfactory system in insects are crucial for recognition of host plants and oviposition sites. General odorant binding proteins (GOBPs) are thought to be involved in detecting odorants released by host plants. Orthaga achatina (Lepidoptera: Pyralidae) is one of the most serious pests of camphor trees, Cinnamomum camphora (L.) Presl, an important urban tree species in southern China. In this study, we study the GOBPs of O. achatina. Firstly, two full-length GOBP genes (OachGOBP1 and OachGOBP2) were successfully cloned according to transcriptome sequencing results, and real-time quantitative PCR measurements showed that both GOBP genes were specifically expressed in the antennae of both sexes, proposing their important roles in olfaction. Then, both GOBP genes were heterologous expressed in Escherichia coli and fluorescence competitive binding assays were conducted. The results showed that OachGOBP1 could bind Farnesol (K_i = 9.49 μM) and Z11-16: OH (K_i = 1.57 μM). OachGOBP2 has a high binding affinity with two camphor plant volatiles (Farnesol, K_i = 7.33 μM; α-Phellandrene, K_i = 8.71 μM) and two sex pheromone components (Z11-16: OAc, K_i = 2.84 μM; Z11-16: OH, K_i = 3.30 μM). These results indicate that OachGOBP1 and OachGOBP2 differ in terms of odorants and other ligands. Furthermore, key amino acid residues that bind to plant volatiles were identified in GOBPs using 3-D structure modeling and ligand molecular docking, predicting the interactions between the GOBPs and the host plant volatiles.

Computational Study of Bis-(1-(Benzoyl)-3-Methyl Thiourea) Platinum (II) Complex Derivatives as Anticancer Candidates

Background

The increasing incidence of cancer every year has resulted in cancer becoming one of the most common causes of death in the world. The most common types of cancer are breast cancer, lung cancer and prostate cancer. Thiourea is one of the compounds that have anticancer effects, and its activity can be increased by structural modifications, one of which involves making a Bis-(1-(benzoyl)-3-methyl thiourea) platinum (II) metal complex.

Purpose

This study aims to obtain platinum (II)-thiourea complex compounds that have a more stable interaction as an anticancer agent compared to cisplatin.

Methods

The methods used are computational studies with molecular docking, simulation of molecular dynamics, and prediction of pharmacokinetics and toxicity.

Results

Based on the molecular docking of the platinum (II)-thiourea complex which has the most stable interaction with lower binding energy than the native ligand and the cisplatin, namely Bis-(3-methyl-1-(naphthalene-2-carbonyl)thiourea)) Platinum (II) against breast cancer receptors (3ERT) and lung cancer (2ITO) and compounds Bis-(1-(3-chlorobenzoyl)-3-methylthiourea) Platinum (II) against prostate cancer receptors (1Z95). The evaluation results of the stability of the interaction using a 50 ns molecular dynamic simulation showed that the Bis-(1-benzoyl-3-methylthiourea) Platinum (II) which binds to the prostate cancer receptor (1Z95) has the most stable interaction. Pharmacokinetic prediction results show that the platinum (II)-thiourea complex has a good pharmacokinetic profile, but there are several compounds that are mutagenic and hepatotoxic.

Conclusion

The Bis-(1-(3,4-dichlorobenzoyl)-3-methyl thiourea) platinum (II) compounds could be a suitable anticancer agent for the lungs.

Immunomodulatory role of metalloproteases in cancers: Current progress and future trends

Metalloproteinases (MPs) is a large family of proteinases with metal ions in their active centers. According to the different domains metalloproteinases can be divided into a variety of subtypes mainly including Matrix Metalloproteinases (MMPs), A Disintegrin and Metalloproteases (ADAMs) and ADAMs with Thrombospondin Motifs (ADAMTS). They have various functions such as protein hydrolysis, cell adhesion and remodeling of extracellular matrix. Metalloproteinases expressed in multiple types of cancers and participate in many pathological processes involving tumor genesis and development, invasion and metastasis by regulating signal transduction and tumor microenvironment. In this review, based on the current research progress, we summarized the structure of MPs, their expression and especially immunomodulatory role and mechanisms in cancers. Additionally, a relevant and timely update of recent advances and future directions were provided for the diagnosis and immunotherapy targeting MPs in cancers.

Computational insights into diverse aspects of glutathione S-transferase gene family in Papaver somniferum

Plant glutathione S-transferases are an ancient protein superfamily having antioxidant activity. These proteins are primarily involved in diverse plant functions such as plant growth and development, secondary metabolism, signaling pathways and defense against biotic and abiotic stresses. The current study aimed to comprehensively identify and characterize the GST gene family in the medicinally important crop Papaver somniferum. A total of 93 GST proteins were identified belonging to eight GST classes and found to be majorly localized in the cytoplasm. All GST genes were found on eleven opium chromosomes. Gene duplication analysis showed segmental duplication as a key factor for opium GST gene family expansion under strong purifying selection. Phylogenetic analysis with gymnosperm, angiosperm and bryophyte revealed the evolution of GSTs earlier than their division into separate groups and also prior to the divergence of monocot and dicot. The secondary structure prediction showed the dominance of α-helices indicative of PsomGSTs as structurally stable and elastic proteins. Gene architecture showed the conservation of number of exons across the classes. MEME analysis revealed only a few class specific and many across class conserved motifs. Ser was found to be the active site residue of tau, phi, theta and zeta class and Cys was catalytic residue of DHAR, lambda and GHR class. Promoter analyses identified many cis-acting regulatory elements related to hormonal, cellular, stress and light response functions. Ser was the key phosphorylation site. Only three glycosylation sites were found across the 93 PsomGSTs. 3D structure prediction was also performed and was validated. Interactome analyses revealed the correlation of PsomGSTs with glutathione metabolizing proteins. Gene enrichment analysis and KEGG pathway analyzed the involvement of PsomGSTs in three major pathways i.e. glutathione metabolism, tyrosine metabolism and ascorbate metabolism. The outcome revealed high model quality of PsomGSTs. The results of the current study will be of potential significance to understand the functional and structural importance of the GST gene family in opium, a medicinally important crop.

Investigation of bioactive compounds from Bacillus sp. against protein homologs CDC42 of Colletotrichum gloeosporioides causing anthracnose disease in cassava by using molecular docking and dynamics studies

Manihot esculenta, commonly called cassava, is an economically valuable crop and important staple food, grown in tropical and subtropical regions of the world. Demand for cassava in the food and fuel industry is growing worldwide. However, anthracnose disease caused by Colletotrichum gloeosporioides severely affects cassava yield and production. The bioactive molecules from Bacillus are widely used to control fungal diseases in several plants. Therefore, in this study, bioactive compounds (erucamide, behenic acid, palmitic acid, phenylacetic acid, and β-sitosterol) from Bacillus megaterium were assessed against CDC42, a key protein for virulence, from C. gloeosporioides. Structure of the CDC42 protein was generated through the comparative homology modeling method. The binding site of the ligands and the stability of the complex were analyzed through docking and molecular dynamics simulation studies, respectively. Furthermore, a protein interaction network was envisaged through the STRING database, followed by enrichment analysis in the WebGestalt tool. From the enrichment analysis, it is apparent that bioactive from B. megaterium chiefly targets the MAP kinase pathway that is essential for filamentous growth and virulence. Further exploration through experimental studies could be advantageous for cassava improvement as well as to combat against C. gloeosporioides pathogen.

Potential of the Liquid Fermentation of Fishery Waste by Paenibacillus elgii for Metalloprotease Production

This study attempted to use fishery processing wastes to produce protease by Paenibacillus elgii TKU051. Of the tested wastes, tuna head powder (THP) was found to be the most effective carbon and nitrogen (C/N) source, and the optimal conditions were as follows: 0.811% THP, 0.052% K₂HPO₄, 0.073% MgSO₄, initial pH of 8.96, incubation temperature of 31.4 °C, and incubation time of 3.092 days to achieve the maximum protease activity of 2.635 ± 0.124 U/mL. A protease with a molecular weight of 29 kDa was purified and biochemically characterized. Liquid chromatography with tandem mass spectrometry analysis revealed an amino acid sequence of STVHYSTR of P. elgii TKU051 protease, suggesting that the enzyme may belong to the M4 family of metalloproteases. The optimal activity of the enzyme was achieved at 60 °C and pH 8. P. elgii TKU051 protease was strongly inhibited by ethylenediaminetetraacetic acid and 1,10-phenanthroline, indicating its precise metalloprotease property. P. elgii TKU051 protease displayed the activity toward casein and raw fishery wastes such as tuna heads, tuna viscera, shrimp heads, and squid pens. Finally, the purified P. elgii TKU051 protease could improve the free-radical scavenging activity of fishery wastes. In short, P. elgii TKU051 has potential application in eco-friendly approaches to efficiently convert fishery wastes to metalloprotease.

Functional differentiation of three pheromone binding proteins in Orthaga achatina using mixed-type sex pheromones

Pheromone-binding proteins (PBPs) play important roles in perception of insect sex pheromones, functioning to recognize and transport pheromone components onto the olfactory receptors of the odorant sensing neurons. Orthaga achatina, a serious pest of camphor trees, uses a mixture of three Type I (Z11-16:OAc, Z11-16:OH and Z11-16:Ald) and one Type II (Z3,Z6,Z9,Z12,Z15-23:H) sex pheromone components in its sex communication, in which Z11-16:OAc is the major component and others are minor components. In this study, we for the first time demonstrated that the three PBPs differentiated in recognition among pheromone components in a moth using mixed-type sex pheromones. First, tissue expression study showed that all three PBPs of O. achatina were expressed only in antennae and highly male-biased, suggesting their involvement in perception of the sex pheromones. Second, the three PBPs were expressed in Escherichia coli and the binding affinities of PBPs to four sex pheromone components and some pheromone analogs were determined by the fluorescence competition binding assays. The results showed that OachPBP1 bound all four sex pheromone components with high binding affinity, while OachPBP2 had high or moderate binding affinity only to three Type I components, and OachPBP3 had high binding affinity only to three minor pheromone components. Furthermore, key amino acid residues that bind to sex pheromone components were identified in three PBPs by 3-D structure modeling and ligand molecular docking, predicting the interactions between PBPs and pheromone components. Our study provides a fundamental insight into the olfactory mechanism in moths that use mixed-type sex pheromones.

Computational Study on Temperature Driven Structure-Function Relationship of Polysaccharide Producing Bacterial Glycosyl Transferase Enzyme

Glycosyltransferase (GTs) is a wide class of enzymes that transfer sugar moiety, playing a key role in the synthesis of bacterial exopolysaccharide (EPS) biopolymer. In recent years, increased demand for bacterial EPSs has been observed in pharmaceutical, food, and other industries. The application of the EPSs largely depends upon their thermal stability, as any industrial application is mainly reliant on slow thermal degradation. Keeping this in context, EPS producing GT enzymes from three different bacterial sources based on growth temperature (mesophile, thermophile, and hyperthermophile) are considered for in silico analysis of the structural–functional relationship. From the present study, it was observed that the structural integrity of GT increases significantly from mesophile to thermophile to hyperthermophile. In contrast, the structural plasticity runs in an opposite direction towards mesophile. This interesting temperature-dependent structural property has directed the GT–UDP-glucose interactions in a way that thermophile has finally demonstrated better binding affinity (−5.57 to −10.70) with an increased number of hydrogen bonds (355) and stabilizing amino acids (Phe, Ala, Glu, Tyr, and Ser). The results from this study may direct utilization of thermophile-origin GT as best for industrial-level bacterial polysaccharide production.

In silico Structural, Functional and Phylogenetic Analyses of cellulase from Ruminococcus albus

Background

Cellulose is the primary component of the plant cell wall and an important source of energy for the ruminant and microbial protein synthesis in the rumen. Cell wall content is digested by anaerobic fermentation activity mainly of bacteria belonging to species Fibrobacter succinogenes, Ruminicoccus albus, Ruminococcus flavefaciens, and Butyrivibrio fibrisolvens. Bacteria belonging to the species Ruminococcus albus contain cellulosomes that enable it to adhere to and digest cellulose, and its genome encodes cellulases and hemicellulases.

This study aimed to perform an in silico comparative characterization and functional analysis of cellulase from Ruminococcus albus to explore physicochemical properties and to estimate primary, secondary, and tertiary structure using various bio-computational tools.

The protein sequences of cellulases belonging to 6 different Ruminococcus albus strains were retrieved using UniProt. In in silico composition of amino acids, basic physicochemical characteristics were analyzed using ProtParam and Protscale. Multiple sequence alignment of retrieved sequences was performed using Clustal Omega and the phylogenetic tree was constructed using Mega X software. Bioinformatics tools are used to better understand and determine the 3D structure of cellulase. The predicted model was refined by ModRefiner. Structure alignment between the best-predicted model and the template is applied to evaluate the similarity between structures.

Results

In this study are demonstrated several physicochemical characteristics of the cellulase enzyme. The instability index values indicate that the proteins are highly stable. Proteins are dominated by random coils and alpha helixes. The aliphatic index was higher than 71 providing information that the proteins are highly thermostable. No transmembrane domain was found in the protein, and the enzyme is extracellular and moderately acidic. The best tertiary structure model of the enzyme was obtained by the use of Raptor X, which was refined by ModRefiner. Raptor X suggested the 6Q1I_A as one of the best homologous templates for the predicted 3D protein structure. Ramachandran plot analysis showed that 90.1% of amino acid residues are within the most favored regions.

Conclusions

This study provides for the first time insights about the physicochemical properties, structure, and function of cellulase, from Ruminococcus albus, that will help for detection and identification of such enzyme in vivo or in silico.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-021-00162-x.