Network Of Electrostatic Interactions Research Articles

Computational biophysical characterization of the effect of gatekeeper mutations on the binding of ponatinib to the FGFR kinase

Fibroblast Growth Factor Receptor (FGFR) is connected to numerous downstream signalling cascades regulating cellular behavior. Any dysregulation leads to a plethora of illnesses, including cancer. Therapeutics are available, but drug resistance driven by gatekeeper mutation impedes the treatment. Ponatinib is an FDA-approved drug against BCR-ABL kinase and has shown effective results against FGFR-mediated carcinogenesis. Herein, we undertake molecular dynamics simulation-based analysis on ponatinib against all the FGFR isoforms having Val to Met gatekeeper mutations. The results suggest that ponatinib is a potent and selective inhibitor for FGFR1, FGFR2, and FGFR4 gatekeeper mutations. The extensive electrostatic and van der Waals interaction network accounts for its high potency. The FGFR3_VM mutation has shown resistance towards ponatinib, which is supported by their lesser binding affinity than wild-type complexes. The disengaged molecular brake and engaged hydrophobic spine were believed to be the driving factors for weak protein-ligand interaction. Taken together, the inhibitory and structural characteristics exhibited by ponatinib may aid in thwarting resistance based on Val-to-Met gatekeeper mutations at an earlier stage of treatment and advance the design and development of other inhibitors targeted at FGFRs harboring gatekeeper mutations.

Structural effects of charge destabilization and amino acid substitutions in amyloid fragments of CsgA

The most studied functional amyloid is the CsgA, major curli subunit protein, which is produced by numerous strains of Enterobacteriaceae. Although CsgA sequences are highly conserved, they exhibit species diversity, which reflects the specific evolutionary and functional adaptability of the major curli subunit. Herein, we performed bioinformatics analyses to uncover the differences in the amyloidogenic properties of the R4 fragments in Escherichia coli and Salmonella enterica and proposed four mutants for more detailed studies: M1, M2, M3, and M4. The mutated sequences were characterized by various experimental techniques, such as circular dichroism, ATR-FTIR, FT-Raman, thioflavin T, transmission electron microscopy and confocal microscopy. Additionally, molecular dynamics simulations were performed to determine the role of buffer ions in the aggregation process. Our results demonstrated that the aggregation kinetics, fibril morphology, and overall structure of the peptide were significantly affected by the positions of charged amino acids within the repeat sequences of CsgA. Notably, substituting glycine with lysine resulted in the formation of distinctive spherically packed globular aggregates. The differences in morphology observed are attributed to the influence of phosphate ions, which disrupt the local electrostatic interaction network of the polypeptide chains. This study provides knowledge on the preferential formation of amyloid fibrils based on charge states within the polypeptide chain.

Antifouling Coatings from Glassy Polyelectrolyte Complex Films.

Coatings that prevent or decrease fouling are sought for many applications, including those that inhibit the attachment of organisms in aquatic environments. To date, antifouling coatings have mostly followed design criteria assembled over decades: surfaces should be well/strongly hydrated, possess low net charge, and maintain a hydrophilic character when exposed to the location of use. Thus, polymers based on ethylene glycol or zwitterionic repeat units have been shown to be highly effective. Unfortunately, hydrated materials can be quite soft, limiting their use in some environments. In a major paradigm shift, this work describes glassy antifouling films made from certain complexes of positive and negative polyelectrolytes. The dense network of electrostatic interactions yields tough materials below the glass transition temperature, Tg, in normal use, while the highly ionic character of these polyelectrolyte complexes ensures strong hydration. The proximity of equal numbers of opposite charges within these complexes mimics zwitterionic structures. Films, assembled layer-by-layer from aqueous solutions, contained sulfonated poly(ether ether ketone), SPEEK, a rigid polyelectrolyte that binds strongly to a selection of quaternary ammonium polycations. Layer-by-layer buildup of SPEEK and polycations was linear, indicating strong complexes between polyelectrolytes. Calorimetry also showed that complex formation was exothermic. Surfaces coated with these films in the 100 nm thickness range completely resisted adhesion of the common flagellate green algae, Chlamydomonas reinhardtii, which were removed from surfaces at a minimum applied flow rate of 0.8 cm s-1. The total surface charge density of adsorbed cations, determined with a sensitive radioisotopic label, was very low, around 10% of a monolayer, which minimized adsorption driven by counterion release from the surface. The viscoelastic properties of the complexes, which were stable even in concentrated salt solutions, were explored using rheology of bulk samples. When fully hydrated, their Tg values were observed to be above 75 °C.

Exploring the Molecular Dynamics of a Lipid-A Vesicle at the Atom Level: Morphology and Permeation Mechanism.

Lipid-A was previously shown to spontaneously aggregate into a vesicle via the hybrid particle field approach. We assess the validity of the proposed vesiculation mechanism by simulating the resulting lipid-A vesicle at the atom level. The spatial confinement imposed by the vesicle geometry on the conformation and packing of lipid-A induces significant heterogeneity of physical properties in the inner and outer leaflets. It also induces tighter molecular packing and lower acyl chain order compared to the lamellar arrangement. Around 5% of water molecules passively permeates the vesicle membrane inward and outward. The permeation is facilitated by interactions with water molecules that are transported across the membrane by a network of electrostatic interactions with the hydrogen bond donors/acceptors in the N-acetylglucosamine ring and upper region of the acyl chains of lipid-A. The permeation process takes place at low rates but still at higher frequencies than observed for the lamellar arrangement of lipid-A. These findings not only substantiate the proposed lipid-A vesiculation mechanism but also reveal the complex structural dynamics of an important nonlamellar arrangement of lipid-A.

Atomistic Insights into the Origin of High-Performance Thermoelectric Response in Hybrid Perovskites.

Due to their tantalizing prospect of heat-electricity interconversion, hybrid organic-inorganic perovskites have sparked considerable research interests recently. Nevertheless, understanding their complex interplay between the macroscopic properties, nonintuitive transport processes, and basic chemical structures still remains far from completion, although it plays a fundamental role in systematic materials development. On the basis of multiscale first-principles calculations, this understanding is herein advanced by establishing a comprehensive picture consisting of atomic and charge dynamics. It is unveiled that the ultralow room-temperature lattice thermal conductivity (≈0.20W m-1 K-1 ) of hybrid perovskites is critical to their decent thermoelectric figure of merit (≈0.34), and such phonon-glass behavior stems from not only the inherent softness but also the strong anharmonicity. It is identified that the 3D electrostatic interaction and hydrogen-bonded networks between the PbI3- cage and embedded cations result in the strongly coupled motions of inorganic framework and cation, giving rise to their high degree of anharmonicity. Furthermore, such coupled motions bring about low-frequency optical vibrational modes, which leads to the dominant role of electron scattering with optical phonons in charge transport. It is expected that these new atomistic-level insights offer a standing point where the performance of thermoelectric perovskites can be further enhanced.

In-silico identification of Tyr232 in AMPKα2 as a dephosphorylation site for the protein tyrosine phosphatase PTP-PEST.

The AMP-activated protein kinase (AMPK) is known to be activated by the protein tyrosine phosphatase non-receptor type 12 (PTP-PEST) under hypoxic conditions. This activation is mediated by tyrosine dephosphorylation of the AMPKα subunit. However, the identity of the phosphotyrosine residues that PTP-PEST dephosphorylates remains unknown. In this study, we first predicted the structure of the complex of the AMPKα2 subunit and PTP-PEST catalytic domain using bioinformatics tools and further confirmed the stability of the complex using molecular dynamics simulations. Evaluation of the protein-protein interfaces indicated that residue Tyr232 is the most likely dephosphorylation site on AMPKα2. In addition, we explored the effect of phosphorylation of PTP-PEST residue Tyr64 on the stability of the complex. Phosphorylation of the highly conserved Tyr64, an interface residue, enhances the stability of the complex via the rearrangement of a network of electrostatic interactions in conjunction with conformational changes in the catalytic WPD loop. We generated a phosphomimetic (PTP-PEST-Y64D) mutant and used co-immunoprecipitation to study the effect of PTP-PEST phosphorylation on AMPKα2 binding. The mutant exhibited an increased affinity for AMPKα2 and corroborated the in-silico predictions. Together, our findings present a plausible structural basis of AMPK regulation by PTP-PEST and show how phosphorylation of PTP-PEST affects its interaction with AMPKα2.

Evolutionary Conservation of PP2A Antagonism and G2/M Cell Cycle Arrest in Maedi-Visna Virus Vif.

The canonical function of lentiviral Vif proteins is to counteract the mutagenic potential of APOBEC3 antiviral restriction factors. However, recent studies have discovered that Vif proteins from diverse HIV-1 and simian immunodeficiency virus (SIV) isolates degrade cellular B56 phosphoregulators to remodel the host phosphoproteome and induce G2/M cell cycle arrest. Here, we evaluate the conservation of this activity among non-primate lentiviral Vif proteins using fluorescence-based degradation assays and demonstrate that maedi-visna virus (MVV) Vif efficiently degrades all five B56 family members. Testing an extensive panel of single amino acid substitution mutants revealed that MVV Vif recognizes B56 proteins through a conserved network of electrostatic interactions. Furthermore, experiments using genetic and pharmacologic approaches demonstrate that degradation of B56 proteins requires the cellular cofactor cyclophilin A. Lastly, MVV Vif-mediated depletion of B56 proteins induces a potent G2/M cell cycle arrest phenotype. Therefore, remodeling of the cellular phosphoproteome and induction of G2/M cell cycle arrest are ancient and conserved functions of lentiviral Vif proteins, which suggests that they are advantageous for lentiviral pathogenesis.

Light- and pH-dependent structural changes in cyanobacteriochrome AnPixJg2

Cyanobacteriochromes (CBCRs) are phytochrome-related photosensory proteins that play an essential role in regulating phototaxis, chromatic acclimation, and cell aggregation in cyanobacteria. Here, we apply solid-state NMR spectroscopy to the red/green GAF2 domain of the CBCR AnPixJ assembled in vitro with a uniformly 13C- and 15N-labeled bilin chromophore, tracking changes in electronic structure, geometry, and structural heterogeneity of the chromophore as well as intimate contacts between the chromophore and protein residues in the photocycle. Our data confirm that the bilin ring D is strongly twisted with respect to the B–C plane in both dark and photoproduct states. We also identify a greater structural heterogeneity of the bilin chromophore in the photoproduct than in the dark state. In addition, the binding pocket is more hydrated in the photoproduct. Observation of interfacial 1H contacts of the photoproduct chromophore, together with quantum mechanics/molecular mechanics (QM/MM)-based structural models for this photoproduct, clearly suggests the presence of a biprotonated (cationic) imidazolium side-chain for a conserved histidine residue (322) at a distance of ~2.7 Å, generalizing the recent theoretical findings that explicitly link the structural heterogeneity of the dark-state chromophore to the protonation of this specific residue. Moreover, we examine pH effects on this in vitro assembled holoprotein, showing a substantially altered electronic structure and protonation of the photoproduct chromophore even with a small pH drop from 7.8 to 7.2. Our studies provide further information regarding the light- and pH-induced changes of the chromophore and the rearrangements of the hydrogen-bonding and electrostatic interaction network around it. Possible correlations between structural heterogeneity of the chromophore, protonation of the histidine residue nearby, and hydration of the pocket in both photostates are discussed.Graphical abstract

Ab initio molecular dynamics simulations on the adsorption of 1-hydroxyethane-1,1-diphosphonic acid on the iron (100) surface

Through IRI, Bader charges and ELF analysis, an electrostatic interaction network exists between HEDP and the iron (100) surface, which is responsible for the stability of the adsorption state.

P53 Transactivation Domain Mediates Binding and Phase Separation with Poly-PR/GR

The most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) is the presence of poly-PR/GR dipeptide repeats, which are encoded by the chromosome 9 open reading frame 72 (C9orf72) gene. Recently, it was shown that poly-PR/GR alters chromatin accessibility, which results in the stabilization and enhancement of transcriptional activity of the tumor suppressor p53 in several neurodegenerative disease models. A reduction in p53 protein levels protects against poly-PR and partially against poly-GR neurotoxicity in cells. Moreover, in model organisms, a reduction of p53 protein levels protects against neurotoxicity of poly-PR. Here, we aimed to study the detailed molecular mechanisms of how p53 contributes to poly-PR/GR-mediated neurodegeneration. Using a combination of biophysical techniques such as nuclear magnetic resonance (NMR) spectroscopy, fluorescence polarization, turbidity assays, and differential interference contrast (DIC) microscopy, we found that p53 physically interacts with poly-PR/GR and triggers liquid–liquid phase separation of p53. We identified the p53 transactivation domain 2 (TAD2) as the main binding site for PR25/GR25 and showed that binding of poly-PR/GR to p53 is mediated by a network of electrostatic and/or hydrophobic interactions. Our findings might help to understand the mechanistic role of p53 in poly-PR/GR-associated neurodegeneration.

Solution Structure and Conformational Flexibility of a Polyketide Synthase Module.

Polyketide synthases (PKSs) are versatile C–C bond-forming enzymes that are broadly distributed in bacteria and fungi. The polyketide compound family includes many clinically useful drugs such as the antibiotic erythromycin, the antineoplastic epothilone, and the cholesterol-lowering lovastatin. Harnessing PKSs for custom compound synthesis remains an open challenge, largely because of the lack of knowledge about key structural properties. Particularly, the domains—well characterized on their own—are poorly understood in their arrangement, conformational dynamics, and interplay in the intricate quaternary structure of modular PKSs. Here, we characterize module 2 from the 6-deoxyerythronolide B synthase by small-angle X-ray scattering and cross-linking mass spectrometry with coarse-grained structural modeling. The results of this hybrid approach shed light on the solution structure of a cis-AT type PKS module as well as its inherent conformational dynamics. Supported by a directed evolution approach, we also find that acyl carrier protein (ACP)-mediated substrate shuttling appears to be steered by a nonspecific electrostatic interaction network.

Design of Bacillus fastidious Uricase Mutants Bearing Long Lagging Phases Before Exponential Decreases of Activities Under Physiological Conditions.

Under physiological conditions, Bacillus fastidious uricase (BFU) activity shows negligible lagging phase before the exponential decrease; mutants are thus designed for long lagging phases before exponential activity decreases. On homodimer surface of BFU (4R8X.pdb), the last fragment ANSEYVAL at the C-terminus forms a loop whose Y319 is H-bonded by the buried D257 in the same monomer. Within 1.5nm from the α-carboxyl group of the last leucine (L322), E30, K26, D257, R258, E311, K312 and E318 from the same monomer plus D126 and K127 from a monomer of the other homodimer generate an electrostatic interaction network. Within 1.5nm from Y319, D307 and R310 in the same monomer interact with ionized residues around the inter-chain β-sheet in the same homodimer. Mutagenesis of Y319R is designed to strengthen the original interactions and concomitantly generate new electrostatic attractions between homodimers. Under physiological conditions, the mutant V144A/Y319R showed an approximately 4week lagging phase before the exponential activity decrease, an apparent half-life of activity nearly three folds of mutant V144A, but comparable activity. The introduction of ionizable residues into the C-terminus contacting the other homodimer for additional and/or stronger electrostatic attractions between homodimers may be a universal approach to thermostable mutants of uricases.

Critical Interactions Between the SARS-CoV-2 Spike Glycoprotein and the Human ACE2 Receptor.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects human cells by binding its spike (S) glycoproteins to angiotensin-converting enzyme 2 (ACE2) receptors and causes the coronavirus disease 2019 (COVID-19). Therapeutic approaches to prevent SARS-CoV-2 infection are mostly focused on blocking S-ACE2 binding, but critical residues that stabilize this interaction are not well understood. By performing all-atom molecular dynamics (MD) simulations, we identified an extended network of salt bridges, hydrophobic and electrostatic interactions, and hydrogen bonds between the receptor-binding domain (RBD) of the S protein and ACE2. Mutagenesis of these residues on the RBD was not sufficient to destabilize binding but reduced the average work to unbind the S protein from ACE2. In particular, the hydrophobic end of RBD serves as the main anchor site and is the last to unbind from ACE2 under force. We propose that blocking the hydrophobic surface of RBD via neutralizing antibodies could prove to be an effective strategy to inhibit S-ACE2 interactions.

Synthesis, in vitro cytotoxicity, and molecular docking study of novel 3,4‐dihydroisoquinolin‐1(2H)‐one based piperlongumine analogues

With the aim of expanding the scope of SAR on piperlongumine (PL), a naturally occurring anticancer molecule, we have designed a novel hybrid molecule bearing 3,4‐dihydroisoquinolin‐1(2H)‐one and trans‐cinnamic acids. The structure, based on hybridization strategy, is used for hybridization of naturally occurring scaffolds. We have synthesized 14 hybrid molecules by coupling 3,4‐dihydroisoquinolin‐1(2H)‐one core with cinnamic acids using the mix anhydride approach. The newly synthesized inhibitors were evaluated for cell viability against breast cancer MCF‐7 and cervical cancer HeLa cell lines. Furthermore, the active compounds were screened for their potential in breast cancer MDA‐MB‐231, cervical cancer C33A cell lines, prostate cancer DU‐145, PC‐3, and normal VERO cells. From the series, compound 10g was seen to inhibit MCF‐7 cell growth significantly with GI50 < 0.1 μM along with growth inhibition in MDA‐MB‐231 (GI50 = 20 μM) and C33A (GI50 = 3.2 μM). While the inhibitor 10i inhibits MCF‐7 breast cancer cell growth GI50 = 3.42 μM along with inhibition of cell growth in MDA‐MB‐231 (GI50 = 30 μM), HeLa (GI50 = 7.67 μM), C33A (GI50 = 13 μM), DU‐145 (GI50 = 6.45 μM), PC‐3 (GI50 = 8.68 μM), and VERO (GI50 = 2.93 μM), respectively. Furthermore, molecular docking study demonstrated these compounds could bind tightly to the colchicine domain of tubulin through a network of favorable steric and electrostatic interactions and thus act as a tubulin polymerization inhibitor.

A new good thermal stability hybrid material based on heptamolybdate cluster: Synthesis and structural characterization

Abstract A new microporous polyoxomolybdate, (C4H7N2)9(NH4)3[Mo7O24]2•4H2O ([C4H7N2]+ = 2-methylimidazolium), was prepared at room temperature by reacting ammonium heptamolybdate tetrahydrate and 2-methylimidazole in a 1:4 molar ratio in water and characterized by infrared, TGA-DTA analysis and single-crystal X-ray diffraction. The compound is a heptamolybdate anion which crystallizes as a mixed salt of both ammonium [NH4]+ and 2-methylimidazolium [C4H7N2]+. The crystal structure of the compound is composed of [Mo7O24]6– clusters connected in a tridimensional network by [C4H7N2]+ and [NH4]+ cations and water molecules via hydrogen-bonding interactions. The relatively high thermal stability may be due to the existence of the network of hydrogen bonds and electrostatic interactions. Obtaining of compound from (NH4)6[Mo7O24]•4H2O shows that the transition from one polyoxomolybdate to another is possible through substitution reactions at room temperature, despite their stability. This material obtained could be used as catalyst and catalyst support.

Parkinson’s disease-related phosphorylation at Tyr39 rearranges α-synuclein amyloid fibril structure revealed by cryo-EM

Posttranslational modifications (PTMs) of α-synuclein (α-syn), e.g., phosphorylation, play an important role in modulating α-syn pathology in Parkinson's disease (PD) and α-synucleinopathies. Accumulation of phosphorylated α-syn fibrils in Lewy bodies and Lewy neurites is the histological hallmark of these diseases. However, it is unclear how phosphorylation relates to α-syn pathology. Here, by combining chemical synthesis and bacterial expression, we obtained homogeneous α-syn fibrils with site-specific phosphorylation at Y39, which exhibits enhanced neuronal pathology in rat primary cortical neurons. We determined the cryo-electron microscopy (cryo-EM) structure of the pY39 α-syn fibril, which reveals a fold of α-syn with pY39 in the center of the fibril core forming an electrostatic interaction network with eight charged residues in the N-terminal region of α-syn. This structure composed of residues 1 to 100 represents the largest α-syn fibril core determined so far. This work provides structural understanding on the pathology of the pY39 α-syn fibril and highlights the importance of PTMs in defining the polymorphism and pathology of amyloid fibrils in neurodegenerative diseases.

Characterization of hydroxymethylpyrimidine phosphate kinase from mesophilic and thermophilic bacteria and structural insights into their differential thermal stability

The hydroxymethylpyrimidine phosphate kinases (HMPPK) encoded by the thiD gene are involved in the thiamine biosynthesis pathway, can perform two consecutive phosphorylations of 4-amino-5-hydroxymethyl-2-methyl pyrimidine (HMP) and are found in thermophilic and mesophilic bacteria, but only a few characterizations of mesophilic enzymes are available. The presence of another homolog enzyme (pyridoxal kinase) that can only catalyze the first phosphorylation of HMP and encoded by pdxK gene, has hampered a precise annotation in this enzyme family. Here we report the kinetic characterization of two HMPPK with structure available, the mesophilic and thermophilic enzyme from Salmonella typhimurium (StHMPPK) and Thermus thermophilus (TtHMPPK), respectively. Also, given their high structural similarity, we have analyzed the structural determinants of protein thermal stability in these enzymes by molecular dynamics simulation. The results show that pyridoxal kinases (PLK) from gram-positive bacteria (PLK/HMPPK-like enzymes) constitute a phylogenetically separate group from the canonical PLK, but closely related to the HMPPK, so the PLK/HMPPK-like and canonical PLK, both encoded by pdxK genes, are different and must be annotated distinctly. The kinetic characterization of StHMPPK and TtHMPPK, shows that they perform double phosphorylation on HMP, both enzymes are specific for HMP, not using pyridoxal-like molecules as substrates and their kinetic mechanism involves the formation of a ternary complex. Molecular dynamics simulation shows that StHMPPK and TtHMPPK have striking differences in their conformational flexibility, which can be correlated with the hydrophobic packing and electrostatic interaction network given mainly by salt bridge bonds, but interestingly not by the number of hydrogen bond interactions as reported for other thermophilic enzymes. EnzymesEC 2.7.1.49, EC 2.7.4.7, EC 2.7.1.35, EC 2.7.1.50

Hierarchical Self-Assembly of Proteins Through Rationally Designed Supramolecular Interfaces.

With the increasing advances in the basic understanding of pathogenesis mechanism and fabrication of advanced biological materials, protein nanomaterials are being developed for their potential bioengineering research and biomedical applications. Among different fabrication strategies, supramolecular self-assembly provides a versatile approach to construct hierarchical nanostructures from polyhedral cages, filaments, tubules, monolayer sheets to even cubic crystals through rationally designed supramolecular interfaces. In this mini review, we will briefly recall recent progress in reconstituting protein interfaces for hierarchical self-assembly and classify by the types of designed protein-protein interactions into receptor-ligand recognition, electrostatic interaction, metal coordination, and non-specific interaction networks. Moreover, some attempts on functionalization of protein superstructures for bioengineering and/or biomedical applications are also shortly discussed. We believe this mini review will outline the stream of hierarchical self-assembly of proteins through rationally designed supramolecular interfaces, which would open minds in visualizing protein-protein recognition and assembly in living cells and organisms, and even constructing multifarious functional bionanomaterials.

Origins of the Increased Affinity of Phosphorothioate-Modified Therapeutic Nucleic Acids for Proteins.

The phosphorothioate backbone modification (PS) is one of the most widely used chemical modifications for enhancing the drug-like properties of nucleic acid-based drugs, including antisense oligonucleotides (ASOs). PS-modified nucleic acid therapeutics show improved metabolic stability from nuclease-mediated degradation and exhibit enhanced interactions with plasma, cell-surface, and intracellular proteins, which facilitates their tissue distribution and cellular uptake in animals. However, little is known about the structural basis of the interactions of PS nucleic acids with proteins. Here, we report a crystal structure of the DNA-binding domain of a model ASO-binding protein PC4, in complex with a full PS 2'-OMe DNA gapmer ASO. To our knowledge this is the first structure of a complex between a protein and fully PS nucleic acid. Each PC4 dimer comprises two DNA-binding interfaces. In the structure one interface binds the 5'-terminal 2'-OMe PS flank of the ASO, while the other interface binds the regular PS DNA central part in the opposite polarity. As a result, the ASO forms a hairpin-like structure. ASO binding also induces the formation of a dimer of dimers of PC4, which is stabilized by base pairing between homologous regions of the ASOs bound by each dimer of PC4. The protein interacts with the PS nucleic acid through a network of electrostatic and hydrophobic interactions, which provides insights into the origins for the enhanced affinity of PS for proteins. The importance of these contacts was further confirmed in a NanoBRET binding assay using a Nano luciferase tagged PC4 acting as the BRET donor, to a fluorescently conjugated ASO acting as the BRET acceptor. Overall, our results provide insights into the molecular forces that govern the interactions of PS ASOs with cellular proteins and provide a potential model for how these interactions can template protein-protein interactions causative of cellular toxicity.

A role of Gln596 in fine-tuning mammalian ALOX15 specificity, protein stability and allosteric properties

His596 of human ALOX12 has been suggested to interact with the COO−-group of arachidonic acid during ALOX catalysis. In mammalian ALOX15 orthologs Gln596 occupies this position and this amino acid exchange might contribute to the functional differences between the two ALOX-isoforms. To explore the role of Gln596 for ALOX15 functionality we mutated this amino acid to different residues in rabbit and human ALOX15 and investigated the impact of these mutations on structural, catalytic and allosteric enzyme properties. To shed light on the molecular basis of the observed functional alterations we performed in silico substrate docking studies and molecular dynamics simulations and also explored the impact of Gln596 exchange on the protein structure. The combined theoretical and experimental data suggest that Gln596 may not directly interact with the COO−-group of arachidonic acid. In contrast, mutations at Gln596 destabilize the secondary and tertiary structure of ALOX15 orthologs, which may be related to a disturbance of the electrostatic interaction network with other amino acids in the immediate surrounding. Moreover, our MD-simulations suggest that the geometry of the dimer interface depends on the structure of substrate bound inside the substrate-binding pocket and that Gln596Ala exchange impairs the allosteric properties of the enzyme. Taken together, these data indicate the structural and functional importance of Gln596 for ALOX15 catalysis.