Lipid composition may significantly affect membrane proteins function, yet its impact on the protein structural determinants is not well understood. Here we present a comparative molecular dynamics (MD) study of the human adenosine receptor type 2A (hA(2A)R) in complex with caffeine--a system of high neuro-pharmacological relevance--within different membrane types. These are POPC, mixed POPC/POPE and cholesterol-rich membranes. 0.8-μs MD simulations unambiguously show that the helical folding of the amphipathic helix 8 depends on membrane contents. Most importantly, the distinct cholesterol binding into the cleft between helix 1 and 2 stabilizes a specific caffeine-binding pose against others visited during the simulation. Hence, cholesterol presence (~33%-50% in synaptic membrane in central nervous system), often neglected in X-ray determination of membrane proteins, affects the population of the ligand binding poses. We conclude that including a correct description of neuronal membranes may be very important for computer-aided design of ligands targeting hA(2A)R and possibly other GPCRs.

Sirtuin 2 (SIRT2) is a NAD+-dependent deacetylase that has been associated with neurodegeneration and cancer. SIRT2 is composed of a central catalytic domain, the structure of which has been solved, and N- and C-terminal extensions that are thought to control SIRT2 function. However structural information of these N- and C-terminal regions is missing. Here, we provide the first full-length molecular models of SIRT2 in the absence and presence of NAD+. We also predict the structural alterations associated with phosphorylation of SIRT2 at S331, a modification that inhibits catalytic activity. Bioinformatics tools and molecular dynamics simulations, complemented by in vitro deacetylation assays, provide a consistent picture based on which the C-terminal region of SIRT2 is suggested to function as an autoinhibitory region. This has the capacity to partially occlude the NAD+ binding pocket or stabilize the NAD+ in a non-productive state. Furthermore, our simulations suggest that the phosphorylation at S331 causes large conformational changes in the C-terminal region that enhance the autoinhibitory activity, consistent with our previous findings that phosphorylation of S331 by cyclin-dependent kinases inhibits SIRT2 catalytic activity. The molecular insight into the role of the C-terminal region in controlling SIRT2 function described in this study may be useful for future design of selective inhibitors targeting SIRT2 for therapeutic applications.

Macrodomains are conserved protein folds associated with ADP-ribose binding and turnover. ADP-ribosylation is a posttranslational modification catalyzed primarily by ARTD (aka PARP) enzymes in cells. ARTDs transfer either single or multiple ADP-ribose units to substrates, resulting in mono- or poly-ADP-ribosylation. TARG1/C6orf130 is a macrodomain protein that hydrolyzes mono-ADP-ribosylation and interacts with poly-ADP-ribose chains. Interactome analyses revealed that TARG1 binds strongly to ribosomes and proteins associated with rRNA processing and ribosomal assembly factors. TARG1 localized to transcriptionally active nucleoli, which occurred independently of ADP-ribose binding. TARG1 shuttled continuously between nucleoli and nucleoplasm. In response to DNA damage, which activates ARTD1/2 (PARP1/2) and promotes synthesis of poly-ADP-ribose chains, TARG1 re-localized to the nucleoplasm. This was dependent on the ability of TARG1 to bind to poly-ADP-ribose. These findings are consistent with the observed ability of TARG1 to competitively interact with RNA and PAR chains. We propose a nucleolar role of TARG1 in ribosome assembly or quality control that is stalled when TARG1 is re-located to sites of DNA damage.

Positron emission tomography (PET) radioligands targeting the human translocator membrane protein (TSPO) are broadly used for the investigations of neuroinflammatory conditions associated with neurological disorders. Structural information on the mammalian protein homodimers-the suggested functional state of the protein-is limited to a solid-state nuclear magnetic resonance (NMR) study and to a model based on the previously-deposited solution NMR structure of the monomeric mouse protein. Computational studies performed here suggest that the NMR-solved structure in the presence of detergents is not prone to dimer formation and is furthermore unstable in its native membrane environment. We, therefore, propose a new model of the functionally-relevant dimeric form of the mouse protein, based on a prokaryotic homologue. The model, fully consistent with solid-state NMR data, is very different from the previous predictions. Hence, it provides, for the first time, structural insights into this pharmaceutically-important target which are fully consistent with experimental data.

Advances in coarse-grained molecular dynamics (CGMD) simulations have extended the use of computational studies on biological macromolecules and their complexes, as well as the interactions of membrane protein and lipid complexes at a reduced level of representation, allowing longer and larger molecular dynamics simulations. Here, we present a computational platform dedicated to the preparation, running, and analysis of CGMD simulations. The platform is built on a completely revisited version of our Martini coarsE gRained MembrAne proteIn Dynamics (MERMAID) web server, and it integrates this with other three dedicated services. In its current version, the platform expands the existing implementation of the Martini force field for membrane proteins to also allow the simulation of soluble proteins using the Martini and the SIRAH force fields. Moreover, it offers an automated protocol for carrying out the backmapping of the coarse-grained description of the system into an atomistic one.

Mitochondrial energy metabolism plays an important role in the pathophysiology of insulin resistance. Recently, a missense N437S variant was identified in the MRPP3 gene, which encodes a mitochondrial RNA processing enzyme within the RNase P complex, with predicted impact on metabolism. We used CRISPR-Cas9 genome editing to introduce this variant into the mouse Mrpp3 gene and show that the variant causes insulin resistance on a high-fat diet. The variant did not influence mitochondrial gene expression markedly, but instead, it reduced mitochondrial calcium that lowered insulin release from the pancreatic islet β cells of the Mrpp3 variant mice. Reduced insulin secretion resulted in lower insulin levels that contributed to imbalanced metabolism and liver steatosis in the Mrpp3 variant mice on a high-fat diet. Our findings reveal that the MRPP3 variant may be a predisposing factor to insulin resistance and metabolic disease in the human population.

The NEET proteins constitute a unique class of [2Fe-2S] proteins. The metal ions bind to three cysteines and one histidine. The proteins' clusters exist in two redox states; the oxidized protein (containing two FeIII ions) can transfer the cluster to apo-acceptor protein(s), while the reduced form (containing one ferrous ion) remains bound to the protein frame. Here, we perform in silico and in vitro studies on human NEET proteins in both reduced and oxidized forms. Quantum chemical calculations on all available human NEET proteins structures suggest that reducing the cluster weakens the Fe-NHis and Fe-SCys bonds, similar to what is seen in other Fe-S proteins (e.g., ferredoxin and Rieske protein). We further show that the extra electron in the [2Fe-2S]+ clusters of one of the NEET proteins (mNT) is localized on the His-bound iron ion, consistently with our previous spectroscopic studies. Kinetic measurements demonstrate that the mNT [2Fe-2S]+ is released only by an increase in temperature. Thus, the reduced state of human NEET proteins [2Fe-2S] cluster is kinetically inert. This previously unrecognized kinetic inertness of the reduced state, along with the reactivity of the oxidized state, is unique across all [2Fe-2S] proteins. Finally, using a coevolutionary analysis, along with molecular dynamics simulations, we provide insight on the observed allostery between the loop L2 and the cluster region. Specifically, we show that W75, R76, K78, K79, F82 and G85 in the latter region share similar allosteric characteristics in both redox states.

Cisplatin is one of the most widely used anticancer drugs. Its efficiency is unfortunately severely hampered by resistance. The High Mobility Group Box (HMGB) proteins may sensitize tumor cells to cisplatin by specifically binding to platinated DNA (PtDNA) lesions. In vivo, the HMGB/PtDNA binding is regulated by multisite post-translational modifications (PTMs). The impact of PTMs on the HMGB/PtDNA complex at atomistic level is here investigated by enhanced sampling molecular simulations. The PTMs turn out to affect the structure of the complex, the mobility of several regions (including the platinated site), and the nature of the protein/PtDNA non-covalent interactions. Overall, the multisite PTMs increase significantly the apparent synchrony of all the contacts between the protein and PtDNA. Consequently, the hydrophobic anchoring of the side chain of F37 between the two cross-linked guanines at the platinated site-a key element of the complexes formation - is more stable than in the complex without PTM. These differences can account for the experimentally measured greater affinity for PtDNA of the protein isoforms with PTMs. The collective behavior of multisite PTMs, as revealed here by the synchrony of contacts, may have a general significance for the modulation of intermolecular recognitions occurring in vivo.

The cyclin-dependent kinase 2 (CDK2) together with its cyclin E and A partners is a central regulator of cell growth and division. Deregulation of CDK2 activity is associated with diseases such as cancer. The analysis of substrates identified S/T-P-X-R/K/H as the CDK2 consensus sequence. The crystal structure of cyclin A/CDK2 with a short model peptide supports this sequence and identifies key interactions. However, CDKs use additional determinants to recognize substrates, including the RXL motif that is read by the cyclin subunits. We were interested to determine whether additional amino acids beyond the minimal consensus sequence of the well-studied substrate and tumor suppressor p27KIP1 were relevant for catalysis.

Human G-protein coupled receptors (GPCRs) are important targets for pharmaceutical intervention against neurological diseases. Here, we use molecular simulation to investigate the key step in ligand recognition governed by the extracellular domains in the neuronal adenosine receptor type 2A (hA2AR), a target for neuroprotective compounds. The ligand is the high-affinity antagonist (4-(2-(7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a][1,3,5]triazin-5-ylamino)ethyl)phenol), embedded in a neuronal membrane mimic environment. Free energy calculations, based on well-tempered metadynamics, reproduce the experimentally measured binding affinity. The results are consistent with the available mutagenesis studies. The calculations identify a vestibular binding site, where lipids molecules can actively participate to stabilize ligand binding. Bioinformatic analyses suggest that such vestibular binding site and, in particular, the second extracellular loop, might drive the ligand toward the orthosteric binding pocket, possibly by allosteric modulation. Taken together, these findings point to a fundamental role of the interaction between extracellular loops and membrane lipids for ligands' molecular recognition and ligand design in hA2AR.

Prion diseases or transmissible spongiform encephalopathies (TSEs) are a group of rare neurodegenerative disorders. TSEs are characterized by the accumulation of prions (PrPSc) that represent pathological isoforms of the physiological cellular prion protein PrPC. Although the conversion of PrPC to PrPSc is still not completely understood, blocking this process may lead to develop new therapies. Here, we have generated a pharmacophore model, based on anti-prion molecules reported in literature to be effective in phenotypic assay. The model was used to conduct a virtual screen of commercial compound databases that selected a small library of ten compounds. These molecules were then screened in mouse neuroblastoma cell line chronically infected with prions (ScN2a) after excluding neurotoxicity. 1 has been identified as the therapeutic hit on the basis of the following evidence: chronic treatments of ScN2a cells using 1 eliminate PrPSc loaded in both Western blotting analysis and Real-Time Quaking-Induced Conversion (RT-QuIC) assay. We also proposed the mechanism of action of 1 by which it has the ability to bind PrPC and consequentially blocks prion conversion. Herein we describe the results of these efforts.

The influence of environmental insults on the onset and progression of mitochondrial diseases is unknown. To evaluate the effects of infection on mitochondrial disease we used a mouse model of Leigh Syndrome, where a missense mutation in the Taco1 gene results in the loss of the translation activator of cytochrome c oxidase subunit I (TACO1) protein. The mutation leads to an isolated complex IV deficiency that mimics the disease pathology observed in human patients with TACO1 mutations. We infected Taco1 mutant and wild-type mice with a murine cytomegalovirus and show that a common viral infection exacerbates the complex IV deficiency in a tissue-specific manner. We identified changes in neuromuscular morphology and tissue-specific regulation of the mammalian target of rapamycin pathway in response to viral infection. Taken together, we report for the first time that a common stress condition, such as viral infection, can exacerbate mitochondrial dysfunction in a genetic model of mitochondrial disease.

Many homodimeric enzymes tune their functions by exploiting either negative or positive cooperativity between subunits. In the SARS-CoV-2 Main protease (Mpro) homodimer, the latter has been suggested by symmetry in most of the 500 reported protease/ligand complex structures solved by macromolecular crystallography (MX). Here we apply the latter to both covalent and noncovalent ligands in complex with Mpro. Strikingly, our experiments show that the occupation of both active sites of the dimer originates from an excess of ligands. Indeed, cocrystals obtained using a 1:1 ligand/protomer stoichiometry lead to single occupation only. The empty binding site exhibits a catalytically inactive geometry in solution, as suggested by molecular dynamics simulations. Thus, Mpro operates through negative cooperativity with the asymmetric activity of the catalytic sites. This allows it to function with a wide range of substrate concentrations, making it resistant to saturation and potentially difficult to shut down, all properties advantageous for the virus' adaptability and resistance.

We provide a molecular-level description of the thermodynamics and mechanistic aspects of drug permeation through the cell membrane. As a case study, we considered the antimalaria FDA approved drug chloroquine. Molecular dynamics simulations of the molecule (in its neutral and protonated form) were performed in the presence of different lipid bilayers, with the aim of uncovering key aspects of the permeation process, a fundamental step for the drug's action. Free energy values obtained by well-tempered metadynamics simulations suggest that the neutral form is the only permeating protomer, consistent with experimental data. H-bond interactions of the drug with water molecules and membrane headgroups play a crucial role for permeation. The presence of the transmembrane potential, investigated here for the first time in a drug permeation study, does not qualitatively affect these conclusions.

Cell-cell communication is mediated by membrane receptors and their ligands, such as the Eph/ephrin system, orchestrating cell migration during development and in diverse cancer types. Epigenetic mechanisms are key for integrating external "signals", e.g., from neighboring cells, into the transcriptome in health and disease. Previously, we reported ephrinA5 to trigger transcriptional changes of lncRNAs and protein-coding genes in cerebellar granule cells, a cell model for medulloblastoma. LncRNAs represent important adaptors for epigenetic writers through which they regulate gene expression. Here, we investigate a lncRNA-mediated targeting of DNMT1 to specific gene loci by the combined power of in silico modeling of RNA/DNA interactions and wet lab approaches, in the context of the clinically relevant use case of ephrinA5-dependent regulation of cellular motility of cerebellar granule cells. We provide evidence that Snhg15, a cancer-related lncRNA, recruits DNMT1 to the Ncam1 promoter through RNA/DNA triplex structure formation and the interaction with DNMT1. This mediates DNA methylation-dependent silencing of Ncam1, being abolished by ephrinA5 stimulation-triggered reduction of Snhg15 expression. Hence, we here propose a triple helix recognition mechanism, underlying cell motility regulation via lncRNA-targeted DNA methylation in a clinically relevant context.

STAT5 is an essential transcription factor in hematopoiesis, which is activated through tyrosine phosphorylation in response to cytokine stimulation. Constitutive activation of STAT5 is a hallmark of myeloid and lymphoblastic leukemia. Using homology modeling and molecular dynamics simulations, a model of the STAT5 phosphotyrosine-SH2 domain interface was generated providing first structural information on the activated STAT5 dimer including a sequence, for which no structural information is available for any of the STAT proteins. We identified a novel intramolecular interaction mediated through F706, adjacent to the phosphotyrosine motif, and a unique hydrophobic interface on the surface of the SH2 domain. Analysis of corresponding STAT5 mutants revealed that this interaction is dispensable for Epo receptor-mediated phosphorylation of STAT5 but essential for dimer formation and subsequent nuclear accumulation. Moreover, the herein presented model clarifies molecular mechanisms of recently discovered leukemic STAT5 mutants and will help to guide future drug development.

Compound repurposing is an important strategy for the identification of effective treatment options against SARS-CoV-2 infection and COVID-19 disease. In this regard, SARS-CoV-2 main protease (3CL-Pro), also termed M-Pro, is an attractive drug target as it plays a central role in viral replication by processing the viral polyproteins pp1a and pp1ab at multiple distinct cleavage sites. We here report the results of a repurposing program involving 8.7 K compounds containing marketed drugs, clinical and preclinical candidates, and small molecules regarded as safe in humans. We confirmed previously reported inhibitors of 3CL-Pro and have identified 62 additional compounds with IC50 values below 1 μM and profiled their selectivity toward chymotrypsin and 3CL-Pro from the Middle East respiratory syndrome virus. A subset of eight inhibitors showed anticytopathic effect in a Vero-E6 cell line, and the compounds thioguanosine and MG-132 were analyzed for their predicted binding characteristics to SARS-CoV-2 3CL-Pro. The X-ray crystal structure of the complex of myricetin and SARS-Cov-2 3CL-Pro was solved at a resolution of 1.77 Å, showing that myricetin is covalently bound to the catalytic Cys145 and therefore inhibiting its enzymatic activity.

Recent studies revealed that molecular events related with the physiology and pathology of αS might be regulated by specific sequence motifs in the primary sequence of αS. The importance of individual residues in these motifs remains an important open avenue of investigation. In this work, we have addressed the structural details related to the amyloid fibril assembly and lipid-binding features of αS through the design of site-directed mutants at position 39 of the protein and their study by in vitro and in vivo assays. We demonstrated that aromaticity at position 39 of αS primary sequence influences strongly the aggregation properties and the membrane-bound conformations of the protein, molecular features that might have important repercussions for the function and dysfunction of αS. Considering that aggregation and membrane damage is an important driver of cellular toxicity in amyloid diseases, future work is needed to link our findings with studies based on toxicity and neuronal cell death. BRIEF STATEMENT OUTLINING SIGNIFICANCE: Modulation by distinct sequential motifs and specific residues of αS on its physiological and pathological states is an active area of research. Here, we demonstrated that aromaticity at position 39 of αS modulates the membrane-bound conformations of the protein, whereas removal of aromatic functionality at position 39 reduces strongly the amyloid assembly in vitro and in vivo. Our study provides new evidence for the modulation of molecular events related with the physiology and pathology of αS.

Glioblastomas strongly invade the brain by infiltrating into the white matter along myelinated nerve fiber tracts even though the myelin protein Nogo-A prevents cell migration by activating inhibitory RhoA signaling. The mechanisms behind this long-known phenomenon remained elusive so far, precluding a targeted therapeutic intervention. This study demonstrates that the prevalent activation of AKT in gliomas increases the ER protein-folding capacity and enables tumor cells to utilize a side effect of RhoA activation: the perturbation of the IRE1α-mediated decay of SPARC mRNA. Once translation is initiated, glioblastoma cells rapidly secrete SPARC to block Nogo-A from inhibiting migration via RhoA. By advanced ultramicroscopy for studying single-cell invasion in whole, undissected mouse brains, we show that gliomas require SPARC for invading into white matter structures. SPARC depletion reduces tumor dissemination that significantly prolongs survival and improves response to cytostatic therapy. Our finding of a novel RhoA-IRE1 axis provides a druggable target for interfering with SPARC production and underscores its therapeutic value.

Human NEET proteins contain two [2Fe-2S] iron-sulfur clusters, bound to three Cys residues and one His residue. They exist in two redox states. Recently, these proteins have revealed themselves as attractive drug targets for mitochondrial dysfunction-related diseases, such as type 2 diabetes, Wolfram syndrome 2, and cancers. Unfortunately, the lack of information and mechanistic understanding of ligands binding to the whole functional, cytoplasmatic domain has limited rational drug design approaches. Here, we use an enhanced sampling technique, volume-based metadynamics, recently developed by a team involving some of us, to predict the poses and affinity of the 2-benzamido-4-(1,2,3,4-tetrahydronaphthalen-2-yl)-thiophene-3-carboxylate ligand to the entire surface of the cytoplasmatic domain of the human NEET protein mitoNEET (mNT) in an aqueous solution. The calculations, based on the recently published X-ray structure of the complex, are consistent with the measured affinity. The calculated free energy landscape revealed that the ligand can bind in multiple sites and with poses other than the one found in the X-ray. This difference is likely to be caused by crystal packing effects that allow the ligand to interact with multiple adjacent NEET protein copies. Such extra contacts are of course absent in the solution; therefore, the X-ray pose is only transient in our calculations, where the binding free energy correlates with the number of contacts. We further evaluated how the reduction and protonation of the Fe-bound histidine, as well as temperature, can affect ligand binding. Both such modifications introduce the possibility for the ligand to bind in an area of the protein other than the one observed in the X-ray, with no or little impact on affinity. Overall, our study can provide insights on the molecular recognition mechanisms of ligand binding to mNT in different oxidative conditions, possibly helping rational drug design of NEET ligands.