Electrostatics of nanosystems: application to microtubules and the ribosome

Rationally Engineering pH Adaptation of Acid-Induced Arginine Decarboxylase from Escherichia coli to Alkaline Environments to Efficiently Biosynthesize Putrescine

Acid-induced arginine decarboxylase AdiA is a typical homo-oligomeric protein biosynthesizing alkaline nylon monomer putrescine. However, upon loss of the AdiA decamer oligomeric state at neutral and alkaline conditions the activity also diminishes, obstructing the whole-cell biosynthesis of alkaline putrescine. Here, a structure cohesion strategy is proposed to change the pH adaptation of AdiA to alkaline environments based on the rational engineering of meridional and latitudinal oligomerization interfaces. After integrating substitutions of E467K at the latitudinal interface and H736E at the meridional channel interface, the structural stability of AdiA decamer and its substrate transport efficiency at neutral and alkaline conditions are improved. Finally, E467K_H736E is well adapted to neutral and alkaline environments (pH 7.0-9.0), and its enzymatic activity is 35-fold higher than that of wild AdiA at pH 8.0. Using E467K_H736E in the putrescine synthesis pathway, the titer of putrescine is up to 128.9 g·L-1 with a conversion of 0.94 mol·mol-1 in whole-cell catalysis. Additionally, the neutral pH adaptation of lysine decarboxylase, with a decamer structure similar to AdiA, is also improved using this cohesion strategy, providing an option for pH-adaptation engineering of other oligomeric decarboxylases.

Localization of four class I glutaredoxins in the cytosol and the secretory pathway and characterization of their biochemical diversification

Class I glutaredoxins (GRXs) are catalytically active oxidoreductases and considered key proteins mediating reversible glutathionylation and deglutathionylation of protein thiols during development and stress responses. To narrow in on putative target proteins, it is mandatory to know the subcellular localization of the respective GRXs and to understand their catalytic activities and putative redundancy between isoforms in the same compartment. We show that in Arabidopsis thaliana, GRXC1 and GRXC2 are cytosolic proteins with GRXC1 being attached to membranes through myristoylation. GRXC3 and GRXC4 are identified as type II membrane proteins along the early secretory pathway with their enzymatic function on the luminal side. Unexpectedly, neither single nor double mutants lacking both GRXs isoforms in the cytosol or the ER show phenotypes that differ from wild-type controls. Analysis of electrostatic surface potentials and clustering of GRXs based on their electrostatic interaction with roGFP2 mirrors the phylogenetic classification of class I GRXs, which clearly separates the cytosolic GRXC1 and GRXC2 from the luminal GRXC3 and GRXC4. Comparison of all four studied GRXs for their oxidoreductase function highlights biochemical diversification with GRXC3 and GRXC4 being better catalysts than GRXC1 and GRXC2 for the reduction of bis(2-hydroxyethyl) disulfide. With oxidized roGFP2 as an alternative substrate, GRXC1 and GRXC2 catalyze the reduction faster than GRXC3 and GRXC4, which suggests that catalytic efficiency of GRXs in reductive reactions depends on the respective substrate. Vice versa, GRXC3 and GRXC4 are faster than GRXC1 and GRXC2 in catalyzing the oxidation of pre-reduced roGFP2 in the reverse reaction.

Does Acinetobacter calcoaceticus glucose dehydrogenase produce self-damaging H2O2?

The soluble glucose dehydrogenase (sGDH) from Acinetobacter calcoaceticus has been widely studied and is used, in biosensors, to detect the presence of glucose, taking advantage of its high turnover and insensitivity to molecular oxygen. This approach, however, presents two drawbacks: the enzyme has broad substrate specificity (leading to imprecise blood glucose measurements) and shows instability over time (inferior to other oxidizing glucose enzymes). We report the characterization of two sGDH mutants: the single mutant Y343F and the double mutant D143E/Y343F. The mutants present enzyme selectivity and specificity of 1.2 (Y343F) and 5.7 (D143E/Y343F) times higher for glucose compared with that of the wild-type. Crystallographic experiments, designed to characterize these mutants, surprisingly revealed that the prosthetic group PQQ (pyrroloquinoline quinone), essential for the enzymatic activity, is in a cleaved form for both wild-type and mutant structures. We provide evidence suggesting that the sGDH produces H2O2, the level of production depending on the mutation. In addition, spectroscopic experiments allowed us to follow the self-degradation of the prosthetic group and the disappearance of sGDH's glucose oxidation activity. These studies suggest that the enzyme is sensitive to its self-production of H2O2. We show that the premature aging of sGDH can be slowed down by adding catalase to consume the H2O2 produced, allowing the design of a more stable biosensor over time. Our research opens questions about the mechanism of H2O2 production and the physiological role of this activity by sGDH.

Designing of xanthine-based DPP-4 inhibitors: a structure-guided alignment dependent Multifacet 3D-QSAR modeling, and molecular dynamics simulation study

The DPP-4 enzyme degrades incretin hormones GLP-1 and GIP. DPP-4 inhibitors are found effective in the prevention of the degradation of incretins. Xanthine scaffold-bearing molecules are reported as potential DPP-4 inhibitors for treating type 2 diabetes mellitus, e.g. the marketed drug linagliptin. In this work, structure-guided alignment-dependent atom- and Gaussian field-based 3D-QSAR have been performed on a dataset of 75 molecules. The robustness and predictive ability of the developed multifacet 3D-QSAR models were validated on different statistical parameters and found to be statistically fit. The favorable and unfavorable pharmacophoric features were mapped for each multifacet 3D-QSAR model based on three alignment sets (1-3). A five-point common pharmacophore hypothesis was generated separately for each set of alignments. The molecular dynamics simulations (up to 100 ns) were performed for the potent molecule from each alignment set (Compounds 12, 40 and 57) compared to reference standard linagliptin to study the binding energy and stability of target-ligand complexes. The MM-PBSA calculations revealed that the binding free energy and stability of compounds 12 (-40.324 ± 17.876 kJ/mol), 40 (-80.543 ± 21.782 kJ/mol) and 57 (-50.202 ± 16.055 kJ/mol) were better than the reference drug linagliptin (-20.390 ± 63.200 kJ/mol). The generated contour maps from structure-guided alignment-dependent multifacet 3D-QSAR models offer information about the structure-activity relationship (SAR) and ligand-target binding energy and stability data from MD simulation may be utilized to design and develop target selective xanthine-based novel DPP-4 inhibitors.Communicated by Ramaswamy H. Sarma.

Discovery of a New Class of Lipophilic Pyrimidine-Biphenyl Herbicides Using an Integrated Experimental-Computational Approach

Herbicides are useful tools for managing weeds and promoting food production and sustainable agriculture. In this study, we report on the development of a novel class of lipophilic pyrimidine-biphenyl (PMB) herbicides. Firstly, three PMBs, Ia, IIa, and IIIa, were rationally designed via a scaffold hopping strategy and were determined to inhibit acetohydroxyacid synthase (AHAS). Computational simulation was carried out to investigate the molecular basis for the efficiency of PMBs against AHAS. With a rational binding mode, and the highest in vitro as well as in vivo potency, Ia was identified as a preferable hit. Furthermore, these integrated analyses guided the design of eighteen new PMBs, which were synthesized via a one-step Suzuki-Miyaura cross-coupling reaction. These new PMBs, Iba-ic, were more effective in post-emergence control of grass weeds compared with Ia. Interestingly, six of the PMBs displayed 98-100% inhibition in the control of grass weeds at 750 g ai/ha. Remarkably, Ica exhibited ≥ 80% control against grass weeds at 187.5 g ai/ha. Overall, our comprehensive and systematic investigation revealed that a structurally distinct class of lipophilic PMB herbicides, which pair excellent herbicidal activities with new interactions with AHAS, represent a noteworthy development in the pursuit of sustainable weed control solutions.

In-depth investigation of the effect of pH on the autofluorescence properties of DPF3b and DPF3a amyloid fibrils

Double PHD fingers 3 (DPF3) protein exists as two splicing variants, DPF3b and DPF3a, the involvement of which in human cancer and neurodegeneration is beginning to be increasingly recognised. Both isoforms have recently been identified as intrinsically disordered proteins able to undergo amyloid fibrillation. Upon their aggregation, DPF3 proteins exhibit an intrinsic fluorescence in the visible range, referred to as deep-blue autofluorescence (dbAF). Comprehension of such phenomenon remaining elusive, we investigated in the present study the influence of pH on the optical properties of DPF3b and DPF3a fibrils. By varying the excitation wavelength and the pH condition, the two isoforms were revealed to display several autofluorescence modes that were defined as violet, deep-blue, and blue-green according to their emission range. Complementarily, analysis of excitation spectra and red edge shift plots allowed to better decipher their photoselection mechanism and to highlight isoform-specific excitation-emission features. Furthermore, the observed violation to Kasha-Vavilov's rule was attributed to red edge excitation shift effects, which were impacted by pH-mediated H-bond disruption, leading to changes in intramolecular charge and proton transfer, or π-electrons delocalisation. Finally, emergence of different autofluorescence emitters was likely related to structurally distinct fibrillar assemblies between isoforms, as well as to discrepancies in the amino acid composition of their aggregation prone regions.

Directional molecular transport in iron redox flow batteries by interfacial electrostatic forces

The mounting global energy demand urges surplus electricity generation. Due to dwindling fossil resources and environmental concerns, shifting from carbon-based fuels to renewables is vital. Though renewables are affordable, their intermittent nature poses supply challenges. In these contexts, aqueous flow batteries (AFBs), are a viable energy storage solution. This study tackles AFBs' energy density and efficiency challenges. Conventional strategies focus on altering molecule's solubility but overlook interface's transport kinetics. We show that triggering electrostatic forces at the interface can significantly enhance the mass transport kinetics of redox active molecules by introducing a powerful electrostatic flux over the diffusional flux, thereby exerting a precise directionality on the molecular transport. This approach of controlling the directionality of molecular flux in an all iron redox flow battery amplifies the current and power rating with approximately 140 % enhancement in the energy density.

Leveraging Bidirectional Nature of Allostery To Inhibit Protein-Protein Interactions (PPIs): A Case Study of PCSK9-LDLR Interaction

The protein PCSK9 (proprotein convertase subtilisin/Kexin type 9) negatively regulates the recycling of LDLR (low-density lipoprotein receptor), leading to an elevated plasma level of LDL. Inhibition of PCSK9-LDLR interaction has emerged as a promising therapeutic strategy to manage hypercholesterolemia. However, the large interaction surface area between PCSK9 and LDLR makes it challenging to identify a small molecule competitive inhibitor. An alternative strategy would be to identify distal cryptic sites as targets for allosteric inhibitors that can remotely modulate PCSK9-LDLR interaction. Using several microseconds long molecular dynamics (MD) simulations, we demonstrate that on binding with LDLR, there is a significant conformational change (population shift) in a distal loop (residues 211-222) region of PCSK9. Consistent with the bidirectional nature of allostery, we establish a clear correlation between the loop conformation and the binding affinity with LDLR. Using a thermodynamic argument, we establish that the loop conformations predominantly present in the apo state of PCSK9 would have lower LDLR binding affinity, and they would be potential targets for designing allosteric inhibitors. We elucidate the molecular origin of the allosteric coupling between this loop and the LDLR binding interface in terms of the population shift in a set of salt bridges and hydrogen bonds. Overall, our work provides a general strategy toward identifying allosteric hotspots: compare the conformational ensemble of the receptor between the apo and bound states of the protein and identify distal conformational changes, if any. The inhibitors should be designed to bind and stabilize the apo-specific conformations.

Structural basis of human NOX5 activation

NADPH oxidase 5 (NOX5) catalyzes the production of superoxide free radicals and regulates physiological processes from sperm motility to cardiac rhythm. Overexpression of NOX5 leads to cancers, diabetes, and cardiovascular diseases. NOX5 is activated by intracellular calcium signaling, but the underlying molecular mechanism of which - in particular, how calcium triggers electron transfer from NADPH to FAD - is still unclear. Here we capture motions of full-length human NOX5 upon calcium binding using single-particle cryogenic electron microscopy (cryo-EM). By combining biochemistry, mutagenesis analyses, and molecular dynamics (MD) simulations, we decode the molecular basis of NOX5 activation and electron transfer. We find that calcium binding to the EF-hand domain increases NADPH dynamics, permitting electron transfer between NADPH and FAD and superoxide production. Our structural findings also uncover a zinc-binding motif that is important for NOX5 stability and enzymatic activity, revealing modulation mechanisms of reactive oxygen species (ROS) production.

Combined effect of SAR-endolysin LysKpV475 with polymyxin B and Salmonella bacteriophage phSE-5

Endolysins are bacteriophage (or phage)-encoded enzymes that catalyse the peptidoglycan breakdown in the bacterial cell wall. The exogenous action of recombinant phage endolysins against Gram-positive organisms has been extensively studied. However, the outer membrane acts as a physical barrier when considering the use of recombinant endolysins to combat Gram-negative bacteria. This study aimed to evaluate the antimicrobial activity of the SAR-endolysin LysKpV475 against Gram-negative bacteria as single or combined therapies, using an outer membrane permeabilizer (polymyxin B) and a phage, free or immobilized in a pullulan matrix. In the first step, the endolysin LysKpV475 in solution, alone and combined with polymyxin B, was tested in vitro and in vivo against ten Gram-negative bacteria, including highly virulent strains and multidrug-resistant isolates. In the second step, the lyophilized LysKpV475 endolysin was combined with the phage phSE-5 and investigated, free or immobilized in a pullulan matrix, against Salmonella enterica subsp. enterica serovar Typhimurium ATCC 13311. The bacteriostatic action of purified LysKpV475 varied between 8.125 μg ml-1 against Pseudomonas aeruginosa ATCC 27853, 16.25 μg ml-1 against S. enterica Typhimurium ATCC 13311, and 32.50 μg ml-1 against Klebsiella pneumoniae ATCC BAA-2146 and Enterobacter cloacae P2224. LysKpV475 showed bactericidal activity only for P. aeruginosa ATCC 27853 (32.50 μg ml-1) and P. aeruginosa P2307 (65.00 μg ml-1) at the tested concentrations. The effect of the LysKpV475 combined with polymyxin B increased against K. pneumoniae ATCC BAA-2146 [fractional inhibitory concentration index (FICI) 0.34; a value lower than 1.0 indicates an additive/combined effect] and S. enterica Typhimurium ATCC 13311 (FICI 0.93). A synergistic effect against S. enterica Typhimurium was also observed when the lyophilized LysKpV475 at ⅔ MIC was combined with the phage phSE-5 (m.o.i. of 100). The lyophilized LysKpV475 immobilized in a pullulan matrix maintained a significant Salmonella reduction of 2 logs after 6 h of treatment. These results demonstrate the potential of SAR-endolysins, alone or in combination with other treatments, in the free form or immobilized in solid matrices, which paves the way for their application in different areas, such as in biocontrol at the food processing stage, biosanitation of food contact surfaces and biopreservation of processed food in active food packing.

Cyanobacterial compound Tolyporphine K as an inhibitor of Apo-PBP (penicillin-binding protein) in A. baumannii and its ADME assessment

Antibiotic-resistant Acinetobacter baumannii, is a common pathogen found in hospital settings and has become nosocomial due to its high infection-causing tendency amongst ICU patients. The present study explores the cyanocompoundswhich were capable to inhibit the Penicillin Binding Protein of A. baumannii through molecular docking, ADMET, and molecular dynamicssimulation strategy. A database having structural and origin details was generated for 85 bioactive compounds in MS Excel. The 3-D structures weredownloaded from the PubChem database and minimized. The receptor protein was minimized and validated for structure correctness. The database was screened against the penicillin-binding protein of A. baumannii through PyRx software. The top 5 compounds including the control molecule werefurther redocked to the receptor molecule through Autodock Vina software. The molecule pose having the highest affinity was further subjected to 100ns MD- simulation and simultaneously the in-vitro activity of the methanol extract and hexane extract was checked through agar well diffusion assay.Docking studies indicate Tolyporphine K to be a lead molecule which was further assessed through Molecular dynamics and MM/PBSA. The in-silicoresults suggested that the protein-ligand complex was found to be stable over the 100 ns trajectory with a binding free energy of -8.56 Kcalmol-1. Theligand did not induce any major structural conformation in the protein moiety and was largely stabilized by hydrophobic interactions. The bioactivityscore and ADME properties of the compounds were also calculated. The in-vitro agar well diffusion assay showed a moderate zone of inhibition of12.33mm. The results indicate that the compound Tolyporphin- K could be a potential inhibitor of penicillin-binding protein in A. baumannii. Yet furtherwork needs to be done to have a more concrete basis for the pathway of inhibition.Communicated by Ramaswamy H. Sarma.

Natural Polymer-Polyphenol Bioadhesive Coacervate with Stable Wet Adhesion, Antibacterial Activity, and On-Demand Detachment

Medical adhesives are emerging as an important clinical tool as adjuvants for sutures and staples in wound closure and healing and in the achievement of hemostasis. However, clinical adhesives combining cytocompatibility, as well as strong and stable adhesion in physiological conditions, are still in demand. Herein, a mussel-inspired strategy is explored to produce adhesive coacervates using tannic acid (TA) and methacrylate pullulan (PUL-MA). TA|PUL-MA coacervates mainly comprise van der Waals forces and hydrophobic interactions. The methacrylic groups in the PUL backbone increase the number of interactions in the adhesives matrix, resulting in enhanced cohesion and adhesion strength (72.7 Jm-2), compared to the non-methacrylated coacervate. The adhesive properties are kept in physiologic-mimetic solutions (72.8 Jm-2) for 72 h. The photopolymerization of TA|PUL-MA enables the on-demand detachment of the adhesive. The poor cytocompatibility associated with the use of phenolic groups is here circumvented by mixing reactive oxygen species-degrading enzyme in the adhesive coacervate. This addition does not hamper the adhesive character of the materials, nor their anti-microbial or hemostatic properties. This affordable and straightforward methodology, together with the tailorable adhesivity even in wet environments, high cytocompatibility, and anti-bacterial activity, enables foresee TA|PUL-MA as a promising ready-to-use bioadhesive for biomedical applications.

Coupling finite and boundary element methods to solve the Poisson-Boltzmann equation for electrostatics in molecular solvation

The Poisson-Boltzmann equation is widely used to model electrostatics in molecular systems. Available software packages solve it using finite difference, finite element, and boundary element methods, where the latter is attractive due to the accurate representation of the molecular surface and partial charges, and exact enforcement of the boundary conditions at infinity. However, the boundary element method is limited to linear equations and piecewise constant variations of the material properties. In this work, we present a scheme that couples finite and boundary elements for the linearised Poisson-Boltzmann equation, where the finite element method is applied in a confined solute region and the boundary element method in the external solvent region. As a proof-of-concept exercise, we use the simplest methods available: Johnson-Nédélec coupling with mass matrix and diagonal preconditioning, implemented using the Bempp-cl and FEniCSx libraries via their Python interfaces. We showcase our implementation by computing the polar component of the solvation free energy of a set of molecules using a constant and a Gaussian-varying permittivity. As validation, we compare against well-established finite difference solvers for an extensive binding energy data set, and with the finite difference code APBS (to 0.5%) for Gaussian permittivities. We also show scaling results from protein G B1 (955 atoms) up to immunoglobulin G (20,148 atoms). For small problems, the coupled method was efficient, outperforming a purely boundary integral approach. For Gaussian-varying permittivities, which are beyond the applicability of boundary elements alone, we were able to run medium to large-sized problems on a single workstation. The development of better preconditioning techniques and the use of distributed memory parallelism for larger systems remains an area for future work. We hope this work will serve as inspiration for future developments that consider space-varying field parameters, and mixed linear-nonlinear schemes for molecular electrostatics with implicit solvent models.

Polyphosphate and tyrosine phosphorylation in the N-terminal domain of the human mitochondrial Lon protease disrupts its functions

Phosphorylation plays a crucial role in the regulation of many fundamental cellular processes. Phosphorylation levels are increased in many cancer cells where they may promote changes in mitochondrial homeostasis. Proteomic studies on various types of cancer identified 17 phosphorylation sites within the human ATP-dependent protease Lon, which degrades misfolded, unassembled and oxidatively damaged proteins in mitochondria. Most of these sites were found in Lon's N-terminal (NTD) and ATPase domains, though little is known about the effects on their function. By combining the biochemical and cryo-electron microscopy studies, we show the effect of Tyr186 and Tyr394 phosphorylations in Lon's NTD, which greatly reduce all Lon activities without affecting its ability to bind substrates or perturbing its tertiary structure. A substantial reduction in Lon's activities is also observed in the presence of polyphosphate, whose amount significantly increases in cancer cells. Our study thus provides an insight into the possible fine-tuning of Lon activities in human diseases, which highlights Lon's importance in maintaining proteostasis in mitochondria.

Studying the ssDNA loaded adeno-associated virus aggregation using coarse-grained molecular dynamics simulations

The aggregation of adeno-associated viral (AAV) capsids in an aqueous environment was investigated via coarse-grained molecular dynamics (CG-MD) simulations. The primary driving force and mechanism of the aggregation were investigated with or without single-strand DNA (ssDNA) loaded at various process temperatures. Capsid aggregation appeared to involve multiple residue interactions (i.e., hydrophobic, polar and charged residues) leading to complex protein aggregation. In addition, two aggregation mechanisms (i.e., the fivefold face-to-face contact and the edge-to-edge contact) were identified from this study. The ssDNA with its asymmetric structure could be the reason for destabilizing protein subunits and enhancing the interaction between the charged residues, and further result in the non-reversible face-to-face contact. At higher temperature, the capsid structure was found to be unstable with the significant size expansion of the loaded ssDNA which could be attributed to reduced number of intramolecular hydrogen bonds, the increased conformational deviations of protein subunits and the higher residue fluctuations. The CG-MD model was further validated with previous experimental and simulation data, including the full capsid size measurement and the capsid internal pressure. Thus, a good understanding of AAV capsid aggregation, instability and the role of ssDNA were revealed by applying the developed computational model.

Extracellular Kir2.1C122Y Mutant Upsets Kir2.1-PIP2 Bonds and Is Arrhythmogenic in Andersen-Tawil Syndrome

BACKGROUND

Andersen-Tawil syndrome type 1 is a rare heritable disease caused by mutations in the gene coding the strong inwardly rectifying K+ channel Kir2.1. The extracellular Cys (cysteine)122-to-Cys154 disulfide bond in the channel structure is crucial for proper folding but has not been associated with correct channel function at the membrane. We evaluated whether a human mutation at the Cys122-to-Cys154 disulfide bridge leads to Kir2.1 channel dysfunction and arrhythmias by reorganizing the overall Kir2.1 channel structure and destabilizing its open state.

METHODS

We identified a Kir2.1 loss-of-function mutation (c.366 A>T; p.Cys122Tyr) in an ATS1 family. To investigate its pathophysiological implications, we generated an AAV9-mediated cardiac-specific mouse model expressing the Kir2.1C122Y variant. We employed a multidisciplinary approach, integrating patch clamping and intracardiac stimulation, molecular biology techniques, molecular dynamics, and bioluminescence resonance energy transfer experiments.

RESULTS

Kir2.1C122Y mice recapitulated the ECG features of ATS1 independently of sex, including corrected QT prolongation, conduction defects, and increased arrhythmia susceptibility. Isolated Kir2.1C122Y cardiomyocytes showed significantly reduced inwardly rectifier K+ (IK1) and inward Na+ (INa) current densities independently of normal trafficking. Molecular dynamics predicted that the C122Y mutation provoked a conformational change over the 2000-ns simulation, characterized by a greater loss of hydrogen bonds between Kir2.1 and phosphatidylinositol 4,5-bisphosphate than wild type (WT). Therefore, the phosphatidylinositol 4,5-bisphosphate-binding pocket was destabilized, resulting in a lower conductance state compared with WT. Accordingly, on inside-out patch clamping, the C122Y mutation significantly blunted Kir2.1 sensitivity to increasing phosphatidylinositol 4,5-bisphosphate concentrations. In addition, the Kir2.1C122Y mutation resulted in channelosome degradation, demonstrating temporal instability of both Kir2.1 and NaV1.5 proteins.

CONCLUSIONS

The extracellular Cys122-to-Cys154 disulfide bond in the tridimensional Kir2.1 channel structure is essential for the channel function. We demonstrate that breaking disulfide bonds in the extracellular domain disrupts phosphatidylinositol 4,5-bisphosphate-dependent regulation, leading to channel dysfunction and defects in Kir2.1 energetic stability. The mutation also alters functional expression of the NaV1.5 channel and ultimately leads to conduction disturbances and life-threatening arrhythmia characteristic of Andersen-Tawil syndrome type 1.

Biomimetic MDSCs membrane coated black phosphorus nanosheets system for photothermal therapy/photodynamic therapy synergized chemotherapy of cancer

Photothermal therapy is favored by cancer researchers due to its advantages such as controllable initiation, direct killing and immune promotion. However, the low enrichment efficiency of photosensitizer in tumor site and the limited effect of single use limits the further development of photothermal therapy. Herein, a photo-responsive multifunctional nanosystem was designed for cancer therapy, in which myeloid-derived suppressor cell (MDSC) membrane vesicle encapsulated decitabine-loaded black phosphorous (BP) nanosheets (BP@ Decitabine @MDSCs, named BDM). The BDM demonstrated excellent biosafety and biochemical characteristics, providing a suitable microenvironment for cancer cell killing. First, the BDM achieves the ability to be highly enriched at tumor sites by inheriting the ability of MDSCs to actively target tumor microenvironment. And then, BP nanosheets achieves hyperthermia and induces mitochondrial damage by its photothermal and photodynamic properties, which enhancing anti-tumor immunity mediated by immunogenic cell death (ICD). Meanwhile, intra-tumoral release of decitabine induced G2/M cell cycle arrest, further promoting tumor cell apoptosis. In vivo, the BMD showed significant inhibition of tumor growth with down-regulation of PCNA expression and increased expression of high mobility group B1 (HMGB1), calreticulin (CRT) and caspase 3. Flow cytometry revealed significantly decreased infiltration of MDSCs and M2-macrophages along with an increased proportion of CD4+, CD8+ T cells as well as CD103+ DCs, suggesting a potentiated anti-tumor immune response. In summary, BDM realizes photothermal therapy/photodynamic therapy synergized chemotherapy for cancer.

Vitamin K Epoxide Reductase Complex-Protein Disulphide Isomerase Assemblies in the Thiol-Disulphide Exchange Reactions: Portrayal of Precursor-to-Successor Complexes

The human Vitamin K Epoxide Reductase Complex (hVKORC1), a key enzyme that converts vitamin K into the form necessary for blood clotting, requires for its activation the reducing equivalents supplied by its redox partner through thiol-disulphide exchange reactions. The functionally related molecular complexes assembled during this process have never been described, except for a proposed de novo model of a 'precursor' complex of hVKORC1 associated with protein disulphide isomerase (PDI). Using numerical approaches (in silico modelling and molecular dynamics simulation), we generated alternative 3D models for each molecular complex bonded either covalently or non-covalently. These models differ in the orientation of the PDI relative to hVKORC1 and in the cysteine residue involved in forming protein-protein disulphide bonds. Based on a comparative analysis of these models' shape, folding, and conformational dynamics, the most probable putative complexes, mimicking the 'precursor', 'intermediate', and 'successor' states, were suggested. In addition, we propose using these complexes to develop the 'allo-network drugs' necessary for treating blood diseases.

Synthesis, topology, molecular docking and dynamics studies of o-phenylenediamine derivative

The N, N'-(1,2-phenylene) bis (1- (4- chlorophenyl) methanimine) (CS4) was synthesized and characterized by infrared (IR), absorption (UV-vis) and NMR (1H and 13C) spectral analyses. The structural parameters, vibrational frequencies, potential energy and the distribution analysis (PED) were calculated by using DFT with the basis set of B3LYP/cc-pVDZ and these spectral values were compared to the experimental values. HOMO and LUMO studied were performed in order to understand the stability and biological activity of the compound. The most reactive sites on the compound were investigated by utilizing MEP energy surface and Fukui function descriptor with the natural population analysis (NPA) of the charges. The study of the natural bond orbitals (NBO) reveals the delocalization of the intramolecular interaction of the charges in the compound. Additionally, topological investigations (ELF, LOL), determination of thermodynamic parameters and noncovalent interaction (NCI) study by using topology (RDG) analysis were also carried out. Finally, the molecular docking and molecular dynamics simulations was carried out by examining against glycosylphosphatidylinositol phospholipase D inhibitor receptor for distinct protein targets (3MZG).

Antifungal plant flavonoids identified in silico with potential to control rice blast disease caused by Magnaporthe oryzae

Rice blast disease, caused by the fungus Magnaporthe oryzae, poses a severe threat to rice production, particularly in Asia where rice is a staple food. Concerns over fungicide resistance and environmental impact have sparked interest in exploring natural fungicides as potential alternatives. This study aimed to identify highly potent natural fungicides against M. oryzae to combat rice blast disease, using advanced molecular dynamics techniques. Four key proteins (CATALASE PEROXIDASES 2, HYBRID PKS-NRPS SYNTHETASE TAS1, MANGANESE LIPOXYGENASE, and PRE-MRNA-SPLICING FACTOR CEF1) involved in M. oryzae's infection process were identified. A list of 30 plant metabolites with documented antifungal properties was compiled for evaluation as potential fungicides. Molecular docking studies revealed that 2-Coumaroylquinic acid, Myricetin, Rosmarinic Acid, and Quercetin exhibited superior binding affinities compared to reference fungicides (Azoxystrobin and Tricyclazole). High throughput molecular dynamics simulations were performed, analyzing parameters like RMSD, RMSF, Rg, SASA, hydrogen bonds, contact analysis, Gibbs free energy, and cluster analysis. The results revealed stable interactions between the selected metabolites and the target proteins, involving important hydrogen bonds and contacts. The SwissADME server analysis indicated that the metabolites possess fungicide properties, making them effective and safe fungicides with low toxicity to the environment and living beings. Additionally, bioactivity assays confirmed their biological activity as nuclear receptor ligands and enzyme inhibitors. Overall, this study offers valuable insights into potential natural fungicides for combating rice blast disease, with 2-Coumaroylquinic acid, Myricetin, Rosmarinic Acid, and Quercetin standing out as promising and environmentally friendly alternatives to conventional fungicides. These findings have significant implications for developing crop protection strategies and enhancing global food security, particularly in rice-dependent regions.

Functional convergence in gastric lysozymes of foregut-fermenting rodents, ruminants, and primates is not attributed to convergent molecular evolution

Convergent evolution is a widespread phenomenon. While there are many examples of convergent evolution at the phenotypic scale, convergence at the molecular level has been more difficult to identify. A classic example of convergent evolution across scales is that of the digestive lysozyme found in ruminants and Colobine monkeys. These herbivorous species rely on foregut fermentation, which has evolved to function more optimally under acidic conditions. Here, we explored if rodents with similar dietary strategies and digestive morphologies have convergently evolved a lysozyme with digestive functions. At the phenotypic level, we find that rodents with bilocular stomach morphologies exhibited a lysozyme that maintained higher relative activities at low pH values, similar to the lysozymes of ruminants and Colobine monkeys. Additionally, the lysozyme of Peromyscus leucopus shared a similar predicted protonation state as that observed in previously identified digestive lysozymes. However, we found limited evidence of positive selection acting on the lysozyme gene in foregut-fermenting species and did not identify patterns of convergent molecular evolution in this gene. This study emphasizes that phenotypic convergence need not be the result of convergent genetic modifications, and we encourage further exploration into the mechanisms regulating convergence across biological scales.

A combined in vitro and molecular dynamics simulation studies unveil the molecular basis of the anticancer potential of piperine targeting AKT1 against prostate cancer

The present study investigates the activity of the natural compound piperine on prostate cancer cell line (PC-3), followed by exploring its mechanistic inhibition on the RAC-alpha serine/threonine-protein kinase (AKT1) protein. The 3-(4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazolium Bromide (MTT) assay showed that after 24 hrs of exposure to piperine (15 µmol/ml), cell viability fell to 50% compared to the standard drug flutamide (SDF) (51 µmol/ml) with a lower IC50 concentration. However, the Dual acridine orange/ethidium bromide (AO/EtBr) staining demonstrated that, as compared to the SDF, piperine caused substantial cellular death in PC-3 cells, presumably by triggering DNA fragmentation. In addition, compared to untreated cells, the proportion of the sub-G0/G1 and G2/M stages population increased considerably in piperine-treated cells. The cell cycle's sub-G0/G1 and G2/M phases were also arrested in piperine-treated cells compared to the SDF in cell cycle analysis. Based on our systems pharmacology and molecular docking studies, AKT1 is predicted as a potential target against piperine. The complementary charge between AKT1 and piperine was emphasized in the transient ligand-protein binding interaction in molecular dynamic modeling over 100 ns, and stable hydrogen bond interaction between Lys268 and Ser205 amino acid residues of the active pocket was hypothesized. Overall, the findings from our in vitro and MD simulations provide insights into the mechanism of piperine targeting AKT1 and offer a possible candidate for future prostate cancer therapeutic development.Communicated by Ramaswamy H. Sarma.

High-throughput screening and molecular dynamics simulations of natural products targeting LuxS/AI-2 system as a novel antibacterial strategy for antibiotic resistance in Helicobacter pylori

The main goal of treating any Helicobacter pylori (H. pylori)-related gastrointestinal disease is completely eradicating infection. Falling eradication efficiency, off-target effects, and patient noncompliance with prolonged and broad spectrums have sparked clinical interest in exploring other effective, safer therapeutic choices. As natural substances are risk-free and privileged with high levels of antibacterial activity, most of these natural chemical's specific modes of action are unknown. With the aid of in silico molecular docking-based virtual screening studies and molecular dynamic simulations, the current study is intended to gather data on numerous such natural chemicals and assess their affinity for the S-ribosyl homocysteine lyase (LuxS) protein of H. pylori. The ligand with the highest binding energy with LuxS, glucoraphanin, catechin gallate and epigallocatechin gallate were rationally selected for further computational analysis. The solution stability of the three compounds' optimal docking postures with LuxS was initially assessed using long-run molecular dynamics simulations. Using molecular dynamics simulation, the epigallocatechin gallate was found to be the most stable molecule with the highest binding free energy, indicating that it might compete with the natural ligand of the inhibitors. According to ADMET analysis, his phytochemical was a promising therapeutic candidate for an antibacterial action since it had a range of physicochemical, pharmacokinetic, and drug-like qualities and had no discernible adverse effects. Additional in vitro, in vivo, and clinical trials are needed to confirm the drug's precise efficacy during H. pylori infection.Communicated by Ramaswamy H. Sarma.

Engineering Human Pancreatic RNase 1 as an Immunotherapeutic Agent for Cancer Therapy Through Computational and Experimental Studies

Most plant and bacterial toxins are highly immunogenic with non-specific toxic effects. Human ribonucleases are thought to provide a promising basis for reducing the toxic agent's immunogenic properties, which are candidates for cancer therapy. In the cell, the ribonuclease inhibitor (RI) protein binds to the ribonuclease enzyme and forms a tight complex. This study aimed to engineer and provide a gene construct encoding an improved version of Human Pancreatic RNase 1 (HP-RNase 1) to reduce connection to RI and modulate the immunogenic effects of immunotoxins. To further characterize the interaction complex of HP-RNase 1 and RI, we established various in silico and in vitro approaches. These methods allowed us to specifically monitor interactions within native and engineered HP-RNase 1/RI complexes. In silico research involved molecular dynamics (MD) simulations of native and mutant HP-RNase 1 in their free form and when bound to RI. For HP-RNase 1 engineering, we designed five mutations (K8A/N72A/N89A/R92D/E112/A) based on literature studies, as this combination proved effective for the intended investigation. Then, the cDNA encoding HP-RNase 1 was generated by RT-PCR from blood and cloned into the pSYN2 expression vector. Consequently, wild-type and the engineered HP-RNase 1 were over-expressed in E. coli TG1 and purified using an IMAC column directed against a poly-his tag. The protein products were detected by SDS-PAGE and Western blot analysis. HP-RNase 1 catalytic activity, in the presence of various concentrations of RI, demonstrated that the mutated version of the protein is able to escape the ribonuclease inhibitor and target the RNA substrate 2.5 folds more than that of the wild type. From these data, we tend to suggest the engineered recombinant HP-RNase 1 potentially as a new immunotherapeutic agent for application in human cancer therapy.

Molecular mechanism of cellulose depolymerization by the two-domain BlCel9A enzyme from the glycoside hydrolase family 9

Carbohydrate-active enzymes from the glycoside hydrolase family 9 (GH9) play a key role in processing lignocellulosic biomass. Although the structural features of some GH9 enzymes are known, the molecular mechanisms that drive their interactions with cellulosic substrates remain unclear. To investigate the molecular mechanisms that the two-domain Bacillus licheniformis BlCel9A enzyme utilizes to depolymerize cellulosic substrates, we used a combination of biochemical assays, X-ray crystallography, small-angle X-ray scattering, and molecular dynamics simulations. The results reveal that BlCel9A breaks down cellulosic substrates, releasing cellobiose and glucose as the major products, but is highly inefficient in cleaving oligosaccharides shorter than cellotetraose. In addition, fungal lytic polysaccharide oxygenase (LPMO) TtLPMO9H enhances depolymerization of crystalline cellulose by BlCel9A, while exhibiting minimal impact on amorphous cellulose. The crystal structures of BlCel9A in both apo form and bound to cellotriose and cellohexaose were elucidated, unveiling the interactions of BlCel9A with the ligands and their contribution to substrate binding and products release. MD simulation analysis reveals that BlCel9A exhibits higher interdomain flexibility under acidic conditions, and SAXS experiments indicate that the enzyme flexibility is induced by pH and/or temperature. Our findings provide new insights into BlCel9A substrate specificity and binding, and synergy with the LPMOs.

Discovery of novel BRD4-BD2 inhibitors via in silico approaches: QSAR techniques, molecular docking, and molecular dynamics simulations

Bromodomain-containing protein 4(BRD4) plays an important role in the occurrence and development of various malignant tumors, which has attracted the attention of scientific research institutions and pharmaceutical companies. The structural modification of most currently available BRD4 inhibitors is relatively simple, but the drug effectiveness is limited. Research has found that the inhibition of BD1 may promote the differentiation of oligodendrocyte progenitor cell; however, the inhibition of BD2 will not cause this outcome. Therefore, newly potential drugs which target BRD4-BD2 need further research. Herein, we initially built QSAR models out of 49 compounds using HQSAR, CoMFA, CoMSIA, and Topomer CoMFA technology. All of the models have shown suitable reliabilities (q2 = 0.778, 0.533, 0.640, 0.702, respectively) and predictive abilities (r2pred = 0.716, 0.6289, 0.6153, 0.7968, respectively) for BRD4-BD2 inhibitors. On the basis of QSAR results and the search of the R-group in the topomer search module, we designed 20 new compounds with high activity that showed appropriate docking score and suitable ADMET. Docking studies and MD simulation were carried out to reveal the amino acid residues (Asn351, Cys347, Tyr350, Pro293, and Asp299) at the active site of BRD4-BD2. Free energy calculations and free energy landscapes verified the stable binding results and indicated stable conformations of the complexes. These theoretical studies provide guidance and theoretical basis for designing and developing novel BRD4-BD2 inhibitors.

Geometry-Based Molecular Generation With Deep Constrained Variational Autoencoder

Finding target molecules with specific chemical properties plays a decisive role in drug development. We proposed GEOM-CVAE, a constrained variational autoencoder based on geometric representation for molecular generation with specific properties, which is protein-context-dependent. In terms of machine learning, it includes continuous feature embedding encoder and molecular generation decoder. Our key contribution is to propose an efficient geometric embedding method, including the spatial structure representations of drug molecule (converting the 3-D coordinates into image) and the geometric graph representations of protein target (modeling the protein surface as a mesh). The 3-D geometric information is vital to successful molecular generation, which is different from previous molecular generative methods based on 1-D or 2-D. Our model framework generates specific molecules in two phases, by first generating special image with molecular 3-D information to learn latent representations and generating molecules with constrained condition based on geometric graph convolution for specific protein and then inputting the generated structural molecules into a parser network for obtaining Simplified Molecular Input Line Entry System (SMILES) strings. Our model achieves competitive performance that implies its potential effectiveness to enable the exploration of the vast chemical space for drug discovery.

Targeting GRP75 with a Chlorpromazine Derivative Inhibits Endometrial Cancer Progression Through GRP75-IP3R-Ca2+-AMPK Axis

Tumors often overexpress glucose-regulated proteins, and agents that interfere with the production or activity of these proteins may represent novel cancer treatments. The chlorpromazine derivative JX57 exhibits promising effects against endometrial cancer with minimal extrapyramidal side effects; however, its mechanisms of action are currently unknown. Here, glucose-regulated protein 75 kD (GRP75) is identified as a direct target of JX57 using activity-based protein profiling and loss-of-function experiments. The findings show that GRP75 is necessary for the biological activity of JX57, as JX57 exhibits moderate anticancer properties in GRP75-deficient cancer cells, both in vitro and in vivo. High GRP75 expression is correlated with poor differentiation and poor survival in patients with endometrial cancer, whereas the knockdown of GRP75 can significantly suppress tumor growth. Mechanistically, the direct binding of JX57 to GRP75 impairs the structure of the mitochondria-associated endoplasmic reticulum membrane and disrupts the endoplasmic reticulum-mitochondrial calcium homeostasis, resulting in a mitochondrial energy crisis and AMP-activated protein kinase activation. Taken together, these findings highlight GRP75 as a potential prognostic biomarker and direct therapeutic target in endometrial cancer and suggest that the chlorpromazine derivative JX57 can potentially be a new therapeutic option for endometrial cancer.

Antiviral peptides inhibiting the main protease of SARS-CoV-2 investigated by computational screening and in vitro protease assay

The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plays an important role in viral replication and transcription and received great attention as a vital target for drug/peptide development. Therapeutic agents such as small-molecule drugs or peptides that interact with the Cys-His present in the catalytic site of Mpro are an efficient way to inhibit the protease. Although several emergency-approved vaccines showed good efficacy and drastically dropped the infection rate, evolving variants are still infecting and killing millions of people globally. While a small-molecule drug (Paxlovid) received emergency approval, small-molecule drugs have low target specificity and higher toxicity. Besides small-molecule drugs, peptide therapeutics are thus gaining increasing popularity as they are easy to synthesize and highly selective and have limited side effects. In this study, we investigated the therapeutic value of 67 peptides targeting Mpro using molecular docking. Subsequently, molecular dynamics (MD) simulations were implemented on eight protein-peptide complexes to obtain molecular-level information on the interaction between these peptides and the Mpro active site, which revealed that temporin L, indolicidin, and lymphocytic choriomeningitis virus (LCMV) GP1 are the best candidates in terms of stability, interaction, and structural compactness. These peptides were synthesized using the solid-phase peptide synthesis protocol, purified by reversed-phase high-performance liquid chromatography (RP-HPLC), and authenticated by mass spectrometry (MS). The in vitro fluorometric Mpro activity assay was used to validate the computational results, where temporin L and indolicidin were observed to be very active against SARS-CoV-2 Mpro with IC50 values of 38.80 and 87.23 μM, respectively. A liquid chromatography-MS (LC-MS) assay was developed, and the IC50 value of temporin L was measured at 23.8 μM. The solution-state nuclear magnetic resonance (NMR) structure of temporin L was determined in the absence of sodium dodecyl sulfate (SDS) micelles and was compared to previous temporin structures. This combined investigation provides critical insights and assists us to further develop peptide inhibitors of SARS-CoV-2 Mpro through structural guided investigation.

Structure and function of Semaphorin-5A glycosaminoglycan interactions

Integration of extracellular signals by neurons is pivotal for brain development, plasticity, and repair. Axon guidance relies on receptor-ligand interactions crosstalking with extracellular matrix components. Semaphorin-5A (Sema5A) is a bifunctional guidance cue exerting attractive and inhibitory effects on neuronal growth through the interaction with heparan sulfate (HS) and chondroitin sulfate (CS) glycosaminoglycans (GAGs), respectively. Sema5A harbors seven thrombospondin type-1 repeats (TSR1-7) important for GAG binding, however the underlying molecular basis and functions in vivo remain enigmatic. Here we dissect the structural basis for Sema5A:GAG specificity and demonstrate the functional significance of this interaction in vivo. Using x-ray crystallography, we reveal a dimeric fold variation for TSR4 that accommodates GAG interactions. TSR4 co-crystal structures identify binding residues validated by site-directed mutagenesis. In vitro and cell-based assays uncover specific GAG epitopes necessary for TSR association. We demonstrate that HS-GAG binding is preferred over CS-GAG and mediates Sema5A oligomerization. In vivo, Sema5A:GAG interactions are necessary for Sema5A function and regulate Plexin-A2 dependent dentate progenitor cell migration. Our study rationalizes Sema5A associated developmental and neurological disorders and provides mechanistic insights into how multifaceted guidance functions of a single transmembrane cue are regulated by proteoglycans.

Deciphering the mannose transfer mechanism of mycobacterial PimE by molecular dynamics simulations

Phosphatidyl-myo-inositol mannosides (PIMs), Lipomannan (LM), and Lipoarabinomannan (LAM) are essential components of the cell envelopes of mycobacteria. At the beginning of the biosynthesis of these compounds, phosphatidylinositol (PI) is mannosylated and acylated by various enzymes to produce Ac1/2PIM4, which is used to synthesize either Ac1/2PIM6 or LM/LAM. The protein PimE, a membrane-bound glycosyltransferase (GT-C), catalyzes the addition of a mannose group to Ac1PIM4 to produce Ac1PIM5, using polyprenolphosphate mannose (PPM) as the mannose donor. PimE-deleted Mycobacterium smegmatis (Msmeg) showed structural deformity and increased antibiotic and copper sensitivity. Despite knowing that the mutation D58A caused inactivity in Msmeg, how PimE catalyzes the transfer of mannose from PPM to Ac1/2PIM4 remains unknown. In this study, analyzing the AlphaFold structure of PimE revealed the presence of a tunnel through the D58 residue with two differently charged gates. Molecular docking suggested PPM binds to the hydrophobic tunnel gate, whereas Ac1PIM4 binds to the positively charged tunnel gate. Molecular dynamics (MD) simulations further demonstrated the critical roles of the residues N55, F87, L89, Y163, Q165, K197, L198, R251, F277, W324, H326, and I375 in binding PPM and Ac1PIM4. The mutation D58A caused a faster release of PPM from the catalytic tunnel, explaining the loss of PimE activity. Along with a hypothetical mechanism of mannose transfer by PimE, we also observe the presence of tunnels through a negatively charged aspartate or glutamate with two differently-charged gates among most GT-C enzymes. Common hydrophobic gates of GT-C enzymes probably harbor sugar donors, whereas, differently-charged tunnel gates accommodate various sugar-acceptors.

Diverging co-translational protein complex assembly pathways are governed by interface energy distribution

Protein-protein interactions are at the heart of all cellular processes, with the ribosome emerging as a platform, orchestrating the nascent-chain interplay dynamics. Here, to study the characteristics governing co-translational protein folding and complex assembly, we combine selective ribosome profiling, imaging, and N-terminomics with all-atoms molecular dynamics. Focusing on conserved N-terminal acetyltransferases (NATs), we uncover diverging co-translational assembly pathways, where highly homologous subunits serve opposite functions. We find that only a few residues serve as "hotspots," initiating co-translational assembly interactions upon exposure at the ribosome exit tunnel. These hotspots are characterized by high binding energy, anchoring the entire interface assembly. Alpha-helices harboring hotspots are highly thermolabile, folding and unfolding during simulations, depending on their partner subunit to avoid misfolding. In vivo hotspot mutations disrupted co-translational complexation, leading to aggregation. Accordingly, conservation analysis reveals that missense NATs variants, causing neurodevelopmental and neurodegenerative diseases, disrupt putative hotspot clusters. Expanding our study to include phosphofructokinase, anthranilate synthase, and nucleoporin subcomplex, we employ AlphaFold-Multimer to model the complexes' complete structures. Computing MD-derived interface energy profiles, we find similar trends. Here, we propose a model based on the distribution of interface energy as a strong predictor of co-translational assembly.

The origin of different bending stiffness between double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The subtle differences in the chemical structures of double-stranded (ds) RNA and DNA lead to significant variations in their biological roles and medical implications, largely due to their distinct biophysical properties, such as bending stiffness. Although it is well known that A-form dsRNA is stiffer than B-form dsDNA under physiological salt conditions, the underlying cause of this difference remains unclear. In this study, we employ high-precision magnetic-tweezer experiments along with molecular dynamics simulations and reveal that the relative bending stiffness between dsRNA and dsDNA is primarily determined by the structure- and salt-concentration-dependent ion distribution around their helical structures. At near-physiological salt conditions, dsRNA shows a sparser ion distribution surrounding its phosphate groups compared to dsDNA, causing its greater stiffness. However, at very high monovalent salt concentrations, phosphate groups in both dsRNA and dsDNA become fully neutralized by excess ions, resulting in a similar intrinsic bending persistence length of approximately 39 nm. This similarity in intrinsic bending stiffness of dsRNA and dsDNA is coupled to the analogous fluctuations in their total groove widths and further coupled to the similar fluctuation of base-pair inclination, despite their distinct A-form and B-form helical structures.

Dissecting the chiral recognition of TLR4/MD2 with Neoseptin-3 enantiomers by molecular dynamics simulations

Toll-like receptor 4 (TLR4) is a pivotal innate immune recognition receptor that regulates intricate signaling pathways within the immune system. Neoseptin-3 (Neo-3), a recently identified small-molecule agonist for mouse TLR4/MD2, exhibits chiral recognition properties. Specifically, the L-enantiomer of Neo-3 (L-Neo-3) effectively activates the TLR4 signaling pathway, while D-Neo-3 fails to induce TLR4 activation. However, the underlying mechanism by which TLR4 enantioselectively recognizes Neo-3 enantiomers remains poorly understood. In this study, in silico simulations were performed to investigate the mechanism of chiral recognition of Neo-3 enantiomers by TLR4/MD2. Two L-Neo-3 molecules stably resided within the cavity of MD2 as a dimer, and the L-Neo-3 binding stabilized the (TLR4/MD2)2 dimerization state. However, the strong electrostatic repulsion between the hydrogen atoms on the chiral carbon of D-Neo-3 molecules caused the relative positions of two D-Neo-3 molecules to continuously shift during the simulation process, thus preventing the formation of D-Neo-3 dimer as well as their stable interactions with the surrounding residues in (TLR4/MD2)2. Considering that L-Neo-3 could not sustain a stable dimeric state in the bulk aqueous environment, it is unlikely that L-Neo-3 entered the cavity of MD2 as a dimeric unit. Umbrella sampling simulations revealed that the second L-Neo-3 molecule entering the cavity of MD2 exhibited a lower binding energy (-25.75 kcal mol-1) than that of the first L-Neo-3 molecule (-14.31 kcal mol-1). These results imply that two L-Neo-3 molecules enter the cavity of MD2 sequentially, with the binding of the first L-Neo-3 molecule facilitating the entry of the second one. This study dissects the binding process of Neo-3 enantiomers, offering a comprehensive understanding of the atomic-level mechanism underlying TLR4's chiral recognition of Neo-3 molecules.

AlphaFold2-Enabled Atomistic Modeling of Structure, Conformational Ensembles, and Binding Energetics of the SARS-CoV-2 Omicron BA.2.86 Spike Protein with ACE2 Host Receptor and Antibodies: Compensatory Functional Effects of Binding Hotspots in Modulating Mechanisms of Receptor Binding and Immune Escape

The latest wave of SARS-CoV-2 Omicron variants displayed a growth advantage and increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with atomistic simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both the structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that the AlphaFold2-predicted structural ensemble of the BA.2.86 spike protein complex with ACE2 can accurately capture the main conformational states of the Omicron variant. Complementary to AlphaFold2 structural predictions, microsecond molecular dynamics simulations reveal the details of the conformational landscape and produced equilibrium ensembles of the BA.2.86 structures that are used to perform mutational scanning of spike residues and characterize structural stability and binding energy hotspots. The ensemble-based mutational profiling of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 revealed a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 convergent mutational hotspots R403K, F486P, and R493Q. To examine the immune evasion properties of BA.2.86 in atomistic detail, we performed structure-based mutational profiling of the spike protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against the BA.2.86 variant. The results revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have evolved to outcompete other Omicron subvariants by improving immune evasion while preserving binding affinity with ACE2 via through a compensatory effect of R493Q and F486P convergent mutational hotspots. This study demonstrated that an integrative approach combining AlphaFold2 predictions with complementary atomistic molecular dynamics simulations and robust ensemble-based mutational profiling of spike residues can enable accurate and comprehensive characterization of structure, dynamics, and binding mechanisms of newly emerging Omicron variants.

Defining the EM-signature of successful cell-transfection

In this report, we describe the architecture of Lipofectamine 2000 and 3000 transfection- reagents, as they appear inside of transfected cells, using classical transmission electron microscopy (EM). We also demonstrate that they provoke consistent structural changes after they have entered cells, changes that not only provide new insights into the mechanism of action of these particular transfection-reagents, but also provide a convenient and robust method for identifying by EM which cells in any culture have been successfully transfected. This also provides clues to the mechanism(s) of their toxic effects, when they are applied in excess. We demonstrate that after being bulk-endocytosed by cells, the cationic spheroids of Lipofectamine remain intact throughout the entire time of culturing, but escape from their endosomes and penetrate directly into the cytoplasm of the cell. In so doing, they provoke a stereotypical recruitment and rearrangement of endoplasmic reticulum (ER), and they ultimately end up escaping into the cytoplasm and forming unique 'inclusion-bodies.' Once free in the cytoplasm, they also invariably develop dense and uniform coatings of cytoplasmic ribosomes on their surfaces, and finally, they become surrounded by 'annulate' lamellae' of the ER. In the end, these annulate-lamellar enclosures become the ultrastructural 'signatures' of these inclusion-bodies, and serve to positively and definitively identify all cells that have been effectively transfected. Importantly, these new EM-observations define several new and unique properties of these classical Lipofectamines, and allow them to be discriminated from other lipoidal or particulate transfection-reagents, which we find do not physically break out of endosomes or end up in inclusion bodies, and in fact, provoke absolutely none of these 'signature' cytoplasmic reactions.

Coherent light scattering from cellular dynamics in living tissues

This review examines the biological physics of intracellular transport probed by the coherent optics of dynamic light scattering from optically thick living tissues. Cells and their constituents are in constant motion, composed of a broad range of speeds spanning many orders of magnitude that reflect the wide array of functions and mechanisms that maintain cellular health. From the organelle scale of tens of nanometers and upward in size, the motion inside living tissue is actively driven rather than thermal, propelled by the hydrolysis of bioenergetic molecules and the forces of molecular motors. Active transport can mimic the random walks of thermal Brownian motion, but mean-squared displacements are far from thermal equilibrium and can display anomalous diffusion through Lévy or fractional Brownian walks. Despite the average isotropic three-dimensional environment of cells and tissues, active cellular or intracellular transport of single light-scattering objects is often pseudo-one-dimensional, for instance as organelle displacement persists along cytoskeletal tracks or as membranes displace along the normal to cell surfaces, albeit isotropically oriented in three dimensions. Coherent light scattering is a natural tool to characterize such tissue dynamics because persistent directed transport induces Doppler shifts in the scattered light. The many frequency-shifted partial waves from the complex and dynamic media interfere to produce dynamic speckle that reveals tissue-scale processes through speckle contrast imaging and fluctuation spectroscopy. Low-coherence interferometry, dynamic optical coherence tomography, diffusing-wave spectroscopy, diffuse-correlation spectroscopy, differential dynamic microscopy and digital holography offer coherent detection methods that shed light on intracellular processes. In health-care applications, altered states of cellular health and disease display altered cellular motions that imprint on the statistical fluctuations of the scattered light. For instance, the efficacy of medical therapeutics can be monitored by measuring the changes they induce in the Doppler spectra of livingex vivocancer biopsies.

Developmental and epileptic encephalopathy 82 (DEE82) with novel compound heterozygous mutations of GOT2 gene

PURPOSE

Developmental and Epileptic Encephalopathies (DEEs) are rare neurological disorders characterized by early-onset medically resistant epileptic seizures, structural brain malformations, and severe developmental delays. These disorders can arise from mutations in genes involved in vital metabolic pathways, including those within the brain. Recent studies have implicated defects in the mitochondrial malate aspartate shuttle (MAS) as potential contributors to the clinical manifestation of infantile epileptic encephalopathy. Although rare, mutations in MDH1, MDH2, AGC1, or GOT2 genes have been reported in patients exhibiting neurological symptoms such as global developmental delay, epilepsy, and progressive microcephaly.

METHOD

In this study, we employed exome data analysis of a patient diagnosed with DEE, focusing on the screening of 1896 epilepsy-related genes listed in the HPO and ClinVar databases. Sanger sequencing was subsequently conducted to validate and assess the inheritance pattern of the identified variants within the family. The evolutionary conservation scores of the mutated residues were evaluated using the ConSurf Database. Furthermore, the impacts of the causative variations on protein stability were analyzed through I-Mutant and MuPro bioinformatic tools. Structural comparisons between wild-type and mutant proteins were performed using PyMOL, and the physicochemical effects of the mutations were assessed using Project Hope.

RESULTS

Exome data analysis unveiled the presence of novel compound heterozygous mutations in the GOT2 gene coding for mitochondrial glutamate aspartate transaminase. Sanger sequencing confirmed the paternal inheritance of the p.Asp257Asn mutation and the maternal inheritance of the p.Arg262Cys mutation. The affected individual exhibited plasma metabolic disturbances, including hyperhomocysteinemia, hyperlactatemia, and reduced levels of methionine and arginine. Detailed bioinformatic analysis indicated that the mutations were located within evolutionarily conserved domains of the enzyme, resulting in disruptions to protein stability and structure.

CONCLUSION

Herein, we describe a case with DEE82 (MIM: # 618721) with pathologic novel biallelic mutations in the GOT2 gene. Early genetic diagnosis of metabolic epilepsies is crucial for long-term neurodevelopmental improvements and seizure control as targeted treatments can be administered based on the affected metabolic pathways.

A novel cellulolytic/xylanolytic SbAA14 from Sordaria brevicollis with a branched chain preference and its synergistic effects with glycoside hydrolases on lignocellulose

In this study, the novel auxiliary activity (AA) family 14 lytic polysaccharide monooxygenase (LPMO) SbAA14 from Sordaria brevicollis was successfully characterized. It was active against heteroxylan, xyloglucan and cellulose in β-cellulose and released native oligosaccharides and corresponding C1- and/or C4-oxidized products. SbAA14 showed a branched chain preference, because partial removal of arabinosyl substituents from heteroxylan led to a decrease in activity. SbAA14 had synergistic effects with the debranching enzyme EpABF62C in an enzyme- and ascorbic acid-dependent manner. SbAA14 had synergistic effects with the GH10 endoxylanase EpXYN1, and the degree of synergy was greater with step-by-step addition than with simultaneous addition. SbAA14 could also synergize with Celluclast® 1.5 L on NaOH-pretreated wheat straw and on NaOH-pretreated and hydrogen peroxide-acetic acid (HPAC)-H2SO4-pretreated bamboo substrates. The greatest synergistic effect between SbAA14 and Celluclast® 1.5 L was observed for HPAC-H2SO4-200 mM pretreated bamboo, in which the degree of synergy reached approximately 1.61. The distinctive substrate preference of SbAA14 indicated that it is a novel AA14 LPMO that may act mainly on heteroxylan with numerous arabinosyl substituents between cellulose fibers rather than on recalcitrant xylan tightly associated with cellulose. These findings broaden the understanding of enigmatic AA14 LPMOs and provide new insights into the substrate specificities and biological functionalities of AA14 LPMOs in fungi.

Exploring preferred binding domains of IgG1 mAbs to multimodal adsorbents using a combined biophysics and simulation approach

In this work, we employ a recently developed biophysical technique that uses diethylpyrocarbonate (DEPC) covalent labeling and mass spectrometry for the identification of mAb binding patches to two multimodal cation exchange resins at different pH. This approach compares the labeling results obtained in the bound and unbound states to identify residues that are sterically shielded and thus located in the mAb binding domains. The results at pH 6 for one mAb (mAb B) indicated that while the complementarity determining region (CDR) had minimal interactions with both resins, the FC domain was actively involved in binding. In contrast, DEPC/MS data with another mAb (mAb C) indicated that both the CDR and FC domains were actively involved in binding. These results corroborated chromatographic retention data with these two mAbs and their fragments and helped to explain the significantly stronger retention of both the intact mAb C and its Fab fragment. In contrast, labeling results with mAb C at pH 7, indicated that only the CDR played a significant role in resin binding, again corroborating chromatographic data. The binding domains identified from the DEPC/MS experiments were also examined using protein surface hydrophobicity maps obtained using a recently developed sparse sampling molecular dynamics (MD) approach in concert with electrostatic potential maps. These results demonstrate that the DEPC covalent labeling/mass spectrometry technique can provide important information about the domain contributions of multidomain proteins such as monoclonal antibodies when interacting with multimodal resins over a range of pH conditions.

PL-PatchSurfer3: Improved Structure-Based Virtual Screening for Structure Variation Using 3D Zernike Descriptors

Structure-based virtual screening (SBVS) is a widely used method in silico drug discovery. It necessitates a receptor structure or binding site to predict the binding pose and fitness of a ligand. Therefore, the performance of the SBVS is affected by the protein conformation. The most frequently used method in SBVS is the protein-ligand docking program, which utilizes atomic distance-based scoring functions. Hence, they are highly prone to sensitivity towards variation in receptor structure, and it is reported that the conformational change significantly drops the performance of the docking program. To address the problem, we have introduced a novel program of SBVS, named PL-PatchSurfer. This program makes use of molecular surface patches and the Zernike descriptor. The surfaces of the pocket and ligand are segmented into several patches by the program. These patches are then mapped with physico-chemical properties such as shape and electrostatic potential before being converted into the Zernike descriptor, which is rotationally invariant. A complementarity between the protein and the ligand is assessed by comparing the descriptors and geometric distribution of the patches in the molecules. A benchmarking study showed that PL-PatchSurfer2 was able to screen active molecules regardless of the receptor structure change with fast speed. However, the program could not achieve high performance for the targets that the hydrogen bonding feature is important such as nuclear hormone receptors. In this paper, we present the newer version of PL-PatchSurfer, PL-PatchSurfer3, which incorporates two new features: a change in the definition of hydrogen bond complementarity and consideration of visibility that contains curvature information of a patch. Our evaluation demonstrates that the new program outperforms its predecessor and other SBVS methods while retaining its characteristic tolerance to receptor structure changes. Interested individuals can access the program at kiharalab.org/plps3.

Mechanism and cellular function of direct membrane binding by the ESCRT and ERES-associated Ca2+-sensor ALG-2

Apoptosis linked Gene-2 (ALG-2) is a multifunctional intracellular Ca2+ sensor and the archetypal member of the penta-EF hand protein family. ALG-2 functions in the repair of damage to both the plasma and lysosome membranes and in COPII-dependent budding at endoplasmic reticulum exit sites (ERES). In the presence of Ca2+, ALG-2 binds to ESCRT-I and ALIX in membrane repair and to SEC31A at ERES. ALG-2 also binds directly to acidic membranes in the presence of Ca2+ by a combination of electrostatic and hydrophobic interactions. By combining giant unilamellar vesicle-based experiments and molecular dynamics simulations, we show that charge-reversed mutants of ALG-2 at these locations disrupt membrane recruitment. ALG-2 membrane binding mutants have reduced or abrogated ERES localization in response to Thapsigargin-induced Ca2+ release but still localize to lysosomes following lysosomal Ca2+ release. In vitro reconstitution shows that the ALG-2 membrane-binding defect can be rescued by binding to ESCRT-I. These data thus reveal the nature of direct Ca2+-dependent membrane binding and its interplay with Ca2+-dependent protein binding in the cellular functions of ALG-2.

Insilico screening to identify novel inhibitors targeting 30S ribosomal protein S12 in meningitis-causing organism 'Elizabethkingia meningoseptica'

The current trend in biomedical research is on prioritizing infections based on multidrug resistance. Elizabethkingia meningoseptica, a nosocomial infection-causing organism emerging from Neonatal Intensive Care Units (NICUs), leads to neonatal meningitis and sepsis resulting in severe illness, and, in some cases, fatal. Finding a solution remains challenging due to limited prior work. Translational S12 ribosomal proteins play a crucial role in decoding the codon-anticodon helix, which is essential for the survival of E. meningoseptica. These proteins do not exhibit significant similarity with humans, making them potential drug targets. An in silico study aims to identify specific inhibitors for E. meningoseptica ribosomal proteins among known bioactive compounds targeting prokaryotic 30S ribosomal protein. A 3D model of the 7JIL_h protein from Flavobacterium johnsoniae, showing 90% sequence similarity with the target protein was generated using SWISS-MODEL software. The model was validated through Molprobity v4.4, VERIFY 3D, Errata, and ProSA analysis, confirming conserved residues of the target protein. Insilico screening of known bioactive compounds and their analogs identified potential ligands for the target protein. Molecular Docking and post-docking analysis assessed the stability of the protein-ligand complexes among the shortlisted compounds. The top two compounds with high Gold fitness scores and low predicted binding energy underwent MD simulation and further estimation of free binding energy using the MM_PBSA module. These computationally shortlisted compounds, namely chEMBL 1323619 and chEMBL 312490 may be considered for future in-vivo studies as potential inhibitors against the modeled 30S ribosomal protein S12 of E. meningoseptica.

A novel AA14 LPMO from Talaromyces rugulosus with bifunctional cellulolytic/hemicellulolytic activity boosted cellulose hydrolysis

BACKGROUND

The recently discovered PcAA14A and B from white-rot basidiomycete Pycnoporus coccineus enriched our understanding of the oxidative degradation of xylan in fungi, however, the unusual mode of action of AA14 LPMOs has sparked controversy. The substrate specificity and functionality of AA14 LPMOs still remain enigmatic and need further investigation.

RESULTS

In this study, a novel AA14 LPMO was characterized from the ascomycete Talaromyces rugulosus. TrAA14A has a broad substrate specificity with strong oxidative activity on pure amorphous cellulose and xyloglucan. It could simultaneously oxidize cellulose, xylan and xyloglucan in natural hemi/cellulosic substrate such as fibrillated eucalyptus pulp, and released native and oxidized cello-oligosaccharides, xylo-oligosaccharides and xyloglucan oligosaccharides from this substrate, but its cellulolytic/hemicellulolytic activity became weaker as the contents of xylan increase in the alkaline-extracted hemi/cellulosic substrates. The dual cellulolytic/hemicellulolytic activity enables TrAA14A to possess a profound boosting effect on cellulose hydrolysis by cellulolytic enzymes. Structure modelling of TrAA14A revealed that it exhibits a relatively flat active-site surface similar to the active-site surfaces in AA9 LPMOs but quite distinct from PcAA14B, despite TrAA14A is strongly clustered together with AA14 LPMOs. Remarkable difference in electrostatic potentials of L2 and L3 surfaces was also observed among TrAA14A, PcAA14B and NcLPMO9F. We speculated that the unique feature in substrate-binding surface might contribute to the cellulolytic/hemicellulolytic activity of TrAA14A.

CONCLUSIONS

The extensive cellulolytic/hemicellulolytic activity on natural hemi/cellulosic substrate indicated that TrAA14A from ascomycete is distinctively different from previously characterized xylan-active AA9 or AA14 LPMOs. It may play as a bifunctional enzyme to decompose some specific network structures formed between cellulose and hemicellulose in the plant cell walls. Our findings shed new insights into the novel substrate specificities and biological functionalities of AA14 LPMOs, and will contribute to developing novel bifunctional LPMOs as the booster in commercial cellulase cocktails to efficiently break down the hemicellulose-cellulose matrix in lignocellulose.

The reversible activation of norovirus by metal ions

Murine norovirus (MNV) undergoes extremely large conformational changes in response to the environment. The T = 3 icosahedral capsid is composed of 180 copies of ~58-kDa VP1 comprised of N-terminus (N), shell (S), and C-terminal protruding (P) domains. At neutral pH, the P domains are loosely tethered to the shell and float ~15 Å above the surface. At low pH or in the presence of bile salts, the P domain drops onto the shell and this movement is accompanied by conformational changes within the P domain that enhance receptor interactions while blocking antibody binding. While previous crystallographic studies identified metal binding sites in the isolated P domain, the ~2.7-Å cryo-electron microscopy structures of MNV in the presence of Mg2+ or Ca2+ presented here show that metal ions can recapitulate the contraction observed at low pH or in the presence of bile. Further, we show that these conformational changes are reversed by dialysis against EDTA. As observed in the P domain crystal structures, metal ions bind to and contract the G'H' loop. This movement is correlated with the lifting of the C'D' loop and rotation of the P domain dimers about each other, exposing the bile salt binding pocket. Isothermal titration calorimetry experiments presented here demonstrate that the activation signals (bile salts, low pH, and metal ions) act in a synergistic manner that, individually, all result in the same activated structure. We present a model whereby these reversible conformational changes represent a uniquely dynamic and tissue-specific structural adaptation to the in vivo environment.IMPORTANCEThe highly mobile protruding domains on the calicivirus capsids are recognized by cell receptor(s) and antibodies. At neutral pH, they float ~15 Å above the shell but at low pH or in the presence of bile salts, they contract onto the surface. Concomitantly, changes within the P domain block antibody binding while enhancing receptor binding. While we previously demonstrated that metals also block antibody binding, it was unknown whether they might also cause similar conformational changes in the virion. Here, we present the near atomic cryo-electron microscopy structures of infectious murine norovirus (MNV) in the presence of calcium or magnesium ions. The metal ions reversibly induce the same P domain contraction as low pH and bile salts and act in a synergistic manner with the other stimuli. We propose that, unlike most other viruses, MNV facilely changes conformations as a unique means to escape immune surveillance as it moves through various tissues.

Mutation of AtPME2, a pH-Dependent Pectin Methylesterase, Affects Cell Wall Structure and Hypocotyl Elongation

Pectin methylesterases (PMEs) modify homogalacturonan's chemistry and play a key role in regulating primary cell wall mechanical properties. Here, we report on Arabidopsis AtPME2, which we found to be highly expressed during lateral root emergence and dark-grown hypocotyl elongation. We showed that dark-grown hypocotyl elongation was reduced in knock-out mutant lines as compared to the control. The latter was related to the decreased total PME activity as well as increased stiffness of the cell wall in the apical part of the hypocotyl. To relate phenotypic analyses to the biochemical specificity of the enzyme, we produced the mature active enzyme using heterologous expression in Pichia pastoris and characterized it through the use of a generic plant PME antiserum. AtPME2 is more active at neutral compared to acidic pH, on pectins with a degree of 55-70% methylesterification. We further showed that the mode of action of AtPME2 can vary according to pH, from high processivity (at pH8) to low processivity (at pH5), and relate these observations to the differences in electrostatic potential of the protein. Our study brings insights into how the pH-dependent regulation by PME activity could affect the pectin structure and associated cell wall mechanical properties.

In silico approaches to investigate enzyme immobilization: a comprehensive systematic review

Enzymes are popular catalysts with many applications, especially in industry. Biocatalyst usage on a large scale is facing some limitations, such as low operational stability, low recyclability, and high enzyme cost. Enzyme immobilization is a beneficial strategy to solve these problems. Bioinformatics tools can often correctly predict immobilization outcomes, resulting in a cost-effective experimental phase with the least time consumed. This study provides an overview of in silico methods predicting immobilization processes via a comprehensive systematic review of published articles till 11 December 2022. It also mentions the strengths and weaknesses of the processes and explains the computational analyses in each method that are required for immobilization assessment. In this regard, Web of Science and Scopus databases were screened to gain relevant publications. After screening the gathered documents (n = 3873), 60 articles were selected for the review. The selected papers have applied in silico procedures including only molecular dynamics (MD) simulations (n = 20), parallel tempering Monte Carlo (PTMC) and MD simulations (n = 3), MD and docking (n = 1), density functional theory (DFT) and MD (n = 1), only docking (n = 11), metal ion binding site prediction (MIB) server and docking (n = 2), docking and DFT (n = 1), docking and analysis of enzyme surfaces (n = 1), only DFT (n = 1), only MIB server (n = 2), analysis of an enzyme structure and surface (n = 12), rational design of immobilized derivatives (RDID) software (n = 3), and dissipative particle dynamics (DPD; n = 2). In most included studies (n = 51), enzyme immobilization was investigated experimentally in addition to in silico evaluation.

Rebamipide and Derivatives are Potent, Selective Inhibitors of Histidine Phosphatase Activity of the Suppressor of T Cell Receptor Signaling Proteins

The suppressor of T cell receptor signaling (Sts) proteins are negative regulators of immune signaling. Genetic inactivation of these proteins leads to significant resistance to infection. From a 590,000 compound high-throughput screen, we identified the 2-(1H)-quinolinone derivative, rebamipide, as a putative inhibitor of Sts phosphatase activity. Rebamipide, and a small library of derivatives, are competitive, selective inhibitors of Sts-1 with IC50 values from low to submicromolar. SAR analysis indicates that the quinolinone, the acid, and the amide moieties are all essential for activity. A crystal structure confirmed the SAR and reveals key interactions between this class of compound and the protein. Although rebamipide has poor cell permeability, we demonstrated that a liposomal preparation can inactivate the phosphatase activity of Sts-1 in cells. These studies demonstrate that Sts-1 enzyme activity can be pharmacologically inactivated and provide foundational tools and insights for the development of immune-enhancing therapies that target the Sts proteins.

Structural and biophysical properties of a [4Fe4S] ferredoxin-like protein from Synechocystis sp. PCC 6803 with a unique two domain structure

Electron carrier proteins (ECPs), binding iron-sulfur clusters, are vital components within the intricate network of metabolic and photosynthetic reactions. They play a crucial role in the distribution of reducing equivalents. In Synechocystis sp. PCC 6803, the ECP network includes at least nine ferredoxins. Previous research, including global expression analyses and protein binding studies, has offered initial insights into the functional roles of individual ferredoxins within this network. This study primarily focuses on Ferredoxin 9 (slr2059). Through sequence analysis and computational modeling, Ferredoxin 9 emerges as a unique ECP with a distinctive two-domain architecture. It consists of a C-terminal iron‑sulfur binding domain and an N-terminal domain with homology to Nil-domain proteins, connected by a structurally rigid 4-amino acid linker. Notably, in contrast to canonical [2Fe2S] ferredoxins exemplified by PetF (ssl0020), which feature highly acidic surfaces facilitating electron transfer with photosystem I reaction centers, models of Ferredoxin 9 reveal a more neutral to basic protein surface. Using a combination of electron paramagnetic resonance spectroscopy and square-wave voltammetry on heterologously produced Ferredoxin 9, this study demonstrates that the protein coordinates 2×[4Fe4S]2+/1+ redox-active and magnetically interacting clusters, with measured redox potentials of -420 ± 9 mV and - 516 ± 10 mV vs SHE. A more in-depth analysis of Fdx9's unique structure and protein sequence suggests that this type of Nil-2[4Fe4S] multi-domain ferredoxin is well conserved in cyanobacteria, bearing structural similarities to proteins involved in homocysteine synthesis in methanogens.

Bothrops snake venom L-amino acid oxidases impair biofilm formation of clinically relevant bacteria

The present work addressed the abilities of two L-amino acid oxidases isolated from Bothrops moojeni (BmooLAAO-I) and Bothrops jararacussu (BjussuLAAO-II) snake venoms to control the growth and prevent the biofilm formation of clinically relevant bacterial pathogens. Upon S. aureus (ATCC BAA44) and S. aureus (clinical isolates), BmooLAAO-I (MIC = 0.12 and 0.24 μg/mL, respectively) and BjussuLAAO-II (MIC = 0.15 μg/mL) showed a potent bacteriostatic effect. Against E. coli (ATCC BAA198) and E. coli (clinical isolates), BmooLAAO-I (MIC = 15.6 and 62.5 μg/mL, respectively) and BjussuLAAO-II (MIC = 4.88 and 9.76 μg/mL, respectively) presented a lower extent effect. Also, BmooLAAO-I (MICB50 = 0.195 μg/mL) and BjussuLAAO-II (MICB50 = 0.39 μg/mL) inhibited the biofilm formation of S. aureus (clinical isolates) in 88% and 89%, respectively, and in 89% and 53% of E. coli (clinical isolates). Moreover, scanning electron microscopy confirmed that the toxins affected bacterial morphology by increasing the roughness of the cell surface and inhibited the biofilm formation. Furthermore, analysis of the tridimensional structures of the toxins showed that the surface-charge distribution presents a remarkable positive region close to the glycosylation motif, which is more pronounced in BmooLAAO-I than BjussuLAAO-II. This region may assist the interaction with bacterial and biofilm surfaces. Collectively, our findings propose that venom-derived antibiofilm agents are promising biotechnological tools which could provide novel strategies for biofilm-associated infections.