Structural and quantum mechanical computations to elucidate the altered binding mechanism of metal and drug with pyrazinamidase from Mycobacterium tuberculosis due to mutagenicity

Hemolytic-Pred: A machine learning-based predictor for hemolytic proteins using position and composition-based features

Objective

The objective of this study is to propose a novel in-silico method called Hemolytic-Pred for identifying hemolytic proteins based on their sequences, using statistical moment-based features, along with position-relative and frequency-relative information.

Methods

Primary sequences were transformed into feature vectors using statistical and position-relative moment-based features. Varying machine learning algorithms were employed for classification. Computational models were rigorously evaluated using four different validation. The Hemolytic-Pred webserver is available for further analysis at http://ec2-54-160-229-10.compute-1.amazonaws.com/.

Results

XGBoost outperformed the other six classifiers with an accuracy value of 0.99, 0.98, 0.97, and 0.98 for self-consistency test, 10-fold cross-validation, Jackknife test, and independent set test, respectively. The proposed method with the XGBoost classifier is a workable and robust solution for predicting hemolytic proteins efficiently and accurately.

Conclusions

The proposed method of Hemolytic-Pred with XGBoost classifier is a reliable tool for the timely identification of hemolytic cells and diagnosis of various related severe disorders. The application of Hemolytic-Pred can yield profound benefits in the medical field.

iGluK-Deep: computational identification of lysine glutarylation sites using deep neural networks with general pseudo amino acid compositions

Lysine glutarylation is a post-translation modification which plays an important regulatory role in a variety of physiological and enzymatic processes including mitochondrial functions and metabolic processes both in eukaryotic and prokaryotic cells. This post-translational modification influences chromatin structure and thereby results in global regulation of transcription, defects in cell-cycle progression, DNA damage repair, and telomere silencing. To better understand the mechanism of lysine glutarylation, its identification in a protein is necessary, however, experimental methods are time-consuming and labor-intensive. Herein, we propose a new computational prediction approach to supplement experimental methods for identification of lysine glutarylation site prediction by deep neural networks and Chou's Pseudo Amino Acid Composition (PseAAC). We employed well-known deep neural networks for feature representation learning and classification of peptide sequences. Our approach opts raw pseudo amino acid compositions and obsoletes the need to separately perform costly and cumbersome feature extraction and selection. Among the developed deep learning-based predictors, the standard neural network-based predictor demonstrated highest scores in terms of accuracy and all other performance evaluation measures and outperforms majority of previously reported predictors without requiring expensive feature extraction process. iGluK-Deep:Computational Identification of lysine glutarylationsites using deep neural networks with general Pseudo Amino Acid Compositions Sheraz Naseer, Rao Faizan Ali, Yaser Daanial Khan, P.D.D DominicCommunicated by Ramaswamy H. Sarma.

Hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) Simulation: A Tool for Structure-based Drug Design and Discovery

Quantum mechanics (QM) is physics based theory which explains the physical properties of nature at the level of atoms and sub-atoms. Molecular mechanics (MM) construct molecular systems through the use of classical mechanics. So, hybrid quantum mechanics and molecular mechanics (QM/MM) when combined together can act as computer-based methods which can be used to calculate structure and property data of molecular structures. Hybrid QM/MM combines the strengths of QM with accuracy and MM with speed. QM/MM simulation can also be applied for the study of chemical process in solutions as well as in the proteins, and has a great scope in structure-based drug design (CADD) and discovery. Hybrid QM/MM also applied to HTS, to derive QSAR models and due to availability of many protein crystal structures; it has a great role in computational chemistry, especially in structure- and fragment-based drug design. Fused QM/MM simulations have been developed as a widespread method to explore chemical reactions in condensed phases. In QM/MM simulations, the quantum chemistry theory is used to treat the space in which the chemical reactions occur; however the rest is defined through molecular mechanics force field (MMFF). In this review, we have extensively reviewed recent literature pertaining to the use and applications of hybrid QM/MM simulations for ligand and structure-based computational methods for the design and discovery of therapeutic agents.

Using CHOU'S 5-Steps Rule to Predict O-Linked Serine Glycosylation Sites by Blending Position Relative Features and Statistical Moment

Glycosylation of proteins in eukaryote cells is an important and complicated post-translation modification due to its pivotal role and association with crucial physiological functions within most of the proteins. Identification of glycosylation sites in a polypeptide chain is not an easy task due to multiple impediments. Analytical identification of these sites is expensive and laborious. There is a dire need to develop a reliable computational method for precise determination of such sites which can help researchers to save time and effort. Herein, we propose a novel predictor namely iGlycoS-PseAAC by integrating the Chou's Pseudo Amino Acid Composition (PseAAC) and relative/absolute position-based features. The self-consistency results show that the accuracy revealed by the model using the benchmark dataset for prediction of O-linked glycosylation having serine sites is 98.8 percent. The overall accuracy of predictor achieved through 10-fold cross validation by combining the positive and negative results is 97.2 percent. The overall accuracy achieved through Jackknife test is 96.195 percent by aggregating of all the prediction results. Thus the proposed predictor can help in predicting the O-linked glycosylated serine sites in an efficient and accurate way. The overall results show that the accuracy of the iGlycoS-PseAAC is higher than the existing tools.

Virtual Screening of Phytochemicals by Targeting HR1 Domain of SARS-CoV-2 S Protein: Molecular Docking, Molecular Dynamics Simulations, and DFT Studies

The recent COVID-19 pandemic has impacted nearly the whole world due to its high morbidity and mortality rate. Thus, scientists around the globe are working to find potent drugs and designing an effective vaccine against COVID-19. Phytochemicals from medicinal plants are known to have a long history for the treatment of various pathogens and infections; thus, keeping this in mind, this study was performed to explore the potential of different phytochemicals as candidate inhibitors of the HR1 domain in SARS-CoV-2 spike protein by using computer-aided drug discovery methods. Initially, the pharmacological assessment was performed to study the drug-likeness properties of the phytochemicals for their safe human administration. Suitable compounds were subjected to molecular docking to screen strongly binding phytochemicals with HR1 while the stability of ligand binding was analyzed using molecular dynamics simulations. Quantum computation-based density functional theory (DFT) analysis was constituted to analyze the reactivity of these compounds with the receptor. Through analysis, 108 phytochemicals passed the pharmacological assessment and upon docking of these 108 phytochemicals, 36 were screened passing a threshold of -8.5 kcal/mol. After analyzing stability and reactivity, 5 phytochemicals, i.e., SilybinC, Isopomiferin, Lycopene, SilydianinB, and Silydianin are identified as novel and potent candidates for the inhibition of HR1 domain in SARS-CoV-2 spike protein. Based on these results, it is concluded that these compounds can play an important role in the design and development of a drug against COVID-19, after an exhaustive in vitro and in vivo examination of these compounds, in future.

iPhosH-PseAAC: Identify Phosphohistidine Sites in Proteins by Blending Statistical Moments and Position Relative Features According to the Chou's 5-Step Rule and General Pseudo Amino Acid Composition

Protein phosphorylation is one of the key mechanism in prokaryotes and eukaryotes and is responsible for various biological functions such as protein degradation, intracellular localization, the multitude of cellular processes, molecular association, cytoskeletal dynamics, and enzymatic inhibition/activation. Phosphohistidine (PhosH) has a key role in a number of biological processes, including central metabolism to signalling in eukaryotes and bacteria. Thus, identification of phosphohistidine sites in a protein sequence is crucial, and experimental identification can be expensive, time-taking, and laborious. To address this problem, here, we propose a novel computational model namely iPhosH-PseAAC for prediction of phosphohistidine sites in a given protein sequence using pseudo amino acid composition (PseAAC), statistical moments, and position relative features. The results of the proposed predictor are validated through self-consistency testing, 10-fold cross-validation, and jackknife testing. The self-consistency validation gave the 100 percent accuracy, whereas, for cross-validation, the accuracy achieved is 94.26 percent. Moreover, jackknife testing gave 97.07 percent accuracy for the proposed model. Thus, the proposed model iPhosH-PseAAC for prediction of iPhosH site has the great ability to predict the PhosH sites in given proteins.

Phytochemicals from Selective Plants Have Promising Potential against SARS-CoV-2: Investigation and Corroboration through Molecular Docking, MD Simulations, and Quantum Computations

Coronaviruses have been reported previously due to their association with the severe acute respiratory syndrome (SARS). After SARS, these viruses were known to be causing Middle East respiratory syndrome (MERS) and caused 35% evanescence amid victims pursuing remedial care. Nowadays, beta coronaviruses, members of Coronaviridae, family order Nidovirales, have become subjects of great importance due to their latest pandemic originating from Wuhan, China. The virus named as human-SARS-like coronavirus-2 contains four structural as well as sixteen nonstructural proteins encoded by single-stranded ribonucleic acid of positive polarity. As there is no vaccine available to treat the infection caused by these viruses, there is a dire need for taking necessary steps against this virus. Herein, we have targeted two nonstructural proteins of SARS-CoV-2, namely, methyltransferase (nsp16) and helicase (nsp13), respectively, due to their substantial activity in viral pathogenesis. A total of 2035 compounds were analyzed for their pharmacokinetics and pharmacological properties. The screened 108 compounds were docked against both targeted proteins and were compared with previously reported known compounds. Compounds with high binding affinity were analyzed for their reactivity through DFT analysis, and binding was analyzed using molecular dynamics simulations. Through the analyses performed in this study, it is concluded that EryvarinM, Silydianin, Osajin, and Raddeanine can be considered potential inhibitors for MTase, while TomentodiplaconeB, Osajin, Sesquiterpene Glycoside, Rhamnetin, and Silydianin for helicase after these compounds are validated thoroughly using in vitro and in vivo protocols.

Insights into Machine Learning-based Approaches for Virtual Screening in Drug Discovery: Existing Strategies and Streamlining Through FP-CADD

BACKGROUND

Machine learning is an active area of research in computer science by the availability of big data collection of all sorts prompting interest in the development of novel tools for data mining. Machine learning methods have wide applications in computer-aided drug discovery methods. Most incredible approaches to machine learning are used in drug designing, which further aid the process of biological modelling in drug discovery. Mainly, two main categories are present which are Ligand-Based Virtual Screening (LBVS) and Structure-Based Virtual Screening (SBVS), however, the machine learning approaches fall mostly in the category of LBVS.

OBJECTIVES

This study exposits the major machine learning approaches being used in LBVS. Moreover, we have introduced a protocol named FP-CADD which depicts a 4-steps rule of thumb for drug discovery, the four protocols of computer-aided drug discovery (FP-CADD). Various important aspects along with SWOT analysis of FP-CADD are also discussed in this article.

CONCLUSION

By this thorough study, we have observed that in LBVS algorithms, Support Vector Machines (SVM) and Random Forest (RF) are those which are widely used due to high accuracy and efficiency. These virtual screening approaches have the potential to revolutionize the drug designing field. Also, we believe that the process flow presented in this study, named FP-CADD, can streamline the whole process of computer-aided drug discovery. By adopting this rule, the studies related to drug discovery can be made homogeneous and this protocol can also be considered as an evaluation criterion in the peer-review process of research articles.

iProtease-PseAAC(2L): A two-layer predictor for identifying proteases and their types using Chou's 5-step-rule and general PseAAC

Proteases are a type of enzymes, which perform the process of proteolysis. Proteolysis normally refers to protein and peptide degradation which is crucial for the survival, growth and wellbeing of a cell. Moreover, proteases have a strong association with therapeutics and drug development. The proteases are classified into five different types according to their nature and physiochemical characteristics. Mostly the methods used to differentiate protease from other proteins and identify their class requires a clinical test which is usually time-consuming and operator dependent. Herein, we report a classifier named iProtease-PseAAC (2L) for identifying proteases and their classes. The predictor is developed employing the flow of 5-step rule, initiating from the collection of benchmark dataset and terminating at the development of predictor. Rigorous verification and validation tests are performed and metrics are collected to calculate the authenticity of the trained model. The self-consistency validation gives the 98.32% accuracy, for cross-validation the accuracy is 90.71% and jackknife gives 96.07% accuracy. The average accuracy for level-2 i.e. protease classification is 95.77%. Based on the above-mentioned results, it is concluded that iProtease-PseAAC (2L) has the great ability to identify the proteases and their classes using a given protein sequence.

Revelation of enzyme activity of mutant pyrazinamidases from Mycobacterium tuberculosis upon binding with various metals using quantum mechanical approach

Pyrazinamide (PZA) is one of the most potent bacteriostatic drug against tuberculosis, a deadliest disease with high mortality and morbidity rate. PZA metabolizes into its active form pyrazinoic acid (POA) with the help of a metalloenzyme, pyrazinamidase (PZase). Mutagenicity and metal substitution in PZase weakens the binding of PZA with PZase and increases the drug resistance in Mycobacterium tuberculosis. The present study aims at the quantum mechanistic analysis of mutant-metal substituted PZase complexes by studying the mechanics of metals and PZA binding at MCS and catalytic site, respectively. A total of 66 complexes are scrutinised in this study to elucidate the effect of mutations on the enzymatic function of PZase. Among the 10 mutations considered in this study, 7 different mutations i.e. Asp₄₉ → Asn, His₅₁ → Arg, Gly₇₈ → Cys, Asp₁₂ → Gly, Asp₁₂ → Ala, Thr₁₃₅ → Pro and Asp₁₃₆ → Gly cause a detrimental effect on the activity of PZase. In addition to this, the substitution of iron with cobalt enhances the enzymatic activity of both wild type and mutant PZase while zinc, magnesium and copper reduce it. Based on these results, it is concluded that upon substitution of iron with zinc, magnesium and copper, PZase cannot function properly. Due to mutations, the reactivity of the drug also reduces as its binding with PZase weakens and this phenomenon enhances the resistance of Mycobacterium tuberculosis against drug.

iMethylK_pseAAC: Improving Accuracy of Lysine Methylation Sites Identification by Incorporating Statistical Moments and Position Relative Features into General PseAAC via Chou's 5-steps Rule

BACKGROUND

Methylation is one of the most important post-translational modifications in the human body which usually arises on lysine among the most intensely modified residues. It performs a dynamic role in numerous biological procedures, such as regulation of gene expression, regulation of protein function and RNA processing. Therefore, to identify lysine methylation sites is an important challenge as some experimental procedures are time-consuming.

OBJECTIVE

Herein, we propose a computational predictor named iMethylK_pseAAC to identify lysine methylation sites.

METHODS

Firstly, we constructed feature vectors based on PseAAC using position and composition rel-ative features and statistical moments. A neural network is trained based on the extracted features. The performance of the proposed method is then validated using cross-validation and jackknife testing.

RESULTS

The objective evaluation of the predictor showed accuracy of 96.7% for self-consistency, 91.61% for 10-fold cross-validation and 93.42% for jackknife testing.

CONCLUSION

It is concluded that iMethylK_pseAAC outperforms the counterparts to identify lysine methylation sites such as iMethyl_pseACC, BPB_pPMS and PMeS.