Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values

The differential impacts of equivalent gating-charge mutations in voltage-gated sodium channels

Voltage-gated sodium (Nav) channels are pivotal for cellular signaling, and mutations in Nav channels can lead to excitability disorders in cardiac, muscular, and neural tissues. A major cluster of pathological mutations localizes in the voltage-sensing domains (VSDs), resulting in either gain-of-function, loss-of-function effects, or both. However, the mechanism behind this functional diversity of mutations at equivalent positions remains elusive. Through hotspot analysis, we identified three gating charges (R1, R2, and R3) as major mutational hotspots in VSDs. The same amino acid substitutions at equivalent gating-charge positions in VSDI and VSDII of the cardiac sodium channel Nav1.5 show differential gating property impacts in electrophysiology measurements. We conducted molecular dynamics (MD) simulations on wild-type channels and six mutants to elucidate the structural basis of their differential impacts. Our 120-µs MD simulations with applied external electric fields captured VSD state transitions and revealed the differential structural dynamics between equivalent R-to-Q mutants. Notably, we observed transient leaky conformations in some mutants during structural transitions, offering a detailed structural explanation for gating-pore currents. Our salt-bridge network analysis uncovered VSD-specific and state-dependent interactions among gating charges, countercharges, and lipids. This detailed analysis revealed how mutations disrupt critical electrostatic interactions, thereby altering VSD permeability and modulating gating properties. By demonstrating the crucial importance of considering the specific structural context of each mutation, our study advances our understanding of structure-function relationships in Nav channels. Our work establishes a robust framework for future investigations into the molecular basis of ion channel-related disorders.

Ligand discrimination in hOR1A1 based on the capacitive response

Odorant discrimination mechanisms are based on the differential interactions between odorant molecules and olfactory receptors (ORs). Biohybrid sensors based on ORs described to date show selectivity towards specific versus non-specific binding of odorants, being unable to distinguish between specific ligands of different affinity. Here we disclose a method that enables odorant discrimination based on the modulation of the capacitive response of the receptor, which allows the differentiation of three high-affinity hOR1A1 agonists. We performed voltammetry and impedance measurements of the hOR1A1 receptor selectively immobilized on a gold electrode in the absence and presence of the agonists. Binding induces a decrease in the capacitive response of the receptor that is proportional to the ligand potency, reaching up to a 40% decrease for the cognate ligand dihydrojasmone, which is attributed to changes in the magnitude and orientation of the electric dipole in the receptor, thereby regulating its response to the applied electric field.

Structural insight into CD20/CD3-bispecific antibodies by molecular modeling

Non-Hodgkin's Lymphoma (NHL) remains a significant challenge in hematology, with chemotherapy and radiation therapy as conventional treatment options, albeit with limitations such as adverse effects. Immunotherapy, particularly bispecific antibodies (BsAbs) T cell engagers (TCEs), has emerged as a promising approach. Despite their potential, TCEs pose challenges, including adverse events like cytokine release syndrome. Understanding the structural details of TCEs and their interactions with target proteins is crucial for optimizing their therapeutic efficacy and toxicity. In this study, we further developed our protocol MCCS-Docker for protein-protein interactions and applied it to investigate the structural intricacies of CD3 interactions with therapeutic antibodies such as OKT3, UCHT1, Mosunetuzumab, Odronextumab, Glofitamab, and Epcoritamab using computational modeling techniques. Our analysis not only approved the effectiveness of our updated MCCS-Docker protocol but also revealed detailed binding interactions between the BsAbs and CD3, elucidating key residues of Tyrosine and Asparagine in the antibodies involved in the binding interface. Molecular dynamics simulations validated the stability of these interactions over time, confirming the reliability of the binding poses generated from docking studies. Overall, our study offered a novel method to predict critical residues in protein-protein interactions and enhanced the understanding of the structural determinants governing BsAb interactions with target proteins, offering valuable insights for designing and optimizing immunotherapeutic agents for NHL and related hematologic malignancies.

Different receptor models show differences in ligand binding strength and location: a computational drug screening for the tick-borne encephalitis virus

The tick-borne encephalitis virus (TBE) is a neurotrophic disease that has spread more rapidly throughout Europe and Asia in the past few years. At the same time, no cure or specific therapy is known to battle the illness apart from vaccination. To find a pharmacologically relevant drug, a computer-aided drug screening was initiated. Such a procedure probes a possible binding of a drug to the RNA Polymerase of TBE. The crystal structure of the receptor, however, includes missing and partially modeled regions, which rendered the structure incomplete and of questionable use for a thorough drug screening procedure. The quality of the receptor model was addressed by studying three putative structures created. We show that the choice of receptor models greatly influences the binding affinity of potential drug molecules and that the binding location could also be significantly impacted. We demonstrate that some drug candidates are unsuitable for one model but show decent results for another. Without any prejudice on the three employed receptor models, the study reveals the imperative need to investigate the receptor structure before drug binding is probed whether experimentally or computationally.

Bat genomes illuminate adaptations to viral tolerance and disease resistance

Zoonoses are infectious diseases transmitted from animals to humans. Bats have been suggested to harbour more zoonotic viruses than any other mammalian order1. Infections in bats are largely asymptomatic2,3, indicating limited tissue-damaging inflammation and immunopathology. To investigate the genomic basis of disease resistance, the Bat1K project generated reference-quality genomes of ten bat species, including potential viral reservoirs. Here we describe a systematic analysis covering 115 mammalian genomes that revealed that signatures of selection in immune genes are more prevalent in bats than in other mammalian orders. We found an excess of immune gene adaptations in the ancestral chiropteran branch and in many descending bat lineages, highlighting viral entry and detection factors, and regulators of antiviral and inflammatory responses. ISG15, which is an antiviral gene contributing to hyperinflammation during COVID-19 (refs. 4,5), exhibits key residue changes in rhinolophid and hipposiderid bats. Cellular infection experiments show species-specific antiviral differences and an essential role of protein conjugation in antiviral function of bat ISG15, separate from its role in secretion and inflammation in humans. Furthermore, in contrast to humans, ISG15 in most rhinolophid and hipposiderid bats has strong anti-SARS-CoV-2 activity. Our work reveals molecular mechanisms that contribute to viral tolerance and disease resistance in bats.

Active Site Engineering of a Glycerol Dehydrogenase Improves its Oxidative Activity and Scope Toward Glycerol Derivatives

Regioselective oxidation of glyceryl alkyl ethers is of utmost importance for the fabrication of substituted hydroxy ketones and enantiopure 1,2-diols as green solvents and pharmaceutical building blocks, respectively. An engineered glycerol dehydrogenase from Bacillus stearothermophilus was described to perform the regioselective oxidation of alkyl glycerol ethers, identifying position 252 as key for accepting larger substrates than glycerol. In this work, we further engineer that position through partial saturation mutagenesis to broaden the substrate scope toward other glycerol derivatives, improving enzyme kinetics and minimizing product inhibition. In particular, the BsGlyDH-L252S variant becomes the most efficient biocatalyst for the deracemization of alkyl glyceryl ethers in a two-step, one-pot immobilized system. The discovery and use of these alternative mutants of GlyDH opens the road to more applications and increases the enzymatic toolbox for the modifications of glyceryl ethers.

Novel inhibitors of the (VIBVN) NAT protein identified through pharmacophore modeling

Arylamine N-acetyltransferases (NATs, E.C. 2.3.1.5) constitute a family of phase II drug metabolizing enzymes. These enzymes catalyze the transfer of acetyl groups from acetyl-CoA to a variety of substrates including arylamines, arylhydrazines, and N-hydroxyarylamines. By facilitating these reactions, NATs play a pivotal role in the detoxification and metabolic processing of a wide range of drugs and carcinogens. NAT in marine V. vulnificus plays a role in the metabolism of drugs, leading to the development of drug resistance in marine V. vulnificus. However, inhibitors targeted marine V. vulnificus NAT [(VIBVN)NAT] remain unclear. Therefore, our research aimed to identify potential hit compounds that target (VIBVN)NAT. We integrated multiple computational approaches to screen for effective inhibitors. From this process, we identified two hit compounds, AK-968-11563024 and AG-205-36710025, with IC50 values of 18.86 µM and 33.27 µM, respectively. Molecular dynamics simulations further elucidated the binding mechanism between (VIBVN)NAT and AK-968-11563024. Our study revealed that AK-968-11563024 forms stable interactions with PHE124, HIS167, and TRP230, which may contribute to its biological activity. Our findings provide a valuable foundation for the future development of drugs targeted therapeutics against (VIBVN)NAT.

Development of a Second-Generation, In Vivo Chemical Probe for PIKfyve

We optimized our highly potent and cell-active chemical probe for phosphatidylinositol-3-phosphate 5-kinase (PIKfyve), SGC-PIKFYVE-1, resulting in compounds with improved potency and demonstrated in vivo stability. Use of an in-cell, kinome-wide selectivity panel allowed for confirmation of excellent in-cell selectivity of our lead compound, 40, and another promising analogue, 46. Evaluation of the pharmacokinetic (PK) profiles of these two compounds revealed that both are well tolerated systemically and orally bioavailable. Coupled with its subnanomolar cellular potency and impressive selectivity in cells, the long half-life of 40 makes it an ideal candidate for the evaluation of the consequences of PIKfyve inhibition in vivo. PIKfyve inhibition has been investigated clinically for indications including rheumatoid arthritis, Crohn's disease, COVID-19, and ALS using a single compound (apilimod), supporting the development of orthogonal PIKfyve inhibitors with in vivo stability.

No Bridge between Us: EXAFS and Computations Confirm Two Distant Iron Ions Comprise the Active Site of Alkane Monooxygenase (AlkB)

Alkane monooxygenase (AlkB) is the dominant enzyme that catalyzes the oxidation of liquid alkanes in the environment. Two recent structural models derived from cryo-electron microscopy (cryo-EM) reveal an unusual active site: a histidine-rich center that binds two iron ions without a bridging ligand. To ensure that potential photoreduction and radiation damage are not responsible for the absence of a bridging ligand in the cryo-EM structures, spectroscopic methods are needed. We present the results of extended X-ray absorption fine structure (EXAFS) experiments collected under conditions where photodamage was avoided. Careful data analysis reveals an active site structure consistent with the cryo-EM structures in which the two iron ions are ligated by nine histidines and separated by at least 5 Å. The EXAFS data were used to inform structural models for molecular dynamics (MD) simulations. The MD simulations corroborate EXAFS observations that neither of the two conserved carboxylate-containing residues (E281 and D190) near the active site are likely candidates for metal ion bridging. Mutagenesis experiments, spectroscopy, and additional MD simulations were used to further explore the role of these carboxylate residues. A variant in which a carboxylate containing residue (E281) was changed to a methyl residue (E281A) showed little change in pre-edge features, consistent with the observation that it is not essential for activity and hence unlikely to serve as a bridging ligand at any point in the catalytic cycle. D190 variants had substantially diminished activity, suggesting an important role in catalysis not yet fully understood.

Structure and function of a β-1,2-galactosidase from Bacteroides xylanisolvens, an intestinal bacterium

Galactosides are major carbohydrates that are found in plant cell walls and various prebiotic oligosaccharides. Studying the detailed biochemical functions of β-galactosidases in degrading these carbohydrates is important. In particular, identifying β-galactosidases with new substrate specificities could help in the production of potentially beneficial oligosaccharides. In this study, we identify a β-galactosidase with novel substrate specificity from Bacteroides xylanisolvens, an intestinal bacterium. The enzyme do not show hydrolytic activity toward natural β-galactosides during the first screening. However, when α-D-galactosyl fluoride (α-GalF) as a donor substrate and galactose or D-fucose as an acceptor substrate are incubated with a nucleophile mutant, reaction products are detected. The galactobiose produced from the α-GalF and galactose is identified as β-1,2-galactobiose using NMR. Kinetic analysis reveals that this enzyme effectively hydrolyzes β-1,2-galactobiose and β-1,2-galactotriose. In the complex structure with methyl β-galactopyranose as a ligand, the ligand is only located at subsite +1. The 2-hydroxy group and the anomeric methyl group of methyl β-galactopyranose faces in the direction of subsite -1 and the solvent, respectively. This observation is consistent with the substrate specificity of the enzyme regarding linkage position and chain length. Overall, we conclude that the enzyme is a β-galactosidase acting on β-1,2-galactooligosaccharides.

Towards New Anti-Inflammatory Agents: Design, Synthesis and Evaluation of Molecules Targeting XIAP-BIR2

The X-chromosome-linked inhibitor of apoptosis protein (XIAP) plays a crucial role in controlling cell survival across multiple regulated cell death pathways and coordinating a range of inflammatory signalling events. The discovery of selective inhibitors for XIAP-BIR2, able to disrupt the direct physical interaction between XIAP and RIPK2, offer promising therapeutic options for NOD2-mediated diseases like Crohn's disease, sarcoidosis, and Blau syndrome. The objective of this study was to design, synthesize, and evaluate small synthetic molecules with binding selectivity to XIAP-BIR2 domain. To achieve this, we applied an interdisciplinary drug design approach and firstly we have synthesized an initial fragment library to achieve a first XIAP inhibition activity. Then using a growing strategy, larger compounds were synthesized and one of them presents a good selectivity for XIAP-BIR2 versus XIAP-BIR3 domain, compound 20 c. The ability of compound 20 c to block the NOD1/2 pathway was confirmed in cell models. These data show that we have synthesized molecules capable of blocking NOD1/2 signalling pathways in cellulo, and ultimately leading to new anti-inflammatory compounds.

QM/MM study reveals novel mechanism of KRAS and KRASG12R catalyzed GTP hydrolysis

As a crucial drug target, KRAS can regulate most cellular processes involving guanosine triphosphate (GTP) hydrolysis. However, the mechanism of GTP hydrolysis has remained controversial over the past decades. Here, several different GTP hydrolysis mechanisms catalyzed by wild-type KRAS (WT-KRAS) and KRASG12R mutants were discussed via four QM/MM calculation models. Based on the computational results, a Mg2+-coordinated H2O-mediated GTP hydrolysis mechanism was proposed. In this mechanism, a Mg2+-coordinated H2O first protonates the fully deprotonated GTP, and then the Mg2+ coordinated hydroxyl anion is generated. The Pγ-O bond is formed via the SN2 attack of the second H2O on the Pγ atom of the GTP, leading to the simultaneous cleavage of the Pγ-O bond. Meanwhile, the hydroxyl anion coordinated to Mg2+ and generated in the first step acts as a proton acceptor from water. This Mg2+ coordinated H2O-involved GTP hydrolysis mechanism may also be suitable for Mg2+-catalyzed ATP hydrolysis. Furthermore, the mechanism of GTP hydrolysis catalyzed by the KRASG12R mutant, whose hydrolysis rate was approximately 40-fold slower than WT-KRAS, was also discussed. Our QM/MM calculations reveal that GTP is easily protonated by the residue R12, and the energy barrier of GTP hydrolysis catalyzed by the KRASG12R mutant is lower than the corresponding barrier for WT-KRAS. Nevertheless, molecular dynamics (MD) simulations reveal that R12, a residue characterized by significant steric hindrance, is positioned at the GTP site where the nucleophilic attack by water occurs during Pγ-O bond formation, thereby strongly impeding the approach of water molecules to GTP. As a result, the GTP hydrolysis rate catalyzed by the KRASG12R mutant was severely impaired. Uncovering the GTP hydrolysis mechanism catalyzed by the WT-KRAS and KRASG12R mutant may also give a reasonable explanation for the relationship between the KRASG12R mutation and the occurrence of cancer. We hope this finding will provide useful guidance for drug discovery that targets KRAS.

Additive CHARMM Force Field for Pterins and Folates

Folates comprise a crucial class of biologically active compounds related to folic acid, playing a vital role in numerous enzymatic reactions. One-carbon metabolism, facilitated by the folate cofactor, supports numerous physiological processes, including biosynthesis, amino acid homeostasis, epigenetic maintenance, and redox defense. Folates share a common pterin heterocyclic ring structure capable of undergoing redox reactions and existing in various protonation states. This study aimed to derive molecular mechanics (MM) parameters compatible with the CHARMM36 all-atom additive force field for pterins and biologically important folates, including pterin, biopterin, and folic acid. Three redox forms were considered: oxidized, dihydrofolate, and tetrahydrofolate states. Across all protonation states, a total of 18 folates were parameterized. Partial charges were derived using the CHARMM force field parametrization protocol, based on targeting reference quantum mechanics monohydrate interactions, electrostatic potential, and dipole moment. Bonded terms were parameterized using one-dimensional adiabatic potential energy surface scans, and two-dimensional scans to parametrize in-ring torsions associated with the puckering states of dihydropterin and tetrahydropterin. The quality of the model was demonstrated through simulations of three protein complexes using optimized and initial parameters. These simulations underscored the significantly enhanced performance of the folate model developed in this study compared to the initial model without optimization in reproducing structural properties of folate-protein complexes. Overall, the presented MM model will be valuable for modeling folates in various redox states and serve as a starting point for parameterizing other folate derivatives.

Exploring the Potential of Malvidin and Echiodinin as Probable Antileishmanial Agents Through In Silico Analysis and In Vitro Efficacy

Leishmaniasis, a neglected tropical disease caused by Leishmania species, presents serious public health challenges due to limited treatment options, toxicity, high costs, and drug resistance. In this study, the in vitro potential of malvidin and echioidinin is examined as antileishmanial agents against L. amazonensis, L. braziliensis, and L. infantum, comparing their effects to amphotericin B (AmpB), a standard drug. Malvidin demonstrated greater potency than echioidinin across all parasite stages and species. Against L. amazonensis, malvidin's IC50 values were 197.71 ± 17.20 µM (stationary amastigotes) and 258.07 ± 17 µM (axenic amastigotes), compared to echioidinin's 272.99 ± 29.90 μM and 335.96 ± 19.35 μM. AmpB was more potent, with IC50 values of 0.06 ± 0.01 µM and 0.10 ± 0.03 µM. Malvidin exhibited lower cytotoxicity (CC50: 2920.31 ± 80.29 µM) than AmpB (1.06 ± 0.12 µM) and a favorable selectivity index. It reduced infection rates by 35.75% in L. amazonensis-infected macrophages. The in silico analysis revealed strong binding between malvidin and Leishmania arginase, with the residues HIS139 and PRO258 playing key roles. Gene expression analysis indicated malvidin's modulation of oxidative stress and DNA repair pathways, involving genes like GLO1 and APEX1. These findings suggest malvidin's potential as a safe, natural antileishmanial compound, warranting further in vivo studies to confirm its therapeutic efficacy and pharmacokinetics in animal models.

Determinants of Product Outcome in Two Sesquiterpene Synthases from the Thermotolerant Bacterium Rubrobacter radiotolerans

Rubrobacter radiotolerans nerolidol synthase (NerS) and trans-α-bergamotene synthase (BerS) are among the first terpene synthases (TPSs) discovered from thermotolerant bacteria, and, despite sharing the same substrate, make terpenoid products with different carbon scaffolds. Here, the potential thermostability of NerS and BerS was investigated, and NerS was found to retain activity up to 55 °C. A library of 22 NerS and BerS variants was designed to probe the differing reaction mechanisms of NerS and BerS, including residues putatively involved in substrate sequestration, cation-π stabilisation of reactive intermediates, and shaping of the active site contour. Two BerS variants showed improved in vivo titres vs the WT enzyme, and also yielded different ratios of the related sesquiterpenoids (E)-β-farnesene and trans-α-bergamotene. BerS-L86F was proposed to encourage substrate isomerisation by cation-π stabilisation of the first cationic intermediate, resulting in a greater proportion of trans-α-bergamotene. By contrast, BerS-S82L significantly preferred (E)-β-farnesene formation, attributed to steric blocking of the isomerisation step, consistent with what has been observed in several plant TPSs. Our work highlights the importance of isomerisation as a key determinant of product outcome in TPSs, and shows how a combined computational and experimental approach can characterise TPSs and variants with improved and altered functionality.

pH-Dependent Assembly and Stability of Toll-Like Receptor 3/dsRNA Signaling Complex: Insights from Constant pH Molecular Dynamics and Metadynamics Simulations

The pH-dependent assembly of Toll-like receptors (TLRs), which triggers a threshold-like response, is a key principle in immune signaling. While crystallography has revealed the intricate structure of these assembly complexes, the mechanisms underlying their pH dependency remain unclear. Herein, constant pH simulations and metadynamics are employed to investigate the pH-dependent assembly and stability of the TLR3/dsRNA signaling complex. The findings demonstrate that system pH regulates complex assembly and stability by modulating the protonation and charge states of histidines. Histidines in TLR3 act as pH-dependent, positively charged binding sites that capture negatively charged dsRNA. Additionally, these histidines form a [H682⁺]-[E626⁻] dipole, facilitating the assembly of two TLR3 molecules into an antisymmetric dimer through dipole-dipole interactions. Surprisingly, TLR3 can shift the pKa values of key histidines from their model pKa of 6.5, increasing protonation likelihood and enhancing ligand binding. Notably, the aromatic residue Phe84, located within the dsRNA binding site [His39⁺-His60⁺-Phe84-His108⁺], alters the pKa of His60 through cation-π interactions with its protonated state. This study offers new insights into the molecular mechanisms underlying pH-dependent immune signaling via higher-order assemblies and suggests potential applications for histidine in self-assembling biomaterials.

Nanobody screening and machine learning guided identification of cross-variant anti-SARS-CoV-2 neutralizing heavy-chain only antibodies

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) continues to persist, demonstrating the risks posed by emerging infectious diseases to national security, public health, and the economy. Development of new vaccines and antibodies for emerging viral threats requires substantial resources and time, and traditional development platforms for vaccines and antibodies are often too slow to combat continuously evolving immunological escape variants, reducing their efficacy over time. Previously, we designed a next-generation synthetic humanized nanobody (Nb) phage display library and demonstrated that this library could be used to rapidly identify highly specific and potent neutralizing heavy chain-only antibodies (HCAbs) with prophylactic and therapeutic efficacy in vivo against the original SARS-CoV-2. In this study, we used a combination of high throughput screening and machine learning (ML) models to identify HCAbs with potent efficacy against SARS-CoV-2 viral variants of interest (VOIs) and concern (VOCs). To start, we screened our highly diverse Nb phage display library against several pre-Omicron VOI and VOC receptor binding domains (RBDs) to identify panels of cross-reactive HCAbs. Using HCAb affinity for SARS-CoV-2 VOI and VOCs (pre-Omicron variants) and model features from other published data, we were able to develop a ML model that successfully identified HCAbs with efficacy against Omicron variants, independent of our experimental biopanning workflow. This biopanning informed ML approach reduced the experimental screening burden by 78% to 90% for the Omicron BA.5 and Omicron BA.1 variants, respectively. The combined approach can be applied to other emerging viruses with pandemic potential to rapidly identify effective therapeutic antibodies against emerging variants.

The structural influence of the oncogenic driver mutation N642H in the STAT5B SH2 domain

The point mutation N642H of the signal transducer and activator of transcription 5B (STAT5B) protein is associated with aggressive and drug-resistant forms of leukemia. This mutation is thought to promote cancer due to hyperactivation of STAT5B caused by increased stability of the active, parallel dimer state. However, the molecular mechanism leading to this stabilization is not well understood as there is currently no structure of the parallel dimer. To investigate the mutation's mechanism of action, we conducted extensive all-atom molecular dynamics simulations of multiple oligomeric forms of both STAT5B and STAT5BN642H, including a model for the parallel dimer. The N642H mutation directly affects the hydrogen bonding network within the phosphotyrosine (pY)-binding pocket of the parallel dimer, enhancing the pY-binding interaction. The simulations indicate that apo STAT5B is highly flexible, exploring a diverse conformational space. In contrast, apo STAT5BN642H accesses two distinct conformational states, one of which resembles the conformation of the parallel dimer. The simulation predictions of the effects of the mutation on structure and dynamics are supported by the results of hydrogen-deuterium exchange (HDX) mass spectrometry measurements carried out on STAT5B and STAT5BN642H in which a phosphopeptide was used to mimic the effects of parallel dimerization on the SH2 domain. The molecular-level information uncovered in this work contributes to our understanding of STAT5B hyperactivation by the N642H mutation and could help pave the way for novel therapeutic strategies targeting this mutation.

Modeling and interaction study of alcohol oxidase and ProteaseA in methylotrophic yeast C. boidinii: insights from In-silico analysis

Flavin adenine nucleotide (FAD)-dependent oxidoreductase enzyme Alcohol oxidase (AOX) facilitates the growth of methylotrophic yeast C. boidinii by catabolizing methanol, producing formaldehyde and hydrogen peroxide. Vacuolar Protease-A (PrA) from C. boidinii is responsible for the proteolytic activity of AOX. However, no experimental structures for these enzymes have been reported. This in-silico study involves modeling and interaction analysis of AOX and PrA. A protein-protein interaction study shows that Thr75, Gly74, Arg72, Tyr73, and Met289 amino acids of PrA have shown interaction with AOX. These residues may be crucial for AOX proteolysis. An in-silico study predicts that serine protease inhibitors bind to specific amino acids, potentially obstructing PrA's degradable activity on AOX. PrA does not interact with the FAD binding sites in AOX. Instead, it interacts with AOX at sites (Ser337, Ala34, and Tyr343) where AOX monomers interact, hindering octamer formation the active form of AOX. During simulation, strong dynamics in PrA were found in the loop regions of the structure, as observed in the complexes. This in-silico work aims to corroborate the experimental research, which lacks structural studies on the proteolysis process.

Mannich Base Derived from Lawsone Inhibits PKM2 and Induces Neoplastic Cell Death

Background/Objectives: Pyruvate kinase M2, a central regulator of cancer cell metabolism, has garnered significant attention as a promising target for disrupting the metabolic adaptability of tumor cells. This study explores the potential of the Mannich base derived from lawsone (MB-6a) to interfere with PKM2 enzymatic activity both in vitro and in silico. Methods: The antiproliferative potential of MB-6a was tested using MTT assay in various cell lines, including SCC-9, Hep-G2, HT-29, B16-F10, and normal human gingival fibroblast (HGF). The inhibition of PKM2 mediated by MB-6a was assessed using an LDH-coupled assay and by measuring ATP production. Docking studies and molecular dynamics calculations were performed using Autodock 4 and GROMACS, respectively, on the tetrameric PKM2 crystallographic structure. Results: The Mannich base 6a demonstrated selective cytotoxicity against all cancer cell lines tested without affecting cell migration, with the highest selectivity index (SI) of 4.63 in SCC-9, followed by B16-F10 (SI = 3.9), Hep-G2 (SI = 3.4), and HT-29 (SI = 2.03). The compound effectively inhibited PKM2 glycolytic activity, leading to a reduction of ATP production both in the enzymatic reaction and in cells treated with this naphthoquinone derivative. MB-6a showed favorable binding to PKM2 in the ATP-bound monomers through docking studies (PDB ID: 4FXF; binding affinity scores ranging from -6.94 to -9.79 kcal/mol) and MD simulations, revealing binding affinities stabilized by key interactions including hydrogen bonds, halogen bonds, and hydrophobic contacts. Conclusions: The findings suggest that MB-6a exerts its antiproliferative activity by disrupting cell glucose metabolism, consequently reducing ATP production and triggering energetic collapse in cancer cells. This study highlights the potential of MB-6a as a lead compound targeting PKM2 and warrants further investigation into its mechanism of action and potential clinical applications.

Crystallographic and Selectivity Studies on the Approved HIF Prolyl Hydroxylase Inhibitors Desidustat and Enarodustat

Prolyl hydroxylase domain-containing proteins 1-3 (PHD1-3) are 2-oxoglutarate (2OG)-dependent oxygenases catalysing C-4 hydroxylation of prolyl residues in α-subunits of the heterodimeric transcription factor hypoxia-inducible factor (HIF), modifications that promote HIF-α degradation via the ubiquitin-proteasome pathway. Pharmacological inhibition of the PHDs induces HIF-α stabilisation, so promoting HIF target gene transcription. PHD inhibitors are used to treat anaemia caused by chronic kidney disease (CKD) due to their ability to stimulate erythropoietin (EPO) production. We report studies on the effects of the approved PHD inhibitors Desidustat and Enarodustat, and the clinical candidate TP0463518, on activities of a representative set of isolated recombinant human 2OG oxygenases. The three molecules manifest selectivity for PHD inhibition over that of the other 2OG oxygenases evaluated. We obtained crystal structures of Desidustat and Enarodustat in complex with the human 2OG oxygenase factor inhibiting hypoxia-inducible factor-α (FIH), which, together with modelling studies, inform on the binding modes of Desidustat and Enarodustat to active site Fe(II) in 2OG oxygenases, including PHD1-3. The results will help in the design of selective inhibitors of both the PHDs and other 2OG oxygenases, which are of medicinal interest due to their involvement inter alia in metabolic regulation, epigenetic signalling, DNA-damage repair, and agrochemical resistance.

Exploring CRISPR-Cas9 HNH-Domain-Catalyzed DNA Cleavage Using Accelerated Quantum Mechanical Molecular Mechanical Free Energy Simulation

The target DNA (tDNA) cleavage catalyzed by the CRISPR Cas9 enzyme is a critical step in the Cas9-based genome editing technologies. Previously, the tDNA cleavage from an active SpyCas9 enzyme conformation was modeled by Palermo and co-workers (Nierzwicki et al., Nat. Catal. 2022 5, 912) using ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulations, where the free energy barrier was found to be more favorable than that from a pseudoactive enzyme conformation. In this work, we performed ai-QM/MM simulations based on another catalytically active conformation (PDB 7Z4J) of the Cas9 HNH domain from cryo-electron microscopy experiments. For the wildtype enzyme, we acquired a free energy profile for the tDNA cleavage that is largely consistent with the previous report. Furthermore, we explored the role of the active-site K866 residue on the catalytic efficiency by modeling the K866A mutant and found that the K866A mutation increased the reaction free energy barrier, which is consistent with the experimentally observed reduction in the enzyme activity.

Accurate Predictions of Molecular Properties of Proteins via Graph Neural Networks and Transfer Learning

Machine learning has emerged as a promising approach for predicting molecular properties of proteins, as it addresses limitations of experimental and traditional computational methods. Here, we introduce GSnet, a graph neural network (GNN) trained to predict physicochemical and geometric properties including solvation free energies, diffusion constants, and hydrodynamic radii, based on three-dimensional protein structures. By leveraging transfer learning, pre-trained GSnet embeddings were adapted to predict solvent-accessible surface area (SASA) and residue-specific pKa values, achieving high accuracy and generalizability. Notably, GSnet outperformed existing protein embeddings for SASA prediction, and a locally charge-aware variant, aLCnet, approached the accuracy of simulation-based and empirical methods for pKa prediction. Our GNN framework demonstrated robustness across diverse datasets, including intrinsically disordered peptides, and scalability for high-throughput applications. These results highlight the potential of GNN-based embeddings and transfer learning to advance protein structure analysis, providing a foundation for integrating predictive models into proteome-wide studies and structural biology pipelines.

Docking heparan sulfate-based ligands as a promising inhibitor for SARS-CoV-2

CONTEXT

Heparan sulfate (HS) linear polysaccharide glycosaminoglycan compound is linked to components from the cell surface and the extracellular matrix. HS mediates SARS-CoV-2 infection through spike protein binding to cell surface receptors and is required to bind ACE2, prompting the need for electronic structure and molecular docking evaluation of this core system to exploit this attachment in developing new derivatives. Therefore, we have studied five molecules based on HS using molecular docking and electronic structure analysis. Non-covalent interaction analysis shows hydrogen bonding and van der Waals interactions in the binding to RBD-ACE2 interface and 3CLpro. SDM3 and SDM1 molecules present the lowest gap, including solvent effect under 154.6 kcal/mol, and exhibit the most reactivity behavior in this group, potentially leading to enhanced interaction in docking studies.

METHODS

Heparan sulfate and four derivatives were optimized using B3LYP functional with two basis sets 6-31 + G(d,p) and def2SVP. Electronic structure was used to explore the main interactions and the reactivity of these molecules, and these optimized structures were used in the molecular docking study against 3CLpro, RBD, and ACE2.

Hypoxia increases methylated histones to prevent histone clipping and heterochromatin redistribution during Raf-induced senescence

Hypoxia enhances histone methylation by inhibiting oxygen- and α-ketoglutarate-dependent demethylases, resulting in increased methylated histones. This study reveals how hypoxia-induced methylation affects histone clipping and the reorganization of heterochromatin into senescence-associated heterochromatin foci (SAHF) during oncogene-induced senescence (OIS) in IMR90 human fibroblasts. Notably, using top-down proteomics, we discovered specific cleavage sites targeted by Cathepsin L (CTSL) in H3, H2B and H4 during Raf activation, identifying novel sites in H2B and H4. Hypoxia counteracts CTSL-mediated histone clipping by promoting methylation without affecting CTSL's activity. This increase in methylation under hypoxia protects against clipping, reshaping the epigenetic landscape and influencing chromatin accessibility, as shown by ATAC-seq analysis. These insights underscore the pivotal role of hypoxia-induced histone methylation in protecting chromatin from significant epigenetic shifts during cellular aging.

Divergent Oxidation Reactions of E- and Z-Allylic Primary Alcohols by an Unspecific Peroxygenase

Unspecific peroxygenases (UPOs) catalyze the selective oxygenation of organic substrates using only hydrogen peroxide as the external oxidant. The PaDa-I variant of the UPO from Agrocybe aegerita catalyses the oxidation of Z- and E-allylic alcohols with complementary selectivity, giving epoxide and carboxylic acid/aldehyde products respectively. Both reactions can be performed on preparative scale with isolated yields up to 80 %, and the epoxidations proceed with excellent enantioselectivity (>99 % ee). The divergent reactions can also be used to transform E/Z mixtures of allylic alcohols, enabling both product series to be isolated from a single reaction. The utility of the epoxidation method is exemplified in the total synthesis of both enantiomers of the insect pheromone disparlure, including a highly enantioselective gram-scale transformation. These reactions provide further evidence for the potential of UPOs as catalysts for the scalable preparation of important oxygenated intermediates.

PRDX6 contributes to selenocysteine metabolism and ferroptosis resistance

Selenocysteine (Sec) metabolism is crucial for cellular function and ferroptosis prevention and begins with the uptake of the Sec carrier, selenoprotein P (SELENOP). Following uptake, Sec released from SELENOP is metabolized via selenocysteine lyase (SCLY), producing selenide, a substrate for selenophosphate synthetase 2 (SEPHS2), which provides the essential selenium donor, selenophosphate (H2SePO3-), for the biosynthesis of the Sec-tRNA. Here, we discovered an alternative pathway in Sec metabolism mediated by peroxiredoxin 6 (PRDX6), independent of SCLY. Mechanistically, we demonstrate that PRDX6 can readily react with selenide and interact with SEPHS2, potentially acting as a selenium delivery system. Moreover, we demonstrate the functional significance of this alternative route in human cancer cells, revealing a notable association between elevated expression of PRDX6 and human MYCN-amplified neuroblastoma subtype. Our study sheds light on a previously unrecognized aspect of Sec metabolism and its implications in ferroptosis, offering further possibilities for therapeutic exploitation.

Role of the αC-β4 loop in protein kinase structure and dynamics

Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, local spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.

3D variability analysis reveals a hidden conformational change controlling ammonia transport in human asparagine synthetase

Advances in X-ray crystallography and cryogenic electron microscopy (cryo-EM) offer the promise of elucidating functionally relevant conformational changes that are not easily studied by other biophysical methods. Here we show that 3D variability analysis (3DVA) of the cryo-EM map for wild-type (WT) human asparagine synthetase (ASNS) identifies a functional role for the Arg-142 side chain and test this hypothesis experimentally by characterizing the R142I variant in which Arg-142 is replaced by isoleucine. Support for Arg-142 playing a role in the intramolecular translocation of ammonia between the active site of the enzyme is provided by the glutamine-dependent synthetase activity of the R142 variant relative to WT ASNS, and MD simulations provide a possible molecular mechanism for these findings. Combining 3DVA with MD simulations is a generally applicable approach to generate testable hypotheses of how conformational changes in buried side chains might regulate function in enzymes.

Enhancing the specific activity of 3α-hydroxysteroid dehydrogenase through cross-regional combinatorial mutagenesis

3α-Hydroxysteroid dehydrogenase (3α-HSD) from Comamonas testosteroni is widely used in clinical settings to measure serum total bile acid levels. However, its low enzymatic activity leads to high operational costs. In this study, we employed a combinatorial mutagenesis approach to systematically identify potential key mutation sites within the enzyme. The enzyme molecule was segmented into distinct regions, and a comprehensive strategy integrating substrate pocket engineering, binding energy calculations, and deep learning techniques was used. Through experimental verification, single-point mutants from the mutation library with enhanced enzymatic activity by at least 1.5-fold were identified. Through iterative combinatorial mutations of them, the optimal mutant H119A/R201G/R216L was obtained. This mutant exhibited a specific activity of 34.18 U/mg towards deoxycholic acid, representing a 6.85-fold increase over the wild-type (WT) enzyme. Additionally, the optimal temperature of the mutant increased from 35 °C to 40 °C, and its turnover number and catalytic efficiency increased by 6.4-fold and 9.4-fold, respectively. Quantum mechanics/molecular mechanics (QM/MM) calculations indicated that the energy barrier of the dehydrogenase reaction was reduced in the H119A/R201G/R216L mutant compared to that of the WT enzyme. Specifically, the R201G mutation significantly reduced the electric field strength along the 3α-hydroxyl group, facilitating its deprotonation. This study provides insights into enhancing enzymatic efficiency through strategic mutagenesis and elucidates mechanistic changes that optimize enzyme performance for clinical and biotechnological applications.

Free energy landscape of the PI3Kα C-terminal activation

The gene PIK3CA, encoding the catalytic subunit p110α of PI3Kα, is the second most frequently mutated gene in cancer, with the highest frequency oncogenic mutants occurring in the C-terminus of the kinase domain. The C-terminus has a dual function in regulating the kinase, playing a putative auto-inhibitory role for kinase activity and being absolutely essential for binding to the cell membrane. However, the molecular mechanisms by which these C-terminal oncogenic mutations cause PI3Kα overactivation remain unclear. To understand how a spectrum of C-terminal mutations of PI3Kα alter kinase activity compared to the WT, we perform unbiased and biased Molecular Dynamics simulations of several C-terminal mutants and report the free energy landscapes for the C-terminal "closed-to-open" transition in the WT, H1047R, G1049R, M1043L and N1068KLKR mutants. Results are consistent with HDX-MS experimental data and provide a molecular explanation why H1047R and G1049R reorient the C-terminus with a different mechanism compared to the WT and M1043L and N1068KLKR mutants. Moreover, we show that in the H1047R mutant, the cavity, where the allosteric ligands STX-478 and RLY-2608 bind, is more accessible contrary to the WT. This study provides insights into the molecular mechanisms underlying activation of oncogenic PI3Kα by C-terminal mutations and represents a valuable resource for continued efforts in the development of mutant selective inhibitors as therapeutics.

A Computational Strategy for the Rapid Identification and Ranking of Patient-Specific T Cell Receptors Bound to Neoantigens

T cell receptor (TCR) recognition of a peptide-major histocompatibility complex (pMHC) is crucial for adaptive immune response. The identification of therapeutically relevant TCR-pMHC protein pairs is a bottleneck in the implementation of TCR-based immunotherapies. The ability to computationally design TCRs to target a specific pMHC requires automated integration of next-generation sequencing, protein-protein structure prediction, molecular dynamics, and TCR ranking. A pipeline to evaluate patient-specific, sequence-based TCRs to a target pMHC is presented. Using the three most frequently expressed TCRs from 16 colorectal cancer patients, the protein-protein structure of the TCRs to the target CEA peptide-MHC is predicted using Modeller and ColabFold. TCR-pMHC structures are compared using automated equilibration and successive analysis. ColabFold generated configurations require an ≈2.5× reduction in equilibration time of TCR-pMHC structures compared to Modeller. The structural differences between Modeller and ColabFold are demonstrated by root mean square deviation (≈0.20 nm) between clusters of equilibrated configurations, which impact the number of hydrogen bonds and Lennard-Jones contacts between the TCR and pMHC. TCR ranking criteria that may prioritize TCRs for evaluation of in vitro immunogenicity are identified, and this ranking is validated by comparing to state-of-the-art machine learning-based methods trained to predict the probability of TCR-pMHC binding.

Correlation Between Antimicrobial Structural Classes and Membrane Partitioning: Role of Emerging Lipid Packing Defects

In this study, a combination of bioinformatics and molecular dynamics simulations is employed to investigate the partitioning behavior of different classes of antimicrobial peptides (AMPs) into model membranes. The main objective is to identify any correlations between the structural characteristics of AMPs and their membrane identification and early-stage partitioning mechanisms. The simulation results reveal distinct membrane interactions among the various structural classes of AMPs, particularly in relation to the generation and subsequent interaction with lipid packing defects. Notably, AMPs with a structure-less coil conformation generate a higher number of deep and shallow defects, which are larger in size compared to other classes of AMPs. AMPs with helical component demonstrated the deepest insertion into the membrane. On the other hand, AMPs with a significant percentage of beta sheets tend to adsorb onto the membrane surface, suggesting a potentially distinct partitioning mechanism attributed to their structural rigidity. These findings highlight the diverse membrane interactions and partitioning mechanisms exhibited by different structural classes of AMPs.

Enhanced and Efficient Predictions of Dynamic Ionization through Constant-pH Adiabatic Free Energy Dynamics

Dynamic or structurally induced ionization is a critical aspect of many physical, chemical, and biological processes. Molecular dynamics (MD) based simulation approaches, specifically constant pH MD methods, have been developed to simulate ionization states of molecules or proteins under experimentally or physiologically relevant conditions. While such approaches are now widely utilized to predict ionization sites of macromolecules or to study physical or biological phenomena, they are often computationally expensive and require long simulation times to converge. In this article, using the principles of adiabatic free energy dynamics, we introduce an efficient technique for performing constant pH MD simulations within the framework of the adiabatic free energy dynamics (AFED) approach. We call the new approach pH-AFED. We show that pH-AFED provides highly accurate predictions of protein residue pKa values, with a MUE of 0.5 pKa units when coupled with driven adiabatic free energy dynamics (d-AFED), while reducing the required simulation times by more than an order of magnitude. In addition, pH-AFED can be easily integrated into most constant pH MD codes or implementations and flexibly adapted to work in conjunction with enhanced sampling algorithms that target collective variables. We demonstrate that our approaches, with both pH-AFED standalone as well as pH-AFED combined with collective variable based enhanced sampling, provide promising predictive accuracy, with a MUE of 0.6 and 0.5 pKa units respectively, on a diverse range of proteins and enzymes, ranging up to 186 residues and 21 titratable sites. Lastly, we demonstrate how this approach can be utilized to understand the in vivo performance engineered antibodies for immunotherapy.

Comprehensive molecular epidemiology of influenza viruses in Brazil: insights from a nationwide analysis

Influenza A and B viruses represent significant global health threats, contributing substantially to morbidity and mortality rates. However, a comprehensive understanding of the molecular epidemiology of these viruses in Brazil, a continental-size country and a crucial hub for the entry, circulation, and dissemination of influenza viruses within South America, still needs to be improved. This study addresses this gap by consolidating data and samples across all Brazilian macroregions, as part of the Center for Viral Surveillance and Serological Assessment project, together with an extensive number of other Brazilian sequences provided by a public database during the epidemic seasons spanning 2021-23. Phylogenetic analysis of the hemagglutinin segment of influenza A/H1N1pdm09, A/H3N2, and influenza B/Victoria-lineage viruses revealed that in 2021 and in the first semester of 2022, the A/H3N2 2a.3 strain was the predominant circulating strain. Subsequently, the A/H3N2 2b became the prevalent strain until October, when it was substituted by A/H1N1pdm09 5a.2a and 5a.2a.1 lineages. This scenario was maintained during the year of 2023. B/Victoria emerged and circulated at low levels between December 2021 and September 2022 and then became coprevalent with A/H1N1pdm09 5a.2a and 5a.2a.1 lineages. The comparison between the vaccine strain A/Darwin/9/2021 and circulating viruses revealed shared mutations to aspartic acid at residues 186 and 225 across all A/H3N2 lineages from 2021 to 2023, altering the charge in the receptor-binding domain. For A/H1N1pdm09, the 2022 consensus of 5a.2a.1 and the vaccine strain A/Victoria/2570/2019 showed 14 amino acid substitutions. Key residues H180, D187, K219, R223, E224, and T133 are involved in hydrogen interactions with sialic acids, while N130, K142, and D222 may contribute to distance interactions based on docking analyses. Importantly, distinct influenza A lineage frequency patterns were observed across Brazil's macroregions, underscoring the regional variations in virus circulation. This study characterizes influenza A and B viruses circulating in Brazil, providing insights into their circulation patterns and dynamics across Brazilian macroregions. These findings hold significant implications for public health interventions, informing strategies to mitigate transmission risks, optimize vaccination efforts, and enhance outbreak control measures.

Structural Insights into Dopamine Receptor-Ligand Interactions: From Agonists to Antagonists

This study explores the intricacies of dopamine receptor-ligand interactions, focusing on the D1R and D5R subtypes. Using molecular modeling techniques, we investigated the binding of the pan-agonist rotigotine, revealing a universal binding mode at the orthosteric binding pocket. Additionally, we analyze the stability of antagonist-receptor complexes with SKF83566 and SCH23390. By examining the impact of specific mutations on ligand-receptor interactions through computational simulations and thermostability assays, we gain insights into binding stability. Our research also delves into the structural and energetic aspects of antagonist binding to D1R and D5R in their inactive states. These findings enhance our understanding of dopamine receptor pharmacology and hold promise for drug development in central nervous system disorders, opening doors to future research and innovation in this field.

An Interdisciplinary Approach Provides Insights into the Pronounced Selectivity of Compound 42 for DPP9

Dipeptidyl peptidase 8 (DPP8) and 9 (DPP9) are proteases gaining significant attention for their role in health and disease. Distinctive studies of these proteases are hampered by their close homology. Furthermore, designing selective compounds is a major challenge due to the highly conserved catalytic site. Here, we provide mechanistic insights underlying the DPP9-over-DPP8 selectivity of the semi-selective inhibitor "Compound 42". We performed enhanced sampling molecular dynamics simulations to investigate the binding pose of "Compound 42", which enabled the design of various DPP9 mutants that were characterized through a combination of biochemical (Ki determinations) and in silico approaches. Our findings show that DPP9 residue F253 is an important selectivity-determining factor. This work marks the discovery and validation of a structural feature that can be exploited for the design of DPP8 or DPP9 selective inhibitors.

Computational insight into the peptide-based inhibition of α-cobratoxin

Snakebite envenoming results in the death of thousands of people each year and has been classified as a neglected tropical disease by the World Health Organization (WHO). The toxins released into the bloodstream of the victim bind to the nicotinic acetylcholine receptor and restrict transmission of nerve impulses leading to paralysis and cardiac arrest. Conventional antibody-based treatments often have adverse side effects or are difficult to perform. Hence, efforts are underway to devise alternative forms of treatment that address these therapeutic shortcomings. Peptide-based inhibitors have recently gained attention due to their high specificity and ease of preparation. Here, we explore the mechanism of a peptide inhibitor of α-cobratoxin using all-atom molecular dynamics (MD) simulations. We also quantify the energetics of the toxin-peptide dissociation process using the non-equilibrium steered MD technique. Our study reveals that the inhibitor migrates close to Loop-II of α-cobratoxin and alters its dimerization tendency. From energy studies, we found that the peptide first binds to one unit of α-cobratoxin in a particular orientation, followed by the binding of a second toxin molecule, which effectively masks the residues that interact with the nicotinic acetylcholine receptor. Our work provides atomic-level insight into the inhibition process that can be utilized in the future design of inhibitors with superior binding capabilities.

Structure and Energetics of PET-Hydrolyzing Enzyme Complexes: A Systematic Comparison from Molecular Dynamics Simulations

Discovered in 2016, the enzyme PETase, secreted by bacterial Ideonella Sakaiensis 201-F6, has an excellent hydrolytic activity toward poly(ethylene terephthalate) (PET) at room temperature, while it decreases at higher temperatures due to the low thermostability. Many variants have been engineered to overcome this limitation, which hinders industrial application. In this work, we systematically compare PETase wild-type (WT) and four mutants (DuraPETase, ThermoPETase, FastPETase, and HotPETase) using standard molecular dynamics (MD) simulations and unbinding free energy calculations. In particular, we analyze the enzymes' structural characteristics and binding to a tetrameric PET chain (PET4) under two temperature conditions: T1─300 K and T2─350 K. Our results indicate that (i) PET4 forms stable complexes with the five enzymes at room temperature (∼300 K) and (ii) most of the interactions are localized close to the active site of the protein, where the W185 and Y87 residues interact with the aromatic rings of the substrate. Specifically, (iii) the W185 side-chain explores different conformations in each variant (a phenomenon known in the literature as "W185 wobbling"). This suggests that the binding pocket retains structural plasticity and flexibility among the variants, facilitating substrate recognition and localization events at moderate temperatures. Moreover, (iv) PET4 establishes aromatic interactions with the catalytic H237 residue, stabilizing the catalytic triad composed of residues S160-H237-D206, and helping the system achieve an effective configuration for the hydrolysis reaction. Conversely, (v) the binding affinity decreases at a higher temperature (∼350 K), retaining moderate interactions only for HotPETase. Finally, (vi) MD simulations of complexes formed with poly(ethylene-2,5-furan dicarboxylate) (PEF) show no persistent interactions, suggesting that these enzymes are not yet optimized for binding this alternative semiaromatic plastic polymer. Our study offers valuable insights into the structural stability of these enzymes and the molecular determinants driving PET binding onto their surfaces, sheds light on the mechanistic steps that precede the onset of hydrolysis, and provides a foundation for future enzyme optimization.

Exploring the Structural and Dynamic Properties of a Chimeric Glycoside Hydrolase Protein in the Presence of Calcium Ions

GH10 xylanases and GH62 Arabinofuranosidases play key roles in the breakdown of arabinoxylans and are important tools in various industrial and biotechnological processes, such as renewable biofuel production, the paper industry, and the production of short-chain xylooligosaccharides (XOS) from plant biomass. However, the use of these enzymes in industrial settings is often limited due to their relatively low thermostability and reduced catalytic efficiency. To overcome these limitations, strategies based on enzymatic chimera construction and the use of metal ions and other cofactors have been proposed to produce new recombinant enzymes with improved catalytic activity and thermostability. Here, we examine the conformational dynamics of a GH10-GH62 chimera at different calcium ion concentrations through molecular dynamics simulations. While experimental data have demonstrated improved activity and thermostability in GH10-GH62 chimera, the mechanistic basis for these enhancements remains unclear. We explored the structural details of the binding subsites of Ca2+ in the parental enzymes GH62 from Aspergillus fumigatus (Afafu62) and a recombinant GH10 from Cryptococcus flavescens (Xyn10cf), as well as their chimeric combination, and how negatively charged electron pairing located at the protein surface affects Ca2+ capture. The results indicate that Ca2+ binding significantly contributes to structural stability and catalytic cavity modulation in the chimera, particularly evident at a concentration of 0.01 M. This effect, not observed in the parental GH10 and GH62 enzymes, highlights how Ca2+ enhances stability in the overall chimeric enzyme, while supporting a larger cavity volume in the chimera GH62 subunit. The increased catalytic site volume and reduced structural flexibility in response to Ca2+ suggest that calcium binding minimizes non-productive conformational states, which could potentially improve catalytic turnover. The findings presented here may aid in the development of more thermostable and efficient catalytic systems.

Accuracy and Reproducibility of Lipari-Szabo Order Parameters From Molecular Dynamics

The Lipari-Szabo generalized order parameter probes the picosecond to nanosecond time scale motions of a protein and is useful for rationalizing a multitude of biological processes such as protein recognition and ligand binding. Although these fast motions are an important and intrinsic property of proteins, it remains unclear what simulation conditions are most suitable to reproduce methyl symmetry axis side chain order parameter data (Saxis2) from molecular dynamics simulations. In this study, we show that, while Saxis2 tends to converge within tens of nanoseconds, it is essential to run 10 to 20 replicas starting from configurations close to the experimental structure to obtain the best agreement with experimental Saxis2 values. Additionally, in a comparison of force fields, AMBER ff14SB outperforms CHARMM36m in accurately capturing these fast time scale motions, and we suggest that the origin of this performance gap is likely attributed to differences in side chain torsional parametrization and not due to differences in the global protein conformations sampled by the force fields. This study provides insight into obtaining accurate and reproducible Saxis2 values from molecular simulations and underscores the necessity of using replica simulations to compute equilibrium properties.

Oligomerization of protein arginine methyltransferase 1 and its functional impact on substrate arginine methylation

Protein arginine methyltransferases (PRMTs) are important posttranslational modifying enzymes in eukaryotic proteins and regulate diverse pathways from gene transcription, RNA splicing, and signal transduction to metabolism. Increasing evidence supports that PRMTs exhibit the capacity to form higher-order oligomeric structures, but the structural basis of PRMT oligomerization and its functional consequence are elusive. Herein, we revealed for the first time different oligomeric structural forms of the predominant arginine methyltransferase PRMT1 using cryo-EM, which included tetramer (dimer of dimers), hexamer (trimer of dimers), octamer (tetramer of dimers), decamer (pentamer of dimers), and also helical filaments. Through a host of biochemical assays, we showed that PRMT1 methyltransferase activity was substantially enhanced as a result of the high-ordered oligomerization. High-ordered oligomerization increased the catalytic turnover and the multimethylation processivity of PRMT1. Presence of a catalytically dead PRMT1 mutant also enhanced the activity of WT PRMT1, pointing out a noncatalytic role of oligomerization. Structural modeling demonstrates that oligomerization enhances substrate retention at the PRMT1 surface through electrostatic force. Our studies offered key insights into PRMT1 oligomerization and established that oligomerization constitutes a novel molecular mechanism that positively regulates the enzymatic activity of PRMTs in biology.

Slow release of a synthetic auxin induces formation of adventitious roots in recalcitrant woody plants

Clonal propagation of plants by induction of adventitious roots (ARs) from stem cuttings is a requisite step in breeding programs. A major barrier exists for propagating valuable plants that naturally have low capacity to form ARs. Due to the central role of auxin in organogenesis, indole-3-butyric acid is often used as part of commercial rooting mixtures, yet many recalcitrant plants do not form ARs in response to this treatment. Here we describe the synthesis and screening of a focused library of synthetic auxin conjugates in Eucalyptus grandis cuttings and identify 4-chlorophenoxyacetic acid-L-tryptophan-OMe as a competent enhancer of adventitious rooting in a number of recalcitrant woody plants, including apple and argan. Comprehensive metabolic and functional analyses reveal that this activity is engendered by prolonged auxin signaling due to initial fast uptake and slow release and clearance of the free auxin 4-chlorophenoxyacetic acid. This work highlights the utility of a slow-release strategy for bioactive compounds for more effective plant growth regulation.

A Dynamic Loop in Halohydrin Dehalogenase HheG Regulates Activity and Enantioselectivity in Epoxide Ring Opening

Halohydrin dehalogenase HheG and its homologues are remarkable enzymes for the efficient ring opening of sterically demanding internal epoxides using a variety of nucleophiles. The enantioselectivity of the respective wild-type enzymes, however, is usually insufficient for application and frequently requires improvement by protein engineering. We herein demonstrate that the highly flexible N-terminal loop of HheG, comprising residues 39 to 47, has a tremendous impact on the activity as well as enantioselectivity of this enzyme in the ring opening of structurally diverse epoxide substrates. Thus, highly active and enantioselective HheG variants could be accessed through targeted engineering of this loop. In this regard, variant M45F displayed almost 10-fold higher specific activity than wild type in the azidolysis of cyclohexene oxide, yielding the corresponding product (1S,2S)-2-azidocyclohexan-1-ol in 96%eeP (in comparison to 49%eeP for HheG wild type). Moreover, this variant was also improved regarding activity and enantioselectivity in the ring opening of cyclohexene oxide with other nucleophiles, demonstrating even inverted enantioselectivity with cyanide and cyanate. In contrast, a complete loop deletion yielded an inactive enzyme. Concomitant computational analyses of HheG M45F in comparison to wild type enzyme revealed that mutation M45F promotes the productive binding of cyclohexene oxide and azide in the active site by establishing noncovalent C-H ··π interactions between epoxide and F45. These interactions further position one of the two carbon atoms of the epoxide ring closer to the azide, resulting in higher enantioselectivity. Additionally, stable and enantioselective cross-linked enzyme crystals of HheG M45F were successfully generated after combination with mutation D114C. Overall, our study highlights that a highly flexible loop in HheG governs the enzyme's activity and selectivity in epoxide ring opening and should thus be considered in future protein engineering campaigns of HheG.

Evidence for the Quercetin Binding Site of Glycogen Phosphorylase as a Target for Liver-Isoform-Selective Inhibitors against Glioblastoma: Investigation of Flavanols Epigallocatechin Gallate and Epigallocatechin

Glycogen phosphorylase (GP) is the rate-determining enzyme in glycogenolysis, and its druggability has been extensively studied over the years for the development of therapeutics against type 2 diabetes (T2D) and, more recently, cancer. However, the conservation of binding sites between the liver and muscle isoforms makes the inhibitor selectivity challenging. Using a combination of kinetic, crystallographic, modeling, and cellular studies, we have probed the binding of dietary flavonoids epigallocatechin gallate (EGCG) and epigallocatechin (EGC) to GP isoforms. The structures of rmGPb-EGCG and rmGPb-EGC complexes were determined by X-ray crystallography, showing binding at the quercetin binding site (QBS) in agreement with kinetic studies that revealed both compounds as noncompetitive inhibitors of GP, with EGCG also causing a significant reduction in cell viability and migration of U87-MG glioblastoma cells. Interestingly, EGCG exhibits different binding modes to GP isoforms, revealing QBS as a promising site for GP targeting, offering new opportunities for the design of liver-selective GP inhibitors.

Effective patchiness from critical points of a coarse-grained protein model with explicit shape and charge anisotropy

Colloidal model systems are successful in rationalizing emergent phenomena like aggregation, rheology and phase behaviour of protein solutions. Colloidal theory in conjunction with isotropic interaction models is often employed to estimate the stability of such solutions. In particular, a universal criterion for the reduced second virial coefficient at the critical point is frequently invoked which is based on the behavior of short-range attractive fluids (Noro-Frenkel rule, ). However, if anisotropic models for the protein-protein interaction are considered, e.g. the Kern-Frenkel (KF) patchy particle model, the value of the criterion is shifted to lower values and explicitly depends on the number of patches. If an explicit shape anisotropy is considered, as e.g. in a coarse-grained protein model, the normalization of becomes ambiguous to some extent, as no unique exclusion volume can be defined anymore. Here, we investigate a low-resolution, coarse-grained model for the globular protein bovine serum albumin (BSA) and study effects of charge-anisotropy on the phase diagram (determined by simulations) at the isoelectric point. We present methods of assigning an "effective patchiness" to our protein model by comparing its critical properties to the KF model. We find that doubling the native charges increases the critical temperature Tc by ≈14% and that our BSA model can be compared to a 3 to 5 patch KF model. Finally, we argue that applying existing criteria from colloidal theory should be done with care, due to multiple, physically plausible ways of how to assign effective diameters to shape-anisotropic models.

NetSci: A Library for High Performance Biomolecular Simulation Network Analysis Computation

We present the NetSci program-an open-source scientific software package designed for estimating mutual information (MI) between data sets using GPU acceleration and a k-nearest-neighbor algorithm. This approach significantly enhances calculation speed, achieving improvements of several orders of magnitude over traditional CPU-based methods, with data set size limits dictated only by available hardware. To validate NetSci, we accurately compute MI for an analytically verifiable two-dimensional Gaussian distribution and replicate the generalized correlation (GC) analysis previously conducted on the B1 domain of protein G. We also apply NetSci to molecular dynamics simulations of the Sarcoendoplasmic Reticulum Calcium-ATPase (SERCA) pump, exploring the allosteric mechanisms and pathways influenced by ATP and 2'-deoxy-ATP (dATP) binding. Our analysis reveals distinct allosteric effects induced by ATP compared to dATP, with predicted information pathways from the bound nucleotide to the calcium-binding domain differing based on the nucleotide involved. NetSci proves to be a valuable tool for estimating MI and GC in various data sets and is particularly effective for analyzing intraprotein communication and information transfer.

Ene-Reductase-Catalyzed Aromatization of Simple Cyclohexanones to Phenols

Direct aromatization of cyclohexanones to synthesize substituted phenols represents a significant challenge in modern synthetic chemistry. Herein, we describe a novel ene-reductase (TsER) catalytic system that converts substituted cyclohexanones into the corresponding phenols. This process involves the successive dehydrogenation of two saturated carbon-carbon bonds within the six-membered ring of cyclohexanones and utilizes molecular oxygen to drive the reaction cycle. It demonstrates a versatile and efficient approach for the synthesis of substituted phenols, providing a valuable complement to existing chemical methodologies.

Electrostatic Preorganization in Three Distinct Heterogeneous Proteasome β-Subunits

The origin of the enzyme's powerful role in accelerating chemical reactions is one of the most critical and still widely discussed questions. It is already accepted that enzymes impose an electrostatic field onto their substrates by adopting complex three-dimensional structures; therefore, the preorganization of electric fields inside protein active sites has been proposed as a crucial contributor to catalytic mechanisms and rate constant enhancement. In this work, we focus on three catalytically active β-subunits of 20S proteasomes with low sequence identity (∼30%) whose active sites, although situated in an electrostatically miscellaneous environment, catalyze the same chemical reaction with similar catalytic efficiency. Our in silico experiments reproduce the experimentally observed equivalent reactivity of the three sites and show that obliteration of the electrostatic potential in all active sites would deprive the enzymes of their catalytic power by slowing down the chemical process by a factor of 1035. To regain enzymatic efficiency, besides catalytic Thr1 and Lys33 residues, the presence of aspartic acid in position 17 and an aqueous solvent is required, proving that the electrostatic potential generated by the remaining residues is insignificant for catalysis. Moreover, it was found that the gradual decay of atomic charges on Asp17 strongly correlates with the enzyme's catalytic rate deterioration as well as with a change in the charge distributions due to introduced mutations. The computational procedure used and described here may help identify key residues for catalysis in other biomolecular systems and consequently may contribute to the process of designing enzyme-like synthetic catalysts.

Charge transfer characteristics in rhodopsin mimics during photoexcitation

To gain insights into the light-harvesting capabilities of the chromophores, it is essential to understand their molecular and electronic structures within their natural chemical or biological contexts. Rhodopsins display varied absorption characteristics due to the interaction between the chromophore retinal and its surrounding protein environments. In this study, we employed a quantum mechanics/molecular mechanics approach to examine a series of artificially designed rhodopsin mimics based on human cellular retinol acid binding protein 2 (hCRABP II). We elucidated the electron transfer within the all-trans protonated Schiff base upon light excitation, and our calculated absorption spectra show well consistency with the experimental result. Furthermore, the interaction mechanisms between the chromophore and the protein were investigated, and the relationship between the blueshifts and redshifts in the absorption spectra was analyzed. Our calculation demonstrates that the blueshifts and redshifts in the rhodopsin mimics correlate well with attractive (such as the hydrogen bonds or electrostatic interactions) and repulsive interactions (such as the steric effects) between the chromophore and the protein environment, respectively. These findings could provide hints for designing rhodopsin with absorption spectra at different wavelengths.