Eigenvector centrality for characterization of protein allosteric pathways

Subtle Changes at the RBD/hACE2 Interface During SARS-CoV-2 Variant Evolution: A Molecular Dynamics Study

The SARS-CoV-2 Omicron variants show different behavior compared to the previous variants, especially with respect to the Delta variant, which promotes a lower morbidity despite being much more contagious. In this perspective, we performed molecular dynamics (MD) simulations of the different spike RBD/hACE2 complexes corresponding to the WT, Delta and four Omicron variants. Carrying out a comprehensive analysis of residue interactions within and between the two partners allowed us to draw the profile of each variant by using complementary methods (PairInt, hydrophobic potential, contact PCA). PairInt calculations highlighted the residues most involved in electrostatic interactions, which make a strong contribution to the binding with highly stable interactions between spike RBD and hACE2. Apolar contacts made a substantial and complementary contribution in Omicron with the detection of two hydrophobic patches. Contact networks and cross-correlation matrices were able to detect subtle changes at point mutations as the S375F mutation occurring in all Omicron variants, which is likely to confer an advantage in binding stability. This study brings new highlights on the dynamic binding of spike RBD to hACE2, which may explain the final persistence of Omicron over Delta.

Functionally important residues from graph analysis of coevolved dynamic couplings

The relationship between protein dynamics and function is essential for understanding biological processes and developing effective therapeutics. Functional sites within proteins are critical for activities such as substrate binding, catalysis, and structural changes. Existing computational methods for the predictions of functional residues are trained on sequence, structural, and experimental data, but they do not explicitly model the influence of evolution on protein dynamics. This overlooked contribution is essential as it is known that evolution can fine-tune protein dynamics through compensatory mutations either to improve the proteins' performance or diversify its function while maintaining the same structural scaffold. To model this critical contribution, we introduce DyNoPy, a computational method that combines residue coevolution analysis with molecular dynamics simulations, revealing hidden correlations between functional sites. DyNoPy constructs a graph model of residue-residue interactions, identifies communities of key residue groups, and annotates critical sites based on their roles. By leveraging the concept of coevolved dynamical couplings-residue pairs with critical dynamical interactions that have been preserved during evolution-DyNoPy offers a powerful method for predicting and analysing protein evolution and dynamics. We demonstrate the effectiveness of DyNoPy on SHV-1 and PDC-3, chromosomally encoded β-lactamases linked to antibiotic resistance, highlighting its potential to inform drug design and address pressing healthcare challenges.

Campesterol and dithymoquinone as a potent inhibitors of SARS cov-2 main proteases-promising drug candidates for targeting its novel variants

The sudden outbreak of the COVID-19 pandemic has currently taken approximately 2.4 million lives, with no specific medication and fast-tracked tested vaccines for prevention. These vaccines have their own adverse effects, which have severely affected the global healthcare system. The discovery of the main protease structure of coronavirus (Mpro/Clpro) has resulted in the identification of compounds having antiviral potential, especially from the herbal system. In this study, the computer-associated drug design tools were utilised to analyze the reported phytoconstituents of Nigella sativa for their antiviral activity against the main protease. Fifty-eight compounds were subjected to pharmacological parameter analysis to determine their lead likeness in comparison to the standard drugs (chloroquine and nirmatrelvir) used in the treatment of SARS-CoV-2. Nearly 31 compounds were docked against five different SARS-CoV-2 main proteases, and all compounds showed better binding affinity and inhibition constant against the proteases. However, dithymoquinone and campesterol displayed the best binding scores and hence were further subjected to dynamics and MMPBSA study for 100 ns. The stability analysis shows that dithymoquinone and campesterol show less variation in fluctuation in residues compared to standard complexes. Moreover, dithymoquinone exhibited higher binding affinity and favorable interaction followed by campesterol as compared to the standard drug. The in silico computational analysis provides a promising hit for regulating the main proteases activity.

MMFuncPhos: A Multi-Modal Learning Framework for Identifying Functional Phosphorylation Sites and Their Regulatory Types

Protein phosphorylation plays a crucial role in regulating a wide range of biological processes, and its dysregulation is strongly linked to various diseases. While many phosphorylation sites have been identified so far, their functionality and regulatory effects are largely unknown. Here, a deep learning model MMFuncPhos, based on a multi-modal deep learning framework, is developed to predict functional phosphorylation sites. MMFuncPhos outperforms existing functional phosphorylation site prediction approaches. EFuncType is further developed based on transfer learning to predict whether phosphorylation of a residue upregulates or downregulates enzyme activity for the first time. The functional phosphorylation sites predicted by MMFuncPhos and the regulatory types predicted by EFuncType align with experimental findings from several newly reported protein phosphorylation studies. The study contributes to the understanding of the functional regulatory mechanism of phosphorylation and provides valuable tools for precision medicine, enzyme engineering, and drug discovery. For user convenience, these two prediction models are integrated into a web server which can be accessed at http://pkumdl.cn:8000/mmfuncphos.

Analyzing Information Exchange in Parkinson's Disease via Eigenvector Centrality: A Source-Level Magnetoencephalography Study

Background: Parkinson's disease (PD) is a progressive neurodegenerative disorder that manifests through motor and non-motor symptoms. Understanding the alterations in brain connectivity associated with PD remains a challenge that is crucial for enhancing diagnosis and clinical management. Methods: This study utilized Magnetoencephalography (MEG) to investigate brain connectivity in PD patients compared to healthy controls (HCs) by applying eigenvector centrality (EC) measures across different frequency bands. Results: Our findings revealed significant differences in EC between PD patients and HCs in the alpha (8-12 Hz) and beta (13-30 Hz) frequency bands. To go into further detail, in the alpha frequency band, PD patients in the frontal lobe showed higher EC values compared to HCs. Additionally, we found statistically significant correlations between EC measures and clinical impairment scores (UPDRS-III). Conclusions: The proposed results suggest that MEG-derived EC measures can reveal important alterations in brain connectivity in PD, potentially serving as biomarkers for disease severity.

Delving into Macrolide Binding Affinities and Associated Structural Modulations in Erythromycin Esterase C: Insights into the Venus Flytrap Mechanism

Since their inception in antibacterial therapy, macrolide-based antibiotics have significantly shaped the evolutionary pathways of pathogenic bacteria, driving them to develop diverse antimicrobial resistance (AMR) mechanisms. Among these, macrolide esterase, commonly referred to as erythromycin esterase, emerged as a critical defense mechanism, enabling bacteria to detoxify macrolides by hydrolyzing the macrolactone ring within the bacterial cell. In this study, we delve into the intricate interactions and conformational dynamics of erythromycin esterase C (EreC), a key member of the Ere enzyme family. We have focused on three FDA-approved and widely prescribed macrolides─erythromycin, clarithromycin, and azithromycin─by employing classical molecular dynamics, absolute binding free energy calculations, and 2D well-tempered metadynamics simulations to explore their interactions with EreC. To estimate the absolute binding free energies, we have used the recently developed and robust "Streamlined Alchemical Free Energy Perturbation (SAFEP)" protocol. The results from our molecular dynamics simulations and advanced analyses portrayed the crucial role of hydrophobic interactions within the macrolide binding cleft of EreC, along with the significant influence of the minor lobe in facilitating overall structural fluctuation. In silico alanine scanning identified top three hydrophobic residues, i.e., PHE248, MET333, and PHE344, responsible for macrolide binding inside that cleft. According to the free energy calculations, azithromycin and clarithromycin showed greater binding affinities toward EreC than the parent macrolide erythromycin. Moreover, 2D metadynamics simulations along with graph theory-based eigenvector centrality analyses revealed a metastable "semiopen" state during the hypothesized "active loop closure" of the EreC protein triggered by subtle conformational changes of an important histidine residue, HIS289, upon macrolide capture, drawing a fascinating parallel to the renowned "Venus flytrap" mechanism.

Deep-HPI-pred: An R-Shiny applet for network-based classification and prediction of Host-Pathogen protein-protein interactions

Host-pathogen interactions (HPIs) are vital in numerous biological activities and are intrinsically linked to the onset and progression of infectious diseases. HPIs are pivotal in the entire lifecycle of diseases: from the onset of pathogen introduction, navigating through the mechanisms that bypass host cellular defenses, to its subsequent proliferation inside the host. At the heart of these stages lies the synergy of proteins from both the host and the pathogen. By understanding these interlinking protein dynamics, we can gain crucial insights into how diseases progress and pave the way for stronger plant defenses and the swift formulation of countermeasures. In the framework of current study, we developed a web-based R/Shiny app, Deep-HPI-pred, that uses network-driven feature learning method to predict the yet unmapped interactions between pathogen and host proteins. Leveraging citrus and CLas bacteria training datasets as case study, we spotlight the effectiveness of Deep-HPI-pred in discerning Protein-protein interaction (PPIs) between them. Deep-HPI-pred use Multilayer Perceptron (MLP) models for HPI prediction, which is based on a comprehensive evaluation of topological features and neural network architectures. When subjected to independent validation datasets, the predicted models consistently surpassed a Matthews correlation coefficient (MCC) of 0.80 in host-pathogen interactions. Remarkably, the use of Eigenvector Centrality as the leading topological feature further enhanced this performance. Further, Deep-HPI-pred also offers relevant gene ontology (GO) term information for each pathogen and host protein within the system. This protein annotation data contributes an additional layer to our understanding of the intricate dynamics within host-pathogen interactions. In the additional benchmarking studies, the Deep-HPI-pred model has proven its robustness by consistently delivering reliable results across different host-pathogen systems, including plant-pathogens (accuracy of 98.4% and 97.9%), human-virus (accuracy of 94.3%), and animal-bacteria (accuracy of 96.6%) interactomes. These results not only demonstrate the model's versatility but also pave the way for gaining comprehensive insights into the molecular underpinnings of complex host-pathogen interactions. Taken together, the Deep-HPI-pred applet offers a unified web service for both identifying and illustrating interaction networks. Deep-HPI-pred applet is freely accessible at its homepage: https://cbi.gxu.edu.cn/shiny-apps/Deep-HPI-pred/ and at github: https://github.com/tahirulqamar/Deep-HPI-pred.

Multi-Epitopic Peptide Vaccine Against Newcastle Disease Virus: Molecular Dynamics Simulation and Experimental Validation

BACKGROUND

Newcastle disease virus (NDV) is a highly contagious and economically devastating pathogen affecting poultry worldwide, leading to significant losses in the poultry industry. Despite existing vaccines, outbreaks continue to occur, highlighting the need for more effective vaccination strategies. Developing a multi-epitopic peptide vaccine offers a promising approach to enhance protection against NDV.

OBJECTIVES

Here, we aimed to design and evaluate a multi-epitopic vaccine against NDV using molecular dynamics (MD) simulation.

METHODOLOGY

We retrieved NDV sequences for the fusion (F) protein and hemagglutinin-neuraminidase (HN) protein. Subsequently, B-cell and T-cell epitopes were predicted. The top potential epitopes were utilized to design the vaccine construct, which was subsequently docked against chicken TLR4 and MHC1 receptors to assess the immunological response. The resulting docked complex underwent a 1 microsecond (1000 ns) MD simulation. For experimental evaluation, the vaccine's efficacy was assessed in mice and chickens using a controlled study design, where animals were randomly divided into groups receiving either a local ND vaccine or the peptide vaccine or a control treatment.

RESULTS

The 40 amino acid peptide vaccine demonstrated strong binding affinity and stability within the TLR4 and MHC1 receptor-peptide complexes. The root mean square deviation of peptide vaccine and TLR4 receptor showed rapid stabilization after an initial repositioning. The root mean square fluctuation revealed relatively low fluctuations (below 3 Å) for the TLR4 receptor, while the peptide exhibited higher fluctuations. The overall binding energy of the peptide vaccine with TLR4 and MHC1 receptors amounted to -15.7 kcal·mol-1 and -36.8 kcal·mol-1, respectively. For experimental evaluations in mice and chicken, the peptide vaccine was synthesized using services of GeneScript Biotech® (Singapore) PTE Limited. Experimental evaluations showed a significant immune response in both mice and chickens, with the vaccine eliciting robust antibody production, as evidenced by increasing HI titers over time. Statistical analysis was performed using an independent t-test with Type-II error to compare the groups, calculating the p-values to determine the significance of the immune response between different groups.

CONCLUSIONS

Multi-epitopic peptide vaccine has demonstrated a good immunological response in natural hosts.

Essential cancer protein identification using graph-based random walk with restart

Protein-protein interaction (PPI) network analysis holds significant promise for cancer diagnosis and drug target identification. This paper introduces a novel random walk-based method called essential cancer protein identification using graph-based random walk with restart (EPI-GBRWR) to address this gap. This proposed method incorporates local and global topological features of proteins, enhancing the accuracy of essential protein identification in PPI networks. Starting with meticulous preprocessing of cancer gene datasets from NCBI, including breast, lung, colorectal, and ovarian cancers, and identifying a core set of common genes. The proposed method constructs PPI networks to capture complex protein interactions from these common cancer genes. Topological analysis, including a centrality measures matrix, is generated to perform the analysis to identify essential nodes. The study revealed that 40 essential proteins among breast, colorectal, lung and ovarian cancer showcase the potency of integrative methodologies in unravelling cancer complexity, signalling a transformative era in cancer research and treatment. The strength of the findings from the study has direct clinical relevance in cancer diseases. It contributes to the field of precision medicine to guide personalized treatment strategies.

Differential network analysis reveals the key role of the ECM-receptor pathway in α-particle-induced malignant transformation

Space particle radiation is a major environmental factor in spaceflight, and it is known to cause body damage and even trigger cancer, but with unknown molecular etiologies. To examine these causes, we developed a systems biology approach by focusing on the co-expression network analysis of transcriptomics profiles obtained from single high-dose (SE) and multiple low-dose (ME) α-particle radiation exposures of BEAS-2B human bronchial epithelial cells. First, the differential network and pathway analysis based on the global network and the core modules showed that genes in the ME group had higher enrichment for the extracellular matrix (ECM)-receptor interaction pathway. Then, collagen gene COL1A1 was screened as an important gene in the ME group assessed by network parameters and an expression study of lung adenocarcinoma samples. COL1A1 was found to promote the emergence of the neoplastic characteristics of BEAS-2B cells by both in vitro experimental analyses and in vivo immunohistochemical staining. These findings suggested that the degree of malignant transformation of cells in the ME group was greater than that of the SE, which may be caused by the dysregulation of the ECM-receptor pathway.

Dynamically driven correlations in elastic net models reveal sequence of events and causality in proteins

An explicit analytic solution is given for the Langevin equation applied to the Gaussian Network Model of a protein subjected to both a random and a deterministic periodic force. Synchronous and asynchronous components of time correlation functions are derived and an expression for phase differences in the time correlations of residue pairs is obtained. The synchronous component enables the determination of dynamic communities within the protein structure. The asynchronous component reveals causality, where the time correlation function between residues i and j differs depending on whether i is observed before j or vice versa, resulting in directional information flow. Driver and driven residues in the allosteric process of cyclophilin A and human NAD-dependent isocitrate dehydrogenase are determined by a perturbation-scanning technique. Factors affecting phase differences between fluctuations of residues, such as network topology, connectivity, and residue centrality, are identified. Within the constraints of the isotropic Gaussian Network Model, our results show that asynchronicity increases with viscosity and distance between residues, decreases with increasing connectivity, and decreases with increasing levels of eigenvector centrality.

Conformational Modulation of a Mobile Loop Controls Catalysis in the (βα)8-Barrel Enzyme of Histidine Biosynthesis HisF

The overall significance of loop motions for enzymatic activity is generally accepted. However, it has largely remained unclear whether and how such motions can control different steps of catalysis. We have studied this problem on the example of the mobile active site β1α1-loop (loop1) of the (βα)8-barrel enzyme HisF, which is the cyclase subunit of imidazole glycerol phosphate synthase. Loop1 variants containing single mutations of conserved amino acids showed drastically reduced rates for the turnover of the substrates N'-[(5'-phosphoribulosyl) formimino]-5-aminoimidazole-4-carboxamide ribonucleotide (PrFAR) and ammonia to the products imidazole glycerol phosphate (ImGP) and 5-aminoimidazole-4-carboxamide-ribotide (AICAR). A comprehensive mechanistic analysis including stopped-flow kinetics, X-ray crystallography, NMR spectroscopy, and molecular dynamics simulations detected three conformations of loop1 (open, detached, closed) whose populations differed between wild-type HisF and functionally affected loop1 variants. Transient stopped-flow kinetic experiments demonstrated that wt-HisF binds PrFAR by an induced-fit mechanism whereas catalytically impaired loop1 variants bind PrFAR by a simple two-state mechanism. Our findings suggest that PrFAR-induced formation of the closed conformation of loop1 brings active site residues in a productive orientation for chemical turnover, which we show to be the rate-limiting step of HisF catalysis. After the cyclase reaction, the closed loop conformation is destabilized, which favors the formation of detached and open conformations and hence facilitates the release of the products ImGP and AICAR. Our data demonstrate how different conformations of active site loops contribute to different catalytic steps, a finding that is presumably of broad relevance for the reaction mechanisms of (βα)8-barrel enzymes and beyond.

Spatiotemporal and Multilayer Trade Network Patterns of the Global Cobalt Cycle

Recent years have witnessed increasing attempts to track trade flows of critical materials across world regions and along the life cycle for renewable energy and the low carbon transition. Previous studies often had limited spatiotemporal coverage, excluded end-use products, and modeled different life cycle stages as single-layer networks. Here, we integrated material flow analysis and complex network analysis into a multilayer framework to characterize the spatiotemporal and multilayer trade network patterns of the global cobalt cycle from 1988 to 2020. We found substantial growth and notable structural changes in global cobalt trade over the past 30 years. China, Germany, and the United States play pivotal roles in different layers and stages of the global cobalt cycle. The interlayer relationships among alloys, batteries, and materials are robust and continually strengthening, indicating a trend toward synergistic trade. However, cobalt ore-exporting countries are highly concentrated and rarely involved in later life cycle stages, resulting in the weakest relationship between the ore layer and other layers. This causes fluctuations and uncertainty in the global cobalt trade. Our model, linking industrial ecology, supply chain analysis, and network analysis, can be extended to other materials that are critical for the future green transition.

Insights into Allosteric Inhibition of the AcrB Efflux Pump: Role of Distinct Binding Pockets, Protomer Preferences, and Crosstalk Disruption

AcrB, a key component in bacterial efflux processes, exhibits distinct binding pockets that influence inhibitor interactions. In addition to the well-known distal binding pocket within the periplasmic domain, a noteworthy pocket amidst the transmembrane (TM) helices serves as an alternate binding site for inhibitors. The bacterial efflux mechanism involves a pivotal functional rotation of the TM protein, inducing conformational changes in each protomer and propelling drugs toward the outer membrane domain. Surprisingly, inhibitors binding to the TM domain display a preference for L protomers over T protomers. Metadynamics simulations elucidate that Lys940 in the TM domain of AcrB can adopt two conformations in L protomers, whereas the energy barrier for such transitions is higher in T protomers. This phenomenon results in stable inhibitor binding in l protomers. Upon a detailed analysis of unbinding pathways using random accelerated molecular dynamics and umbrella sampling, we have identified three distinct routes for ligand exit from the allosteric site, specifically involving regions within the TM domains─TM4, TM5, and TM10. To explore allosteric crosstalk, we focused on the following key residues: Val452 from the TM domain and Ala831 from the porter domain. Surprisingly, our findings reveal that inhibitor binding disrupts this communication. The shortest path connecting Val452 and Ala831 increases upon inhibitor binding, suggesting sabotage of the natural interdomain communication dynamics. This result highlights the intricate interplay between inhibitor binding and allosteric signaling within our studied system.

Exploring transmembrane allostery in the MexB: DB08385 variant as a promising inhibitor-like candidate against Pseudomonas aeruginosa antibiotic resistance: a computational study

Pseudomonas aeruginosa, a formidable pathogen renowned for its antimicrobial resistance, poses a significant threat to immunocompromised individuals. In this regard, the MexAB-OprM efflux pump acts as a pivotal line of defense by extruding antimicrobials from bacterial cells. The inner membrane homotrimeric protein MexB captures antibiotics and translocates them into the outer membrane OprM channel protein connected through the MexA adaptor protein. Despite extensive efforts, competitive inhibitors targeting the tight (T) protomer of the MexB protein have not received FDA approval for medical use. Over the past few years, allosteric inhibitors have become popular as alternatives to the classical competitive inhibitor-based approach because of their higher specificity, lower dosage, and reduced toxicological effects. Hence, in this study, we unveiled the existence of a transmembrane allosteric binding pocket of MexB inspired by the recent discovery of an important allosteric inhibitor, BDM88855, for the homolog AcrB protein. While repurposing BDM88855 proved ineffective in controlling the MexB loose (L) protomer, our investigation identified a promising alternative: a chlorine-containing variant of DB08385 (2-Cl DB08385 or Variant 1). Molecular dynamics simulations, including binding free energy estimation coupled with heterogeneous dielectric implicit membrane model (implicit-membrane MM/PBSA), interaction entropy (IE) analysis and potential of mean force (PMF) calculation, demonstrated Variant 1's superior binding affinity to the transmembrane pocket, displaying the highest energy barrier in the ligand unbinding process. To elucidate the allosteric crosstalk between the transmembrane and porter domain of MexB, we employed the 'eigenvector centrality' measure in the linear mutual information obtained from the protein correlation network. Notably, this study confirmed the presence of an allosteric transmembrane site in the MexB L protomer. In addition to this, Variant 1 emerged as a potent regulator of allosteric crosstalk, inducing an 'O-L intermediate state' in the MexB L protomer. This induced state might hold the potential to diminish substrate intake into the access pocket, leading to the ineffective efflux of antibiotics.

Topological and spectral properties of random digraphs

We investigate some topological and spectral properties of Erdős-Rényi (ER) random digraphs of size n and connection probability p, D(n,p). In terms of topological properties, our primary focus lies in analyzing the number of nonisolated vertices V_{x}(D) as well as two vertex-degree-based topological indices: the Randić index R(D) and sum-connectivity index χ(D). First, by performing a scaling analysis, we show that the average degree 〈k〉 serves as a scaling parameter for the average values of V_{x}(D), R(D), and χ(D). Then, we also state expressions relating the number of arcs, largest eigenvalue, and closed walks of length 2 to (n,p), the parameters of ER random digraphs. Concerning spectral properties, we observe that the eigenvalue distribution converges to a circle of radius sqrt[np(1-p)]. Subsequently, we compute six different invariants related to the eigenvalues of D(n,p) and observe that these quantities also scale with sqrt[np(1-p)]. Additionally, we reformulate a set of bounds previously reported in the literature for these invariants as a function (n,p). Finally, we phenomenologically state relations between invariants that allow us to extend previously known bounds.

The synergy between the PscC subunits for electron transfer to the P840 special pair in Chlorobaculum tepidum

The primary photochemical reaction of photosynthesis in green sulfur bacteria occurs in the homodimer PscA core proteins by a special chlorophyll pair. The light induced excited state of the special pair producing P840+ is rapidly reduced by electron transfer from one of the two PscC subunits. Molecular dynamics (MD) simulations are combined with bioinformatic tools herein to provide structural and dynamic insight into the complex between the two PscA core proteins and the two PscC subunits. The microscopic dynamic model involves extensive sampling at atomic resolution and at a cumulative time-scale of 22µs and reveals well defined protein-protein interactions. The membrane complex is composed of the two PscA and the two PscC subunits and macroscopic connections are revealed within a putative electron transfer pathway from the PscC subunit to the special pair P840 located within the PscA subunits. Our results provide a structural basis for understanding the electron transport to the homodimer RC of the green sulfur bacteria. The MD based approach can provide the basis to further probe the PscA-PscC complex dynamics and observe electron transfer therein at the quantum level. Furthermore, the transmembrane helices of the different PscC subunits exert distinct dynamics in the complex.

Probing allosteric communication with combined molecular dynamics simulations and network analysis

Understanding the allosteric mechanisms within biomolecules involved in diseases is of paramount importance for drug discovery. Indeed, characterizing communication pathways and critical hotspots in signal transduction can guide a rational approach to leverage allosteric modulation for therapeutic purposes. While the atomistic signatures of allosteric processes are difficult to determine experimentally, computational methods can be a remarkable resource. Network analysis built on Molecular Dynamics simulation data is particularly suited in this respect and is gradually becoming of routine use. Herein, we collect the recent literature in the field, discussing different aspects and available options for network construction and analysis. We further highlight interesting refinements and extensions, eventually providing our perspective on this topic.

scDrugPrio: a framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases

BACKGROUND

Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs.

METHODS

Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs.

RESULTS

scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment.

CONCLUSIONS

We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package ( https://github.com/SDTC-CPMed/scDrugPrio ).

What is allosteric regulation? Exploring the exceptions that prove the rule!

"Allosteric" was first introduced to mean the other site (i.e., a site distinct from the active or orthosteric site), an adjective for "regulation" to imply a regulatory outcome resulting from ligand binding at another site. That original idea outlines a system with two ligand-binding events at two distinct locations on a macromolecule (originally a protein system), which defines a four-state energy cycle. An allosteric energy cycle provides a quantifiable allosteric coupling constant and focuses our attention on the unique properties of the four equilibrated protein complexes that constitute the energy cycle. Because many observed phenomena have been referenced as "allosteric regulation" in the literature, the goal of this work is to use literature examples to explore which systems are and are not consistent with the two-ligand thermodynamic energy cycle-based definition of allosteric regulation. We emphasize the need for consistent language so comparisons can be made among the ever-increasing number of allosteric systems. Building on the mutually exclusive natures of an energy cycle definition of allosteric regulation versus classic two-state models, we conclude our discussion by outlining how the often-proposed Rube-Goldberg-like mechanisms are likely inconsistent with an energy cycle definition of allosteric regulation.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study, we employ this paradigm to answer a basic question: in enzyme superfamilies, where the catalytic mechanism, active sites, and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as representatives of the conserved protein tyrosine phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of the catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

CardiOmics signatures reveal therapeutically actionable targets and drugs for cardiovascular diseases

Cardiovascular diseases are the leading cause of death worldwide, with heart failure being a complex condition that affects millions of individuals. Single-nucleus RNA sequencing has recently emerged as a powerful tool for unraveling the molecular mechanisms behind cardiovascular diseases. This cutting-edge technology enables the identification of molecular signatures, intracellular networks, and spatial relationships among cardiac cells, including cardiomyocytes, mast cells, lymphocytes, macrophages, lymphatic endothelial cells, endocardial cells, endothelial cells, epicardial cells, adipocytes, fibroblasts, neuronal cells, pericytes, and vascular smooth muscle cells. Despite these advancements, the discovery of essential therapeutic targets and drugs for precision cardiology remains a challenge. To bridge this gap, we conducted comprehensive in silico analyses of single-nucleus RNA sequencing data, functional enrichment, protein interactome network, and identification of the shortest pathways to physiological phenotypes. This integrated multi-omics analysis generated CardiOmics signatures, which allowed us to pinpoint three therapeutically actionable targets (ADRA1A1, PPARG, and ROCK2) and 15 effective drugs, including adrenergic receptor agonists, adrenergic receptor antagonists, norepinephrine precursors, PPAR receptor agonists, and Rho-associated kinase inhibitors, involved in late-stage cardiovascular disease clinical trials.

Enhancing Exchange-Traded Fund Price Predictions: Insights from Information-Theoretic Networks and Node Embeddings

This study presents a novel approach to predicting price fluctuations for U.S. sector index ETFs. By leveraging information-theoretic measures like mutual information and transfer entropy, we constructed threshold networks highlighting nonlinear dependencies between log returns and trading volume rate changes. We derived centrality measures and node embeddings from these networks, offering unique insights into the ETFs' dynamics. By integrating these features into gradient-boosting algorithm-based models, we significantly enhanced the predictive accuracy. Our approach offers improved forecast performance for U.S. sector index futures and adds a layer of explainability to the existing literature.

High connectivity and human movement limits the impact of travel time on infectious disease transmission

The speed of spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during the coronavirus disease 2019 (COVID-19) pandemic highlights the importance of understanding how infections are transmitted in a highly connected world. Prior to vaccination, changes in human mobility patterns were used as non-pharmaceutical interventions to eliminate or suppress viral transmission. The rapid spread of respiratory viruses, various intervention approaches, and the global dissemination of SARS-CoV-2 underscore the necessity for epidemiological models that incorporate mobility to comprehend the spread of the virus. Here, we introduce a metapopulation susceptible-exposed-infectious-recovered model parametrized with human movement data from 340 cities in China. Our model replicates the early-case trajectory in the COVID-19 pandemic. We then use machine learning algorithms to determine which network properties best predict spread between cities and find travel time to be most important, followed by the human movement-weighted personalized PageRank. However, we show that travel time is most influential locally, after which the high connectivity between cities reduces the impact of travel time between individual cities on transmission speed. Additionally, we demonstrate that only significantly reduced movement substantially impacts infection spread times throughout the network.

Catalytic Effects of Active Site Conformational Change in the Allosteric Activation of Imidazole Glycerol Phosphate Synthase

Imidazole glycerol phosphate synthase (IGPS) is a class-I glutamine amidotransferase (GAT) that hydrolyzes glutamine. Ammonia is produced and transferred to a second active site, where it reacts with N1-(5'-phosphoribosyl)-formimino-5-aminoimidazole-4-carboxamide ribonucleotide (PrFAR) to form precursors to purine and histidine biosynthesis. Binding of PrFAR over 25 Å away from the active site increases glutaminase efficiency by ∼4500-fold, primarily altering the glutamine turnover number. IGPS has been the focus of many studies on allosteric communication; however, atomic details for how the glutamine hydrolysis rate increases in the presence of PrFAR are lacking. We present a density functional theory study on 237-atom active site cluster models of IGPS based on crystallized structures representing the inactive and allosterically active conformations and investigate the multistep reaction leading to thioester formation and ammonia production. The proposed mechanism is supported by similar, well-studied enzyme mechanisms, and the corresponding energy profile is consistent with steady-state kinetic studies of PrFAR + IGPS. Additional active site models are constructed to examine the relationship between active site structural change and transition-state stabilization via energy decomposition schemes. The results reveal that the inactive IGPS conformation does not provide an adequately formed oxyanion hole structure and that repositioning of the oxyanion strand relative to the substrate is vital for a catalysis-competent oxyanion hole, with or without the hVal51 dihedral flip. These findings are valuable for future endeavors in modeling the IGPS allosteric mechanism by providing insight into the atomistic changes required for rate enhancement that can inform suitable reaction coordinates for subsequent investigations.

Second quantization of many-body dispersion interactions for chemical and biological systems

The many-body dispersion (MBD) framework is a successful approach for modeling the long-range electronic correlation energy and optical response of systems with thousands of atoms. Inspired by field theory, here we develop a second-quantized MBD formalism (SQ-MBD) that recasts a system of atomic quantum Drude oscillators in a Fock-space representation. SQ-MBD provides: (i) tools for projecting observables (interaction energy, transition multipoles, polarizability tensors) on coarse-grained representations of the atomistic system ranging from single atoms to large structural motifs, (ii) a quantum-information framework to analyze correlations and (non)separability among fragments in a given molecular complex, and (iii) a path toward the applicability of the MBD framework to molecular complexes with even larger number of atoms. The SQ-MBD approach offers conceptual insights into quantum fluctuations in molecular systems and enables direct coupling of collective plasmon-like MBD degrees of freedom with arbitrary environments, providing a tractable computational framework to treat dispersion interactions and polarization response in intricate systems.

Insight into the conserved structural dynamics of the C-terminus of mammal PrPC identifies structural core and possible structural role of pharmacological chaperones

Misfolding of the prion protein is central to prion disease aetiology. Although understanding the dynamics of the native fold helps to decipher the conformational conversion mechanism, a complete depiction of distal but coupled prion protein sites common across species is lacking. To fill this gap, we used normal mode analysis and network analysis to examine a collection of prion protein structures deposited on the protein data bank. Our study identified a core of conserved residues that sustains the connectivity across the C-terminus of the prion protein. We propose how a well-characterized pharmacological chaperone may stabilize the fold. Also, we provide insight into the effect on the native fold of initial misfolding pathways identified by others using kinetics studies.

scDrugPrio: A framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases

Background

Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs.

Methods

Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs.

Results

scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment.

Conclusion

We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package (https://github.com/SDTC-CPMed/scDrugPrio).

A Proofreading Mutation with an Allosteric Effect Allows a Cluster of SARS-CoV-2 Viruses to Rapidly Evolve

The RNA-dependent RNA polymerase of the severe acute respiratory syndrome coronavirus 2 virus is error prone, with errors being corrected by the exonuclease (NSP14) proofreading mechanism. However, the mutagenesis and subsequent evolutionary trajectory of the virus is mediated by the delicate interplay of replicase fidelity and environmental pressures. Here, we have shown that a single, distal mutation (F60S) in NSP14 can have a profound impact upon proofreading with an increased accumulation of mutations and elevated evolutionary rate being observed. Understanding the implications of these changes is crucial, as these underlying mutational processes may have important implications for understanding the population-wide evolution of the virus. This study underscores the urgent need for continued research into the replicative mechanisms of this virus to combat its continued impact on global health, through the re-emergence of immuno-evasive variants.

Connected Component Analysis of Dynamical Perturbation Contact Networks

Describing protein dynamical networks through amino acid contacts is a powerful way to analyze complex biomolecular systems. However, due to the size of the systems, identifying the relevant features of protein-weighted graphs can be a difficult task. To address this issue, we present the connected component analysis (CCA) approach that allows for fast, robust, and unbiased analysis of dynamical perturbation contact networks (DPCNs). We first illustrate the CCA method as applied to a prototypical allosteric enzyme, the imidazoleglycerol phosphate synthase (IGPS) enzyme from Thermotoga maritima bacteria. This approach was shown to outperform the clustering methods applied to DPCNs, which could not capture the propagation of the allosteric signal within the protein graph. On the other hand, CCA reduced the DPCN size, providing connected components that nicely describe the allosteric propagation of the signal from the effector to the active sites of the protein. By applying the CCA to the IGPS enzyme in different conditions, i.e., at high temperature and from another organism (yeast IGPS), and to a different enzyme, i.e., a protein kinase, we demonstrated how CCA of DPCNs is an effective and transferable tool that facilitates the analysis of protein-weighted networks.

Examining the origins of observed terahertz modes from an optically pumped atomistic model protein in aqueous solution

The microscopic origins of terahertz (THz) vibrational modes in biological systems are an active and open area of current research. Recent experiments [Phys Rev X. 8, 031061 (2018)] have revealed the presence of a pronounced mode at ∼0.3 THz in fluorophore-decorated bovine serum albumin (BSA) protein in aqueous solution under nonequilibrium conditions induced by optical pumping. This result was heuristically interpreted as a collective elastic fluctuation originating from the activation of a low-frequency phonon mode. In this work, we show that the sub-THz spectroscopic response emerges in a statistically significant manner (> 2 σ ) from such collective behavior, illustrating how photoexcitation can alter specific THz vibrational modes. We revisit the theoretical analysis with proof-of-concept molecular dynamics that introduce optical excitations into the simulations. Using information theory techniques, we show that these excitations can give rise to a multiscale response involving two optically excited chromophores (tryptophans), other amino acids in the protein, ions, and water. Our results motivate new experiments and fully nonequilibrium simulations to probe these phenomena, as well as the refinement of atomistic models of Fröhlich condensates that are fundamentally determined by nonlinear interactions in biology.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study we employ this paradigm to answer a basic question: in enzyme superfamilies where the catalytic mechanism, active sites and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as the representatives of the conserved Protein Tyrosine Phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

Gene co-expression network construction and analysis for identification of genetic biomarkers associated with glioblastoma multiforme using topological findings

BACKGROUND

Glioblastoma multiforme (GBM) is one of the most malignant types of central nervous system tumors. GBM patients usually have a poor prognosis. Identification of genes associated with the progression of the disease is essential to explain the mechanisms or improve the prognosis of GBM by catering to targeted therapy. It is crucial to develop a methodology for constructing a biological network and analyze it to identify potential biomarkers associated with disease progression.

METHODS

Gene expression datasets are obtained from TCGA data repository to carry out this study. A survival analysis is performed to identify survival associated genes of GBM patient. A gene co-expression network is constructed based on Pearson correlation between the gene's expressions. Various topological measures along with set operations from graph theory are applied to identify most influential genes linked with the progression of the GBM.

RESULTS

Ten key genes are identified as a potential biomarkers associated with GBM based on centrality measures applied to the disease network. These genes are SEMA3B, APS, SLC44A2, MARK2, PITPNM2, SFRP1, PRLH, DIP2C, CTSZ, and KRTAP4.2. Higher expression values of two genes, SLC44A2 and KRTAP4.2 are found to be associated with progression and lower expression values of seven gens SEMA3B, APS, MARK2, PITPNM2, SFRP1, PRLH, DIP2C, and CTSZ are linked with the progression of the GBM.

CONCLUSIONS

The proposed methodology employing a network topological approach to identify genetic biomarkers associated with cancer.

Mapping N- to C-terminal allosteric coupling through disruption of a putative CD74 activation site in D-dopachrome tautomerase

The macrophage migration inhibitory factor (MIF) protein family consists of MIF and D-dopachrome tautomerase (also known as MIF-2). These homologs share 34% sequence identity while maintaining nearly indistinguishable tertiary and quaternary structure, which is likely a major contributor to their overlapping functions, including the binding and activation of the cluster of differentiation 74 (CD74) receptor to mediate inflammation. Previously, we investigated a novel allosteric site, Tyr99, that modulated N-terminal catalytic activity in MIF through a "pathway" of dynamically coupled residues. In a comparative study, we revealed an analogous allosteric pathway in MIF-2 despite its unique primary sequence. Disruptions of the MIF and MIF-2 N termini also diminished CD74 activation at the C terminus, though the receptor activation site is not fully defined in MIF-2. In this study, we use site-directed mutagenesis, NMR spectroscopy, molecular simulations, in vitro and in vivo biochemistry to explore the putative CD74 activation region of MIF-2 based on homology to MIF. We also confirm its reciprocal structural coupling to the MIF-2 allosteric site and N-terminal enzymatic site. Thus, we provide further insight into the CD74 activation site of MIF-2 and its allosteric coupling for immunoregulation.

Visualizing the Residue Interaction Landscape of Proteins by Temporal Network Embedding

Characterizing the structural dynamics of proteins with heterogeneous conformational landscapes is crucial to understanding complex biomolecular processes. To this end, dimensionality reduction algorithms are used to produce low-dimensional embeddings of the high-dimensional conformational phase space. However, identifying a compact and informative set of input features for the embedding remains an ongoing challenge. Here, we propose to harness the power of Residue Interaction Networks (RINs) and their centrality measures, established tools to provide a graph theoretical view on molecular structure. Specifically, we combine the closeness centrality, which captures global features of the protein conformation at residue-wise resolution, with EncoderMap, a hybrid neural-network autoencoder/multidimensional-scaling like dimensionality reduction algorithm. We find that the resulting low-dimensional embedding is a meaningful visualization of the residue interaction landscape that resolves structural details of the protein behavior while retaining global interpretability. This feature-based graph embedding of temporal protein graphs makes it possible to apply the general descriptive power of RIN formalisms to the analysis of protein simulations of complex processes such as protein folding and multidomain interactions requiring no protein-specific input. We demonstrate this on simulations of the fast folding protein Trp-Cage and the multidomain signaling protein FAT10. Due to its generality and modularity, the presented approach can easily be transferred to other protein systems.

Turning up the heat mimics allosteric signaling in imidazole-glycerol phosphate synthase

Allosteric drugs have the potential to revolutionize biomedicine due to their enhanced selectivity and protection against overdosage. However, we need to better understand allosteric mechanisms in order to fully harness their potential in drug discovery. In this study, molecular dynamics simulations and nuclear magnetic resonance spectroscopy are used to investigate how increases in temperature affect allostery in imidazole glycerol phosphate synthase. Results demonstrate that temperature increase triggers a cascade of local amino acid-to-amino acid dynamics that remarkably resembles the allosteric activation that takes place upon effector binding. The differences in the allosteric response elicited by temperature increase as opposed to effector binding are conditional to the alterations of collective motions induced by either mode of activation. This work provides an atomistic picture of temperature-dependent allostery, which could be harnessed to more precisely control enzyme function.

Identification of region-specific amino acid signatures for doxorubicin-induced chemo brain

Doxorubicin (DOX) is a cornerstone of chemotherapy for solid tumors and leukemias. DOX-induced cognitive impairment, termed chemo brain, has been reported in cancer survivors, whereas its mechanism remains poorly understood. Here we initially evaluated the cognitive impairments of mice treated with clinically relevant, long-term, low-dosage of DOX. Using HILIC-MS/MS-based targeted metabolomics, we presented the changes of 21 amino acids across six anatomical brain regions of mice with DOX-induced chemo brain. By mapping the altered amino acids to the human metabolic network, we constructed an amino acid-based network module for each brain region. We identified phenylalanine, tyrosine, methionine, and γ-aminobutyric acid as putative signatures of three regions (hippocampus, prefrontal cortex, and neocortex) highly associated with cognition. Relying on the reported mouse brain metabolome atlas, we found that DOX might perturb the amino acid homeostasis in multiple brain regions, similar to the changes in the aging brain. Correlation analysis suggested the possible indirect neurotoxicity of DOX that altered the brain levels of phenylalanine, tyrosine, and methionine by causing metabolic disorders in the liver and kidney. In summary, we revealed the region-specific amino acid signatures as actionable targets for DOX-induced chemo brain, which might provide safer treatment and improve the quality of life among cancer survivors.

Heterogeneity in winged helix-turn-helix and substrate DNA interactions: Insights from theory and experiments

Specific interactions between transcription factors (TFs) and substrate DNA constitute the fundamental basis of gene expression. Unlike in TFs like basic helix-loop-helix or basic leucine zippers, prediction of substrate DNA is extremely challenging for helix-turn-helix (HTH). Experimental techniques like chromatin immunoprecipitation combined with massively parallel DNA sequencing remains a viable option. We characterize the molecular basis of heterogeneity in HTH-DNA interaction using in silico tools and thence validate them experimentally. Given the profound functional diversity in HTH, we focus primarily on winged-HTH (wHTH). We consider 180 wHTH TFs, whose experimental three-dimensional structures are available in DNA bound/unbound conformations. Starting with PDB-wide scanning and curation of data, we construct a phylogenetic tree, which distributes 180 wHTH sequences under multiple sub-groups. Structure-sequence alignment followed by detailed intra/intergroup analysis, covariation studies and extensive network theory analysis help us to gain deep insight into heterogeneous wHTH-substrate DNA interactions. A central aim of this study is to find a consensus to predict the substrate DNA sequence for wHTH, amidst heterogeneity. The strength of our exhaustive theoretical investigations including molecular docking are successfully tested through experimental characterization of wHTH TF from Sulfurimonas denitrificans.

The "violin model": Looking at community networks for dynamic allostery

Allosteric regulation of proteins continues to be an engaging research topic for the scientific community. Models describing allosteric communication have evolved from focusing on conformation-based descriptors of protein structural changes to appreciating the role of internal protein dynamics as a mediator of allostery. Here, we explain a "violin model" for allostery as a contemporary method for approaching the Cooper-Dryden model based on redistribution of protein thermal fluctuations. Based on graph theory, the violin model makes use of community network analysis to functionally cluster correlated protein motions obtained from molecular dynamics simulations. This Review provides the theory and workflow of the methodology and explains the application of violin model to unravel the workings of protein kinase A.

Computational analysis of the sequence-structure relation in SARS-CoV-2 spike protein using protein contact networks

The structure of proteins impacts directly on the function they perform. Mutations in the primary sequence can provoke structural changes with consequent modification of functional properties. SARS-CoV-2 proteins have been extensively studied during the pandemic. This wide dataset, related to sequence and structure, has enabled joint sequence-structure analysis. In this work, we focus on the SARS-CoV-2 S (Spike) protein and the relations between sequence mutations and structure variations, in order to shed light on the structural changes stemming from the position of mutated amino acid residues in three different SARS-CoV-2 strains. We propose the use of protein contact network (PCN) formalism to: (i) obtain a global metric space and compare various molecular entities, (ii) give a structural explanation of the observed phenotype, and (iii) provide context dependent descriptors of single mutations. PCNs have been used to compare sequence and structure of the Alpha, Delta, and Omicron SARS-CoV-2 variants, and we found that omicron has a unique mutational pattern leading to different structural consequences from mutations of other strains. The non-random distribution of changes in network centrality along the chain has allowed to shed light on the structural (and functional) consequences of mutations.

Calling for justice with #JusticeforBreonnaTaylor: a case study of hashtag activism in the evolution of the black lives matter movement

Taking a stage-based approach, before and after the release of the 15-h audio recording files of the grand jury's inquiry on the Breonna Taylor case on October 2, 2020, this study examined the #JusticeforBreonnaTaylor Twitter networks. By employing multimethodology, including natural language processing, social network analysis, and qualitative textual analysis, I examined keys connectors of the two Twitter networks and investigated major themes conducting thematic analysis of network discourses and highly associated hashtags with the hashtag #JusticeforBreonnaTaylor. In both networks, several key stakeholders, such as Benjamin Crump, Danial Cameron, and Black women activists were identified as key connectors along with social activists and ordinary participants. Demanding justice to the case was the core agenda of the hashtag activism. The findings of the study revealed that the participants not only shared breaking news and important information but also organized protests and routinely tagged people to spread messages about the Taylor's case on Twitter. The participants conversed major issues about the Taylor case and set the agendas for the next action, such as encouraging to take part in voting for the 2020 presidential election. The thematic analysis concurrently demonstrated that the network participants strongly demanded legal prosecution to the three Louisville cops that involved in the act of killing Breonna Taylor during the botched raid in her apartment.

Time Series Analysis of SARS-CoV-2 Genomes and Correlations among Highly Prevalent Mutations

The efforts of the scientific community to tame the recent pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) seem to have been diluted by the emergence of new viral strains. Therefore, it is imperative to understand the effect of mutations on viral evolution. We performed a time series analysis on 59,541 SARS-CoV-2 genomic sequences from around the world to gain insights into the kinetics of the mutations arising in the viral genomes. These 59,541 genomes were grouped according to month (January 2020 to March 2021) based on the collection date. Meta-analysis of these data led us to identify significant mutations in viral genomes. Pearson correlation of these mutations led us to the identification of 16 comutations. Among these comutations, some of the individual mutations have been shown to contribute to viral replication and fitness, suggesting a possible role of other unexplored mutations in viral evolution. We observed that the mutations 241C>T in the 5' untranslated region (UTR), 3037C>T in nsp3, 14408C>T in the RNA-dependent RNA polymerase (RdRp), and 23403A>G in spike are correlated with each other and were grouped in a single cluster by hierarchical clustering. These mutations have replaced the wild-type nucleotides in SARS-CoV-2 sequences. Additionally, we employed a suite of computational tools to investigate the effects of T85I (1059C>T), P323L (14408C>T), and Q57H (25563G>T) mutations in nsp2, RdRp, and the ORF3a protein of SARS-CoV-2, respectively. We observed that the mutations T85I and Q57H tend to be deleterious and destabilize the respective wild-type protein, whereas P323L in RdRp tends to be neutral and has a stabilizing effect. IMPORTANCE We performed a meta-analysis on SARS-CoV-2 genomes categorized by collection month and identified several significant mutations. Pearson correlation analysis of these significant mutations identified 16 comutations having absolute correlation coefficients of >0.4 and a frequency of >30% in the genomes used in this study. The correlation results were further validated by another statistical tool called hierarchical clustering, where mutations were grouped in clusters on the basis of their similarity. We identified several positive and negative correlations among comutations in SARS-CoV-2 isolates from around the world which might contribute to viral pathogenesis. The negative correlations among some of the mutations in SARS-CoV-2 identified in this study warrant further investigations. Further analysis of mutations such as T85I in nsp2 and Q57H in ORF3a protein revealed that these mutations tend to destabilize the protein relative to the wild type, whereas P323L in RdRp is neutral and has a stabilizing effect. Thus, we have identified several comutations which can be further characterized to gain insights into SARS-CoV-2 evolution.

Most influential countries in the international medical device trade: Network-based analysis

Since the outbreak of the coronavirus disease 2019 (COVID-19) pandemic, the international medical device trade has received extensive attention. To maintain the domestic supply of medical devices, some countries have sought multilateral trade cooperation or simply implemented export restrictions, which has exacerbated the instability and fragility of the global medical device market. It is crucial for government policymakers to identify the most influential countries in the international medical device trade and nip exports in the bud. However, few efforts have been made in previous studies to explore various countries' influence on the international medical device trade in light of their intricate trade relationships. To fill these research gaps, this study constructs a global medical device trade network (GMDTN) and explores the criticality of various countries from a network-based perspective. The evolution patterns and geographical distribution of influence among countries in the GMDTN are revealed. Details on the ways in which the influence of some crucial countries has formed are provided. The results show that the global medical device trade market is export oriented. The formation of some countries' strong influence may be due to their large number of trading partners or the deep dependence of some of those trading partners on that country (namely, breadth- or depth-based patterns). It is worth noting that the US has a dominant position in the international medical device trade in terms of both breadth and depth. In addition, some countries play a critical role as intermediate points in the influence formation process of other countries, although these countries are not critical direct trading partners. The findings of this study provide implications for policymakers seeking to understand the influence of countries on the international medical device trade and to proactively prepare responses to unexpected changes in this trade.

Smoking by altering the peri-implant microbial community structure compromises the responsiveness to treatment

Smoking is an essential risk factor for peri-implant diseases. It also hampers the clinical outcomes of peri-implant therapies. Nonetheless, the effect of smoking can go undetected until the emergence of clinical signs. Bacterial-induced inflammation is responsible for the initiation and progression of peri-implant diseases. We hypothesize that smoking impacts the peri-implant microbiome even in status of clinical health, putting it into a sub-healthy condition that responds poorly to peri-implant treatments. To validate this, peri-implant plaque samples from 18 participants including 10 smokers (S) and 8 non-smokers (NS), who had received implant prostheses were analyzed using metagenomic shotgun sequencing. The results showed that in addition to taxonomical and functional differences, the local stability in the S group was also shown to be much higher than that in the NS group, indicating greater stubbornness of the peri-implant microbiome associated with smoking. Besides, the topological structures were also distinct between the two groups. The highly connected species interacted more preferentially with each other in the S group (eigenvector centralization, 0.0273 in S and 0.0183 in NS), resulting in a greater tendency of forming small-world modules (modularity, 0.714 in S and 0.582 in NS). While in the NS group, inter-species correlations were more evenly distributed (clustering coefficient, 0.532 in S and 0.666 in NS). These alterations overall explained the greater stubbornness of the peri-implant microbiome associated with smoking, which may cause poor responsiveness to peri-implant therapies. From a microbial perspective, this may be a potential reason why smoking impacts negatively on the outcome of peri-implant treatments.

Excavating the social representations and perceived barriers of organ donation in China over the past decade: A hybrid text analysis approach

Background

Organ donation has been claimed as a prosocial behavior to prolong the recipient's life and deliver great love. However, the supply-demand ratio of organs in China is highly unbalanced. Being entangled with multiple factors derived from individual and supra-individual levels, organ donation in China is important but sensitive. Previous scholars usually depended on obtrusive approaches to explore the facilitators and hindrances of organ donation, which is hard to discover genuine perceptions toward organ donation. Besides, relatively limited scholarly attention has been paid to what hampers organ donation in China.

Objective

We intended to excavate the diversified social representations and perceived barriers to organ donation in China over the past decade.

Method

Two kinds of text analysis methods-semantic network analysis and conventional content analysis, were applied to 120,172 posts from ordinary users on the Sina Weibo platform to address the research questions.

Results

Regarding social representations, the "hope, understanding, and acceptance" of organ donation was the most pronounced one (34% of the whole semantic network), followed by "family story" (26%), "the procedure of organ donation in NGOs" (15%), "the practical value of organ donation" (14%), and "organ donation in the medical context" (11%). Regarding perceived barriers, a four-layer framework was constructed, including (1) the individual level, mainly about the fear of death and postmortem autopsy; (2) the familial level, which refers to the opposition from family members; (3) the societal level, which alludes to distrust toward medical institutions and the general society; (4) the cultural level, which covers religious-cultural concerns about fatalism.

Conclusion

In concordance with prior works on social representations regarding organ donation, the current study also uncovered the coexistence of antithetical representations about organ donation-the longing for survival and the fear of death. This representation pair serves as the foundation of Chinese people's ambivalence. Besides, family-related narratives were dispersed over various representations, demonstrating the critical position of family support in organ donation. Moreover, the four-layer framework concerning donation barriers affords a reference for future empirical studies. The practical implications of this work are further discussed.

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct "enhanced network models", describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. "Ab-initio network models" combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.

MptpA Kinetics Enhanced by Allosteric Control of an Active Conformation

Understanding allostery in the Mycobacterium tuberculosis low molecular weight protein tyrosine phosphatase (MptpA) is a subject of great interest since MptpA is one of two protein tyrosine phosphatases (PTPs) from the pathogenic organism Mycobacterium tuberculosis expressed during host cell infection. Here, we combine computational modeling with solution NMR spectroscopy and we find that Q75 is an allosteric site. Removal of the polar side chain of Q75 by mutation to leucine results in a cascade of events that reposition the acid loop over the active site and relocates the catalytic aspartic acid (D126) at an optimal position for proton donation to the leaving aryl group of the substrate and for subsequent hydrolysis of the thiophosphoryl intermediate. The computational analysis is consistent with kinetic data, and NMR spectroscopy, showing that the Q75L mutant exhibits enhanced reaction kinetics with similar substrate binding affinity. We anticipate that our findings will motivate further studies on the possibility that MptpA remains passivated during the chronic state of infection and increases its activity as part of the pathogenic life cycle of M. tuberculosis possibly via allosteric means.

Network-based integration of multi-omics data for clinical outcome prediction in neuroblastoma

Multi-omics data are increasingly being gathered for investigations of complex diseases such as cancer. However, high dimensionality, small sample size, and heterogeneity of different omics types pose huge challenges to integrated analysis. In this paper, we evaluate two network-based approaches for integration of multi-omics data in an application of clinical outcome prediction of neuroblastoma. We derive Patient Similarity Networks (PSN) as the first step for individual omics data by computing distances among patients from omics features. The fusion of different omics can be investigated in two ways: the network-level fusion is achieved using Similarity Network Fusion algorithm for fusing the PSNs derived for individual omics types; and the feature-level fusion is achieved by fusing the network features obtained from individual PSNs. We demonstrate our methods on two high-risk neuroblastoma datasets from SEQC project and TARGET project. We propose Deep Neural Network and Machine Learning methods with Recursive Feature Elimination as the predictor of survival status of neuroblastoma patients. Our results indicate that network-level fusion outperformed feature-level fusion for integration of different omics data whereas feature-level fusion is more suitable incorporating different feature types derived from same omics type. We conclude that the network-based methods are capable of handling heterogeneity and high dimensionality well in the integration of multi-omics.

Identification of key genes with prognostic value in gastric cancer by bioinformatics analysis

Background: Gastric cancer (GC) is a digestive system tumor with high morbidity and mortality. It is urgently required to identify genes to elucidate the underlying molecular mechanisms. The aim of this study is to identify the key genes which may affect the prognosis of GC patients and be a therapeutic strategy for GC patients by bioinformatic analysis.

Methods: The significant prognostic differentially expressed genes (DEGs) were screened out from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) datasets. The protein–protein interaction (PPI) network was established by STRING and screening key genes by MCODE and CytoNCA plug-ins in Cytoscape. Functional enrichment analysis, construction of a prognostic risk model, and nomograms verify key genes as potential therapeutic targets.

Results: In total, 997 genes and 805 genes were related to the prognosis of GC in the GSE84437 and TCGA datasets, respectively. We define the 128 genes shared by the two datasets as prognostic DEGs (P-DEGs). Then, the first four genes (MYLK, MYL9, LUM, and CAV1) with great node importance in the PPI network of P-DEGs were identified as key genes. Independent prognostic risk analysis found that patients with high key gene expression had a poor prognosis, excluding their age, gender, and TNM stage. GO and KEGG enrichment analyses showed that key genes may exert influence through the PI3K-Akt pathway, in which extracellular matrix organization and focal adhesion may play important roles in key genes influencing the prognosis of GC patients.

Conclusion: We found that MYLK, MYL9, LUM, and CAV1 are potential and reliable prognostic key genes that affect the invasion and migration of gastric cancer.

A QUBO formulation for top-τ eigencentrality nodes

The efficient calculation of the centrality or "hierarchy" of nodes in a network has gained great relevance in recent years due to the generation of large amounts of data. The eigenvector centrality (aka eigencentrality) is quickly becoming a good metric for centrality due to both its simplicity and fidelity. In this work we lay the foundations for solving the eigencentrality problem of ranking the importance of the nodes of a network with scores from the eigenvector of the network, using quantum computational paradigms such as quantum annealing and gate-based quantum computing. The problem is reformulated as a quadratic unconstrained binary optimization (QUBO) that can be solved on both quantum architectures. The results focus on correctly identifying a given number of the most important nodes in numerous networks given by the sparse vector solution of our QUBO formulation of the problem of identifying the top-τ highest eigencentrality nodes in a network on both the D-Wave and IBM quantum computers.