Multiple Structure Alignment Research Articles

Improving the inhibitory resistance of xylanase FgXyn11C from Fusarium graminearum to SyXIP-I by site-directed mutagenesis

The aim of this study was to improve the inhibitory resistance of xylanase FgXyn11C from Fusarium graminearum to XIP in cereal flour. Site saturation mutagenesis was performed using computer-aided redesign. Firstly, based on multiple primary structure alignments, the amino acid residues in the active site architecture were identified, and specific residue T144 in the thumb region of FgXyn11C was selected for site-saturation mutagenesis. After screening, FgXyn11CT144F was selected as the best mutant, as it displayed the highest enzymatic activity and resistance simultaneously compared to other mutants. The specific activity of FgXyn11CT144F was 208.8 U/mg and it exhibited complete resistance to SyXIP-I. Compared with the wild-type, FgXyn11CT144F displayed similar activity and the most resistant against SyXIP-I. The optimal temperature and pH of the wild-type and purified FgXyn11CT144F were similar at pH 5.0 and 30 °C. Our findings provided preliminary insight into how the specific residue at position 144 in the thumb region of FgXyn11C influenced the enzymatic properties and interacted with SyXIP-I. The inhibition sensitivity of FgXyn11C was reduced through directed evolution, leading to creation of the mutant enzyme FgXyn11CT144F. The FgXyn11CT144F resistance to SyXIP-I has potential application and can also provide references for engineering other resistant xylanases of the GHF11.

The Complex GH32 Enzyme Orchestra from Priestia megaterium Holds the Key to Better Discriminate Sucrose-6-phosphate Hydrolases from Other β-Fructofuranosidases in Bacteria.

Inulin is widely used as a prebiotic and emerging as a priming compound to counteract plant diseases. We isolated inulin-degrading strains from the lettuce phyllosphere, identified as Bacillus subtilis and Priestia megaterium, species hosting well-known biocontrol organisms. To better understand their varying inulin degradation strategies, three intracellular β-fructofuranosidases from P. megaterium NBRC15308 were characterized after expression in Escherichia coli: a predicted sucrose-6-phosphate (Suc6P) hydrolase (SacAP1, supported by molecular docking), an exofructanase (SacAP2), and an invertase (SacAP3). Based on protein multiple sequence and structure alignments of bacterial glycoside hydrolase family 32 enzymes, we identified conserved residues predicted to be involved in binding phosphorylated (Suc6P hydrolases) or nonphosphorylated substrates (invertases and fructanases). Suc6P hydrolases feature positively charged residues near the structural catalytic pocket (histidine, arginine, or lysine), whereas other β-fructofuranosidases contain tryptophans. This correlates with our phylogenetic tree, grouping all predicted Suc6P hydrolases in a clan associated with genomic regions coding for transporters involved in substrate phosphorylation. These results will help to discriminate between Suc6P hydrolases and other β-fructofuranosidases in future studies and to better understand the interaction of B. subtilis and P. megaterium endophytes with sucrose and/or fructans, sugars naturally present in plants or exogenously applied in the context of defense priming.

RCSB Protein Data Bank: visualizing groups of experimentally determined PDB structures alongside computed structure models of proteins.

Recent advances in Artificial Intelligence and Machine Learning (e.g., AlphaFold, RosettaFold, and ESMFold) enable prediction of three-dimensional (3D) protein structures from amino acid sequences alone at accuracies comparable to lower-resolution experimental methods. These tools have been employed to predict structures across entire proteomes and the results of large-scale metagenomic sequence studies, yielding an exponential increase in available biomolecular 3D structural information. Given the enormous volume of this newly computed biostructure data, there is an urgent need for robust tools to manage, search, cluster, and visualize large collections of structures. Equally important is the capability to efficiently summarize and visualize metadata, biological/biochemical annotations, and structural features, particularly when working with vast numbers of protein structures of both experimental origin from the Protein Data Bank (PDB) and computationally-predicted models. Moreover, researchers require advanced visualization techniques that support interactive exploration of multiple sequences and structural alignments. This paper introduces a suite of tools provided on the RCSB PDB research-focused web portal RCSB. org, tailor-made for efficient management, search, organization, and visualization of this burgeoning corpus of 3D macromolecular structure data.

PC_ali: A tool for improved multiple alignments and evolutionary inference based on a hybrid protein sequence and structure similarity score.

Evolutionary inference depends crucially on the quality of multiple sequence alignments (MSA), which is problematic for distantly related proteins. Since protein structure is more conserved than sequence, it seems natural to use structure alignments for distant homologs. However, structure alignments may not be suitable for inferring evolutionary relationships. Here we examined four protein similarity measures that depend on sequence and structure (fraction of aligned residues, sequence identity, fraction of superimposed residues and contact overlap), finding that they are intimately correlated but none of them provides a complete and unbiased picture of conservation in proteins. Therefore, we propose the new hybrid protein sequence and structure similarity score PC_sim based on their main Principal Component. The corresponding divergence measure PC_div shows the strongest correlation with divergences obtained from individual similarities, suggesting that it infers accurate evolutionary divergences. We developed the program PC_ali that constructs protein MSAs either de novo or modifying an input MSA, using a similarity matrix based on PC_sim. tructs a starting MSA based on the maximal cliques of the graph of these PAs and it refines it through progressive alignments along the tree reconstructed with PC_div. Compared with eight state-of-the-art multiple structure or sequence alignment tools, PC_ali achieve higher or equal aligned fraction and structural scores, sequence identity higher than structure aligners although lower than sequence aligners, highest score PC_sim and highest similarity with the MSAs produced by other tools and with the reference MSA Balibase. https://github.com/ugobas/PC_ali. Supplementary data are available at Bioinformatics online.

Structure-guided engineering for improving the thermal stability of zearalenone hydrolase

Zearalenone is one of the most widely polluted Fusarium toxins in the world, seriously endangering livestock and human health. Zearalenone hydrolase (ZHD) derived from Clonostachys rosea can effectively degrade zearalenone. However, the high temperature environment in feed processing hampers the application of this enzyme. Structure-based rational design may provide guidance for engineering the thermal stability of enzymes. In this paper, we used the multiple structure alignment (MSTA) to screen the structural flexibility regions of ZHD. Subsequently, a candidate mutation library was constructed by sequence conservation scoring and conformational free energy calculation, from which 9 single point mutations based on residues 136 and 220 were obtained. The experiments showed that the thermal melting temperature (Tm) of the 9 mutants increased by 0.4-5.6 ℃. The S220R and S220W mutants showed the best thermal stability, the Tm of which increased by 5.6 ℃ and 4.0 ℃ compared to that of the wild type. Moreover, the thermal half-inactivation time at 45 ℃ were 15.4 times and 3.1 times longer, and the relative activities were 70.6% and 57.3% of the wild type. Molecular dynamics simulation analysis showed that the interaction force at and around the mutation site was enhanced, contributing to the improved thermal stability of ZHD. The probability of 220-K130 hydrogen bond of the mutants S220R and S220W increased by 37.1% and 19.3%, and the probability of K130-D223 salt bridge increased by 30.1% and 12.5%, respectively. This work demonstrated the feasibility of thermal stability engineering strategy where the structural and sequence alignment as well as free energy calculation of natural enzymes were integrated, and obtained ZHD variants with enhanced thermal stability, which may facilitate the industrial application of ZHD.

Structural Homology-Based Drug Repurposing Approach for Targeting NSP12 SARS-CoV-2.

The severe acute respiratory syndrome coronavirus 2, also known as SARS-CoV-2, is the causative agent of the COVID-19 global pandemic. SARS-CoV-2 has a highly conserved non-structural protein 12 (NSP-12) involved in RNA-dependent RNA polymerase (RdRp) activity. For the identification of potential inhibitors for NSP-12, computational approaches such as the identification of homologous proteins that have been previously targeted by FDA-approved antivirals can be employed. Herein, homologous proteins of NSP-12 were retrieved from Protein DataBank (PDB) and the evolutionary conserved sequence and structure similarity of the active site of the RdRp domain of NSP-12 was characterized. The identified homologous structures of NSP-12 belonged to four viral families: Coronaviridae, Flaviviridae, Picornaviridae, and Caliciviridae, and shared evolutionary conserved relationships. The multiple sequences and structural alignment of homologous structures showed highly conserved amino acid residues that were located at the active site of the RdRp domain of NSP-12. The conserved active site of the RdRp domain of NSP-12 was evaluated for binding affinity with the FDA-approved antivirals, i.e., Sofosbuvir and Dasabuvir in a molecular docking study. The molecular docking of Sofosbuvir and Dasabuvir with the active site that contains conserved motifs (motif A-G) of the RdRp domain of NSP-12 revealed significant binding affinity. Furthermore, MD simulation also inferred the potency of Sofosbuvir and Dasabuvir. In conclusion, targeting the active site of the RdRp domain of NSP-12 with Dasabuvir and Sofosbuvir might reduce viral replication and pathogenicity and could be further studied for the treatment of SARS-CoV-2.

US-align: universal structure alignments of proteins, nucleic acids, and macromolecular complexes.

Structure comparison and alignment are of fundamental importance in structural biology studies. We developed the first universal platform, US-align, to uniformly align monomer and complex structures of different macromolecules-proteins, RNAs and DNAs. The pipeline is built on a uniform TM-score objective function coupled with a heuristic alignment searching algorithm. Large-scale benchmarks demonstrated consistent advantages of US-align over state-of-the-art methods in pairwise and multiple structure alignments of different molecules. Detailed analyses showed that the main advantage of US-align lies in the extensive optimization of the unified objective function powered by efficient heuristic search iterations, which substantially improve the accuracy and speed of the structural alignment process. Meanwhile, the universal protocol fusing different molecular and structural types helps facilitate the heterogeneous oligomer structure comparison and template-based protein-protein and protein-RNA/DNA docking.

Multiple alignment of structures using center of proteins

Multiple structure alignment (MStA) is a fundamental tool for correlating the structural similarity of proteins with their functional similarity. A number of algorithms are extant, MUSTANG, POSA, M...

PmTM-align: scalable pairwise and multiple structure alignment with Apache Spark and OpenMP

BackgroundStructure comparison can provide useful information to identify functional and evolutionary relationship between proteins. With the dramatic increase of protein structure data in the Protein Data Bank, computation time quickly becomes the bottleneck for large scale structure comparisons. To more efficiently deal with informative multiple structure alignment tasks, we propose pmTM-align, a parallel protein structure alignment approach based on mTM-align/TM-align. pmTM-align contains two stages to handle pairwise structure alignments with Spark and the phylogenetic tree-based multiple structure alignment task on a single computer with OpenMP.ResultsExperiments with the SABmark dataset showed that parallelization along with data structure optimization provided considerable speedup for mTM-align. The Spark-based structure alignments achieved near ideal scalability with large datasets, and the OpenMP-based construction of the phylogenetic tree accelerated the incremental alignment of multiple structures and metrics computation by a factor of about 2–5.ConclusionspmTM-align enables scalable pairwise and multiple structure alignment computing and offers more timely responses for medium to large-sized input data than existing alignment tools such as mTM-align.

In Vitro and in Vivo Activity of mTOR Kinase and PI3K Inhibitors Against Leishmania donovani and Trypanosoma brucei.

Kinetoplastid parasites, including Leishmania and Trypanosoma spp., are life threatening pathogens with a worldwide distribution. Next-generation therapeutics for treatment are needed as current treatments have limitations, such as toxicity and drug resistance. In this study, we examined the activities of established mammalian target of rapamycin (mTOR)/phosphoinositide 3-kinase (PI3K) inhibitors against these tropical diseases. High-throughput screening of a library of 1742 bioactive compounds against intracellular L. donovani was performed, and seven mTOR/PI3K inhibitors were identified. Dose-dilution assays revealed that these inhibitors had half maximal effective concentration (EC50) values ranging from 0.14 to 13.44 μM for L. donovani amastigotes and from 0.00005 to 8.16 μM for T. brucei. The results of a visceral leishmaniasis mouse model indicated that treatment with Torin2, dactolisib, or NVP-BGT226 resulted in reductions of 35%, 53%, and 54%, respectively, in the numbers of liver parasites. In an acute T. brucei mouse model using NVP-BGT226 parasite numbers were reduced to under the limits of detection by five consecutive days of treatment. Multiple sequence and structural alignment results indicated high similarities between mTOR and kinetoplastid TORs; the inhibitors are predicted to bind in a similar manner. Taken together, these results indicated that the TOR pathways of parasites have potential for the discovery of novel targets and new potent inhibitors.

Insights about the structure of farnesyl diphosphate synthase (FPPS) and the activity of bisphosphonates on the proliferation and ultrastructure of Leishmania and Giardia

BackgroundThe enzyme farnesyl diphosphate synthase (FPPS) is positioned in the intersection of different sterol biosynthesis pathways such as those producing isoprenoids, dolichols and ergosterol. FPPS is ubiquitous in eukaryotes and is inhibited by nitrogen-containing bisphosphonates (N-BP). N-BP activity and the mechanisms of cell death as well as damage to the ultrastructure due to N-BP has not yet been investigated in Leishmania infantum and Giardia. Thus, we evaluated the effect of N-BP on cell viability and ultrastructure and then performed structural modelling and phylogenetic analysis on the FPPS enzymes of Leishmania and Giardia.MethodsWe performed multiple sequence alignment with MAFFT, phylogenetic analysis with MEGA7, and 3D structural modelling for FPPS with Modeller 9.18 and on I-Tasser server. We performed concentration curves with N-BP in Leishmania promastigotes and Giardia trophozoites to estimate the IC50via the MTS/PMS viability method. The ultrastructure was evaluated by transmission electron microscopy, and the mechanism of cell death by flow cytometry.ResultsThe nitrogen-containing bisphosphonate risedronate had stronger anti-proliferative activity in Leishmania compared to other N-BPs with an IC50 of 13.8 µM, followed by ibandronate and alendronate with IC50 values of 85.1 µM and 112.2 µM, respectively. The effect of N-BPs was much lower on trophozoites of Giardia than Leishmania (IC50 of 311 µM for risedronate). Giardia treated with N-BP displayed concentric membranes around the nucleus and nuclear pyknosis. Leishmania had mitochondrial swelling, myelin figures, double membranes, and plasma membrane blebbing. The same population labelled with annexin-V and 7-AAD had a loss of membrane potential (TMRE), indicative of apoptosis. Multiple sequence alignments and structural alignments of FPPS proteins showed that Giardia and Leishmania FPPS display low amino acid identity but possess the conserved aspartate-rich motifs.ConclusionsGiardia and Leishmania FPPS enzymes are phylogenetically distant but display conserved protein signatures. The N-BPs effect on FPPS was more pronounced in Leishmania than Giardia. This might be due to general differences in metabolism and differences in the FPPS catalytic site.

An Implementation of Three-level Multi-objective ABC Algorithm for RNA Multiple Structural Alignment

Background: Structural alignment of ribose nucleic acid (RNA) is one of the most challenging multi-objective real world problems from the field of bioinformatics. Objective: RNA molecules are less stable; hence they require inclusion of most stable secondary structure during their alignment. Therefore, the structural alignment requires the consideration of similarity score and structure score, as two objectives. Trade-off between these two objectives exists since obtaining optimum similarity score at concurrent optimum structure score is not possible. This paper presents artificial bee colony algorithm based three level multi-objective approach for performing structural alignment of RNA sequences, namely MO-3LABC. Methods: Algorithm firstly builds the secondary structure of all sequences at minimum free energy (MFE). Then sequence alignment is performed in level one at average percent sequence identity (APSI) score based gap length, optimized by ABC algorithm. Level two now builds the secondary structure of these aligned sequences based on base-pair probability and co-variation. Now the scores of level one and level two move towards level three for multi-objective optimization at Pareto optimality criteria with few additional strategies. Results: The results of MO-3LABC are compared with an already established efficient strategy MO-TLPSO; multi-objective two level strategy based on particle swarm optimization. The outputs are compared for pairwise and multiple sequence alignment datasets at prediction accuracy and solution quality criteria. Conclusion: MO-3LABC is found significantly better than MO-TLPSO at all the four evaluation criteria for both the datasets.

Dimer Interface Organization is a Main Determinant of Intermonomeric Interactions and Correlates with Evolutionary Relationships of Retroviral and Retroviral-Like Ddi1 and Ddi2 Proteases.

The life cycles of retroviruses rely on the limited proteolysis catalyzed by the viral protease. Numerous eukaryotic organisms also express endogenously such proteases, which originate from retrotransposons or retroviruses, including DNA damage-inducible 1 and 2 (Ddi1 and Ddi2, respectively) proteins. In this study, we performed a comparative analysis based on the structural data currently available in Protein Data Bank (PDB) and Structural summaries of PDB entries (PDBsum) databases, with a special emphasis on the regions involved in dimerization of retroviral and retroviral-like Ddi proteases. In addition to Ddi1 and Ddi2, at least one member of all seven genera of the Retroviridae family was included in this comparison. We found that the studied retroviral and non-viral proteases show differences in the mode of dimerization and density of intermonomeric contacts, and distribution of the structural characteristics is in agreement with their evolutionary relationships. Multiple sequence and structure alignments revealed that the interactions between the subunits depend mainly on the overall organization of the dimer interface. We think that better understanding of the general and specific features of proteases may support the characterization of retroviral-like proteases.

Caretta – A multiple protein structure alignment and feature extraction suite

The vast number of protein structures currently available opens exciting opportunities for machine learning on proteins, aimed at predicting and understanding functional properties. In particular, in combination with homology modelling, it is now possible to not only use sequence features as input for machine learning, but also structure features. However, in order to do so, robust multiple structure alignments are imperative.Here we present Caretta, a multiple structure alignment suite meant for homologous but sequentially divergent protein families which consistently returns accurate alignments with a higher coverage than current state-of-the-art tools. Caretta is available as a GUI and command-line application and additionally outputs an aligned structure feature matrix for a given set of input structures, which can readily be used in downstream steps for supervised or unsupervised machine learning. We show Caretta’s performance on two benchmark datasets, and present an example application of Caretta in predicting the conformational state of cyclin-dependent kinases.

ParMATT: parallel multiple alignment of protein 3D-structures with translations and twists for distributed-memory systems.

Accurate structural alignment of proteins is crucial at studying structure-function relationship in evolutionarily distant homologues. Various software tools were proposed to align multiple protein 3D-structures utilizing one CPU and thus are of limited productivity at large-scale analysis of protein families/superfamilies. The parMATT is a hybrid MPI/pthreads/OpenMP parallel re-implementation of the MATT algorithm to align multiple protein 3D-structures by allowing translations and twists. The parMATT can be faster than MATT on a single multi-core CPU, and provides a much greater speedup when executed on distributed-memory systems, i.e. computing clusters and supercomputers hosting memory-independent computing nodes. The most computationally demanding steps of the MATT algorithm-the initial construction of pairwise alignments between all input structures and further iterative progression of the multiple alignment-were parallelized using MPI and pthreads, and the concluding refinement step was optimized by introducing the OpenMP support. The parMATT can significantly accelerate the time-consuming process of building a multiple structural alignment from a large set of 3D-records of homologous proteins. The source code is available at https://biokinet.belozersky.msu.ru/parMATT. Supplementary data are available at Bioinformatics online.

In silico comparative study of structural homology modeling on envelope (E) glycoprotein of dengue viruses (Denv)

Dengue infection has become global epidemic health problem. The specific drugs for treatment and effective vaccination for prevention of dengue infection against four types of DENV serotypes are currently unknown. The envelope (E) glycoprotein is responsible in viral attachment to cell surface receptors. Homology modelling on E glycoprotein of dengue viruses (DENV) Type 1–4 was constructed using SWISS-MODEL database. The target sequence retrieved from Clustal Omega. The homology model of three dimensional structures was evaluated and checked for quality estimation using PROCHECK, VERIFY 3D and ERRAT server. Multiple structure alignment was analysed using PDBeFold. DENV1 has the highest quality structure while DEN2 has the lowest quality structure. Multiple structural alignments showed that DENV2 and DENV3 have the highest sequence similarity. The structural alignment showed that DENV2 and DENV3 had the lowest RMSD score (1.325), and the highest Q score (0.4433). The homology model structure of DENV2 was the most similar and closely related to the other three homology model structures of DENV serotypes. In conclusion, the homology model of three-dimensional structures on E glycoprotein of each DENV serotype was successfully obtained and can be used for further investigations on E glycoprotein of dengue viruses (DENV) Type 1–4.

MTM-align: a server for fast protein structure database search and multiple protein structure alignment.

With the rapid increase of the number of protein structures in the Protein Data Bank, it becomes urgent to develop algorithms for efficient protein structure comparisons. In this article, we present the mTM-align server, which consists of two closely related modules: one for structure database search and the other for multiple structure alignment. The database search is speeded up based on a heuristic algorithm and a hierarchical organization of the structures in the database. The multiple structure alignment is performed using the recently developed algorithm mTM-align. Benchmark tests demonstrate that our algorithms outperform other peering methods for both modules, in terms of speed and accuracy. One of the unique features for the server is the interplay between database search and multiple structure alignment. The server provides service not only for performing fast database search, but also for making accurate multiple structure alignment with the structures found by the search. For the database search, it takes about 2–5 min for a structure of a medium size (∼300 residues). For the multiple structure alignment, it takes a few seconds for ∼10 structures of medium sizes. The server is freely available at: http://yanglab.nankai.edu.cn/mTM-align/.

The Interplay between Structural Stability and Plasticity Determines Mutation Profiles and Chaperone Dependence in Protein Kinases.

We present a novel comparative analysis of representative protein kinases to characterize the main dynamic and energetic determinants of functional regulation shared among different families. The relationships between stability and plasticity are also used to rationalize kinase tendencies to interact with the molecular chaperone Hsp90. These questions are tackled through newly developed molecular-dynamics-based methods of analysis of internal energy and dynamics applied to a total of 37 different systems, which represent wild-type and mutated proteins, including active and inactive states. Energetic decomposition analysis is coupled to multiple structural alignments and dynamic decomposition methods and identifies, across different families, common elements that underlie fold stabilization and conformational regulation. This analysis also exposes which substructures play a key role in determining chaperone dependence. Overall, the results highlight common interaction networks that underpin kinase stabilization, are modulated by mutations (even if located at a distance), and underlie their tendencies to act as clients or nonclients of Hsp90.

MTM-align: an algorithm for fast and accurate multiple protein structure alignment.

As protein structure is more conserved than sequence during evolution, multiple structure alignment can be more informative than multiple sequence alignment, especially for distantly related proteins. With the rapid increase of the number of protein structures in the Protein Data Bank, it becomes urgent to develop efficient algorithms for multiple structure alignment. A new multiple structure alignment algorithm (mTM-align) was proposed, which is an extension of the highly efficient pairwise structure alignment program TM-align. The algorithm was benchmarked on four widely used datasets, HOMSTRAD, SABmark_sup, SABmark_twi and SISY-multiple, showing that mTM-align consistently outperforms other algorithms. In addition, the comparison with the manually curated alignments in the HOMSTRAD database shows that the automated alignments built by mTM-align are in general more accurate. Therefore, mTM-align may be used as a reliable complement to construct multiple structure alignments for real-world applications. http://yanglab.nankai.edu.cn/mTM-align. zhng@umich.edu or yangjy@nankai.edu.cn. Supplementary data are available at Bioinformatics online.

RNAscClust: clustering RNA sequences using structure conservation and graph based motifs.

MotivationClustering RNA sequences with common secondary structure is an essential step towards studying RNA function. Whereas structural RNA alignment strategies typically identify common structure for orthologous structured RNAs, clustering seeks to group paralogous RNAs based on structural similarities. However, existing approaches for clustering paralogous RNAs, do not take the compensatory base pair changes obtained from structure conservation in orthologous sequences into account.ResultsHere, we present RNAscClust, the implementation of a new algorithm to cluster a set of structured RNAs taking their respective structural conservation into account. For a set of multiple structural alignments of RNA sequences, each containing a paralog sequence included in a structural alignment of its orthologs, RNAscClust computes minimum free-energy structures for each sequence using conserved base pairs as prior information for the folding. The paralogs are then clustered using a graph kernel-based strategy, which identifies common structural features. We show that the clustering accuracy clearly benefits from an increasing degree of compensatory base pair changes in the alignments.Availability and Implementation RNAscClust is available at http://www.bioinf.uni-freiburg.de/Software/RNAscClust.Supplementary information Supplementary data are available at Bioinformatics online.