Inactivation of polyphenol oxidase by low intensity DC field: Experiment and mechanism analysis via molecular dynamics simulation and molecular docking.

In this study, inactivation of mushroom polyphenol oxidase (PPO) by low intensity direct current (DC) electric field and its molecular mechanism were investigated. In the experiments under 3 V/cm, 5 V/cm, 7 V/cm and 9 V/cm electric fields, PPOs were all completely inactivated after different exposure times. Under 1 V/cm, a residual activity of 11.88 % remained. The inactivation kinetics confirms to Weibull model. Under 1-7 V/cm, n value closes to a constant about 1.3. The structural analysis of PPO under 3 V/cm and 5 V/cm by fluorescence emission spectroscopy and molecular dynamics (MD) simulation showed that the tertiary structure was slightly changed with increased radius of gyration, higher potential energy and rate of C-alpha fluctuation. After exposure to the electric field, most of the hydrophobic tryptophan (TRP) residues turned to the hydrophilic surface, resulting the fluorescence red-shifted and quenched. Molecular docking indicated that the receptor binding domain of catechol in PPO was changed. PPO under electric field was MD simulated the first time, revealing the changing mechanism of the electric field itself on PPO, a binuclear copper enzyme, which has a metallic center. All these suggest that the low intensity DC electric field would be a promising option for enzymatic browning inhibition or even enzyme activity inactivation.

Hairpin protein partitioning from the ER to lipid droplets involves major structural rearrangements.

Lipid droplet (LD) function relies on proteins partitioning between the endoplasmic reticulum (ER) phospholipid bilayer and the LD monolayer membrane to control cellular adaptation to metabolic changes. It has been proposed that these hairpin proteins integrate into both membranes in a similar monotopic topology, enabling their passive lateral diffusion during LD emergence at the ER. Here, we combine biochemical solvent-accessibility assays, electron paramagnetic resonance spectroscopy and intra-molecular crosslinking experiments with molecular dynamics simulations, and determine distinct intramembrane positionings of the ER/LD protein UBXD8 in ER bilayer and LD monolayer membranes. UBXD8 is deeply inserted into the ER bilayer with a V-shaped topology and adopts an open-shallow conformation in the LD monolayer. Major structural rearrangements are required to enable ER-to-LD partitioning. Free energy calculations suggest that such structural transition is unlikely spontaneous, indicating that ER-to-LD protein partitioning relies on more complex mechanisms than anticipated and providing regulatory means for this trans-organelle protein trafficking.

Network pharmacology, molecular docking, and dynamics analyses to predict the antiviral activity of ginger constituents against coronavirus infection.

COVID-19 is a global pandemic that caused a dramatic loss of human life worldwide, leading to accelerated research for antiviral drug discovery. Herbal medicine is one of the most commonly used alternative medicine for the prevention and treatment of many conditions including respiratory system diseases. In this study, a computational pipeline was employed, including network pharmacology, molecular docking simulations, and molecular dynamics simulations, to analyze the common phytochemicals of ginger rhizomes and identify candidate constituents as viral inhibitors. Furthermore, experimental assays were performed to analyze the volatile and non-volatile compounds of ginger and to assess the antiviral activity of ginger oil and hydroalcoholic extract. Network pharmacology analysis showed that ginger compounds target human genes that are involved in related cellular processes to the viral infection. Docking analysis highlighted five pungent compounds and zingiberenol as potential inhibitors for the main protease (Mpro), spike receptor-binding domain (RBD), and human angiotensin-converting enzyme 2 (ACE2). Then, (6)-gingerdiacetate was selected for molecular dynamics (MD) simulations as it exhibited the best binding interactions and free energies over the three target proteins. Trajectories analysis of the three complexes showed that RBD and ACE2 complexes with the ligand preserved similar patterns of root mean square deviation (RMSD) and radius of gyration (Rg) values to their respective native structures. Finally, experimental validation of the ginger hydroalcoholic extract confirmed the existence of (6)-gingerdiacetate and revealed the strong antiviral activity of the hydroalcoholic extract with IC 50 of 2.727 µ g / ml . Our study provides insights into the potential antiviral activity of (6)-gingerdiacetate that may enhance the host immune response and block RBD binding to ACE2, thereby, inhibiting SARS-CoV-2 infection.

Functional and structural dissection of glycosyltransferases underlying the glycodiversity of wolfberry-derived bioactive ingredients lycibarbarspermidines.

Lycibarbarspermidines are unusual phenolamide glycosides characterized by a dicaffeoylspermidine core with multiple glycosyl substitutions, and serve as a major class of bioactive ingredients in the wolfberry. So far, little is known about the enzymatic basis of the glycosylation of phenolamides including dicaffeoylspermidine. Here, we identify five lycibarbarspermidine glycosyltransferases, LbUGT1-5, which are the first phenolamide-type glycosyltransferases and catalyze regioselective glycosylation of dicaffeoylspermidines to form structurally diverse lycibarbarspermidines in wolfberry. Notably, LbUGT3 acts as a distinctive enzyme that catalyzes a tandem sugar transfer to the ortho-dihydroxy group on the caffeoyl moiety to form the unusual ortho-diglucosylated product, while LbUGT1 accurately discriminates caffeoyl and dihydrocaffeoyl groups to catalyze a site-selective sugar transfer. Crystal structure analysis of the complexes of LbUGT1 and LbUGT3 with UDP, combined with molecular dynamics simulations, revealed the structural basis of the difference in glycosylation selectivity between LbUGT1 and LbUGT3. Site-directed mutagenesis illuminates a conserved tyrosine residue (Y389 in LbUGT1 and Y390 in LbUGT3) in PSPG box that plays a crucial role in regulating the regioselectivity of LbUGT1 and LbUGT3. Our study thus sheds light on the enzymatic underpinnings of the chemical diversity of lycibarbarspermidines in wolfberry, and expands the repertoire of glycosyltransferases in nature.

Identification of small molecular inhibitors of SIRT3 by computational and biochemical approaches a potential target of breast cancer.

Sirtuin 3 (SIRT3) belongs to the Sirtuin protein family, which consists of NAD+-dependent lysine deacylase, involved in the regulation of various cellular activities. Dysregulation of SIRT3 activity has been linked to several types of cancer, including breast cancer. Because of its ability to stimulate adaptive metabolic pathways, it can aid in the survival and proliferation of breast cancer cells. Finding new chemical compounds targeted towards SIRT3 was the primary goal of the current investigation. Virtual screening of ~ 800 compounds using molecular docking techniques yielded 8 active hits with favorable binding affinities and poses. Docking studies verified that the final eight compounds formed stable contacts with the catalytic domain of SIRT3. Those compounds have good pharmacokinetic/dynamic properties and gastrointestinal absorption. Based on excellent pharmacokinetic and pharmacodynamic properties, two compounds (MI-44 and MI-217) were subjected to MD simulation. Upon drug interaction, molecular dynamics simulations demonstrate mild alterations in the structure of proteins and stability. Binding free energy calculations revealed that compounds MI-44 (- 45.61 ± 0.064 kcal/mol) and MI-217 (- 41.65 ± 0.089 kcal/mol) showed the maximum energy, suggesting an intense preference for the SIRT3 catalytic site for attachment. The in-vitro MTT assay on breast cancer cell line (MDA-MB-231) and an apoptotic assay for these potential compounds (MI-44/MI-217) was also performed, with flow cytometry to determine the compound's ability to cause apoptosis in breast cancer cells. The percentage of apoptotic cells (including early and late apoptotic cells) increased from 1.94% in control to 79.37% for MI-44 and 85.37% for MI-217 at 15 µM. Apoptotic cell death was effectively induced by these two compounds in a flow cytometry assay indicating them as a good inhibitor of human SIRT3. Based on our findings, MI-44 and MI-217 merit additional investigation as possible breast cancer therapeutics.

Core proteome mediated subtractive approach for the identification of potential therapeutic drug target against the honeybee pathogen Paenibacillus larvae.

Background & Objectives: American foulbrood (AFB), caused by the highly virulent, spore-forming bacterium Paenibacillus larvae, poses a significant threat to honey bee brood. The widespread use of antibiotics not only fails to effectively combat the disease but also raises concerns regarding honey safety. The current computational study was attempted to identify a novel therapeutic drug target against P. larvae, a causative agent of American foulbrood disease in honey bee. Methods: We investigated effective novel drug targets through a comprehensive in silico pan-proteome and hierarchal subtractive sequence analysis. In total, 14 strains of P. larvae genomes were used to identify core genes. Subsequently, the core proteome was systematically narrowed down to a single protein predicted as the potential drug target. Alphafold software was then employed to predict the 3D structure of the potential drug target. Structural docking was carried out between a library of phytochemicals derived from traditional Chinese flora (n > 36,000) and the potential receptor using Autodock tool 1.5.6. Finally, molecular dynamics (MD) simulation study was conducted using GROMACS to assess the stability of the best-docked ligand. Results: Proteome mining led to the identification of Ketoacyl-ACP synthase III as a highly promising therapeutic target, making it a prime candidate for inhibitor screening. The subsequent virtual screening and MD simulation analyses further affirmed the selection of ZINC95910054 as a potent inhibitor, with the lowest binding energy. This finding presents significant promise in the battle against P. larvae. Conclusions: Computer aided drug design provides a novel approach for managing American foulbrood in honey bee populations, potentially mitigating its detrimental effects on both bee colonies and the honey industry.

Influence of point mutations on PR65 conformational adaptability: Insights from molecular simulations and nanoaperture optical tweezers.

PR65 is the HEAT repeat scaffold subunit of the heterotrimeric protein phosphatase 2A (PP2A) and an archetypal tandem repeat protein. Its conformational mechanics plays a crucial role in PP2A function by opening/closing substrate binding/catalysis interface. Using in silico saturation mutagenesis, we identified PR65 "hinge" residues whose substitutions could alter its conformational adaptability and thereby PP2A function, and selected six mutations that were verified to be expressed and soluble. Molecular simulations and nanoaperture optical tweezers revealed consistent results on the specific effects of the mutations on the structure and dynamics of PR65. Two mutants observed in simulations to stabilize extended/open conformations exhibited higher corner frequencies and lower translational scattering in experiments, indicating a shift toward extended conformations, whereas another displayed the opposite features, confirmed by both simulations and experiments. The study highlights the power of single-molecule nanoaperture-based tweezers integrated with in silico approaches for exploring the effect of mutations on protein structure and dynamics.

On the role of native contact cooperativity in protein folding.

The consistency of energy landscape theory predictions with available experimental data, as well as direct evidence from molecular simulations, have shown that protein folding mechanisms are largely determined by the contacts present in the native structure. As expected, native contacts are generally energetically favorable. However, there are usually at least as many energetically favorable nonnative pairs owing to the greater number of possible nonnative interactions. This apparent frustration must therefore be reduced by the greater cooperativity of native interactions. In this work, we analyze the statistics of contacts in the unbiased all-atom folding trajectories obtained by Shaw and coworkers, focusing on the unfolded state. By computing mutual cooperativities between contacts formed in the unfolded state, we show that native contacts form the most cooperative pairs, while cooperativities among nonnative or between native and nonnative contacts are typically much less favorable or even anticooperative. Furthermore, we show that the largest network of cooperative interactions observed in the unfolded state consists mainly of native contacts, suggesting that this set of mutually reinforcing interactions has evolved to stabilize the native state.

Assessing the Quality of Cotranscriptional Folding Simulations.

Structural changes in RNAs are an important contributor to controlling gene expression not only at the posttranscriptional stage but also during transcription. A subclass of riboswitches and RNA thermometers located in the 5' region of the primary transcript regulates the downstream functional unit - usually an ORF - through premature termination of transcription. Not only such elements occur naturally, but they are also attractive devices in synthetic biology. The possibility to design such riboswitches or RNA thermometers is thus of considerable practical interest. Since these functional RNA elements act already during transcription, it is important to model and understand the dynamics of folding and, in particular, the formation of intermediate structures concurrently with transcription. Cotranscriptional folding simulations are therefore an important step to verify the functionality of design constructs before conducting expensive and labor-intensive wet lab experiments. For RNAs, full-fledged molecular dynamics simulations are far beyond practical reach because of both the size of the molecules and the timescales of interest. Even at the simplified level of secondary structures, further approximations are necessary. The BarMap approach is based on representing the secondary structure landscape for each individual transcription step by a coarse-grained representation that only retains a small set of low-energy local minima and the energy barriers between them. The folding dynamics between two transcriptional elongation steps is modeled as a Markov process on this representation. Maps between pairs of consecutive coarse-grained landscapes make it possible to follow the folding process as it changes in response to transcription elongation. In its original implementation, the BarMap software provides a general framework to investigate RNA folding dynamics on temporally changing landscapes. It is, however, difficult to use in particular for specific scenarios such as cotranscriptional folding. To overcome this limitation, we developed the user-friendly BarMap-QA pipeline described in detail in this contribution. It is illustrated here by an elaborate example that emphasizes the careful monitoring of several quality measures. Using an iterative workflow, a reliable and complete kinetics simulation of a synthetic, transcription-regulating riboswitch is obtained using minimal computational resources. All programs and scripts used in this contribution are free software and available for download as a source distribution for Linux® or as a platform-independent Docker® image including support for Apple macOS® and Microsoft Windows®.

Fast collective motions of backbone in transmembrane α helices are critical to water transfer of aquaporin.

Fast collective motions are widely present in biomolecules, but their functional relevance remains unclear. Herein, we reveal that fast collective motions of backbone are critical to the water transfer of aquaporin Z (AqpZ) by using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and molecular dynamics (MD) simulations. A total of 212 residue site-specific dipolar order parameters and 158 15N spin relaxation rates of the backbone are measured by combining the 13C- and 1H-detected multidimensional ssNMR spectra. Analysis of these experimental data by theoretic models suggests that the small-amplitude (~10°) collective motions of the transmembrane α helices on the nanosecond-to-microsecond timescales are dominant for the dynamics of AqpZ. The MD simulations demonstrate that these collective motions are critical to the water transfer efficiency of AqpZ by facilitating the opening of the channel and accelerating the water-residue hydrogen bonds renewing in the selectivity filter region.

In Silico Comparison of Bioactive Compounds Characterized from Azadirachta indica with an FDA-Approved Drug against Schistosomal Agents: New Insight into Schistosomiasis Treatment.

The burden of human schistosomiasis, a known but neglected tropical disease in Sub-Saharan Africa, has been worrisome in recent years. It is becoming increasingly difficult to tackle schistosomiasis with praziquantel, a drug known to be effective against all Schistosoma species, due to reports of reduced efficacy and resistance. Therefore, this study seeks to investigate the antischistosomal potential of phytochemicals from Azadirachta indica against proteins that have been implicated as druggable targets for the treatment of schistosomiasis using computational techniques. In this study, sixty-three (63) previously isolated and characterized phytochemicals from A. indica were identified from the literature and retrieved from the PubChem database. In silico screening was conducted to assess the inhibitory potential of these phytochemicals against three receptors (Schistosoma mansoni Thioredoxin glutathione reductase, dihydroorotate dehydrogenase, and Arginase) that may serve as therapeutic targets for schistosomiasis treatment. Molecular docking, ADMET prediction, ligand interaction, MMGBSA, and molecular dynamics simulation of the hit compounds were conducted using the Schrodinger molecular drug discovery suite. The results show that Andrographolide possesses a satisfactory pharmacokinetic profile, does not violate the Lipinski rule of five, binds with favourable affinity with the receptors, and interacts with key amino acids at the active site. Importantly, its interaction with dihydroorotate dehydrogenase, an enzyme responsible for the catalysis of the de novo pyrimidine nucleotide biosynthetic pathway rate-limiting step, shows a glide score and MMGBSA of -10.19 and -45.75 Kcal/mol, respectively. In addition, the MD simulation shows its stability at the active site of the receptor. Overall, this study revealed that Andrographolide from Azadirachta indica could serve as a potential lead compound for the development of an anti-schistosomal drug.

Identifying the Multitarget Pharmacological Mechanism of Action of Genistein on Lung Cancer by Integrating Network Pharmacology and Molecular Dynamic Simulation.

Food supplements have become beneficial as adjuvant therapies for many chronic disorders, including cancer. Genistein, a natural isoflavone enriched in soybeans, has gained potential interest as an anticancer agent for various cancers, primarily by modulating apoptosis, the cell cycle, and angiogenesis and inhibiting metastasis. However, in lung cancer, the exact impact and mechanism of action of genistein still require clarification. To provide more insight into the mechanism of action of genistein, network pharmacology was employed to identify the key targets and their roles in lung cancer pathogenesis. Based on the degree score, the hub genes AKT1, CASP3, EGFR, STAT3, ESR1, SRC, PTGS2, MMP9, PRAG, and AR were significantly correlated with genistein treatment. AKT1, EGFR, and STAT3 were enriched in the non-small cell lung cancer (NSCLC) pathway according to Kyoto Encyclopedia of Genes and Genomes analysis, indicating a significant connection to lung cancer development. Moreover, the binding affinity of genistein to NSCLC target proteins was further verified by molecular docking and molecular dynamics simulations. Genistein exhibited potential binding to AKT1, which is involved in apoptosis, cell migration, and metastasis, thus holding promise for modulating AKT1 function. Therefore, this study aimed to investigate the mechanism of action of genistein and its therapeutic potential for the treatment of NSCLC.

Integrating Explicit and Implicit Fullerene Models into UNRES Force Field for Protein Interaction Studies.

Fullerenes, particularly C60, exhibit unique properties that make them promising candidates for various applications, including drug delivery and nanomedicine. However, their interactions with biomolecules, especially proteins, remain not fully understood. This study implements both explicit and implicit C60 models into the UNRES coarse-grained force field, enabling the investigation of fullerene-protein interactions without the need for restraints to stabilize protein structures. The UNRES force field offers computational efficiency, allowing for longer timescale simulations while maintaining accuracy. Five model proteins were studied: FK506 binding protein, HIV-1 protease, intestinal fatty acid binding protein, PCB-binding protein, and hen egg-white lysozyme. Molecular dynamics simulations were performed with and without C60 to assess protein stability and investigate the impact of fullerene interactions. Analysis of contact probabilities reveals distinct interaction patterns for each protein. FK506 binding protein (1FKF) shows specific binding sites, while intestinal fatty acid binding protein (1ICN) and uteroglobin (1UTR) exhibit more generalized interactions. The explicit C60 model shows good agreement with all-atom simulations in predicting protein flexibility, the position of C60 in the binding pocket, and the estimation of effective binding energies. The integration of explicit and implicit C60 models into the UNRES force field, coupled with recent advances in coarse-grained modeling and multiscale approaches, provides a powerful framework for investigating protein-nanoparticle interactions at biologically relevant scales without the need to use restraints stabilizing the protein, thus allowing for large conformational changes to occur. These computational tools, in synergy with experimental techniques, can aid in understanding the mechanisms and consequences of nanoparticle-biomolecule interactions, guiding the design of nanomaterials for biomedical applications.

Synthesis, In Vivo Anticonvulsant Activity Evaluation and In Silico Studies of Some Quinazolin-4(3H)-One Derivatives.

Two series, "a" and "b", each consisting of nine chemical compounds, with 2,3-disubstituted quinazolin-4(3H)-one scaffold, were synthesized and evaluated for their anticonvulsant activity. They were investigated as dual potential positive allosteric modulators of the GABAA receptor at the benzodiazepine binding site and inhibitors of carbonic anhydrase II. Quinazolin-4(3H)-one derivatives were evaluated in vivo (D1-3 = 50, 100, 150 mg/kg, administered intraperitoneally) using the pentylenetetrazole (PTZ)-induced seizure model in mice, with phenobarbital and diazepam, as reference anticonvulsant agents. The in silico studies suggested the compounds act as anticonvulsants by binding on the allosteric site of GABAA receptor and not by inhibiting the carbonic anhydrase II, because the ligands-carbonic anhydrase II predicted complexes were unstable in the molecular dynamics simulations. The mechanism targeting GABAA receptor was confirmed through the in vivo flumazenil antagonism assay. The pentylenetetrazole experimental anticonvulsant model indicated that the tested compounds, 1a-9a and 1b-9b, present a potential anticonvulsant activity. The evaluation, considering the percentage of protection against PTZ, latency until the onset of the first seizure, and reduction in the number of seizures, revealed more favorable results for the "b" series, particularly for compound 8b.

Ensemble-Based Virtual Screening Led to the Discovery of Novel Lead Molecules as Potential NMBAs.

Neuromuscular blocking agents (NMBAs) are routinely used during anesthesia to relax skeletal muscle. Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels; NMBAs can induce muscle paralysis by preventing the neurotransmitter acetylcholine (ACh) from binding to nAChRs situated on the postsynaptic membranes. Despite widespread efforts, it is still a great challenge to find new NMBAs since the introduction of cisatracurium in 1995. In this work, an effective ensemble-based virtual screening method, including molecular property filters, 3D pharmacophore model, and molecular docking, was applied to discover potential NMBAs from the ZINC15 database. The results showed that screened hit compounds had better docking scores than the reference compound d-tubocurarine. In order to further investigate the binding modes between the hit compounds and nAChRs at simulated physiological conditions, the molecular dynamics simulation was performed. Deep analysis of the simulation results revealed that ZINC257459695 can stably bind to nAChRs' active sites and interact with the key residue Asp165. The binding free energies were also calculated for the obtained hits using the MM/GBSA method. In silico ADMET calculations were performed to assess the pharmacokinetic properties of hit compounds in the human body. Overall, the identified ZINC257459695 may be a promising lead compound for developing new NMBAs as an adjunct to general anesthesia, necessitating further investigations.

Reactive Force Field Molecular Dynamics Investigation of NH3 Generation Mechanism during Protein Pyrolysis Process.

The mechanism of ammonia formation during the pyrolysis of proteins in biomass is currently unclear. To further investigate this issue, this study employed the AMS 2023.104 software to select proteins (actual proteins) as the model compounds and the amino acids contained within them (assembled amino acids) as the comparative models. ReaxFF molecular dynamics simulations were conducted to explore the nitrogen transformation and NH3 generation mechanisms in three-phase products (char, tar, and gas) during protein pyrolysis. The research results revealed several key findings. Regardless of whether the model compounds are actual proteins or assembled amino acids, NH3 is the primary nitrogen-containing product during pyrolysis. However, as the temperature rises to higher levels, such as 2000 K and 2500 K, the amount of NH3 decreases significantly in the later stages of pyrolysis, indicating that it is being converted into other nitrogen-bearing species, such as HCN and N2. Simultaneously, we also observed significant differences between the pyrolysis processes of actual proteins and assembled amino acids. Notably, at 2000 K, the amount of NH3 generated from the pyrolysis of assembled amino acids was twice that of actual proteins. This discrepancy mainly stems from the inherent structural differences between proteins and amino acids. In proteins, nitrogen is predominantly present in a network-like structure (NH-N), which shields it from direct external exposure, thus requiring more energy for nitrogen to participate in pyrolysis reactions, making it more difficult for NH3 to form. Conversely, assembled amino acids can release NH3 through a simpler deamination process, leading to a significant increase in NH3 production during their pyrolysis.

Green and Efficient Extraction of Phenolic Components from Plants with Supramolecular Solvents: Experimental and Theoretical Studies.

The supramolecular solvent (SUPRAS) has garnered significant attention as an innovative, efficient, and environmentally friendly solvent for the effective extraction and separation of bioactive compounds from natural resources. However, research on the use of a SUPRAS for the extraction of phenolic compounds from plants, which are highly valued in food products due to their exceptional antioxidant properties, remains scarce. The present study developed a green, ultra-sound-assisted SUPRAS method for the simultaneous determination of three phenolic acids in Prunella vulgaris using high-performance liquid chromatography (HPLC). The experimental parameters were meticulously optimized. The efficiency and antioxidant properties of the phenolic compounds obtained using different extraction methods were also compared. Under optimal conditions, the extraction efficiency of the SUPRAS, prepared with octanoic acid reverse micelles dispersed in ethanol-water, significantly exceeded that of conventional organic solvents. Moreover, the SUPRAS method demonstrated greater antioxidant capacity. Confocal laser scanning microscopy (CLSM) images revealed the spherical droplet structure of the SUPRAS, characterized by a well-defined circular fluorescence position, which coincided with the position of the phenolic acids. The phenolic acids were encapsulated within the SUPRAS droplets, indicating their efficient extraction capacity. Furthermore, molecular dynamics simulations combined with CLSM supported the proposed method's mechanism and theoretically demonstrated the superior extraction performance of the SUPRAS. In contrast to conventional methods, the higher extraction efficiency of the SUPRAS can be attributed to the larger solvent contact surface area, the formation of more types of hydrogen bonds between the extractants and the supramolecular solvents, and stronger, more stable interaction forces. The results of the theoretical studies corroborate the experimental outcomes.

Structural insights into drug transport by an aquaglyceroporin.

Pentamidine and melarsoprol are primary drugs used to treat the lethal human sleeping sickness caused by the parasite Trypanosoma brucei. Cross-resistance to these two drugs has recently been linked to aquaglyceroporin 2 of the trypanosome (TbAQP2). TbAQP2 is the first member of the aquaporin family described as capable of drug transport; however, the underlying mechanism remains unclear. Here, we present cryo-electron microscopy structures of TbAQP2 bound to pentamidine or melarsoprol. Our structural studies, together with the molecular dynamic simulations, reveal the mechanisms shaping substrate specificity and drug permeation. Multiple amino acids in TbAQP2, near the extracellular entrance and inside the pore, create an expanded conducting tunnel, sterically and energetically allowing the permeation of pentamidine and melarsoprol. Our study elucidates the mechanism of drug transport by TbAQP2, providing valuable insights to inform the design of drugs against trypanosomiasis.

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.

Novel anthraquinone amide derivatives as potential glyoxalase-I inhibitors.

This study aimed to identify novel Glyoxalase-I (Glo-I) inhibitors with potential anticancer properties, focusing on anthraquinone amide-based derivatives. We synthesized a series of these derivatives and conducted in silico docking studies to predict their binding interactions with Glo-I. In vitro assessments were performed to evaluate the anti-Glo-I activity of the synthesized compounds. A comprehensive structure-activity relationship (SAR) analysis identified key features responsible for specific binding affinities of anthraquinone amide-based derivatives to Glo-I. Additionally, a 100 ns molecular dynamics simulation assessed the stability of the most potent compound compared to a co-crystallized ligand. Compound MQ3 demonstrated a remarkable inhibitory effect against Glo-I, with an IC50 concentration of 1.45 µM. The inhibitory potency of MQ3 may be attributed to the catechol ring, amide functional group, and anthraquinone moiety, collectively contributing to a strong binding affinity with Glo-I. Anthraquinone amide-based derivatives exhibit substantial potential as Glo-I inhibitors with prospective anticancer activity. The exceptional inhibitory efficacy of compound MQ3 indicates its potential as an effective anticancer agent. These findings underscore the significance of anthraquinone amide-based derivatives as a novel class of compounds for cancer therapy, supporting further research and advancements in targeting the Glo-I enzyme to combat cancer.