Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS

Impacts of Nutlin-3a and exercise on murine double minute 2-enriched glioma treatment

JOURNAL/nrgr/04.03/01300535-202504000-00029/figure1/v/2024-07-06T104127Z/r/image-tiff Recent research has demonstrated the impact of physical activity on the prognosis of glioma patients, with evidence suggesting exercise may reduce mortality risks and aid neural regeneration. The role of the small ubiquitin-like modifier (SUMO) protein, especially post-exercise, in cancer progression, is gaining attention, as are the potential anti-cancer effects of SUMOylation. We used machine learning to create the exercise and SUMO-related gene signature (ESLRS). This signature shows how physical activity might help improve the outlook for low-grade glioma and other cancers. We demonstrated the prognostic and immunotherapeutic significance of ESLRS markers, specifically highlighting how murine double minute 2 (MDM2), a component of the ESLRS, can be targeted by nutlin-3. This underscores the intricate relationship between natural compounds such as nutlin-3 and immune regulation. Using comprehensive CRISPR screening, we validated the effects of specific ESLRS genes on low-grade glioma progression. We also revealed insights into the effectiveness of Nutlin-3a as a potent MDM2 inhibitor through molecular docking and dynamic simulation. Nutlin-3a inhibited glioma cell proliferation and activated the p53 pathway. Its efficacy decreased with MDM2 overexpression, and this was reversed by Nutlin-3a or exercise. Experiments using a low-grade glioma mouse model highlighted the effect of physical activity on oxidative stress and molecular pathway regulation. Notably, both physical exercise and Nutlin-3a administration improved physical function in mice bearing tumors derived from MDM2-overexpressing cells. These results suggest the potential for Nutlin-3a, an MDM2 inhibitor, with physical exercise as a therapeutic approach for glioma management. Our research also supports the use of natural products for therapy and sheds light on the interaction of exercise, natural products, and immune regulation in cancer treatment.

Discovery of novel and potent dual-targeting AXL/HDAC2 inhibitors for colorectal cancer treatment via structure-based pharmacophore modelling, virtual screening, and molecular docking, molecular dynamics simulation studies, and biological evaluation

Colorectal cancer (CRC) is one of the most common cancers worldwide. Nowadays, owing to the complex mechanism of tumorigenesis, simultaneous inhibition of multiple targets is an important anticancer strategy. Recent studies have demonstrated receptor tyrosine kinase AXL (AXL) and histone deacetylase 2 (HDAC2) are closely associated with colorectal cancer. Herein, we identified five hit compounds concurrently targeting AXL and HDAC2 using virtual screening. Inhibitory experiments revealed these hit compounds potently inhibited AXL and HDAC2 in the nanomolar range. Among them, Hit-3 showed the strongest inhibitory effects which were better than that of the positive control groups. Additionally, MD assays showed that Hit-3 could bind stably to the AXL and HDAC2 active pockets. Further MTT assays demonstrated that Hit-3 showed potent anti-proliferative activity. Most importantly, Hit-3 exhibited significant in vivo antitumor efficacy in xenograft models. Collectively, this study is the first discovery of dual-targeting AXL/HDAC2 inhibitors for colorectal cancer treatment.

Binding affinity between coronavirus spike protein and human ACE2 receptor

Coronaviruses (CoVs) pose a major risk to global public health due to their ability to infect diverse animal species and potential for emergence in humans. The CoV spike protein mediates viral entry into the cell and plays a crucial role in determining the binding affinity to host cell receptors. With particular emphasis on α- and β-coronaviruses that infect humans and domestic animals, current research on CoV receptor use suggests that the exploitation of the angiotensin-converting enzyme 2 (ACE2) receptor poses a significant threat for viral emergence with pandemic potential. This review summarizes the approaches used to study binding interactions between CoV spike proteins and the human ACE2 (hACE2) receptor. Solid-phase enzyme immunoassays and cell binding assays allow qualitative assessment of binding but lack quantitative evaluation of affinity. Surface plasmon resonance, Bio-layer interferometry, and Microscale Thermophoresis on the other hand, provide accurate affinity measurement through equilibrium dissociation constants (KD). In silico modeling predicts affinity through binding structure modeling, protein-protein docking simulations, and binding energy calculations but reveals inconsistent results due to the lack of a standardized approach. Machine learning and deep learning models utilize simulated and experimental protein-protein interaction data to elucidate the critical residues associated with CoV binding affinity to hACE2. Further optimization and standardization of existing approaches for studying binding affinity could aid pandemic preparedness. Specifically, prioritizing surveillance of CoVs that can bind to human receptors stands to mitigate the risk of zoonotic spillover.

The protective effect of Lonicera japonica Thunb. against lipopolysaccharide-induced acute lung injury in mice: Modulation of inflammation, oxidative stress, and ferroptosis

ETHNOPHARMACOLOGICAL RELEVANCE

Various components of Lonicera japonica Thunb. (LJT) exhibit pharmacological activities, including anti-inflammatory and antioxidant effects. Nevertheless, the relationship between LJT and ferroptosis remains largely unexplored.

AIM OF THE STUDY

The purpose of this research was to look into the role of LJT in regulating LPS-induced ferroptosis in ALI and to compare the effects of different parts of LJT.

MATERIALS AND METHODS

We established a mice ALI model by treating with LPS. Administered mice with different doses of Lonicerae Japonicae Flos (LJF), Lonicera Japonica Leaves (LJL) and Lonicerae Caulis (LRC) extracts, respectively. The levels of IL-6, IL-1β, TNF-α, IL-4, IL-10, and PGE2 in bronchoalveolar lavage fluid (BALF) were measured using enzyme-linked immunosorbent assay. Furthermore, the concentrations of superoxide dismutase (SOD), malondialdehyde (MDA), glutathione (GSH), reactive oxygen species (ROS), and total ferrous ions (Fe2+) in lung tissues were evaluated. Hematoxylin and eosin staining was conducted to examine the morphological structure of lung tissues. Transmission electron microscopy was used to investigate the ultrastructural morphology of mitochondria. Furthermore, the effects of LJT were evaluated via immunohistochemical staining, western blotting, and quantitative real-time polymerase chain reaction analyses. Finally, employing molecular docking and molecular dynamics research techniques, we aimed to identify crucial components in LJT that might inhibit ferroptosis by targeting nuclear factor erythroid 2-related factor 2 (Nrf2) and glutathione peroxidase 4 (GPX4).

RESULTS

We observed that pretreatment with LJT significantly mitigated LPS-induced lung injury and suppressed ferroptosis. This was supported by reduced accumulation of pro-inflammatory cytokines, ROS, MDA, and Fe2+, along with increased levels of anti-inflammatory cytokines, SOD, GSH, Nrf2, and GPX4 in the lung tissues of ALI mice. Luteolin-7-O-rutinoside, apigenin-7-O-rutinoside, and amentoflavone in LJT exhibit excellent docking effects with key targets of ferroptosis, Nrf2 and GPX4.

CONCLUSIONS

Pretreatment with LJT may alleviate LPS-induced ALI, possibly by suppressing ferroptosis. Our initial results indicate that LJT activates the Nrf2/GPX4 axis, providing protection against ferroptosis in ALI. This finding offers a promising therapeutic candidate for ALI treatment.

Structure and dynamics of cholesterol-mediated aquaporin-0 arrays and implications for lipid rafts

Aquaporin-0 (AQP0) tetramers form square arrays in lens membranes through a yet unknown mechanism, but lens membranes are enriched in sphingomyelin and cholesterol. Here, we determined electron crystallographic structures of AQP0 in sphingomyelin/cholesterol membranes and performed molecular dynamics (MD) simulations to establish that the observed cholesterol positions represent those seen around an isolated AQP0 tetramer and that the AQP0 tetramer largely defines the location and orientation of most of its associated cholesterol molecules. At a high concentration, cholesterol increases the hydrophobic thickness of the annular lipid shell around AQP0 tetramers, which may thus cluster to mitigate the resulting hydrophobic mismatch. Moreover, neighboring AQP0 tetramers sandwich a cholesterol deep in the center of the membrane. MD simulations show that the association of two AQP0 tetramers is necessary to maintain the deep cholesterol in its position and that the deep cholesterol increases the force required to laterally detach two AQP0 tetramers, not only due to protein-protein contacts but also due to increased lipid-protein complementarity. Since each tetramer interacts with four such 'glue' cholesterols, avidity effects may stabilize larger arrays. The principles proposed to drive AQP0 array formation could also underlie protein clustering in lipid rafts.

Revealing the material basis and mechanism for the inhibition of intestinal peristalsis by Zingiber officinale Roscoe through integrated metabolomics, serum pharmacochemistry, and network pharmacology

Abnormal relaxation and contraction of intestinal smooth muscle can cause various intestinal diseases. Diarrhea is a common and important public health problem worldwide in epidemiology. Zingiber officinale Roscoe (fresh ginger) has been found to treat diarrhea, but the material basis and mechanism of action that inhibits intestinal peristalsis remain unclear. Metabolomics and serum pharmacology were used to identify differential metabolites, metabolic pathways, and pharmacodynamic substances, and were then combined with network pharmacology to explore the potential targets of ginger that inhibit intestinal peristalsis during diarrhea treatment, and the targets identified were verified using molecular docking and molecular dynamic simulation. We found that 25 active components of ginger (the six most relevant components), 35 potential key targets (three core targets), 40 differential metabolites (four key metabolites), and four major metabolic pathways were involved in the process by which ginger inhibits intestinal peristalsis during diarrhea treatment. This study reveals the complex mechanism of action and pharmacodynamic material basis of ginger in the inhibition of intestinal peristalsis, and this information helps in the development of new Chinese medicine to treat diarrhea and lays the foundation for the clinical application of ginger.

Tumor-agnostic cancer therapy using antibodies targeting oncofetal chondroitin sulfate

Molecular similarities between embryonic and malignant cells can be exploited to target tumors through specific signatures absent in healthy adult tissues. One such embryonic signature tumors express is oncofetal chondroitin sulfate (ofCS), which supports disease progression and dissemination in cancer. Here, we report the identification and characterization of phage display-derived antibody fragments recognizing two distinct ofCS epitopes. These antibody fragments show binding affinity to ofCS in the low nanomolar range across a broad selection of solid tumor types in vitro and in vivo with minimal binding to normal, inflamed, or benign tumor tissues. Anti-ofCS antibody drug conjugates and bispecific immune cell engagers based on these targeting moieties disrupt tumor progression in animal models of human and murine cancers. Thus, anti-ofCS antibody fragments hold promise for the development of broadly effective therapeutic and diagnostic applications targeting human malignancies.

Thermodynamic Details of Pinholin S2168 Activation Revealed Using Alchemical Free Energy Simulations

Pinholin S2168 is a viral integral membrane protein whose function is to form nanoscopic "pinholes" in bacterial cell membranes to induce cell lysis as part of the viral replication cycle. Pinholin can transition from an inactive to an active conformation by exposing a transmembrane domain (TMD1) to the extracellular fluid. Upon activation, several copies of the protein assemble via interactions among a second transmembrane domain (TMD2) to form a single pore, thus hastening cell lysis and viral escape. The following experiments provide conformational descriptors of pinholin in active and inactive states and elucidate the molecular driving forces that control pinholin activity. In the present study, molecular dynamics (MD) simulations have been used to refine experimentally derived conformational descriptors into an atomistically detailed model of irsS2168, an antiholin mutant. To provide additional details about the thermodynamics of pinholin activation and to overcome large intrinsic kinetic barriers to activation, alchemical free energy simulations have been conducted. Alchemical mutations reveal the change in folding free energy upon mutation. The results suggest that alchemical mutations are an effective tool to rationalize experimental observations and predict the effects of site mutations on conformational states for proteins integrated into lipid bilayers. S16F, A17Q, A17Q+G21Q, and A17Q+G21Q+G14Q mutants reveal how changes in hydrophilicity and disruption of the glycine zipper motif influence pinholin's thermodynamic equilibrium, favoring the active conformation. These findings align with experimental observations from DEER spectroscopy, demonstrating that mutations increasing the hydrophilicity of TMD1 promote activation by making TMD1 more likely to exit the membrane and enter the extracellular fluid.

Structure-Based Computational Scanning of Chemical Modification Sites in Biologics

To address the challenges of short half-life, immunogenicity, and nonspecific distribution, chemical modifications of peptide and protein-based drugs have emerged as a versatile strategy for improving their therapeutic efficacy. One such modification involves the derivatization of peptides and proteins with fatty acids, which can protract their half-life, modify their biodistribution, and potentially enable targeted delivery to specific tissues or disease sites of interest. However, the present strategies for the synthesis of such synthetically modified biologics require numerous rounds of experimental testing and often yield unstable, inactive, or heterogeneous products. To address the inefficiencies in designing modified biologics, we developed a hybrid computational workflow that integrates RosettaMatch from the Rosetta suite of protein modeling tools with molecular dynamics (MD) simulations. This approach not only reduces the number of amino acid positions that need to be experimentally tested by targeting only the most promising candidates for modification but also expedites the design of chemically modified biologics with the desired properties, ensuring a rapid and cost-effective development cycle. Although we demonstrate the utility of our method on a peptide therapeutic, GLP-1, with different fatty acid derivatizations, this straightforward approach has the potential to streamline the design process of a diverse range of chemically modified therapeutics, enabling tailored enhancements to their pharmacokinetic properties.

Assessment of Four Theoretical Approaches to Predict Protein Flexibility in the Crystal Phase and Solution

In this paper, we evaluated the ability of four coarse-grained methods to predict protein flexible regions with potential biological importance, UNRES-flex, UNRES-DSSP-flex (based on the united residue model of polypeptide chains without and with secondary structure restraints, respectively), CABS-flex (based on the C-α, C-β, and side chain model), and nonlinear rigid block normal mode analysis (NOLB) with a set of 100 protein structures determined by NMR spectroscopy or X-ray crystallography, with all secondary structure types. End regions with high fluctuations were excluded from analysis. The Pearson and Spearman correlation coefficients were used to quantify the conformity between the calculated and experimental fluctuation profiles, the latter determined from NMR ensembles and X-ray B-factors, respectively. For X-ray structures (corresponding to proteins in a crowded environment), NOLB resulted in the best agreement between the predicted and experimental fluctuation profiles, while for NMR structures (corresponding to proteins in solution), the ranking of performance is CABS-flex > UNRES-DSSP-flex > UNRES-flex > NOLB; however, CABS-flex sometimes exaggerated the extent of small fluctuations, as opposed to UNRES-DSSP-flex.

Single-residue effects on the behavior of a nascent polypeptide chain inside the ribosome exit tunnel

Nascent polypeptide chains (NCs) are extruded from the ribosome through an exit tunnel (ET) traversing the large ribosomal subunit. The ET's irregular and chemically complex wall allows for various NC-ET interactions. Translational arrest peptides (APs) bind in the ET to induce translational arrest, a property that can be exploited to study NC-ET interactions by Force Profile Analysis (FPA). We employed FPA and molecular dynamics (MD) simulations to investigate how individual residues placed in a glycine-serine repeat segment within an AP-stalled NC interact with the ET to exert a pulling force on the AP and release stalling. Our results indicate that large and hydrophobic residues generate a pulling force on the NC when placed ≳10 residues away from the peptidyl transfer center (PTC). Moreover, an asparagine placed 12 residues from the PTC makes a specific stabilizing interaction with the tip of ribosomal protein uL22 that reduces the pulling force on the NC, while a lysine or leucine residue in the same position increases the pulling force. Finally, the MD simulations suggest how the Mannheimia succiniproducens SecM AP interacts with the ET to promote translational stalling.

Dynamics and thermochemistry of the negatively charged clusters in a 2-hydroxyethylhydrazinium nitrate ionic liquid system

The formation and fragmentation of negatively charged 2-hydroxyethylhydrazinium nitrate ([HOCH2CH2NH2NH2]+NO3-, HEHN) ionic liquid clusters were examined using a guided-ion beam tandem mass spectrometer furnished with collision-induced dissociation of selected ions with Xe atoms. Measurements included the compositions of cluster ions formed in the ionization source, and the dissociation products, cross sections, and 0 K threshold energies for individually selected cluster ions. To identify the structures of the main cluster ion series [(HEHN)n(HNO3)0-1NO3]- formed, molecular dynamics simulations were employed to create initial geometry guesses, followed by optimization at the ωB97XD/6-31+G(d,p) level of theory, from which global minimum structures were identified for reaction thermodynamics analyses. A comparison was made between the cluster formation and fragmentation in the negatively charged 2-hydroxyethylhydrazinium nitrate with those in the positive mode (reported by W. Zhou et al., Phys. Chem. Chem. Phys., 2023, 25, 17370). In both modes, the cluster ions were predominantly composed of m/z below 350; loss of a neutral 2-hydroxyethylhydrazinium nitrate ion pair represents the most important cluster fragmentation pathway, followed by intra-ion pair proton transfer-mediated 2-hydroxyethylhydrazine and HNO3 elimination; and all clusters started to dissociate at threshold energies less than 1.5 eV. The overwhelming similarities in the formation and fragmentation chemistry of positively vs. negatively charged 2-hydroxyethylhydrazinium nitrate clusters may be attributed to their inherent ionic nature and high electric conductivities.

Interfacial rheology of lanthanide binding peptide surfactants at the air-water interface

Peptide surfactants (PEPS) are studied to capture and retain rare earth elements (REEs) at air-water interfaces to enable REE separations. Peptide sequences, designed to selectively bind REEs, depend crucially on the position of ligands within their binding loop domain. These ligands form a coordination sphere that wraps and retains the cation. We study variants of lanthanide binding tags (LBTs) designed to complex strongly with Tb3+. The peptide LBT5- (with net charge -5) is known to bind Tb3+ and adsorb with more REE cations than peptide molecules, suggesting that undesired non-specific coulombic interactions occur. Rheological characterization of interfaces of LBT5- and Tb3+ solutions reveal the formation of an interfacial gel. To probe whether this gelation reflects chelation among intact adsorbed LBT5-:Tb3+ complexes or destruction of the binding loop, we study a variant, LBT3-, designed to form net neutral LBT3-:Tb3+ complexes. Solutions of LBT3- and Tb3+ form purely viscous layers in the presence of excess Tb3+, indicating that each peptide binds a single REE in an intact coordination sphere. We introduce the variant RR-LBT3- with net charge -3 and anionic ligands outside of the coordination sphere. We find that such exposed ligands promote interfacial gelation. Thus, a nuanced requirement for interfacial selectivity of PEPS is proposed: that anionic ligands outside of the coordination sphere must be avoided to prevent the non-selective recruitment of REE cations. This view is supported by simulation, including interfacial molecular dynamics simulations, and interfacial metadynamics simulations of the free energy landscape of the binding loop conformational space.

Influence of molecular parameters on the representativeness of interfacial properties of simple fluids

New parameterizations for the Lennard-Jones 12/6 potential capable of reproducing the vapor pressure and surface tension with sufficient precision, but not the liquid-vapor equilibrium densities for the case of simple fluids that include Ar, Kr, Xe, Ne, and CH4 are presented in this work. These results are compared with those derived from the family of Mie(n, 6) potentials, which adequately reproduce the coexistence curve and the vapor pressure, leaving aside the surface tension. In addition, a detailed analysis is presented on different parameterizations and methodologies, which have been developed in recent decades to estimate the interfacial properties of interest here for simple fluids, such as argon, which is a molecule that is, in principle, "simple" to study but that clearly reveals the enormous discrepancy between the results reported in the literature throughout these years. These facts undoubtedly reveal one of the fundamental problems in the context of molecular thermodynamics of fluids: reproducing different thermodynamic properties with sufficient precision from a single set of free parameters for some interaction potential. In order to show the scope of the parameterizations presented for the Lennard-Jones model, they were successfully applied to the case of binary mixtures, which included Ar-Kr, Ar-CH4, and Xe-Kr. Finally, and with the aim of showing a possible solution to the problem posed in this research, results of the same interfacial properties above mentioned for argon and methane were presented in this work by using a set of molecular interactions, called ANC2s, whose flexibility allowed to reproduce the experimental evidence with just one parameterization. The results reported in this work were generated using molecular dynamics simulations.

Under the Influence of Water: Molecular Recognition of Organic Hydrophilic Molecules in Water with a Prismarene Host Driven by Hydration Effects

A water-soluble prism[5]arene host can form endo-cavity complexes with hydrophilic organic substances in water by displacing frustrated water molecules from its deep cavity. Water molecules structured at both rims of the prismarene host can mediate hydrogen bonding interactions with the guest. Water-mediated hydrogen bonding interactions were invoked here to elucidate the elevated binding affinities and selectivity of the prismarene host toward hydrophilic organic guests. We show that water at the interface of a host-guest complex can act as an extension of the host structure, facilitating the accommodation of neutral guests within the binding site. This study highlights the crucial role of water in facilitating supramolecular interactions between a deep-cavity prismarene host and organic hydrophilic guests in aqueous medium.

Design of Cryptococcus neoformans multi-epitope vaccine based on immunoinformatics method

Cryptococcus neoformans is a widely distributed opportunistic pathogenic fungus. While C. neoformans commonly infects immunocompromised individuals, it can also affect those who are immunocompetent. Transmission of C. neoformans primarily occurs through the respiratory tract, leading to the development of meningitis. The mortality rate of Cryptococcal meningitis is high, and treatment options are limited. Cryptococcus neoformans infections pose a significant public health threat and currently lack targeted and effective response strategies. This study aimed to screen T lymphocyte (cytotoxic T lymphocyte and helper T lymphocyte) and B lymphocyte epitopes derived from four C. neoformans antigens and develop two multi-epitope vaccines by combining them with various adjuvants. Molecular docking results demonstrated that the vaccines bind stably to Toll-like receptor 4 ( and induce innate immunity. The credibility of the molecular docking results was validated through subsequent molecular dynamics simulations. Furthermore, the results of immune simulation analyses underscored the multi-epitope vaccine's capability to effectively induce robust humoral and cellular immune responses within the host organism. These two vaccines have demonstrated theoretical efficacy against C. neoformans infection as indicated by computer analysis. Nevertheless, additional experimental validation is essential to substantiate the protective efficacy of the vaccines.

Notch1 Phase Separation Coupled Percolation facilitates target gene expression and enhancer looping

The Notch receptor is a pleiotropic signaling protein that translates intercellular ligand interactions into changes in gene expression via the nuclear localization of the Notch intracellular Domain (NICD). Using a combination of immunohistochemistry, RNA in situ, Optogenetics and super-resolution live imaging of transcription in human cells, we show that the N1ICD can form condensates that positively facilitate Notch target gene expression. We determined that N1ICD undergoes Phase Separation Coupled Percolation (PSCP) into transcriptional condensates, which recruit, enrich, and encapsulate a broad set of core transcriptional proteins. We show that the capacity for condensation is due to the intrinsically disordered transcriptional activation domain of the N1ICD. In addition, the formation of such transcriptional condensates acts to promote Notch-mediated super enhancer-looping and concomitant activation of the MYC protooncogene expression. Overall, we introduce a novel mechanism of Notch1 activity in which discrete changes in nuclear N1ICD abundance are translated into the assembly of transcriptional condensates that facilitate gene expression by enriching essential transcriptional machineries at target genomic loci.

Application of the molecular dynamics simulation GROMACS in food science

Food comprises proteins, lipids, sugars and various other molecules that constitute a multicomponent biological system. It is challenging to investigate microscopic changes in food systems solely by performing conventional experiments. Molecular dynamics (MD) simulation serves as a crucial bridge in addressing this research gap. The Groningen Machine for Chemical Simulations (GROMACS) is an open-source, high-performing molecular dynamics simulation software that plays a significant role in food science research owing to its high flexibility and powerful functionality; it has been used to explore the molecular conformations and the mechanisms of interaction between food molecules at the microcosmic level and to analyze their properties and functions. This review presents the workflow of the GROMACS software and emphasizes the recent developments and achievements in its applications in food science research, thus providing important theoretical guidance and technical support for obtaining an in-depth understanding of the properties and functions of food.

Semi-rational engineering of D-allulose 3-epimerase for simultaneously improving the catalytic activity and thermostability based on D-allulose biosensor

BACKGROUND

D-Allulose is one of the most well-known rare sugars widely used in food, cosmetics, and pharmaceutical industries. The most popular method for D-allulose production is the conversion from D-fructose catalyzed by D-allulose 3-epimerase (DAEase). To address the general problem of low catalytic efficiency and poor thermostability of wild-type DAEase, D-allulose biosensor was adopted in this study to develop a convenient and efficient method for high-throughput screening of DAEase variants.

RESULTS

The catalytic activity and thermostability of DAEase from Caballeronia insecticola were simultaneously improved by semi-rational molecular modification. Compared with the wild-type enzyme, DAEaseS37N/F157Y variant exhibited 14.7% improvement in the catalytic activity and the half-time value (t1/2) at 65°C increased from 1.60 to 27.56 h by 17.23-fold. To our delight, the conversion rate of D-allulose was 33.6% from 500-g L-1 D-fructose in 1 h by Bacillus subtilis WB800 whole cells expressing this DAEase variant. Furthermore, the practicability of cell immobilization was evaluated and more than 80% relative activity of the immobilized cells was maintained from the second to seventh cycle.

CONCLUSION

All these results indicated that the DAEaseS37N/F157Y variant would be a potential candidate for the industrial production of D-allulose.

Molecular dynamics simulations reveal how vinculin refolds partially unfolded talin rod helices to stabilize them against mechanical force

Vinculin binds to specific sites of mechanically unfolded talin rod domains to reinforce the coupling of the cell's exterior to its force generation machinery. Force-dependent vinculin-talin complexation and dissociation was previously observed as contraction or extension of the unfolded talin domains respectively using magnetic tweezers. However, the structural mechanism underlying vinculin recognition of unfolded vinculin binding sites (VBSs) in talin remains unknown. Using molecular dynamics simulations, we demonstrate that a VBS dynamically refolds under force, and that vinculin can recognize and bind to partially unfolded VBS states. Vinculin binding enables refolding of the mechanically strained VBS and stabilizes its folded α-helical conformation, providing resistance against mechanical stress. Together, these results provide an understanding of a recognition mechanism of proteins unfolded by force and insight into the initial moments of how vinculin binds unfolded talin rod domains during the assembly of this mechanosensing meshwork.

A Novel AHAS-Inhibiting Herbicide Candidate for Controlling Leptochloa chinensis: A Devastating Weedy Grass in Rice Fields

Given the prevalence of the malignant weed Chinese Sprangletop (Leptochloa chinensis (L.) Nees) in rice fields, the development of novel herbicides against this weed has aroused wide interest. Here, we report a novel diphenyl ether-pyrimidine hybrid, DEP-5, serving as a systematic pre/postemergence herbicide candidate for broad-spectrum weed control in rice fields, specifically for L. chinensis. Notably, DEP-5 exhibits over 80% herbicidal activity against the resistant biotypes even at 37.5 g a.i./ha under greenhouse conditions and has complete control of L. chinensis at 150 g a.i./ha in the rice fields. We uncover that DEP-5 acts as a noncompetitive inhibitor of acetohydroxyacid synthase (AHAS) with an inhibition constant (Ki) of 39.4 μM. We propose that DEP-5 binds to AHAS in two hydrophobic-driven binding modes that differ from commercial AHAS inhibitors. Overall, these findings demonstrate that DEP-5 has great potential to be developed into a herbicide for L. chinensis control and inspire fresh concepts for novel AHAS-inhibiting herbicide design.

Pantothenate-encapsulated liposomes combined with exercise for effective inhibition of CRM1-mediated PKM2 translocation in Alzheimer's therapy

Alzheimer's disease (AD) is a complex neurodegenerative condition characterized by metabolic imbalances and neuroinflammation, posing a formidable challenge in medicine due to the lack of effective treatments. Despite considerable research efforts, a cure for AD remains elusive, with current therapies primarily focused on symptom management rather than addressing the disease's underlying causes. This study initially discerned, through Mendelian randomization analysis that elevating pantothenate levels significantly contributes to the prophylaxis of Alzheimer's disease. We explore the therapeutic potential of pantothenate encapsulated in liposomes (Pan@TRF@Liposome NPs), targeting the modulation of CRM1-mediated PKM2 nuclear translocation, a critical mechanism in AD pathology. Additionally, we investigate the synergistic effects of exercise, proposing a combined approach to AD treatment. Exercise-induced metabolic alterations share significant similarities with those associated with dementia, suggesting a potential complementary effect. The Pan@TRF@Liposome NPs exhibit notable biocompatibility, showing no liver or kidney toxicity in vivo, while demonstrating stability and effectiveness in modulating CRM1-mediated PKM2 nuclear translocation, thereby reducing neuroinflammation and neuronal apoptosis. The combined treatment of exercise and Pan@TRF@Liposome NP administration in an AD animal model leads to improved neurofunctional outcomes and cognitive performance. These findings highlight the nanoparticles' role as effective modulators of CRM1-mediated PKM2 nuclear translocation, with significant implications for mitigating neuroinflammation and neuronal apoptosis. Together with exercise, this dual-modality approach could offer new avenues for enhancing cognitive performance and neurofunctional outcomes in AD, marking a promising step forward in developing treatment strategies for this challenging disorder.

Structural insights into the molecular effects of the anthelmintics monepantel and betaine on the Caenorhabditis elegans acetylcholine receptor ACR-23

Anthelmintics are drugs used for controlling pathogenic helminths in animals and plants. The natural compound betaine and the recently developed synthetic compound monepantel are both anthelmintics that target the acetylcholine receptor ACR-23 and its homologs in nematodes. Here, we present cryo-electron microscopy structures of ACR-23 in apo, betaine-bound, and betaine- and monepantel-bound states. We show that ACR-23 forms a homo-pentameric channel, similar to some other pentameric ligand-gated ion channels (pLGICs). While betaine molecules are bound to the classical neurotransmitter sites in the inter-subunit interfaces in the extracellular domain, monepantel molecules are bound to allosteric sites formed in the inter-subunit interfaces in the transmembrane domain of the receptor. Although the pore remains closed in betaine-bound state, monepantel binding results in an open channel by wedging into the cleft between the transmembrane domains of two neighboring subunits, which causes dilation of the ion conduction pore. By combining structural analyses with site-directed mutagenesis, electrophysiology and in vivo locomotion assays, we provide insights into the mechanism of action of the anthelmintics monepantel and betaine.

MiMiC: A high-performance framework for multiscale molecular dynamics simulations

MiMiC is a framework for performing multiscale simulations in which loosely coupled external programs describe individual subsystems at different resolutions and levels of theory. To make it highly efficient and flexible, we adopt an interoperable approach based on a multiple-program multiple-data (MPMD) paradigm, serving as an intermediary responsible for fast data exchange and interactions between the subsystems. The main goal of MiMiC is to avoid interfering with the underlying parallelization of the external programs, including the operability on hybrid architectures (e.g., CPU/GPU), and keep their setup and execution as close as possible to the original. At the moment, MiMiC offers an efficient implementation of electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) that has demonstrated unprecedented parallel scaling in simulations of large biomolecules using CPMD and GROMACS as QM and MM engines, respectively. However, as it is designed for high flexibility with general multiscale models in mind, it can be straightforwardly extended beyond QM/MM. In this article, we illustrate the software design and the features of the framework, which make it a compelling choice for multiscale simulations in the upcoming era of exascale high-performance computing.

Tuning the Emission of Bis-ethylenedioxythiophene-thiophenes upon Aggregation

The ability of small lipophilic molecules to penetrate the blood-brain barrier through transmembrane diffusion has enabled researchers to explore new diagnostics and therapies for brain disorders. Until now, therapies targeting the brain have mainly relied on biochemical mechanisms, while electrical treatments such as deep brain stimulation often require invasive procedures. An alternative to implanting deep brain stimulation probes could involve administering small molecule precursors intravenously, capable of crossing the blood-brain barrier, and initiating the formation of conductive polymer networks in the brain through in vivo polymerization. This study examines the aggregation behavior of five water-soluble conducting polymer precursors sharing the same conjugate core but differing in side chains, using spectroscopy and various computational chemistry tools. Our findings highlight the significant impact of side chain composition on both aggregation and spectroscopic response.

Simulating the Skin Permeation Process of Ionizable Molecules

It is commonly assumed that ionizable molecules, such as drugs, permeate through the skin barrier in their neutral form. By using molecular dynamics simulations of the charged and neutral states separately, we can study the dynamic protonation behavior during the permeation process. We have studied three weak acids and three weak bases and conclude that the acids are ionized to a larger extent than the bases, when passing through the headgroup region of the lipid barrier structure, at pH values close to their pKa. It can also be observed that even if these dynamic protonation simulations are informative, in the cases studied herein they are not necessary for the calculation of permeability coefficients. It is sufficient to base the calculations only on the neutral form, as is commonly done.

Structural details of helix-mediated TDP-43 C-terminal domain multimerization

The primarily disordered C-terminal domain (CTD) of TAR DNA binding protein-43 (TDP-43), a key nuclear protein in RNA metabolism, forms neuronal inclusions in several neurodegenerative diseases. A conserved region (CR, spanning residues 319-341) in CTD forms transient helix-helix contacts important for its higher-order oligomerization and function that are disrupted by ALS-associated mutations. However, the structural details of CR assembly and the explanation for several ALS-associated variants' impact on phase separation and function remain unclear due to challenges in analyzing the dynamic association of TDP-43 CTD using traditional structural biology approaches. By employing an integrative approach, combining biophysical experiments, biochemical assays, AlphaFold2-Multimer (AF2-Multimer), and atomistic simulations, we generated structural models of helical oligomerization of TDP-43 CR. Using NMR, we first established that the native state of TDP-43 CR under physiological conditions is α-helical. Next, alanine scanning mutagenesis revealed that while hydrophobic residues in the CR are important for CR assembly, phase separation and TDP-43 nuclear retention function, polar residues down regulate these processes. Finally, pairing AF2-Multimer modeling with AAMD simulations indicated that dynamic, oligomeric assemblies of TDP-43 that are stabilized by a methionine-rich core with specific contributions from a tryptophan/leucine pair. In conclusion, our results advance the structural understanding of the mechanisms driving TDP-43 function and provide a window into the initial stages of its conversion into pathogenic aggregates.

Molecular Mechanism of Cynodon dactylon Phytosterols Targeting MAPK3 and PARP1 to Combat Epithelial Ovarian Cancer: A Multifaceted Computational Approach

Epithelial Ovarian Cancer (EOC) presents a global health concern, necessitating the development of innovative therapeutic strategies to combat its impact. This study was employed to investigate the unexplored therapeutic efficacy of Cynodon dactylon phytochemicals against EOC using a multifaceted computational approach. A total of 19 out of 89 rigorously curated phytochemicals were assessed as potential drug targets via ADMET profiling, while protein-protein interaction analysis scrutinized the top 20 hub genes among 264 disease targets, revealing their involvement in cancer-related pathways and underscoring their significance in EOC pathogenesis. In molecular docking, Stigmasterol acetate showed the highest binding affinity (-10.9 kcal/mol) with Poly [ADP-ribose] polymerase-1 (PDB: 1UK1), while Arundoin and Beta-Sitosterol exhibited strong affinities (-10.4 kcal/mol and -10.1 kcal/mol, respectively); additionally, Beta-Sitosterol interacting with Mitogen-activated protein kinase 3 (PDB: 4QTB) showed a binding affinity of -10.1 kcal/mol, forming 2 hydrogen bonds and a total of 10 bonds with 10 residues. Molecular dynamics simulations exhibited the significant structural stability of the Beta-Sitosterol-4QTB complex with superior binding free energy (-36.61 kcal/mol) among the three complexes. This study identified C. dactylon phytosterols, particularly Beta-Sitosterol, as effective in targeting MAPK3 and PARP1 to combat EOC, laying the groundwork for further experimental validation and drug development efforts.

Efficient and novel biosynthesis of myricetin α-triglucoside with improved solubility using amylosucrase from Deinococcus deserti

Although myricetin (3,3',4',5,5',7-hexahydroxyflavone, MYR) has a high antioxidant capacity and health functions, its use as a functional food material is limited owing to its low stability and water solubility. Amylosucrase (ASase) is capable of biosynthesizing flavonol α-glycoside using flavonols as acceptor molecules and sucrose as a donor molecule. Here, ASase from Deinococcus deserti (DdAS) efficiently biosynthesizes a novel MYR α-triglucoside (MYRαG3) using MYR as the acceptor molecule. Comparative homology analysis and computational simulation revealed that DdAS has a different active pocket for the transglycosylation reaction. DdAS produced MYRαG3 with a conversion efficiency of 67.4 % using 10 mM MYR and 50 mM sucrose as acceptor and donor molecules, respectively. The structure of MYRαG3 was identified as MYR 4'-O-4″,6″-tri-O-α-D-glucopyranoside using NMR and LC-MS. In silico analysis confirmed that DdAS has a distinct active pocket compared to other ASases. In addition, molecular docking simulations predicted the synthetic sequence of MYRαG3. Furthermore, MYRαG3 showed a similar DPPH radical scavenging activity of 49 %, comparable to MYR, but with significantly higher water solubility, which increased from 0.03 μg/mL to 511.5 mg/mL. In conclusion, this study demonstrated the efficient biosynthesis of a novel MYRαG3 using DdAS and highlighted the potential of MYRαG3 as a functional material.

GENESIS 2.1: High-Performance Molecular Dynamics Software for Enhanced Sampling and Free-Energy Calculations for Atomistic, Coarse-Grained, and Quantum Mechanics/Molecular Mechanics Models

GENeralized-Ensemble SImulation System (GENESIS) is a molecular dynamics (MD) software developed to simulate the conformational dynamics of a single biomolecule, as well as molecular interactions in large biomolecular assemblies and between multiple biomolecules in cellular environments. To achieve the latter purpose, the earlier versions of GENESIS emphasized high performance in atomistic MD simulations on massively parallel supercomputers, with or without graphics processing units (GPUs). Here, we implemented multiscale MD simulations that include atomistic, coarse-grained, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. They demonstrate high performance and are integrated with enhanced conformational sampling algorithms and free-energy calculations without using external programs except for the QM programs. In this article, we review new functions, molecular models, and other essential features in GENESIS version 2.1 and discuss ongoing developments for future releases.

Towards new bioactive fluorine-containing 1,3,4-oxadiazole-amide derivatives: synthesis, antibacterial activity, molecular docking and molecular dynamics simulation study

Through the approach of molecular hybridization, this study rationally designed and synthesized new trifluoromethyl-1,3,4-oxadiazole amide derivatives, denoted as 1a-1n. The findings reveal that these novel molecules exhibit potent inhibitory effects against various bacterial strains. Thereinto, compounds 1c, 1d, 1i, 1j and 1n, demonstrate relatively superior antimicrobial performance against B. cereus FM314, with a minimum inhibitory concentration (MIC) of 0.03907 μg/mL. Molecular docking analysis suggests the potential importance of the Ser57 and Thr125 amino acid residues (PDB ID: 4EI9) in contributing to the inhibitory activity against B. cereus. The consistency of these results was further corroborated through subsequent molecular dynamics simulations and MMPBSA validations. The insights gained from this study serve to facilitate the rational design and efficient development of novel eco-friendly antimicrobial inhibitors based on the trifluoromethyl-1,3,4-oxadiazole amide scaffold.

Identification of potential modulators for human GPD1 by docking-based virtual screening, molecular dynamics simulations, binding free energy calculations, and DeLA-drug analysis

Cytosolic Glycerol-3-phosphate dehydrogenase 1 (GPD1, EC 1.1.1.8) plays a pivotal role in regulating the Embden-Meyerhof glucose glycolysis pathway (E-M pathway), as well as in conditions such as Huntington's disease, cancer, and its potential role as a specific marker for Dormant Glioma Stem Cells. In this study, we conducted virtual screening using the ZINC database ( http://zinc.docking.org/ ) and the GPD1 structure to identify potential GPD1 modulators. The investigation involved screening active candidate ligands using ADMET (Absorption, Distribution, Metabolism, Excretion, Toxicity) parameters, combined with molecular docking, pose analysis, and interaction analysis based on Lipinski and Veber criteria. Subsequently, the top 10 ligands were subjected to 200 ns all-atom molecular dynamics (M.D.) simulations, and binding free energies were calculated. The findings revealed that specific residues, namely TRP14, PRO94, LYS120, ASN151, THR264, ASP260, and GLN298, played a crucial role in ensuring system stability. Furthermore, through a comprehensive analysis involving molecular docking, molecular M.D., and DeLA-Drug, we identified 10 promising small molecules. These molecules represent potential lead compounds for developing effective therapeutics targeting GPD1-associated diseases, thereby contributing to a deeper understanding of GPD1-associated mechanisms. This study's significance lies in identifying key residues associated with GPD1 and discovering valuable small molecules, providing a foundation for further research and development.

Recent advances in the development of antimicrobial peptides against ESKAPE pathogens

Multi-drug resistant ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) are a global health threat. The severity of the problem lies in its impact on mortality, therapeutic limitations, the threat to public health, and the costs associated with managing infections caused by these resistant strains. Effectively addressing this challenge requires innovative approaches to research, the development of new antimicrobials, and more responsible antibiotic use practices globally. Antimicrobial peptides (AMPs) are a part of the innate immune system of all higher organisms. They are short, cationic and amphipathic molecules with broad-spectrum activity. AMPs interact with the negatively charged bacterial membrane. In recent years, AMPs have attracted considerable interest as potential antibiotics. However, AMPs have low bioavailability and short half-lives, which may be circumvented by chemical modification. This review presents recent in vitro and in silico strategies for the modification of AMPs to improve their stability and application in preclinical experiments.

Molecular Dynamics Simulation of the Interaction within the Carboxymethyl Cellulose-Modified Montmorillonite Lamellae at Aggressive CuCl2 Concentrations

To elucidate the degradation mechanism of the CMC-modified MMT composite at aggressive Cu2+ concentrations, large scale molecular dynamics simulation was conducted for CuCl2 concentrations ranging from 0 to 800 mM. Both CMC and MMT followed the Langmuir isotherm for Cu2+ adsorption, and the adsorption capacity of CMC (8.75 mmol/g) was much higher than that of MMT (0.83 mmol/g). Despite the CMC mass ratio being only 4.1%, it adsorbed up to 34.3% of the total adsorbed Cu2+. The Cu2+ attraction ability hierarchy of oxygen-containing functional groups in the CMC is as follows: carboxylic oxygens > alcoholic oxygens > carbinolic oxygens > bridging oxygens > glucose oxygens. Carboxyls were the most effective in chelating and complexing with Cu2+, and they could be intentionally added in artificially synthesized polymer-MMT composites for Cu2+ containment. Formation of the Cu2+ cation bridge between CMC and MMT at aggressive CuCl2 concentrations contributed to the transition of CMC density distribution from unimodality to bimodality and enhanced resistance of polymer elution. As the CuCl2 concentration increased, the stoichiometric ratio between the chelated Cu2+ and carboxylic oxygens increased from 1:2 to 1:1, suggesting the evolution of the Cu2+ chelation mechanism.

Structure and Ligand Based Virtual Screening and MPO Topological Analysis of Triazolo Thiadiazepine-fused Coumarin Derivatives as Anti-Parkinson Drug Candidates

Parkinson's disease (PD) is a debilitating condition that can cause locomotor problems in affected patients, such as tremors and body rigidity. PD therapy often includes the use of monoamine oxidase B (MAOB) inhibitors, particularly phenylhalogen compounds and coumarin-based semi-synthetic compounds. The objective of this study was to analyze the structural, pharmacokinetic, and pharmacodynamic profile of a series of Triazolo Thiadiazepine-fused Coumarin Derivatives (TDCDs) against MAOB, in comparison with the inhibitor safinamide. To achieve this goal, we utilized structure-based virtual screening techniques, including target prediction and absorption, distribution, metabolism, and excretion (ADME) prediction based on multi-parameter optimization (MPO) topological analysis, as well as ligand-based virtual screening techniques, such as docking and molecular dynamics. The findings indicate that the TDCDs exhibit structural similarity to other bioactive compounds containing coumarin and MAOB-binding azoles, which are present in the ChEMBL database. The topological analyses suggest that TDCD3 has the best ADME profile, particularly due to the alignment between low lipophilicity and high polarity. The coumarin and triazole portions make a strong contribution to this profile, resulting in a permeability with Papp estimated at 2.15 × 10-5 cm/s, indicating high cell viability. The substance is predicted to be metabolically stable. It is important to note that this is an objective evaluation based on the available data. Molecular docking simulations showed that the ligand has an affinity energy of - 8.075 kcal/mol with MAOB and interacts with biological substrate residues such as Pro102 and Phe103. The results suggest that the compound has a safe profile in relation to the MAOB model, making it a promising active ingredient for the treatment of PD.

From GPUs to AI and quantum: three waves of acceleration in bioinformatics

The enormous growth in the amount of data generated by the life sciences is continuously shifting the field from model-driven science towards data-driven science. The need for efficient processing has led to the adoption of massively parallel accelerators such as graphics processing units (GPUs). Consequently, the development of bioinformatics methods nowadays often heavily depends on the effective use of these powerful technologies. Furthermore, progress in computational techniques and architectures continues to be highly dynamic, involving novel deep neural network models and artificial intelligence (AI) accelerators, and potentially quantum processing units in the future. These are expected to be disruptive for the life sciences as a whole and for drug discovery in particular. Here, we identify three waves of acceleration and their applications in a bioinformatics context: (i) GPU computing, (ii) AI and (iii) next-generation quantum computers.

Structural basis for selectivity and antagonism in extracellular GPCR-nanobodies

G protein-coupled receptors (GPCRs) are pivotal therapeutic targets, but their complex structure poses challenges for effective drug design. Nanobodies, or single-domain antibodies, have emerged as a promising therapeutic strategy to target GPCRs, offering advantages over traditional small molecules and antibodies. However, an incomplete understanding of the structural features enabling GPCR-nanobody interactions has limited their development. In this study, we investigate VUN701, a nanobody antagonist targeting the atypical chemokine receptor 3 (ACKR3). We determine that an extended CDR3 loop is required for ACKR3 binding. Uncommon in most nanobodies, an extended CDR3 is prevalent in GPCR-targeting nanobodies. Combining experimental and computational approaches, we map an inhibitory ACKR3-VUN701 interface and define a distinct conformational mechanism for GPCR inactivation. Our results provide insights into class A GPCR-nanobody selectivity and suggest a strategy for the development of these new therapeutic tools.

The NOMAD mini-apps: A suite of kernels from ab initio electronic structure codes enabling co-design in high-performance computing

This article introduces a suite of mini-applications (mini-apps) designed to optimise computational kernels in ab initio electronic structure codes. The suite is developed from flagship applications participating in the NOMAD Center of Excellence, such as the ELPA eigensolver library and the GW implementations of the exciting, Abinit, and FHI-aims codes. The mini-apps were identified by targeting functions that significantly contribute to the total execution time in the parent applications. This strategic selection allows for concentrated optimisation efforts. The suite is designed for easy deployment on various High-Performance Computing (HPC) systems, supported by an integrated CMake build system for straightforward compilation and execution. The aim is to harness the capabilities of emerging (post)exascale systems, which necessitate concurrent hardware and software development - a concept known as co-design. The mini-app suite serves as a tool for profiling and benchmarking, providing insights that can guide both software optimisation and hardware design. Ultimately, these developments will enable more accurate and efficient simulations of novel materials, leveraging the full potential of exascale computing in material science research.

Imidazole[1,5-a]pyridine derivatives as EGFR tyrosine kinase inhibitors unraveled by umbrella sampling and steered molecular dynamics simulations

Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.

Intra-condensate demixing of TDP-43 inside stress granules generates pathological aggregates

Cytosolic aggregation of the nuclear protein TDP-43 is associated with many neurodegenerative diseases, but the triggers for TDP-43 aggregation are still debated. Here, we demonstrate that TDP-43 aggregation requires a double event. One is up-concentration in stress granules beyond a threshold, and the other is oxidative stress. These two events collectively induce intra-condensate demixing, giving rise to a dynamic TDP-43 enriched phase within stress granules, which subsequently transitions into pathological aggregates. Mechanistically, intra-condensate demixing is triggered by local unfolding of the RRM1 domain for intermolecular disulfide bond formation and by increased hydrophobic patch interactions in the C-terminal domain. By engineering TDP-43 variants resistant to intra-condensate demixing, we successfully eliminate pathological TDP-43 aggregates in cells. We conclude that up-concentration inside condensates and simultaneous exposure to environmental stress could be a general pathway for protein aggregation, with intra-condensate demixing constituting a key intermediate step.

Solutes unmask differences in clustering versus phase separation of FET proteins

Phase separation and percolation contribute to phase transitions of multivalent macromolecules. Contributions of percolation are evident through the viscoelasticity of condensates and through the formation of heterogeneous distributions of nano- and mesoscale pre-percolation clusters in sub-saturated solutions. Here, we show that clusters formed in sub-saturated solutions of FET (FUS-EWSR1-TAF15) proteins are affected differently by glutamate versus chloride. These differences on the nanoscale, gleaned using a suite of methods deployed across a wide range of protein concentrations, are prevalent and can be unmasked even though the driving forces for phase separation remain unchanged in glutamate versus chloride. Strikingly, differences in anion-mediated interactions that drive clustering saturate on the micron-scale. Beyond this length scale the system separates into coexisting phases. Overall, we find that sequence-encoded interactions, mediated by solution components, make synergistic and distinct contributions to the formation of pre-percolation clusters in sub-saturated solutions, and to the driving forces for phase separation.

Deep profiling of potential substrate atlas of porcine epidemic diarrhea virus 3C-like protease

Coronavirus (CoV) 3C-like protease (3CLpro) is essential for viral replication and is involved in immune escape by proteolyzing host proteins. Deep profiling the 3CLpro substrates in the host proteome extends our understanding of viral pathogenesis and facilitates antiviral drug discovery. Here, 3CLpro from porcine epidemic diarrhea virus (PEDV), an enteropathogenic CoV, was used as a model which to identify the potential 3CLpro cleavage motifs in all porcine proteins. We characterized the selectivity of PEDV 3CLpro at sites P5-P4'. We then compiled the 3CLpro substrate preferences into a position-specific scoring matrix and developed a 3CLpro profiling strategy to delineate the protein substrate landscape of CoV 3CLpro. We identified 1,398 potential targets in the porcine proteome containing at least one putative cleavage site and experimentally validated the reliability of the substrate degradome. The PEDV 3CLpro-targeted pathways are involved in mRNA processing, translation, and key effectors of autophagy and the immune system. We also demonstrated that PEDV 3CLpro suppresses the type 1 interferon (IFN-I) cascade via the proteolysis of multiple signaling adaptors in the retinoic acid-inducible gene I (RIG-I) signaling pathway. Our composite method is reproducible and accurate, with an unprecedented depth of coverage for substrate motifs. The 3CLpro substrate degradome establishes a comprehensive substrate atlas that will accelerate the investigation of CoV pathogenicity and the development of anti-CoV drugs.IMPORTANCECoronaviruses (CoVs) are major pathogens that infect humans and animals. The 3C-like protease (3CLpro) encoded by CoV not only cleaves the CoV polyproteins but also degrades host proteins and is considered an attractive target for the development of anti-CoV drugs. However, the comprehensive characterization of an atlas of CoV 3CLpro substrates is a long-standing challenge. Using porcine epidemic diarrhea virus (PEDV) 3CLpro as a model, we developed a method that accurately predicts the substrates of 3CLpro and comprehensively maps the substrate degradome of PEDV 3CLpro. Interestingly, we found that 3CLpro may simultaneously degrade multiple molecules responsible for a specific function. For instance, it cleaves at least four adaptors in the RIG-I signaling pathway to suppress type 1 interferon production. These findings highlight the complexity of the 3CLpro substrate degradome and provide new insights to facilitate the development of anti-CoV drugs.

Treatment of flexibility of protein backbone in simulations of protein-ligand interactions using steered molecular dynamics

To ensure that an external force can break the interaction between a protein and a ligand, the steered molecular dynamics simulation requires a harmonic restrained potential applied to the protein backbone. A usual practice is that all or a certain number of protein's heavy atoms or Cα atoms are fixed, being restrained by a small force. This present study reveals that while fixing both either all heavy atoms and or all Cα atoms is not a good approach, while fixing a too small number of few atoms sometimes cannot prevent the protein from rotating under the influence of the bulk water layer, and the pulled molecule may smack into the wall of the active site. We found that restraining the Cα atoms under certain conditions is more relevant. Thus, we would propose an alternative solution in which only the Cα atoms of the protein at a distance larger than 1.2 nm from the ligand are restrained. A more flexible, but not too flexible, protein will be expected to lead to a more natural release of the ligand.

Transient interdomain interactions modulate the monomeric structural ensemble and oligomerization landscape of Huntingtin Exon 1

Polyglutamine expansion (≥ 36 residues) within the N-terminal exon-1 of Huntingtin (Httex1) leads to Huntington's disease, a neurogenerative condition marked by the presence of intranuclear Htt inclusions. Notably, the polyglutamine tract in Httex1 is flanked by an N-terminal coiled-coil domain - N17 (17 amino acids), which undergoes self-association to promote the formation of soluble Httex1 oligomers and brings the aggregation-prone polyQ tracts in close spatial proximity. However, the mechanisms underlying the subsequent conversion of soluble oligomers into insoluble β-rich aggregates with increasing polyQ length, remain unclear. Current knowledge suggests that expansion of the polyQ tract increases its helicity, and this favors its oligomerization and aggregation. In addition, studies utilizing conformation-specific antibodies and a stable coiled-coil heterotetrametric system fused to polyQ indicate that domain "cross-talk" (i.e., interdomain interactions) may be necessary to efficiently promote the emergence of toxic conformations (in monomers and oligomers) and fibrillar aggregation. Here, we performed extensive atomistic molecular dynamics (MD) simulations (aggregate time ∼ 0.7 ms) of N17-polyQ fragments to uncover the interplay between structural transformation and domain "cross-talk" on the monomeric structural ensemble and oligomerization landscape of Httex1. Our simulation ensembles of N17-polyQ monomers validated against 13 C NMR chemical shifts indicated that in addition to elevated α-helicity, polyQ expansion also favors transient, interdomain (N17-polyQ) interactions which result in the emergence of β-conformations. Further, interdomain interactions decreased the overall stability of N17-mediated dimers by counteracting the stabilizing effect of increased α-helicity and promoted a heterogenous oligomerization landscape on the sub-microsecond timescale. Overall, our study uncovers the significance of domain "cross-talk" in modulating the monomeric conformational ensemble and oligomerization landscape of Httex1 to favor the formation of amyloid aggregates.

Deciphering the mechanisms and interactions of the endocrine disruptor bisphenol A and its analogs with the androgen receptor

Bisphenol A (BPA) and its various forms used as BPA alternatives in industries are recognized toxic compounds and antiandrogenic endocrine disruptors. These chemicals are widespread in the environment and frequently detected in biological samples. Concerns exist about their impact on hormones, disrupting natural biological processes in humans, together with their negative impacts on the environment and biotic life. This study aims to characterize the interaction between BPA analogs and the androgen receptor (AR) and the effect on the receptor's normal activity. To achieve this goal, molecular docking was conducted with BPA and its analogs and dihydrotestosterone (DHT) as a reference ligand. Four BPA analogs exhibited higher affinity (-10.2 to -8.7 kcal/mol) for AR compared to BPA (-8.6 kcal/mol), displaying distinct interaction patterns. Interestingly, DHT (-11.0 kcal/mol) shared a binding pattern with BPA. ADMET analysis of the top 10 compounds, followed by molecular dynamics simulations, revealed toxicity and dynamic behavior. Experimental studies demonstrated that only BPA disrupts DHT-induced AR dimerization, thereby affecting AR's function due to its binding nature. This similarity to DHT was observed during computational analysis. These findings emphasize the importance of targeted strategies to mitigate BPA toxicity, offering crucial insights for interventions in human health and environmental well-being.

Exploring the neurotoxicity of chiral dinotefuran towards nicotinic acetylcholine receptors: Enantioselective insights into species selectivity

Dinotefuran is a chiral neonicotinoid that is widely distributed in environmental matrices, but its health risks to different organisms are poorly understood. This study investigated the neurotoxic responses of honeybee/cotton aphid nicotinic acetylcholine receptors (nAChRs) to chiral dinotefuran at the enantiomeric scale and demonstrated the microscopic mechanism of species selectivity in nAChR-mediated enantioselective neurotoxicity. The findings indicated that (S)-dinotefuran had a higher affinity for honeybee nAChR than (R)-dinotefuran whereas both enantiomers exhibited similar bioactivity toward cotton aphid nAChR. The results of dynamic neurotoxic processes indicated the association of conformational changes induced by chiral dinotefuran with its macroscopic neurotoxicity, and (R)-dinotefuran, which exhibit low toxicity to honeybee, was found to induce significant conformational changes in the enantioselective neurotoxic reaction, as supported by the average root-mean-square fluctuation (0.35 nm). Energy decomposition results indicated that electrostatic contribution (ΔGele) is the critical energy term that leads to substantial enantioselectivity, and both Trp-51 (－2.57 kcal mol-1) and Arg-75 (－4.86 kcal mol-1), which form a hydrogen-bond network, are crucial residues in mediating the species selectivity for enantioselective neurotoxic responses. Clearly, this study provides experimental evidence for a comprehensive assessment of the health hazards of chiral dinotefuran.

Identification of Helicobacter pylori-carcinogenic TNF-alpha-inducing protein inhibitors via daidzein derivatives through computational approaches

Gastric cancer is considered a class 1 carcinogen that is closely linked to infection with Helicobacter pylori (H. pylori), which affects over 1 million people each year. However, the major challenge to fight against H. pylori and its associated gastric cancer due to drug resistance. This research gap had led our research team to investigate a potential drug candidate targeting the Helicobacter pylori-carcinogenic TNF-alpha-inducing protein. In this study, a total of 45 daidzein derivatives were investigated and the best 10 molecules were comprehensively investigated using in silico approaches for drug development, namely pass prediction, quantum calculations, molecular docking, molecular dynamics simulations, Lipinski rule evaluation, and prediction of pharmacokinetics. The molecular docking study was performed to evaluate the binding affinity between the target protein and the ligands. In addition, the stability of ligand-protein complexes was investigated by molecular dynamics simulations. Various parameters were analysed, including root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), hydrogen bond analysis, principal component analysis (PCA) and dynamic cross-correlation matrix (DCCM). The results has confirmed that the ligand-protein complex CID: 129661094 (07) and 129664277 (08) formed stable interactions with the target protein. It was also found that CID: 129661094 (07) has greater hydrogen bond occupancy and stability, while the ligand-protein complex CID 129664277 (08) has greater conformational flexibility. Principal component analysis revealed that the ligand-protein complex CID: 129661094 (07) is more compact and stable. Hydrogen bond analysis revealed favourable interactions with the reported amino acid residues. Overall, this study suggests that daidzein derivatives in particular show promise as potential inhibitors of H. pylori.

In silico design of high-affinity antigenic peptides for HLA-B44

Cancer cell-killing by CD8+ T cells demands effective tumor antigen presentation by human leukocyte antigen class I (HLA-I) molecules. Screening and designing highly immunogenic neoantigens require quantitative computations to reliably predict HLA-peptide binding affinities. Here, with all-atom molecular dynamics (MD) simulations and free energy perturbation (FEP) methods, we design a collection of antigenic peptide candidates through in silico mutagenesis studies on immunogenic neoantigens, yielding enhanced binding affinities to HLA-B*44:02. In-depth structural dissection shows that introducing positively charged residues such as arginine to position 6 or lysine to position 7 of the candidates triggers conformational shifts in both peptides and the antigen-binding groove of the HLA, following the "induced-fit" mechanism. Enhancement in binding affinities compared to the wild-type was found in three out of five mutated candidates. The HLA pocket, capable of accommodating positively charged residues in positions from 5 to 7, is designated as the "dynamic pocket". Taken together, we showcase an effective structure-based binding affinity optimization framework for antigenic peptides of HLA-B*44:02 and underscore the importance of dynamic nature of the antigen-binding groove in concert with the anchoring motifs. This work provides structural insights for rational design of favorable HLA-peptide bindings and future developments in neoantigen-based therapeutics.

High-activity recombinant human carboxypeptidase B expression in Pichia pastoris through rational protein engineering and enhancing secretion from the Golgi apparatus to the plasma membrane

Human carboxypeptidase B1 (hCPB1) is vital for recombinant insulin production, holding substantial value in the pharmaceutical industry. Current challenges include limited hCPB1 enzyme activity. In this study, recombinant hCPB1 efficient expression in Pichia pastoris was achieved. To enhance hCPB1 secretion, we conducted signal peptides screening and deleted the Vps10 sortilin domain, reducing vacuolar mis-sorting. Overexpression of Sec4p increased the fusion of secretory vesicles with the plasma membrane and improved hCPB1 secretion by 20%. Rational protein engineering generated twenty-two single-mutation mutants and identified the A178L mutation resulted in a 30% increase in hCPB1 specific activity. However, all combinational mutations that increased specific activities decreased protein expression levels. Therefore, computer-aided global protein design with PROSS was employed for the aim of improving specific activities and preserving good protein expression. Among the six designed mutants, hCPB1-P6 showed a remarkable 114% increase in the catalytic rate constant (kcat), a 137% decrease in the Michaelis constant (Km), and a 490% increase in catalytic efficiency. Most mutations occurred on the surface of hCPB1-P6, with eight sites mutated to proline. In a 5 L fermenter, hCPB1-P6 was produced by the secretion-enhanced P. pastoris chassis to 199.6 ± 20 mg L-1 with a specific activity of 96 ± 0.32 U mg-1, resulting in a total enzyme activity of 19137 ± 1131 U L-1, demonstrating significant potential for industrial applications.

PaCS-Toolkit: Optimized Software Utilities for Parallel Cascade Selection Molecular Dynamics (PaCS-MD) Simulations and Subsequent Analyses

Parallel cascade selection molecular dynamics (PaCS-MD) is an enhanced conformational sampling method conducted as a "repetition of time leaps in parallel worlds", comprising cycles of multiple molecular dynamics (MD) simulations performed in parallel and selection of the initial structures of MDs for the next cycle. We developed PaCS-Toolkit, an optimized software utility enabling the use of different MD software and trajectory analysis tools to facilitate the execution of the PaCS-MD simulation and analyze the obtained trajectories, including the preparation for the subsequent construction of the Markov state model. PaCS-Toolkit is coded with Python, is compatible with various computing environments, and allows for easy customization by editing the configuration file and specifying the MD software and analysis tools to be used. We present the software design of PaCS-Toolkit and demonstrate applications of PaCS-MD variations: original targeted PaCS-MD to peptide folding; rmsdPaCS-MD to protein domain motion; and dissociation PaCS-MD to ligand dissociation from adenosine A2A receptor.