Interacting motif networks located in hotspots associated with RNA release are conserved in Enterovirus capsids

The Structural Framework and Opening Appearance of the VP1-Pocket of Enteroviruses Correlated with Viral Thermostability

Enteroviruses (EVs and RVs) are prevalent worldwide and cause various diseases in humans, of which the VP1-pocket is a target of antivirals, with a lipid molecule as a pocket factor to stabilize the virion. However, the characterization of the structure of the VP1-pocket in EVs is poor. Here, we compared the published capsid crystals of EVs and RVs and proposed a structural framework for the VP1-pocket: Frame 1-4, which is located at the CD loop, GH loop, and C-terminus, presenting with an outward opening appearance or not. The non-outward viral strains-CVB3, Echo 11, RV-A81, and RV-B70-are more thermally stable, with a breakpoint temperature (B.T.) of 51~62 °C for genome releasing, which is 4~10 °C higher than its outward temperature of 41~47 °C, and infectivity preservation when treated at 50 °C for 3 min. Its outward versus non-outward opening is correlated significantly with the B.T. for genome release (r = -0.90; p = 0.0004) and infectivity (r = -0.82, p = 0.0039). The energy of Frames 1, 2, and 4, including Van der Waals attractive and repulsive interactions and hydrogen bonds, showed significant correlations with the B.T. (r = -0.67, 0.75, and -0.8; p = 0.034, 0.013, and 0.006, respectively). These characters of the VP1-pocket could be predictors for virion thermostability and aid in the development of vaccines or antivirals.

Identification of Critical Amino Acids of Coxsackievirus A10 Associated with Cell Tropism and Viral RNA Release during Uncoating

Coxsackievirus A10 (CV-A10) is a prevailing causative agent of hand-foot-mouth disease, necessitating the isolation and adaptation of appropriate strains in cells allowed for human vaccine development. In this study, amino acid sequences of CV-A10 strains with different cell tropism on RD and Vero cells were compared. Various amino acids on the structural and non-structural proteins related to cell tropism were identified. The reverse genetic systems of several CV-A10 strains with RD+/Vero- and RD+/Vero+ cell tropism were developed, and a set of CV-A10 recombinants were produced. The binding, entry, uncoating, and proliferation steps in the life cycle of these viruses were evaluated. P1 replacement of CV-A10 strains with different cell tropism revealed the pivotal role of the structural proteins in cell tropism. Further, seven amino acid substitutions in VP2 and VP1 were introduced to further investigate their roles played in cell tropism. These mutations cooperated in the growth of CV-A10 in Vero cells. Particularly, the valine to isoleucine mutation at the position VP1-236 (V1236I) was found to significantly restrict viral uncoating in Vero cells. Co-immunoprecipitation assays showed that the release of viral RNA from the KREMEN1 receptor-binding virions was restricted in r0195-V1236I compared with the parental strain r0195 (a RD+/Vero+ strain). Overall, this study highlights the dominant effect of structural proteins in CV-A10 adaption in Vero cells and the importance of V1236 in viral uncoating, providing a foundation for the mechanism study of CV-A10 cell tropism, and facilitating the development of vaccine candidates.

Structural and Functional Annotation of Hypothetical Proteins from the Microsporidia Species Vittaforma corneae ATCC 50505 Using in silico Approaches

Microsporidia are spore-forming eukaryotes that are related to fungi but have unique traits that set them apart. They have compact genomes as a result of evolutionary gene loss associated with their complete dependency on hosts for survival. Despite having a relatively small number of genes, a disproportionately high percentage of the genes in microsporidia genomes code for proteins whose functions remain unknown (hypothetical proteins-HPs). Computational annotation of HPs has become a more efficient and cost-effective alternative to experimental investigation. This research developed a robust bioinformatics annotation pipeline of HPs from Vittaforma corneae, a clinically important microsporidian that causes ocular infections in immunocompromised individuals. Here, we describe various steps to retrieve sequences and homologs and to carry out physicochemical characterization, protein family classification, identification of motifs and domains, protein-protein interaction network analysis, and homology modelling using a variety of online resources. Classification of protein families produced consistent findings across platforms, demonstrating the accuracy of annotation utilizing in silico methods. A total of 162 out of 2034 HPs were fully annotated, with the bulk of them categorized as binding proteins, enzymes, or regulatory proteins. The protein functions of several HPs from Vittaforma corneae were accurately inferred. This improved our understanding of microsporidian HPs despite challenges related to the obligate nature of microsporidia, the absence of fully characterized genes, and the lack of homologous genes in other systems.

Deciphering Isoniazid Drug Resistance Mechanisms on Dimeric Mycobacterium tuberculosis KatG via Post-molecular Dynamics Analyses Including Combined Dynamic Residue Network Metrics

Resistance mutations in Mycobacterium tuberculosis (Mtb) catalase peroxidase protein (KatG), an essential enzyme in isoniazid (INH) activation, reduce the sensitivity of Mtb to first-line drugs, hence presenting challenges in tuberculosis (TB) management. Thus, understanding the mutational imposed resistance mechanisms remains of utmost importance in the quest to reduce the TB burden. Herein, effects of 11 high confidence mutations in the KatG structure and residue network communication patterns were determined using extensive computational approaches. Combined traditional post-molecular dynamics analysis and comparative essential dynamics revealed that the mutant proteins have significant loop flexibility around the heme binding pocket and enhanced asymmetric protomer behavior with respect to wild-type (WT) protein. Heme contact analysis between WT and mutant proteins identified a reduction to no contact between heme and residue His270, a covalent bond vital for the heme-enabled KatG catalytic activity. Betweenness centrality calculations showed large hub ensembles with new hubs especially around the binding cavity and expanded to the dimerization domain via interface in the mutant systems, providing possible compensatory allosteric communication paths for the active site as a result of the mutations which may destabilize the heme binding pocket and the loops in its vicinity. Additionally, an interesting observation came from Eigencentrality hubs, most of which are located in the C-terminal domain, indicating relevance of the domain in the protease functionality. Overall, our results provide insight toward the mechanisms involved in KatG-INH resistance in addition to identifying key regions in the enzyme functionality, which can be used for future drug design.

Characterisation of plasmodial transketolases and identification of potential inhibitors: an in silico study

BACKGROUND

Plasmodial transketolase (PTKT) enzyme is one of the novel pharmacological targets being explored as potential anti-malarial drug target due to its functional role and low sequence identity to the human enzyme. Despite this, features contributing to such have not been exploited for anti-malarial drug design. Additionally, there are no anti-malarial drugs targeting PTKTs whereas the broad activity of these inhibitors against PTKTs from other Plasmodium spp. is yet to be reported. This study characterises different PTKTs [Plasmodium falciparum (PfTKT), Plasmodium vivax (PvTKT), Plasmodium ovale (PoTKT), Plasmodium malariae (PmTKT) and Plasmodium knowlesi (PkTKT) and the human homolog (HsTKT)] to identify key sequence and structural based differences as well as the identification of selective potential inhibitors against PTKTs.

METHODS

A sequence-based study was carried out using multiple sequence alignment, phylogenetic tree calculations and motif discovery analysis. Additionally, TKT models of PfTKT, PmTKT, PoTKT, PmTKT and PkTKT were modelled using the Saccharomyces cerevisiae TKT structure as template. Based on the modelled structures, molecular docking using 623 South African natural compounds was done. The stability, conformational changes and detailed interactions of selected compounds were accessed viz all-atom molecular dynamics (MD) simulations and binding free energy (BFE) calculations.

RESULTS

Sequence alignment, evolutionary and motif analyses revealed key differences between plasmodial and the human TKTs. High quality homodimeric three-dimensional PTKTs structures were constructed. Molecular docking results identified three compounds (SANC00107, SANC00411 and SANC00620) which selectively bind in the active site of all PTKTs with the lowest (better) binding affinity ≤ - 8.5 kcal/mol. MD simulations of ligand-bound systems showed stable fluctuations upon ligand binding. In all systems, ligands bind stably throughout the simulation and form crucial interactions with key active site residues. Simulations of selected compounds in complex with human TKT showed that ligands exited their binding sites at different time steps. BFE of protein-ligand complexes showed key residues involved in binding.

CONCLUSIONS

This study highlights significant differences between plasmodial and human TKTs and may provide valuable information for the development of novel anti-malarial inhibitors. Identified compounds may provide a starting point in the rational design of PTKT inhibitors and analogues based on these scaffolds.

Alpha-Carbonic Anhydrases from Hydrothermal Vent Sources as Potential Carbon Dioxide Sequestration Agents: In Silico Sequence, Structure and Dynamics Analyses

With the increase in CO2 emissions worldwide and its dire effects, there is a need to reduce CO2 concentrations in the atmosphere. Alpha-carbonic anhydrases (α-CAs) have been identified as suitable sequestration agents. This study reports the sequence and structural analysis of 15 α-CAs from bacteria, originating from hydrothermal vent systems. Structural analysis of the multimers enabled the identification of hotspot and interface residues. Molecular dynamics simulations of the homo-multimers were performed at 300 K, 363 K, 393 K and 423 K to unearth potentially thermostable α-CAs. Average betweenness centrality (BC) calculations confirmed the relevance of some hotspot and interface residues. The key residues responsible for dimer thermostability were identified by comparing fluctuating interfaces with stable ones, and were part of conserved motifs. Crucial long-lived hydrogen bond networks were observed around residues with high BC values. Dynamic cross correlation fortified the relevance of oligomerization of these proteins, thus the importance of simulating them in their multimeric forms. A consensus of the simulation analyses used in this study suggested high thermostability for the α-CA from Nitratiruptor tergarcus. Overall, our novel findings enhance the potential of biotechnology applications through the discovery of alternative thermostable CO2 sequestration agents and their potential protein design.

Probing the Structural Dynamics of the Plasmodium falciparum Tunneling-Fold Enzyme 6-Pyruvoyl Tetrahydropterin Synthase to Reveal Allosteric Drug Targeting Sites

The de novo folate synthesis pathway is a well-established drug target in the treatment of many infectious diseases. Antimalarial antifolate drugs have proven to be effective against malaria, however, rapid drug resistance has emerged on the two primary targeted enzymes: dihydrofolate reductase and dihydroptoreate synthase. The need to identify alternative antifolate drugs and novel metabolic targets is of imminent importance. The 6-pyruvol tetrahydropterin synthase (PTPS) enzyme belongs to the tunneling fold protein superfamily which is characterized by a distinct central tunnel/cavity. The enzyme catalyzes the second reaction step of the parasite’s de novo folate synthesis pathway and is responsible for the conversion of 7,8-dihydroneopterin to 6-pyruvoyl-tetrahydropterin. In this study, we examine the structural dynamics of Plasmodium falciparum PTPS using the anisotropic network model, to elucidate the collective motions that drive the function of the enzyme and identify potential sites for allosteric modulation of its binding properties. Based on our modal analysis, we identified key sites in the N-terminal domains and central helices which control the accessibility to the active site. Notably, the N-terminal domains were shown to regulate the open-to-closed transition of the tunnel, via a distinctive wringing motion that deformed the core of the protein. We, further, combined the dynamic analysis with motif discovery which revealed highly conserved motifs that are unique to the Plasmodium species and are located in the N-terminal domains and central helices. This provides essential structural information for the efficient design of drugs such as allosteric modulators that would have high specificity and low toxicity as they do not target the PTPS active site that is highly conserved in humans.

The In Silico Prediction of Hotspot Residues that Contribute to the Structural Stability of Subunit Interfaces of a Picornavirus Capsid

The assembly of picornavirus capsids proceeds through the stepwise oligomerization of capsid protein subunits and depends on interactions between critical residues known as hotspots. Few studies have described the identification of hotspot residues at the protein subunit interfaces of the picornavirus capsid, some of which could represent novel drug targets. Using a combination of accessible web servers for hotspot prediction, we performed a comprehensive bioinformatic analysis of the hotspot residues at the intraprotomer, interprotomer and interpentamer interfaces of the Theiler’s murine encephalomyelitis virus (TMEV) capsid. Significantly, many of the predicted hotspot residues were found to be conserved in representative viruses from different genera, suggesting that the molecular determinants of capsid assembly are conserved across the family. The analysis presented here can be applied to any icosahedral structure and provides a platform for in vitro mutagenesis studies to further investigate the significance of these hotspots in critical stages of the virus life cycle with a view to identify potential targets for antiviral drug design.

Mechanism of Action of Non-Synonymous Single Nucleotide Variations Associated with α-Carbonic Anhydrase II Deficiency

Human carbonic anhydrase II (CA-II) is a Zinc (Zn 2 + ) metalloenzyme responsible for maintenance of acid-base balance within the body through the reversible hydration of CO 2 to produce protons (H + ) and bicarbonate (BCT). Due to its importance, alterations to the amino acid sequence of the protein as a result of single nucleotide variations (nsSNVs) have detrimental effects on homeostasis. Six pathogenic CA-II nsSNVs, K18E, K18Q, H107Y, P236H, P236R and N252D were identified, and variant protein models calculated using homology modeling. The effect of each nsSNV was analyzed using motif analysis, molecular dynamics (MD) simulations, principal component (PCA) and dynamic residue network (DRN) analysis. Motif analysis identified 11 functionally important motifs in CA-II. RMSD data indicated subtle SNV effects, while PCA analysis revealed that the presence of BCT results in greater conformational sampling and free energy in proteins. DRN analysis showed variant allosteric effects, and the average betweenness centrality (BC) calculations identified Glu117 as the most important residue for communication in CA-II. The presence of BCT was associated with a reduction to Glu117 usage in all variants, suggesting implications for Zn 2 + dissociation from the CA-II active site. In addition, reductions to Glu117 usage are associated with increases in the usage of the primary and secondary Zn 2 + ligands; His94, His96, His119 and Asn243 highlighting potential compensatory mechanisms to maintain Zn 2 + within the active site. Compared to traditional MD simulation investigation, DRN analysis provided greater insights into SNV mechanism of action, indicating its importance for the study of missense mutation effects in proteins and, in broader terms, precision medicine related research.

Multifunctionality of structural proteins in the enterovirus life cycle

Members of the genus Enterovirus have a significant effect on human health, especially in infants and children. Since the viral genome has limited coding capacity, Enteroviruses subvert a range of cellular processes for viral infection via the interaction of viral proteins and numerous cellular factors. Intriguingly, the capsid-receptor interaction plays a crucial role in viral entry and has significant implications in viral pathogenesis. Moreover, interactions between structural proteins and host factors occur directly or indirectly in multiple steps of viral replication. In this review, we focus on the current understanding of the multifunctionality of structural proteins in the viral life cycle, which may constitute valuable targets for antiviral and therapeutic interventions.