Unraveling the Motions behind Enterovirus 71 Uncoating

Structural shifts in TolC facilitate Efflux-Mediated β-lactam resistance

Efflux-mediated β-lactam resistance is a major public health concern, reducing the effectiveness of β-lactam antibiotics against many bacteria. Structural analyses show the efflux protein TolC in Gram-negative bacteria acts as a channel for antibiotics, impacting bacterial susceptibility and virulence. This study examines β-lactam drug efflux mediated by TolC using experimental and computational methods. Molecular dynamics simulations of drug-free TolC reveal essential movements and key residues involved in TolC opening. A whole-gene-saturation mutagenesis assay, mutating each TolC residue and measuring fitness effects under β-lactam selection, is performed. Here we show the TolC-mediated efflux of three antibiotics: oxacillin, piperacillin, and carbenicillin. Steered molecular dynamics simulations identify general and drug-specific efflux mechanisms, revealing key positions at TolC's periplasmic entry affecting efflux motions. Our findings provide insights into TolC's structural dynamics, aiding the design of new antibiotics to overcome bacterial efflux mechanisms.

Identification and characterization of molecular entities differentially expressed in bacteria genome upon treatment with glyphosate shock

Antimetabolites developed from enzymes in the shikimate pathway are appealing targets. There are, however, certain unidentified molecular entities that show bacterial sensitivity to glyphosate shock. This study aims to identify the expression pattern of such entities following treatment with glyphosate shock and to characterize them structurally and functionally. Understanding such entities' catalytic structure and modulatory role guides the design and development of novel antibiotics. This study's functional profiling of 16S rRNA sequencing data and transcriptome analysis of glyphosate-exposedE. coli revealed that two genes were upregulated and twenty-eight were downregulated after glyphosate shock. We discovered the differential expression of some processes based on functional gene analysis, such as global and overview maps (4.2195 on average), carbohydrate metabolism (0.6858 on average), amino acid metabolism (0.5032 on average), and co-factor and vitamin metabolism (0.5032 on average) (0.2876 on average). After examining the two data sets, we discovered that some unidentified proteins were strongly expressed after glyphosate treatment. After examining the two datasets, we discovered a protein with no unique features expressed when treated with glyphosate. The Ecs2020 model looks to be the most stable in structural modeling investigations, while the catalytic residues sought in drug development are anticipated. Furthermore, biological processes and cellular component enrichment analysis revealed that the differentially expressed genes were strongly related to the trehalose manufacturing process and represented the cell membrane's outer membrane component. To estimate the functional gene content of soil sample metagenomics based on 16S rRNA, predictive functional analysis was done with R using the Tax4Fun2 package. On the other hand, transcriptome analysis was carried out using the R tool GEO2R. The results could be a good starting point for making new antibiotic medicines.

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

Identification of the Intrinsic Motions and Related Key Residues Responsible for the Twofold Channel Opening of Poliovirus Capsid by Using an Elastic Network Model Combined with an Internal Coordinate

Poliovirus (PV) is an infectious virus that causes poliomyelitis, which seriously threatens the health of children. The release of viral RNA is a key step of PV in host cell infection, and multiple lines of evidence have demonstrated that RNA release is initiated by the opening of the twofold channels of the PV capsid. However, the mechanism that controls the twofold channel opening is still not well understood. In addition, the channel opening motion of the recombinant PV capsid leads to the destruction of predominant neutralizing epitopes and thus hinders the capsid as a vaccine immunogen. Therefore, it is important to identify the intrinsic motions and the related key residues controlling the twofold channel opening for understanding the virus infection mechanism and developing capsid-based vaccines. In the present work, the width of the channel was selected as an internal coordinate directly related to the channel opening, and then the elastic network model (ENM) combined with the group theory were employed to extract the intrinsic motion modes that mostly contribute to the opening of the twofold channels. Our results show that the channel opening predominately induced by the breathing motion and the overall rotation of each protomer in the capsid. Then, an internal coordinate-based perturbation method was used to identify the key residues regulating the twofold channel opening of PV. The calculation results showed that the predicted key residues are mainly located at the twofold axes, the bottom of the canyons and the quasi threefold axes. Our study is helpful for better understanding the twofold channel opening mechanism and provides a potential target for preventing the opening of the channels, which is of great significance for PV vaccine design. The source code of this study is available at https://github.com/SJGLAB/CapsidKeyRes.git.

Heat-induced transitions of an empty minute virus of mice capsid in explicit water: all-atom MD simulation

The capsid-like structure of the virus-based protein nanoparticles (NPs) can serve as bionanomaterials, with applications in biomedicines and nanotechnology. Release of packaged material from these nanocontainers is associated with subtle conformational changes of the NP structure, which in vitro, is readily accomplished by heating. Characterizing the structural changes as a function of temperature may provide fresh insights into nanomaterial/antiviral strategies. Here, we have calculated heat induced changes in the properties of an empty minute virus of mice particle using large-scale ≈ 3.0 × 10⁶ all-atom molecular dynamics simulations. We focus on two heat induced structural changes of the NP, namely, dynamical transition (DT) and breathing transition (BT), both characterized by sudden and sharp change of measured parameters at temperatures, T_DT and T_BT, respectively. While DT is assessed by mean-square fluctuation of hydrogen atoms of the NP, BT is monitored through internal volume and permeation rate of water molecules through the NP. Both the transitions, resulting primarily from collective atomistic motion, are found to occur at temperatures widely separated from one another (T_BT>T_DT). The breathing motions, responsible for the translocation events of the packaged materials through the NP to kick off, are further probed by computing atomic resolution stresses from NVE simulations. Distribution of equilibrium atomistic stresses on the NP reveals a largely asymmetric nature and suggests structural breathing may actually represent large dynamic changes in the hotspot regions, far from the NP pores, which is in remarkable resemblance with recently conducted hydrogen-deuterium exchange coupled to mass spectrometry experiment. Communicated by Ramaswamy H. Sarma.

Insight into the Enterovirus A71: A review

Enterovirus A71 is a major causative pathogen of hand, foot and mouth disease. It has become a global public health threat, and is especially important for infants and young children in the Asian-Pacific countries. The enterovirus A71 is a non-enveloped virus of the Picornaviridae family having a single-stranded positive-sense RNA genome of about 7.4 kb which encodes the structural and nonstructural proteins. Currently there are no US FDA-approved vaccines or antiviral therapy available against enterovirus A71 infection. Although enterovirus A71 vaccines have been licenced in China, clinically approved vaccines for widespread vaccination programs are lacking. Substantial progress has recently been achieved on understanding the structure and function of enterovirus A71 proteins together with information on the viral genetic diversity and geographic distribution. The present review is intended to provide an overview on our current understanding of the molecular biology and epidemiology of enterovirus A71 which will aid the development of vaccines, therapeutics and other control strategies so as to bolster the preparedness for future enterovirus A71 outbreaks.

GTP Cyclohydrolase I as a Potential Drug Target: New Insights into Its Allosteric Modulation via Normal Mode Analysis

Guanosine triphosphate (GTP) cyclohydrolase I (GCH1) catalyzes the conversion of GTP into dihydroneopterin triphosphate (DHNP). DHNP is the first intermediate of the folate de novo biosynthesis pathway in prokaryotic and lower eukaryotic microorganisms and the tetrahydrobiopterin (BH4) biosynthesis pathway in higher eukaryotes. The de novo folate biosynthesis provides essential cofactors for DNA replication, cell division, and synthesis of key amino acids in rapidly replicating pathogen cells, such as Plasmodium falciparum (P. falciparum), a causative agent of malaria. In eukaryotes, the product of the BH4 biosynthesis pathway is essential for the production of nitric oxide and several neurotransmitter precursors. An increased copy number of the malaria parasite P. falciparum GCH1 gene has been reported to influence antimalarial antifolate drug resistance evolution, whereas mutations in the human GCH1 are associated with neuropathic and inflammatory pain disorders. Thus, GCH1 stands as an important and attractive drug target for developing therapeutics. The GCH1 intrinsic dynamics that modulate its activity remains unclear, and key sites that exert allosteric effects across the structure are yet to be elucidated. This study employed the anisotropic network model to analyze the intrinsic motions of the GCH1 structure alone and in complex with its regulatory partner protein. We showed that the GCH1 tunnel-gating mechanism is regulated by a global shear motion and an outward expansion of the central five-helix bundle. We further identified hotspot residues within sites of structural significance for the GCH1 intrinsic allosteric modulation. The obtained results can provide a solid starting point to design novel antineuropathic treatments for humans and novel antimalarial drugs against the malaria parasite P. falciparum GCH1 enzyme.

MDM-TASK-web: MD-TASK and MODE-TASK web server for analyzing protein dynamics

The web server, MDM-TASK-web, combines the MD-TASK and MODE-TASK software suites, which are aimed at the coarse-grained analysis of static and all-atom MD-simulated proteins, using a variety of non-conventional approaches, such as dynamic residue network analysis, perturbation-response scanning, dynamic cross-correlation, essential dynamics and normal mode analysis. Altogether, these tools allow for the exploration of protein dynamics at various levels of detail, spanning single residue perturbations and weighted contact network representations, to global residue centrality measurements and the investigation of global protein motion. Typically, following molecular dynamic simulations designed to investigate intrinsic and extrinsic protein perturbations (for instance induced by allosteric and orthosteric ligands, protein binding, temperature, pH and mutations), this selection of tools can be used to further describe protein dynamics. This may lead to the discovery of key residues involved in biological processes, such as drug resistance. The server simplifies the set-up required for running these tools and visualizing their results. Several scripts from the tool suites were updated and new ones were also added and integrated with 2D/3D visualization via the web interface. An embedded work-flow, integrated documentation and visualization tools shorten the number of steps to follow, starting from calculations to result visualization. The Django-powered web server (available at https://mdmtaskweb.rubi.ru.ac.za/) is compatible with all major web browsers. All scripts implemented in the web platform are freely available at https://github.com/RUBi-ZA/MD-TASK/tree/mdm-task-web and https://github.com/RUBi-ZA/MODE-TASK/tree/mdm-task-web.

The functional motions and related key residues behind the uncoating of coxsackievirus A16

The coxsackievirus A16 (CVA16) is a highly contagious virus that causes the hand, foot, and mouth disease, which seriously threatens the health of children. At present, there are still no available antiviral drugs or effective treatments against the infection of CVA16, and thus it is of great significance to develop anti-CVA16 vaccines. However, the intrinsic uncoating property of the capsid may destroy the neutralizing epitopes and influence its immunogenicity, which hinders the vaccine developments. In the present work, the functional-quantity-based elastic network model analysis method developed by our group was extended to combine with group theory to investigate the uncoating motions of the CVA16 capsid, and then the functionally key residues controlling the uncoating motions were identified by our functional-quantity-based perturbation method. Several motion modes encoded in the topological structure of the capsid were revealed to be responsible for the uncoating of CVA16 particle. These modes predominantly contribute to the fluctuation of the gyration radius of the capsid. Then, by using the perturbation method, four clusters of key sites involved in the uncoating motions were identified, whose perturbations induce significant changes in the fluctuation of the gyration radius. These key residues are mainly located at the 2-fold channels, the quasi 3-fold channels, the bottom of the canyons, and the inter-subunit interfaces around the 3-fold axes. Our studies are helpful for better understanding the uncoating mechanism of the CVA16 capsid and provide potential target sites to prevent the uncoating motions, which is valuable for the vaccine design against CVA16.

Normal Mode Analysis of KRas4B Reveals Partner Specific Dynamics

Ras GTPase interacts with its regulators and downstream effectors for its critical function in cellular signaling. Targeting the disrupted mechanisms in Ras-related human cancers requires understanding the distinct dynamics of these protein-protein interactions. We performed normal mode analysis (NMA) of KRas4B in wild-type or mutant monomeric and neurofibromin-1 (NF1), Son of Sevenless 1 (SOS1) or Raf-1 bound dimeric conformational states to reveal partner-specific dynamics of the protein. Gaussian network model (GNM) analysis showed that the known KRas4B lobes further partition into subdomains upon binding to its partners. Furthermore, KRas4B interactions with different partners suppress the flexibility in not only their binding sites but also distant residues in the allosteric lobe in a partner-specific way. The conformational changes can be driven by intrinsic residue fluctuations of the open state KRas4B-GDP, as we illustrated with anisotropic network model (ANM) analysis. The allosteric paths connecting the nucleotide binding residues to the allosteric site at α3-L7 portray differences in the inactive and active states. These findings help in understanding the partner-specific KRas4B dynamics, which could be utilized for therapeutic targeting.

A Coarse-Grained Methodology Identifies Intrinsic Mechanisms That Dissociate Interacting Protein Pairs

We address the problem of triggering dissociation events between proteins that have formed a complex. We have collected a set of 25 non-redundant, functionally diverse protein complexes having high-resolution three-dimensional structures in both the unbound and bound forms. We unify elastic network models with perturbation response scanning (PRS) methodology as an efficient approach for predicting residues that have the propensity to trigger dissociation of an interacting protein pair, using the three-dimensional structures of the bound and unbound proteins as input. PRS reveals that while for a group of protein pairs, residues involved in the conformational shifts are confined to regions with large motions, there are others where they originate from parts of the protein unaffected structurally by binding. Strikingly, only a few of the complexes have interface residues responsible for dissociation. We find two main modes of response: In one mode, remote control of dissociation in which disruption of the electrostatic potential distribution along protein surfaces play the major role; in the alternative mode, mechanical control of dissociation by remote residues prevail. In the former, dissociation is triggered by changes in the local environment of the protein, e.g., pH or ionic strength, while in the latter, specific perturbations arriving at the controlling residues, e.g., via binding to a third interacting partner is required for decomplexation. We resolve the observations by relying on an electromechanical coupling model which reduces to the usual elastic network result in the limit of the lack of coupling. We validate the approach by illustrating the biological significance of top residues selected by PRS on select cases where we show that the residues whose perturbation leads to the observed conformational changes correspond to either functionally important or highly conserved residues in the complex.

Integrated Computational Approaches and Tools forAllosteric Drug Discovery

Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.

Role of protein-protein interactions in allosteric drug design for DNA methyltransferases

DNA methyltransferases (DNMTs) not only play key roles in epigenetic gene regulation, but also serve as emerging targets for several diseases, especially for cancers. Due to the multi-domains of DNMT structures, targeting allosteric sites of protein-protein interactions (PPIs) is becoming an attractive strategy in epigenetic drug discovery. This chapter aims to review the major contemporary approaches utilized for the drug discovery based on PPIs in different dimensions, from the enumeration of allosteric mechanism to the identification of allosteric pockets. These include the construction of protein structure networks (PSNs) based on molecular dynamics (MD) simulations, performing elastic network models (ENMs) and perturbation response scanning (PRS) calculation, the sequence-based conservation and coupling analysis, and the allosteric pockets identification. Furthermore, we complement this methodology by highlighting the role of computational approaches in promising practical applications for the computer-aided drug design, with special focus on two DNMTs, namely, DNMT1 and DNMT3A.

MODE-TASK: large-scale protein motion tools

Summary

MODE-TASK, a novel and versatile software suite, comprises Principal Component Analysis, Multidimensional Scaling, and t-Distributed Stochastic Neighbor Embedding techniques using Molecular Dynamics trajectories. MODE-TASK also includes a Normal Mode Analysis tool based on Anisotropic Network Model so as to provide a variety of ways to analyse and compare large-scale motions of protein complexes for which long MD simulations are prohibitive. Beside the command line function, a GUI has been developed as a PyMOL plugin.

Availability and implementation

MODE-TASK is open source, and available for download from https://github.com/RUBi-ZA/MODE-TASK. It is implemented in Python and C++. It is compatible with Python 2.x and Python 3.x and can be installed by Conda.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.