Random sampling of the Protein Data Bank: RaSPDB

Interplay between hydrogen and chalcogen bonds in cysteine

Protein structures are stabilized by several types of chemical interactions between amino acids, which can compete with each other. This is the case of chalcogen and hydrogen bonds formed by the thiol group of cysteine, which can form three hydrogen bonds with one hydrogen acceptor and two hydrogen donors and a chalcogen bond with a nucleophile along the extension of the CS bond. A survey of the Protein Data Bank shows that hydrogen bonds are about 40-50 more common than chalcogen bonds, suggesting that they are stronger and, consequently, prevail, though not always. It is also observed that frequently a thiol group that forms a chalcogen bond is also involved, as a hydrogen donor, in a hydrogen bond.

Survey of the Intermolecular Disulfide Bonds Observed in Protein Crystal Structures Deposited in the Protein Data Bank

About 5% of the disulfide bonds (DBs) observed in the Protein Data Bank bridge two protein chains. Several of their features were comprehensively analyzed, resulting in a structural atlas of the intermolecular DBs. The analysis was performed on a very large set of data extracted from the Protein Data Bank, according to the RaSPDB procedure. It was observed that the two chains tend to have different sequences and belong to the same structural class. Intermolecular DBs tend to be more solvent accessible and less distorted from the most stable conformation than intermolecular DBs while showing similar B-factors. They tend to occur in beta strands and in mainly-beta structures. These and other data should prove useful in protein modelling and design.

Developing a bioinformatics pipeline for comparative protein classification analysis

BACKGROUND

Protein classification is a task of paramount importance in various fields of biology. Despite the great momentum of modern implementation of protein classification, machine learning techniques such as Random Forest and Neural Network could not always be used for several reasons: data collection, unbalanced classification or labelling of the data.As an alternative, I propose the use of a bioinformatics pipeline to search for and classify information from protein databases. Hence, to evaluate the efficiency and accuracy of the pipeline, I focused on the carotenoid biosynthetic genes and developed a filtering approach to retrieve orthologs clusters in two well-studied plants that belong to the Brassicaceae family: Arabidopsis thaliana and Brassica rapa Pekinensis group. The result obtained has been compared with previous studies on carotenoid biosynthetic genes in B. rapa where phylogenetic analysis was conducted.

RESULTS

The developed bioinformatics pipeline relies on commercial software and multiple databeses including the use of phylogeny, Gene Ontology terms (GOs) and Protein Families (Pfams) at a protein level. Furthermore, the phylogeny is coupled with "population analysis" to evaluate the potential orthologs. All the steps taken together give a final table of potential orthologs. The phylogenetic tree gives a result of 43 putative orthologs conserved in B. rapa Pekinensis group. Different A. thaliana proteins have more than one syntenic ortholog as also shown in a previous finding (Li et al., BMC Genomics 16(1):1-11, 2015).

CONCLUSIONS

This study demonstrates that, when the biological features of proteins of interest are not specific, I can rely on a computational approach in filtering steps for classification purposes. The comparison of the results obtained here for the carotenoid biosynthetic genes with previous research confirmed the accuracy of the developed pipeline which can therefore be applied for filtering different types of datasets.

Network Pharmacology and Molecular Docking to Elucidate the Potential Mechanism of Ligusticum Chuanxiong Against Osteoarthritis

Background: Osteoarthritis (OA) is a degenerative disease which serious affects patients. Ligusticum chuanxiong (CX) has been shown to have a certain curative effect on osteoarthritis in traditional Chinese medicine therapy. This study is based on network pharmacology and molecular docking technology to explore the potential mechanism of CX.

Methods: Components of CX to treat osteoarthritis were screened in the TCMSP database and targets were predicted by the PharmMapper database, the osteoarthritis targets were collected from the GeneCards database, and intersection genes were found to be the possible targets of CX anti-OA. The STRING database and Cytoscape software were utilized for protein-protein interaction analysis and further screening of core targets. The Metascape database was used for KEGG and GO enrichment analyses. Then, the top 10 pathways were selected to construct “drug-compound-target-pathway-disease” network analysis. Finally, molecular docking was used to analyze the binding affinity of seven compounds with core targets and TNF-α.

Results: Seven compounds with 253 non-repetitive targets of CX were screened from the TCMSP database and 60 potential intersection targets of CX anti-OA were found. PPI network analysis showed that the core targets were ALB, AKT1, IGF1, CASP3, MAPK1, ANXA5, and MAPK14, while GO and KEGG pathway enrichment analyses showed that the relevant biological processes involved in the treatment of osteoarthritis by CX might include the MAPK cascade and reactive oxygen species metabolic process. The KEGG pathway analysis result was mainly associated with the MAPK signaling pathway and PI3K-AKT signaling pathway. We further docked seven ingredients with MAPK1 and MAPK14 enriched in the MAPK pathway, and TNF-α as the typical inflammatory cytokine. The results also showed good binding affinity, especially FA, which may be the most important component of CX anti-OA.

Conclusion: Our research revealed the potential mechanism of CX in the treatment of OA, and our findings can also pave the way for subsequent basic experimental verification and a new research direction.