Validation and quality assessment of macromolecular structures using complex network analysis

Albertson (Alb) spectral radii and Albertson (Alb) energies of graph operation

The sum of the absolute eigenvalues of the adjacency matrix make up graph energy. The greatest absolute eigenvalue of the adjacency matrix is represented by the spectral radius of the graph. Both molecular computing and computer science have uses for graph energies and spectral radii. The Albertson (Alb) energies and spectral radii of generalized splitting and shadow graphs constructed on any regular graph is the main focus of this study. The only thing that may be disputed is the comparison of the (Alb) energies and (Alb) spectral radii of the newly formed graphs to those of the base graph. By concentrating on splitting and shadow graph, we compute new correlations between the Alb energies and spectral radius of the new graph and the prior graph.

Computational formulation of a multiepitope vaccine unveils an exceptional prophylactic candidate against Merkel cell polyomavirus

Merkel cell carcinoma (MCC) is a rare neuroendocrine skin malignancy caused by human Merkel cell polyomavirus (MCV), leading to the most aggressive skin cancer in humans. MCV has been identified in approximately 43%-100% of MCC cases, contributing to the highly aggressive nature of primary cutaneous carcinoma and leading to a notable mortality rate. Currently, no existing vaccines or drug candidates have shown efficacy in addressing the ailment caused by this specific pathogen. Therefore, this study aimed to design a novel multiepitope vaccine candidate against the virus using integrated immunoinformatics and vaccinomics approaches. Initially, the highest antigenic, immunogenic, and non-allergenic epitopes of cytotoxic T lymphocytes, helper T lymphocytes, and linear B lymphocytes corresponding to the virus whole protein sequences were identified and retrieved for vaccine construction. Subsequently, the selected epitopes were linked with appropriate linkers and added an adjuvant in front of the construct to enhance the immunogenicity of the vaccine candidates. Additionally, molecular docking and dynamics simulations identified strong and stable binding interactions between vaccine candidates and human Toll-like receptor 4. Furthermore, computer-aided immune simulation found the real-life-like immune response of vaccine candidates upon administration to the human body. Finally, codon optimization was conducted on the vaccine candidates to facilitate the in silico cloning of the vaccine into the pET28+(a) cloning vector. In conclusion, the vaccine candidate developed in this study is anticipated to augment the immune response in humans and effectively combat the virus. Nevertheless, it is imperative to conduct in vitro and in vivo assays to evaluate the efficacy of these vaccine candidates thoroughly. These evaluations will provide critical insights into the vaccine's effectiveness and potential for further development.

Genome-wide identification, comprehensive characterization of transcription factors, cis-regulatory elements, protein homology, and protein interaction network of DREB gene family in Solanum lycopersicum

Tomato is a drought-sensitive crop which has high susceptibility to adverse climatic changes. Dehydration-responsive element-binding (DREB) are significant plant transcription factors that have a vital role in regulating plant abiotic stress tolerance by networking with DRE/CRT cis-regulatory elements in response to stresses. In this study, bioinformatics analysis was performed to conduct the genome-wide identification and characterization of DREB genes and promoter elements in Solanum lycopersicum. In genome-wide coverage, 58 SlDREB genes were discovered on 12 chromosomes that justified the criteria of the presence of AP2 domain as conserved motifs. Intron-exon organization and motif analysis showed consistency with phylogenetic analysis and confirmed the absence of the A3 class, thus dividing the SlDREB genes into five categories. Gene expansion was observed through tandem duplication and segmental duplication gene events in SlDREB genes. Ka/Ks values were calculated in ortholog pairs that indicated divergence time and occurrence of purification selection during the evolutionary period. Synteny analysis demonstrated that 32 out of 58 and 47 out of 58 SlDREB genes were orthologs to Arabidopsis and Solanum tuberosum, respectively. Subcellular localization predicted that SlDREB genes were present in the nucleus and performed primary functions in DNA binding to regulate the transcriptional processes according to gene ontology. Cis-acting regulatory element analysis revealed the presence of 103 motifs in 2.5-kbp upstream promoter sequences of 58 SlDREB genes. Five representative SlDREB proteins were selected from the resultant DREB subgroups for 3D protein modeling through the Phyre2 server. All models confirmed about 90% residues in the favorable region through Ramachandran plot analysis. Moreover, active catalytic sites and occurrence in disorder regions indicated the structural and functional flexibility of SlDREB proteins. Protein association networks through STRING software suggested the potential interactors that belong to different gene families and are involved in regulating similar functional and biological processes. Transcriptome data analysis has revealed that the SlDREB gene family is engaged in defense response against drought and heat stress conditions in tomato. Overall, this comprehensive research reveals the identification and characterization of SlDREB genes that provide potential knowledge for improving abiotic stress tolerance in tomato.

Maximum degree and minimum degree spectral radii of some graph operations

New results relating to the maximum and minimum degree spectral radii of generalized splitting and shadow graphs have been constructed on the basis of any regular graph, referred as base graph. In particular, we establish the relations of extreme degree spectral radii of generalized splitting and shadow graphs of any regular graph.

Structure and dynamics of membrane protein in SARS-CoV-2

SARS-CoV-2 membrane (M) protein performs a variety of critical functions in virus infection cycle. However, the expression and purification of membrane protein structure is difficult despite tremendous progress. In this study, the 3 D structure is modeled followed by intensive validation and molecular dynamics simulation. The lack of suitable homologous templates (>30% sequence identities) leads us to construct the membrane protein models using template-free modeling (de novo or ab initio) approach with Robetta and trRosetta servers. Comparing with other model structures, it is evident that trRosetta (TM-score: 0.64; TM region RMSD: 2 Å) can provide the best model than Robetta (TM-score: 0.61; TM region RMSD: 3.3 Å) and I-TASSER (TM-score: 0.45; TM region RMSD: 6.5 Å). 100 ns molecular dynamics simulations are performed on the model structures by incorporating membrane environment. Moreover, secondary structure elements and principal component analysis (PCA) have also been performed on MD simulation data. Finally, trRosetta model is utilized for interpretation and visualization of interacting residues during protein-protein interactions. The common interacting residues including Phe103, Arg107, Met109, Trp110, Arg131, and Glu135 in the C-terminal domain of M protein are identified in membrane-spike and membrane-nucleocapsid protein complexes. The active site residues are also predicted for potential drug and peptide binding. Overall, this study might be helpful to design drugs and peptides against the modeled membrane protein of SARS-CoV-2 to accelerate further investigation. Communicated by Ramaswamy H. Sarma.

Some New Bounds for the Inverse Sum Indeg Energy of Graphs

Let G be a (molecular) graph with n vertices, and di be the degree of its i-th vertex. Then, the inverse sum indeg matrix of G is the n×n matrix C(G) with entries cij=didjdi+dj, if the i-th and the j-th vertices are adjacent and 0 otherwise. Let μ1≥μ2≥…≥μn be the eigenvalues of C arranged in order. The inverse sum indeg energy of G, εisi(G) can be represented as ∑j=1n|μi|. In this paper, we establish several novel upper and lower sharp bounds on μ1 and εisi(G) via some other graph parameters, and describe the structures of the extremal graphs.

Randic and reciprocal randic spectral radii and energies of some graph operations

The largest absolute eigenvalue of a matrix A associated to the graph G is called the A-Spectral Radius of the graph G, and A-energy of the graph G is defined as the absolute sum of all its eigenvalues. In the present article, we compute Randic energies, Reciprocal Randic energies, Randic spectral radii and Reciprocal Randic radii of s-shadow and s-splitting graph of G. We actually relate these energies and Spectral Radii of new graphs with the energies and Spectral Radii of original graphs.

Genome-wide promoter analysis, homology modeling and protein interaction network of Dehydration Responsive Element Binding (DREB) gene family in Solanum tuberosum

Dehydration Responsive Element Binding (DREB) regulates the expression of numerous stress-responsive genes, and hence plays a pivotal role in abiotic stress responses and tolerance in plants. The study aimed to develop a complete overview of the cis-acting regulatory elements (CAREs) present in S. tuberosum DREB gene promoters. A total of one hundred and four (104) cis-regulatory elements (CREs) were identified from 2.5kbp upstream of the start codon (ATG). The in-silico promoter analysis revealed variable sets of cis-elements and functional diversity with the predominance of light-responsive (30%), development-related (20%), abiotic stress-responsive (14%), and hormone-responsive (12%) elements in StDREBs. Among them, two light-responsive elements (Box-4 and G-box) were predicted in 64 and 61 StDREB genes, respectively. Two development-related motifs (AAGAA-motif and as-1) were abundant in StDREB gene promoters. Most of the DREB genes contained one or more Myeloblastosis (MYB) and Myelocytometosis (MYC) elements associated with abiotic stress responses. Hormone-responsive element i.e. ABRE was found in 59 out of 66 StDREB genes, which implied their role in dehydration and salinity stress. Moreover, six proteins were chosen corresponding to A1-A6 StDREB subgroups for secondary structure analysis and three-dimensional protein modeling followed by model validation through PROCHECK server by Ramachandran Plot. The predicted models demonstrated >90% of the residues in the favorable region, which further ensured their reliability. The present study also anticipated pocket binding sites and disordered regions (DRs) to gain insights into the structural flexibility and functional annotation of StDREB proteins. The protein association network determined the interaction of six selected StDREB proteins with potato proteins encoded by other gene families such as MYB and NAC, suggesting their similar functional roles in biological and molecular pathways. Overall, our results provide fundamental information for future functional analysis to understand the precise molecular mechanisms of the DREB gene family in S. tuberosum.

Escherichia coli expression and characterization of α-amylase from Geobacillus thermodenitrificans DSM-465

Alpha amylase, catalyzing the hydrolysis of starch is a ubiquitous enzyme with tremendous industrial applications. A 1698 bp gene coding for 565 amino acid amylase was PCR amplified from Geobacillus thermodenitrificans DSM-465, cloned in pET21a (+) plasmid, expressed in BL21 (DE3) strain of E. coli and characterized. The recombinant enzyme exhibited molecular weight of 63 kDa, optimum pH 8, optimum temperature 70°C, and KM value of 157.7µM. On pilot scale, the purified enzyme efficiently removed up to 95% starch from the cotton fabric indicating its desizing ability at high temperature. 3D model of enzyme built by Raptor-X and validated by Ramachandran plot appeared as a monomer having 31% α-helices, 15% β-sheets, and 52% loops. Docking studies have shown the best binding affinity of enzyme with amylopectin (∆G -10.59). According to our results, Asp 232, Glu274, Arg448, Glu385, Asp34, Asn276, and Arg175 constitute the potential active site of enzyme.

Protein conformational switch discerned via network centrality properties

Network analysis has emerged as a powerful tool for examining structural biology systems. The spatial organization of the components of a biomolecular structure has been rendered as a graph representation and analyses have been performed to deduce the biophysical and mechanistic properties of these components. For proteins, the analysis of protein structure networks (PSNs), especially via network centrality measurements and cluster coefficients, has led to identifying amino acid residues that play key functional roles and classifying amino acid residues in general. Whether these network properties examined in various studies are sensitive to subtle (yet biologically significant) conformational changes remained to be addressed. Here, we focused on four types of network centrality properties (betweenness, closeness, degree, and eigenvector centralities) for conformational changes upon ligand binding of a sensor protein (constitutive androstane receptor) and an allosteric enzyme (ribonucleotide reductase). We found that eigenvector centrality is sensitive and can distinguish salient structural features between protein conformational states while other centrality measures, especially closeness centrality, are less sensitive and rather generic with respect to the structural specificity. We also demonstrated that an ensemble-informed, modified PSN with static edges removed (which we term PSN*) has enhanced sensitivity at discerning structural changes.

The loops of the N-SH2 binding cleft do not serve as allosteric switch in SHP2 activation

The Src-homology-2 domain-containing phosphatase SHP2 is a critical regulator of signal transduction, being implicated in cell growth and differentiation. Activating mutations cause developmental disorders and act as oncogenic drivers in hematologic cancers. SHP2 is activated by phosphopeptide binding to the N-SH2 domain, triggering the release of N-SH2 from the catalytic PTP domain. Based on early crystallographic data, it has been widely accepted that opening of the binding cleft of N-SH2 serves as the key "allosteric switch" driving SHP2 activation. To test the putative coupling between binding cleft opening and SHP2 activation as assumed by the allosteric switch model, we critically reviewed structural data of SHP2, and we used extensive molecular dynamics (MD) simulation and free energy calculations of isolated N-SH2 in solution, SHP2 in solution, and SHP2 in a crystal environment. Our results demonstrate that the binding cleft in N-SH2 is constitutively flexible and open in solution and that a closed cleft found in certain structures is a consequence of crystal contacts. The degree of opening of the binding cleft has only a negligible effect on the free energy of SHP2 activation. Instead, SHP2 activation is greatly favored by the opening of the central β-sheet of N-SH2. We conclude that opening of the N-SH2 binding cleft is not the key allosteric switch triggering SHP2 activation.

An Educational Bioinformatics Project to Improve Genome Annotation

Scientific advancement is hindered without proper genome annotation because biologists lack a complete understanding of cellular protein functions. In bacterial cells, hypothetical proteins (HPs) are open reading frames with unknown functions. HPs result from either an outdated database or insufficient experimental evidence (i.e., indeterminate annotation). While automated annotation reviews help keep genome annotation up to date, often manual reviews are needed to verify proper annotation. Students can provide the manual review necessary to improve genome annotation. This paper outlines an innovative classroom project that determines if HPs have outdated or indeterminate annotation. The Hypothetical Protein Characterization Project uses multiple well-documented, freely available, web-based, bioinformatics resources that analyze an amino acid sequence to (1) detect sequence similarities to other proteins, (2) identify domains, (3) predict tertiary structure including active site characterization and potential binding ligands, and (4) determine cellular location. Enough evidence can be generated from these analyses to support re-annotation of HPs or prioritize HPs for experimental examinations such as structural determination via X-ray crystallography. Additionally, this paper details several approaches for selecting HPs to characterize using the Hypothetical Protein Characterization Project. These approaches include student- and instructor-directed random selection, selection using differential gene expression from mRNA expression data, and selection based on phylogenetic relations. This paper also provides additional resources to support instructional use of the Hypothetical Protein Characterization Project, such as example assignment instructions with grading rubrics, links to training videos in YouTube, and several step-by-step example projects to demonstrate and interpret the range of achievable results that students might encounter. Educational use of the Hypothetical Protein Characterization Project provides students with an opportunity to learn and apply knowledge of bioinformatic programs to address scientific questions. The project is highly customizable in that HP selection and analysis can be specifically formulated based on the scope and purpose of each student's investigations. Programs used for HP analysis can be easily adapted to course learning objectives. The project can be used in both online and in-seat instruction for a wide variety of undergraduate and graduate classes as well as undergraduate capstone, honor's, and experiential learning projects.

Domain sliding of two Staphylococcus aureus N-acetylglucosaminidases enables their substrate-binding prior to its catalysis

To achieve productive binding, enzymes and substrates must align their geometries to complement each other along an entire substrate binding site, which may require enzyme flexibility. In pursuit of novel drug targets for the human pathogen S. aureus, we studied peptidoglycan N-acetylglucosaminidases, whose structures are composed of two domains forming a V-shaped active site cleft. Combined insights from crystal structures supported by site-directed mutagenesis, modeling, and molecular dynamics enabled us to elucidate the substrate binding mechanism of SagB and AtlA-gl. This mechanism requires domain sliding from the open form observed in their crystal structures, leading to polysaccharide substrate binding in the closed form, which can enzymatically process the bound substrate. We suggest that these two hydrolases must exhibit unusual extents of flexibility to cleave the rigid structure of a bacterial cell wall.

Small Conformational Changes Underlie Evolution of Resistance to NNRTI in HIV Reverse Transcriptase

Despite achieving considerable success in reducing the number of fatalities due to acquired immunodeficiency syndrome, emergence of resistance against the reverse transcriptase (RT) inhibitor drugs remains one of the biggest challenges of the human immunodeficiency virus antiretroviral therapy (ART). Non-nucleoside reverse transcriptase inhibitors (NNRTIs) form a large class of drugs and a crucial component of ART. In NNRTIs, even a single resistance mutation is known to make the drugs completely ineffective. Additionally, several inhibitor-bound RTs with single resistance mutations do not exhibit any significant variations in their three-dimensional structures compared with the inhibitor-bound RT but completely nullify their inhibitory functions. This makes understanding the structural mechanism of these resistance mutations crucial for drug development. Here, we study several single resistance mutations in the allosteric inhibitor (nevirapine)-bound RT to analyze the mechanism of small structural changes leading to these large functional effects. In this study, we have shown that in absence of significant conformational variations in the inhibitor-bound wild-type RT and RT with single resistance mutations, the protein contact network analysis of their static structures, along with molecular dynamics simulations, can be a useful approach to understand the functional effect of small local conformational variations. The simple network analysis exposes the localized contact changes that lead to global rearrangement in the communication pattern within RT. Furthermore, these conformational changes have implications on the overall dynamics of RT. Using various measures, we show that a single resistance mutation can change the network structure and dynamics of RT to behave more like unbound RT, even in the presence of the inhibitor. This combined coarse-grained contact network and molecular dynamics approach promises to be a useful tool to analyze structure-function studies of proteins that show large functional changes with negligible variations in their overall conformation.