M Michael Gromiha

Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant of concern (VoC) Omicron (B.1.1.529) has rapidly spread around the world, presenting a new threat to global public human health. Due to the large number of mutations accumulated by SARS-CoV-2 Omicron, concerns have emerged over potentially reduced diagnostic accuracy of reverse-transcription polymerase chain reaction (RT-qPCR), the gold standard diagnostic test for diagnosing coronavirus disease 2019 (COVID-19). Thus, we aimed to assess the impact of the currently endemic Omicron sublineages BA.4 and BA.5 on the integrity and sensitivity of RT-qPCR assays used for coronavirus disease 2019 (COVID-19) diagnosis via in silico analysis. We employed whole genome sequencing data and evaluated the potential for false negatives or test failure due to mismatches between primers/probes and the Omicron VoC viral genome. Methods: In silico sensitivity of 12 RT-qPCR tests (containing 30 primers and probe sets) developed for detection of SARS-CoV-2 reported by the World Health Organization (WHO) or available in the literature, was assessed for specifically detecting SARS-CoV-2 Omicron BA.4 and BA.5 sublineages, obtained after removing redundancy from publicly available genomes from National Center for Biotechnology Information (NCBI) and Global Initiative on Sharing Avian Influenza Data (GISAID) databases. Mismatches between amplicon regions of SARS-CoV-2 Omicron VoC and primers and probe sets were evaluated, and clustering analysis of corresponding amplicon sequences was carried out. Results: From the 1164 representative SARS-CoV-2 Omicron VoC BA.4 sublineage genomes analyzed, a substitution in the first five nucleotides (C to T) of the amplicon's 3′-end was observed in all samples resulting in 0% sensitivity for assays HKUnivRdRp/Hel (mismatch in reverse primer) and CoremCharite N (mismatch in both forward and reverse primers). Due to a mismatch in the forward primer's 5′-end (3-nucleotide substitution, GGG to AAC), the sensitivity of the ChinaCDC N assay was at 0.69%. The 10 nucleotide mismatches in the reverse primer resulted in 0.09% sensitivity for Omicron sublineage BA.4 for Thai N assay. Of the 1926 BA.5 sublineage genomes, HKUnivRdRp/Hel assay also had 0% sensitivity. A sensitivity of 3.06% was observed for the ChinaCDC N assay because of a mismatch in the forward primer's 5′-end (3-nucleotide substitution, GGG to AAC). Similarly, due to the 10 nucleotide mismatches in the reverse primer, the Thai N assay's sensitivity was low at 0.21% for sublineage BA.5. Further, eight assays for BA.4 sublineage retained high sensitivity (more than 97%) and 9 assays for BA.5 sublineage retained more than 99% sensitivity. Conclusion: We observed four assays (HKUnivRdRp/Hel, ChinaCDC N, Thai N, CoremCharite N) that could potentially result in false negative results for SARS-CoV-2 Omicron VoCs BA.4 and BA.5 sublineages. Interestingly, CoremCharite N had 0% sensitivity for Omicron Voc BA.4 but 99.53% sensitivity for BA.5. In addition, 66.67% of the assays for BA.4 sublineage and 75% of the assays for BA.5 sublineage retained high sensitivity. Further, amplicon clustering and additional substitution analysis along with sensitivity analysis could be used for the modification and development of RT-qPCR assays for detecting SARS-CoV-2 Omicron VoC sublineages.

Background: Protein-carbohydrate interactions play vital roles in several biological processes in living organisms. The comparative analysis of binding site residues along with stabilizing residues in protein-carbohydrate complexes provides ample insights to understand the structure, function and recognition mechanism. Objective: The main objective of this study is to identify and analyze the residues, which are involved in both folding and binding of the protein-carbohydrate complexes. Methods: We have identified the stabilizing residues using the knowledge of hydrophobicity, longrange interactions and conservation, as well as binding site residues using a distance cutoff of 3.5Å between any heavy atoms in protein and ligand. Residues, which are common in stabilizing and binding, are termed as key residues. These key resides are analyzed with various sequence and structure based parameters such as frequency of occurrence, surrounding hydrophobicity, longrange order and conservation score. Results: In this work, we have identified 2.45% binding site residues in a non-redundant dataset of 1130 complexes using distance-based criteria and 7.07% stabilizing residues using the concepts of hydrophobicity, long-range interactions and conservation of residues. Further, 5.9% of binding and 2.04% of stabilizing residues are common to each other, which are termed as key residues. The key residues have been analysed based on protein classes, carbohydrate types, gene ontology functional classifications, amino acid preference and structure-based parameters. We found that all-β, α+β and α/β have more key residues than other protein classes and most of the KRs are present in β-strands, which shows their importance in stability and binding of complexes. On the ligand side, Lsaccharide has the highest number of key residues and it has a high percentage of KRs in SRs and BRs than other carbohydrate types. Further, polar and charged residues have a high tendency to serve as key residues. Classifications based on gene ontology terms revealed that Lys is preferred in all the three groups: molecular functions, biological processes and cellular components. Key residues have 6 to 9 contacts within the protein and make only one contact with the carbohydrate ligand. These contacts are dominant to form polar-nonpolar contacts followed by the contacts between charged atoms. Further, the influence of sequence and structural parameters such as surrounding hydrophobicity, solvent accessibility, secondary structure, long-range order and conservation score has been discussed. Conclusion: The results obtained in the present work provide deep insights for understanding the interplay between stability and binding in protein-carbohydrate complexes. © 2018 Bentham Science Publishers.

Protein-nucleic acid interactions are involved in various biological processes such as gene expression, replication, transcription, translation and packaging. The binding affinities of protein-DNA and protein-RNA complexes are important for elucidating the mechanism of protein-nucleic acid recognition. Although experimental data on binding affinity are reported abundantly in the literature, no well-curated database is currently available for protein-nucleic acid binding affinity. We have developed a database, ProNAB, which contains more than 20 000 experimental data for the binding affinities of protein-DNA and protein-RNA complexes. Each entry provides comprehensive information on sequence and structural features of a protein, nucleic acid and its complex, experimental conditions, thermodynamic parameters such as dissociation constant (Kd), binding free energy (ΔG) and change in binding free energy upon mutation (ΔΔG), and literature information. ProNAB is cross-linked with GenBank, UniProt, PDB, ProThermDB, PROSITE, DisProt and Pubmed. It provides a user-friendly web interface with options for search, display, sorting, visualization, download and upload the data. ProNAB is freely available at https://web.iitm.ac.in/bioinfo2/pronab/ and it has potential applications such as understanding the factors influencing the affinity, development of prediction tools, binding affinity change upon mutation and design complexes with the desired affinity.

Protein-DNA recognition plays an essential role in the regulation of gene expression. Understanding the recognition mechanism of protein-DNA complexes is a challenging task in molecular and computational biology. In this work, a scoring function based approach has been developed for identifying the binding sites and delineating the important residues for binding in protein-DNA complexes. This approach considers both the repulsive interactions and the effect of distance between atoms in protein and DNA. The results showed that positively charged, polar, and aromatic residues are important for binding. These residues influence the formation of electrostatic, hydrogen bonding, and stacking interactions. Our observation has been verified with experimental binding specificity of protein-DNA complexes and found to be in good agreement with experiments. The comparison of protein-RNA and protein-DNA complexes reveals that the contribution of phosphate atoms in DNA is twice as large as in protein-RNA complexes. Furthermore, we observed that the positively charged, polar, and aromatic residues serve as hotspot residues in protein-RNA complexes, whereas other residues also altered the binding specificity in protein-DNA complexes. Based on the results obtained in the present study and related reports, a plausible mechanism has been proposed for the recognition of protein-DNA complexes. © 2011 American Chemical Society.

Several prediction algorithms and tools have been developed in the last two decades to predict protein and peptide aggregation. These in silico tools aid to predict the aggregation propensity and amyloidogenicity as well as the identification of aggregation-prone regions. Despite the immense interest in the field, it is of prime importance to systematically compare these algorithms for their performance. In this review, we have provided a rigorous performance analysis of nine prediction tools using a variety of assessments. The assessments were carried out on several non-redundant datasets ranging from hexapeptides to protein sequences as well as amyloidogenic antibody light chains to soluble protein sequences. Our analysis reveals the robustness of the current prediction tools and the scope for improvement in their predictive performances. Insights gained from this work provide critical guidance to the scientific community on advantages and limitations of different aggregation prediction methods and make informed decisions about their research needs.

Introduction: Alzheimer's disease (AD) is a complex and heterogeneous disease that affects neuronal cells over time and it is prevalent among all neurodegenerative diseases. Next Generation Sequencing (NGS) techniques are widely used for developing high-throughput screening methods to identify biomarkers and variants, which help early diagnosis and treatments. Objective: The primary purpose of this study is to develop a classification model using machine learning for predicting the deleterious effect of variants with respect to AD. Methods: We have constructed a set of 20,401 deleterious and 37,452 control variants from Genome-Wide Association Study (GWAS) and Genotype-Tissue Expression (GTEx) portals, respectively. Recursive feature elimination using cross-validation (RFECV) followed by a forward feature selection method was utilized to select the important features and a random forest classifier was used for distinguishing between deleterious and neutral variants. Results: Our method showed an accuracy of 81.21% on 10-fold cross-validation and 70.63% on a test set of 5785 variants. The same test set was used to compare the performance of CADD and FATHMM and their accuracies are in the range of 54%–62%. Conclusion: Our model is freely available as the Variant Effect Predictor for Alzheimer's Disease (VEPAD) at http://web.iitm.ac.in/bioinfo2/vepad/. VEPAD can be used to predict the effect of new variants associated with AD.

Why do patients suffering from neurodegenerative diseases generate autoantibodies that selectively bind soluble aggregates of amyloidogenic proteins? Presently, molecular basis of interactions between the soluble aggregates and human immune system is unknown. By analyzing sequences of experimentally validated T-cell autoimmune epitopes, aggregating peptides, amyloidogenic proteins and randomly generated peptides, here we report overlapping regions that likely drive aggregation as well as generate autoantibodies against the aggregates. Sequence features, that make short peptides susceptible to aggregation, increase their incidence in human T-cell autoimmune epitopes by 4-6 times. Many epitopes are predicted to be significantly aggregation prone (aggregation propensities ≥10%) and the ones containing experimentally validated aggregating regions are enriched in hydrophobicity by 10-20%. Aggregate morphologies also influence Human Leukocyte Antigen (HLA)-types recognized by the aggregating regions containing epitopes. Most (88%) epitopes that contain amyloid fibril forming regions bind HLA-DR, while majority (63%) of those containing amorphous β-Aggregating regions bind HLA-DQ. More than two-Thirds (70%) of human amyloidogenic proteins contain overlapping regions that are simultaneously aggregation prone and auto-immunogenic. Such regions help clear soluble aggregates by generating selective autoantibodies against them. This can be harnessed for early diagnosis of proteinopathies and for drug/vaccine design against them.

Protein-protein interactions are intrinsic to virtually every cellular process. Predicting the binding affinity of protein-protein complexes is one of the challenging problems in computational and molecular biology. In this work, we related sequence features of protein-protein complexes with their binding affinities using machine learning approaches. We set up a database of 185 protein-protein complexes for which the interacting pairs are heterodimers and their experimental binding affinities are available. On the other hand, we have developed a set of 610 features from the sequences of protein complexes and utilized Ranker search method, which is the combination of Attribute evaluator and Ranker method for selecting specific features. We have analyzed several machine learning algorithms to discriminate protein-protein complexes into high and low affinity groups based on their Kd values. Our results showed a 10-fold cross-validation accuracy of 76.1% with the combination of nine features using support vector machines. Further, we observed accuracy of 83.3% on an independent test set of 30 complexes. We suggest that our method would serve as an effective tool for identifying the interacting partners in protein-protein interaction networks and human-pathogen interactions based on the strength of interactions. © 2014 Wiley Periodicals, Inc.

Protein-carbohydrate interactions play an important role in several biological processes. The mutation of amino acid residues in carbohydrate-binding proteins may alter the binding affinity, affect the functions and lead to diseases. Elucidating the factors influencing the binding affinity change (ΔΔG) of protein-carbohydrate complexes upon mutation is a challenging task. In this work, we have collected the experimental data for the binding affinity change of 318 unique mutants and related with sequence and structural features of amino acid residues at the mutant sites. We found that accessible surface area, secondary structure, mutation preference, conservation score, hydrophobicity and contact energies are important to understand the binding affinity change upon mutation. We have developed multiple regression equations for predicting the binding affinity change upon mutation and our method showed an average correlation of 0.74 and a mean absolute error of 0.70 kcal/mol between experimental and predicted ΔΔG on a 10-fold cross-validation. Further, we have validated our method using an independent test data set of 124 (62 unique) mutations, which showed a correlation and MAE of 0.79 and 0.56 kcal/mol, respectively. We have developed a web server PCA-MutPred, Protein-CArbohydrate complex Mutation affinity Predictor, for predicting the change in binding affinity of protein–carbohydrate complexes and it is freely accessible at https://web.iitm.ac.in/bioinfo2/pcamutpred. We suggest that the method could be a useful resource for designing protein-carbohydrate complexes with desired affinities.

Studying protein–protein interactions leads to a better understanding of the underlying principles of several biological pathways. Cost and labor-intensive experimental techniques suggest the need for computational methods to complement them. Several such state-of-the-art methods have been reported for analyzing diverse aspects such as predicting binding partners, interface residues, and binding affinity for protein–protein complexes with reliable performance. However, there are specific drawbacks for different methods that indicate the need for their improvement. This review highlights various available computational algorithms for analyzing diverse aspects of protein–protein interactions and endorses the necessity for developing new robust methods for gaining deep insights about protein–protein interactions.