In Silico Structural Modeling and Analysis of Elongation Factor-1 Alpha and Elongation Factor-like Protein

In Silico Structural Modeling and Binding Site Analysis of Cerebroside Sulfotransferase (CST): A Therapeutic Target for Developing Substrate Reduction Therapy for Metachromatic Leukodystrophy

Cerebroside sulfotransferase (CST) is emerging as an important therapeutic target to develop substrate reduction therapy (SRT) for metachromatic leukodystrophy (MLD), a rare neurodegenerative lysosomal storage disorder. MLD develops with progressive impairment and destruction of the myelin sheath as a result of accumulation of sulfatide around the nerve cells in the absence of its recycling mechanism with deficiency of arylsulfatase A (ARSA). Sulfatide is the product of the catalytic action of cerebroside sulfotransferase (CST), which needs to be regulated under pathophysiological conditions by inhibitor development. To carry out in silico-based preliminary drug screening or for designing new drug candidates, a high-quality three-dimensional (3D) structure is needed in the absence of an experimentally derived three-dimensional crystal structure. In this study, a 3D model of the protein was developed using a primary sequence with the SWISS-MODEL server by applying the top four GMEQ score-based templates belonging to the sulfotransferase family as a reference. The 3D model of CST highlights the features of the protein responsible for its catalytic action. The CST model comprises five β-strands, which are flanked by ten α-helices from both sides as well as form the upside cover of the catalytic pocket of CST. CST has two catalytic regions: PAPS (-sulfo donor) binding and galactosylceramide (-sulfo acceptor) binding. The catalytic action of CST was proposed via molecular docking and molecular dynamic (MD) simulation with PAPS, galactosylceramide (GC), PAPS-galactosylceramide, and PAP. The stability of the model and its catalytic action were confirmed using molecular dynamic simulation-based trajectory analysis. CST response against the inhibition potential of the experimentally reported competitive inhibitor of CST was confirmed via molecular docking and molecular dynamics simulation, which suggested the suitability of the CST model for future drug discovery to strengthen substrate reduction therapy for MLD.

Proteomics analysis of resistance mechanism of Trichoderma harzianum under U(VI) stress

Trichoderma harzianum has a certain resistance to Hexavalent Uranium (U(VI)), but its resistance mechanism is unknown. Based on proteomics sequencing using DIA mode, differentially expressed proteins (DEPs) of Trichoderma harzianum under U(VI) stress were identified. GO enrichment, KEGG annotation analysis and DEPs annotation were performed. The results showed that 8 DEPs, 8 DEPs and 15 DEPs were obtained in the low-dose, medium-dose and high-dose groups, respectively. The functional classification of GO demonstrated that DEPs were associated with 17 molecular functions, 5 biological processes, and 5 cellular components. Furthermore, DEPs were enriched in transport and catabolism, energy metabolism, translation, and signal transduction. These findings showed that Trichoderma harzianum was significantly changed in protein expression and signaling pathway after U(VI) exposure. Therefore, these results have provided Trichoderma harzianum with a theoretical background that can be applied to environmental cleanup.

In-silico Structural Modeling of Human Immunodeficiency Virus Proteins

Human immunodeficiency virus (HIV) is an infectious virus that depletes the CD4⁺ T lymphocytes of the immune system and causes a chronic life-treating disease-acquired immunodeficiency syndrome (AIDS). The HIV genome encodes different structural and accessory proteins involved in viral entry and life cycle. Determining the 3D structure of HIV proteins is essential for new target position finding, structure-based drug designing, and future planning for computational and laboratory experimentations. Hence, the study aims to predict the 3D structures of all the HIV structural and accessory proteins using computational homology modeling to understand better the structural basis of HIV proteins interacting with host cells and viral replication. The sequences of HIV capsid, matrix, nucleocapsid, p6, reverse transcriptase, invertase, protease, gp120, gp41, virus protein r, viral infectivity factor, virus protein unique, RNA splicing regulator, transactivator protein, negative regulating factor, and virus protein x proteins were retrieved from UniProt. The primary and secondary structures of HIV proteins were predicted by Expasy ProtParam and SOPMA web servers. For the homology modeling, the MODELLER predicted the 3D structures of HIV proteins using templates. Then, the modeled structures were validated by the Ramachandran plot, local and global quality estimation scores, QMEAN scores, and Z-scores. Most of the amino acid residues of HIV proteins were present in the most favored and generously allowed regions in the Ramachandran plots. The local and global quality scores and Z-scores of the HIV proteins confirmed the good quality of modeled structures. The 3D modeled structures of HIV proteins might help further investigate the possible treatment.

Shape-based Machine Learning Models for the Potential Novel COVID-19 Protease Inhibitors Assisted by Molecular Dynamics Simulation

BACKGROUND

The vast geographical expansion of novel coronavirus and an increasing number of COVID-19 affected cases have overwhelmed health and public health services. Artificial Intelligence (AI) and Machine Learning (ML) algorithms have extended their major role in tracking disease patterns, and in identifying possible treatments.

OBJECTIVE

This study aims to identify potential COVID-19 protease inhibitors through shape-based Machine Learning assisted by Molecular Docking and Molecular Dynamics simulations.

METHODS

31 Repurposed compounds have been selected targeting the main coronavirus protease (6LU7) and a machine learning approach was employed to generate shape-based molecules starting from the 3D shape to the pharmacophoric features of their seed compound. Ligand-Receptor Docking was performed with Optimized Potential for Liquid Simulations (OPLS) algorithms to identify highaffinity compounds from the list of selected candidates for 6LU7, which were subjected to Molecular Dynamic Simulations followed by ADMET studies and other analyses.

RESULTS

Shape-based Machine learning reported remdesivir, valrubicin, aprepitant, and fulvestrant as the best therapeutic agents with the highest affinity for the target protein. Among the best shape-based compounds, a novel compound identified was not indexed in any chemical databases (PubChem, Zinc, or ChEMBL). Hence, the novel compound was named 'nCorv-EMBS'. Further, toxicity analysis showed nCorv-EMBS to be suitable for further consideration as the main protease inhibitor in COVID-19.

CONCLUSION

Effective ACE-II, GAK, AAK1, and protease 3C blockers can serve as a novel therapeutic approach to block the binding and attachment of the main COVID-19 protease (PDB ID: 6LU7) to the host cell and thus inhibit the infection at AT2 receptors in the lung. The novel compound nCorv- EMBS herein proposed stands as a promising inhibitor to be evaluated further for COVID-19 treatment.