Search

Abstract 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. Graphical Abstract

The ability to alter the genomes of living cells is key to understanding how genes influence the functions of organisms and will be critical to modify living systems for useful purposes. Here, the authors use computational design to discover Cas9 enzymes with increased activity.

Marjault et al. report two crystal structures of a drug molecule bound to NEET proteins mNT and NAF-1 and the effects of this molecule on diabetic mice, suggesting an approach for reducing mitochondrial iron and ROS levels in diabetic patients.

The RNA-binding protein human antigen R (HuR) regulates stability, translation, and nucleus-to-cytoplasm shuttling of its target mRNAs. This protein has been progressively recognized as a relevant therapeutic target for several pathologies, like cancer, neurodegeneration, as well as inflammation. Inhibitors of mRNA binding to HuR might thus be beneficial against a variety of diseases. Here, we present the rational identification of structurally novel HuR inhibitors. In particular, by combining chemoinformatic approaches, high-throughput virtual screening, and RNA–protein pulldown assays, we demonstrate that the 4-(2-(2,4,6-trioxotetrahydropyrimidin-5(2H)-ylidene)hydrazineyl)benzoate ligand exhibits a dose-dependent HuR inhibition effect in binding experiments. Importantly, the chemical scaffold is new with respect to the currently known HuR inhibitors, opening up a new avenue for the design of pharmaceutical agents targeting this important protein.

Polyglutamine (polyQ) diseases are characterized by an expansion of cytosine-adenine-guanine (CAG) trinucleotide repeats encoding for an uninterrupted prolonged polyQ tract. We previously identified TRMT2A as a strong modifier of polyQ-induced toxicity in an unbiased large-scale screen in Drosophila melanogaster. This work aimed at identifying and validating pharmacological TRMT2A inhibitors as treatment opportunities for polyQ diseases in humans. Computer-aided drug discovery was implemented to identify human TRMT2A inhibitors. Additionally, the crystal structure of one protein domain, the RNA recognition motif (RRM), was determined, and Biacore experiments with the RRM were performed. The identified molecules were validated for their potency to reduce polyQ aggregation and polyQ-induced cell death in human HEK293T cells and patient derived fibroblasts. Our work provides a first step towards pharmacological inhibition of this enzyme and indicates TRMT2A as a viable drug target for polyQ diseases.

The disordered nature of Intrinsically Disordered Proteins (IDPs) makes their structural ensembles particularly susceptible to changes in chemical environmental conditions, often leading to an alteration of their normal functions. A Radial Distribution Function (RDF) is considered a standard method for characterizing the chemical environment surrounding particles during atomistic simulations, commonly averaged over an entire or part of a trajectory. Given their high structural variability, such averaged information might not be reliable for IDPs. We introduce the Time-Resolved Radial Distribution Function (TRRDF), implemented in our open-source Python package SPEADI, which is able to characterize dynamic environments around IDPs. We use SPEADI to characterize the dynamic distribution of ions around the IDPs Alpha-Synuclein (AS) and Humanin (HN) from Molecular Dynamics (MD) simulations, and some of their selected mutants, showing that local ion–residue interactions play an important role in the structures and behaviors of IDPs.

Human NEET proteins, such as NAF-1 and mitoNEET, are homodimeric, redox iron-sulfur proteins characterized by triple cysteine and one histidine-coordinated [2Fe-2S] cluster. They exist in an oxidized and reduced state. Abnormal release of the cluster is implicated in a variety of diseases, including cancer and neurodegeneration. The computer-aided and structure-based design of ligands affecting cluster release is of paramount importance from a pharmaceutical perspective. Unfortunately, experimental structural information so far is limited to only one ligand/protein complex. This is the X-ray structure of furosemide bound to oxidized mitoNEET. Here we employ an enhanced sampling approach, Localized Volume-based Metadynamics, developed by some of us, to identify binding poses of furosemide to human mitoNEET protein in solution. The binding modes show a high variability within the same shallow binding pocket on the protein surface identified in the X-ray structure. Among the different binding conformations, one of them is in agreement with the crystal structure’s one. This conformation might have been overstabilized in the latter because of the presence of crystal packing interactions, absent in solution. The calculated binding affinity is compatible with experimental data. Our protocol can be used in a straightforward manner in drug design campaigns targeting this pharmaceutically important family of proteins.

The long history of local anesthetics (LAs) starts out in the late 19th century when the content of coca plant leaves was discovered to alleviate pain. Soon after, cocaine was established and headed off to an infamous career as a substance causing addiction. Today, LAs and related substances—in modified form—are indispensable in our clinical everyday life for pain relief during and after minor and major surgery, and dental practices. In this review, we elucidate on the interaction of modern LAs with their main target, the voltage-gated sodium channel (Navs), in the light of the recently published channel structures. Knowledge of the 3D interaction sites of the drug with the protein will allow to mechanistically substantiate the comprehensive data available on LA gating modification. In the 1970s it was suggested that LAs can enter the channel pore from the lipid phase, which was quite prospective at that time. Today we know from cryo-electron microscopy structures and mutagenesis experiments, that indeed Navs have side fenestrations facing the membrane, which are likely the entrance for LAs to induce tonic block. In this review, we will focus on the effects of LA binding on fast inactivation and use-dependent inhibition in the light of the proposed new allosteric mechanism of fast inactivation. We will elaborate on subtype and species specificity and provide insights into modelling approaches that will help identify the exact molecular binding orientation, access pathways and pharmacokinetics. With this comprehensive overview, we will provide new perspectives in the use of the drug, both clinically and as a tool for basic ion channel research.

Coarse-grained (CG) modelling with the Martini force field has come of age. By combining a variety of bead types and sizes with a new mapping approach, the newest version of the model is able to accurately simulate large biomolecular complexes at millisecond timescales. In this perspective, we discuss possible applications of the Martini 3 model in drug discovery and development pipelines and highlight areas for future development. Owing to its high simulation efficiency and extended chemical space, Martini 3 has great potential in the area of drug design and delivery. However, several aspects of the model should be improved before Martini 3 CG simulations can be routinely employed in academic and industrial settings. These include the development of automatic parameterisation protocols for a variety of molecule types, the improvement of backmapping procedures, the description of protein flexibility and the development of methodologies enabling efficient sampling. We illustrate our view with examples on key areas where Martini could give important contributions such as drugs targeting membrane proteins, cryptic pockets and protein–protein interactions and the development of soft drug delivery systems.

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.

Abstract 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.

Summary: The regulation of mitochondrial RNA life cycles and their roles in ribosome biogenesis and energy metabolism are not fully understood. We used CRISPR/Cas9 to generate heart- and skeletal-muscle-specific knockout mice of the pentatricopeptide repeat domain protein 1, PTCD1, and show that its loss leads to severe cardiomyopathy and premature death. Our detailed transcriptome-wide and functional analyses of these mice enabled us to identify the molecular role of PTCD1 as a 16S rRNA-binding protein essential for its stability, pseudouridylation, and correct biogenesis of the mitochondrial large ribosomal subunit. We show that impaired mitoribosome biogenesis can have retrograde signaling effects on nuclear gene expression through the transcriptional activation of the mTOR pathway and upregulation of cytoplasmic protein synthesis and pro-survival factors in the absence of mitochondrial translation. Taken together, our data show that impaired assembly of the mitoribosome exerts its consequences via differential regulation of mitochondrial and cytoplasmic protein synthesis. : Perks et al. engineered intron-exon boundaries using CRISPR/Cas9 to conditionally knock out Ptcd1 in mice. The RNA-binding protein PTCD1 is essential for heart function and regulates the stability and maturation of the 16S rRNA. PTCD1 is required for mitoribosome biogenesis, mitochondrial function, and coordinated nuclear transcription. Keywords: RNA, cardiomyopathy, regulatory RNAs, RNA-seq, mitochondrial ribosome, RNA-binding proteins, mitochondria, mitochondrial gene expression, ribosome biogenesis

After almost two years from its first evidence, the COVID-19 pandemic continues to afflict people worldwide, highlighting the need for multiple antiviral strategies. SARS-CoV-2 main protease (Mpro/3CLpro) is a recognized promising target for the development of effective drugs. Because single target inhibition might not be sufficient to block SARS-CoV-2 infection and replication, multi enzymatic-based therapies may provide a better strategy. Here we present a structural and biochemical characterization of the binding mode of MG-132 to both the main protease of SARS-CoV-2, and to the human Cathepsin-L, suggesting thus an interesting scaffold for the development of double-inhibitors. X-ray diffraction data show that MG-132 well fits into the Mpro active site, forming a covalent bond with Cys145 independently from reducing agents and crystallization conditions. Docking of MG-132 into Cathepsin-L well-matches with a covalent binding to the catalytic cysteine. Accordingly, MG-132 inhibits Cathepsin-L with nanomolar potency and reversibly inhibits Mpro with micromolar potency, but with a prolonged residency time. We compared the apo and MG-132-inhibited structures of Mpro solved in different space groups and we identified a new apo structure that features several similarities with the inhibited ones, offering interesting perspectives for future drug design and in silico efforts.

Aberrant RNA–protein complexes are formed in a variety of diseases. Identifying the ligands that interfere with their formation is a valuable therapeutic strategy. Molecular simulation, validated against experimental data, has recently emerged as a powerful tool to predict both the pose and energetics of such ligands. Thus, the use of molecular simulation may provide insight into aberrant molecular interactions in diseases and, from a drug design perspective, may allow for the employment of less wet lab resources than traditional in vitro compound screening approaches. With regard to basic research questions, molecular simulation can support the understanding of the exact molecular interaction and binding mode. Here, we focus on examples targeting RNA–protein complexes in neurodegenerative diseases and viral infections. These examples illustrate that the strategy is rather general and could be applied to different pharmacologically relevant approaches. We close this study by outlining one of these approaches, namely the light-controllable association of small molecules with RNA, as an emerging approach in RNA-targeting therapy.

PIN1 is a promising therapeutic target for cancer treatment. In this study, the authors identify a covalent inhibitor of PIN1 with anti-tumour and anti-metastatic properties thanks to PIN1 inactivation and to the release, after binding to PIN1, of a quinone-mimicking compound that elicits reactive oxygen generation and causes DNA damage.

Flavonoids are plant bioactives that are recognized as hormone-like polyphenols because of their similarity to the endogenous sex steroids 17β-estradiol and testosterone, and to their estrogen- and androgen-like activity. Most efforts to verify flavonoid binding to nuclear receptors (NRs) and explain their action have been focused on ERα, while less attention has been paid to other nuclear and non-nuclear membrane androgen and estrogen receptors. Here, we investigate six flavonoids (apigenin, genistein, luteolin, naringenin, quercetin, and resveratrol) that are widely present in fruits and vegetables, and often used as replacement therapy in menopause. We performed comparative computational docking simulations to predict their capability of binding nuclear receptors ERα, ERβ, ERRβ, ERRγ, androgen receptor (AR), and its variant ART877A and membrane receptors for androgens, i.e., ZIP9, GPRC6A, OXER1, TRPM8, and estrogens, i.e., G Protein-Coupled Estrogen Receptor (GPER). In agreement with data reported in literature, our results suggest that these flavonoids show a relevant degree of complementarity with both estrogen and androgen NR binding sites, likely triggering genomic-mediated effects. It is noteworthy that reliable protein–ligand complexes and estimated interaction energies were also obtained for some suggested estrogen and androgen membrane receptors, indicating that flavonoids could also exert non-genomic actions. Further investigations are needed to clarify flavonoid multiple genomic and non-genomic effects. Caution in their administration could be necessary, until the safe assumption of these natural molecules that are largely present in food is assured.

The 3CL-Protease appears to be a very promising medicinal target to develop anti-SARS-CoV-2 agents. The availability of resolved structures allows structure-based computational approaches to be carried out even though the lack of known inhibitors prevents a proper validation of the performed simulations. The innovative idea of the study is to exploit known inhibitors of SARS-CoV 3CL-Pro as a training set to perform and validate multiple virtual screening campaigns. Docking simulations using four different programs (Fred, Glide, LiGen, and PLANTS) were performed investigating the role of both multiple binding modes (by binding space) and multiple isomers/states (by developing the corresponding isomeric space). The computed docking scores were used to develop consensus models, which allow an in-depth comparison of the resulting performances. On average, the reached performances revealed the different sensitivity to isomeric differences and multiple binding modes between the four docking engines. In detail, Glide and LiGen are the tools that best benefit from isomeric and binding space, respectively, while Fred is the most insensitive program. The obtained results emphasize the fruitful role of combining various docking tools to optimize the predictive performances. Taken together, the performed simulations allowed the rational development of highly performing virtual screening workflows, which could be further optimized by considering different 3CL-Pro structures and, more importantly, by including true SARS-CoV-2 3CL-Pro inhibitors (as learning set) when available.

The translocator protein (TSPO) is a 18kDa transmembrane protein, ubiquitously present in human mitochondria. It is overexpressed in tumor cells and at the sites of neuroinflammation, thus representing an important biomarker, as well as a promising drug target. In mammalian TSPO, there are cholesterol–binding motifs, as well as a binding cavity able to accommodate different chemical compounds. Given the lack of structural information for the human protein, we built a model of human (h) TSPO in the apo state and in complex with PK11195, a molecule routinely used in positron emission tomography (PET) for imaging of neuroinflammatory sites. To better understand the interactions of PK11195 and cholesterol with this pharmacologically relevant protein, we ran molecular dynamics simulations of the apo and holo proteins embedded in a model membrane. We found that: (i) PK11195 stabilizes hTSPO structural fold; (ii) PK11195 might enter in the binding site through transmembrane helices I and II of hTSPO; (iii) PK11195 reduces the frequency of cholesterol binding to the lower, N–terminal part of hTSPO in the inner membrane leaflet, while this impact is less pronounced for the upper, C–terminal part in the outer membrane leaflet, where the ligand binding site is located; (iv) very interestingly, cholesterol most frequently binds simultaneously to the so-called CRAC and CARC regions in TM V in the free form (residues L150–X–Y152–X(3)–R156 and R135–X(2)–Y138–X(2)–L141, respectively). However, when the protein is in complex with PK11195, cholesterol binds equally frequently to the CRAC–resembling motif that we observed in TM I (residues L17–X(2)–F20–X(3)–R24) and to CRAC in TM V. We expect that the CRAC–like motif in TM I will be of interest in future experimental investigations. Thus, our MD simulations provide insight into the structural features of hTSPO and the previously unknown interplay between PK11195 and cholesterol interactions with this pharmacologically relevant protein.

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.

The translocator protein (TSPO) is a transmembrane protein present across the three domains of life. Its functional quaternary structure consists of one or more subunits. In mice, the dimer-to-monomer equilibrium is shifted in vitro towards the monomer by adding cholesterol, a natural component of mammalian membranes. Here, we present a coarse-grained molecular dynamics study on the mouse protein in the presence of a physiological content and of an excess of cholesterol. The latter turns out to weaken the interfaces of the dimer by clusterizing mostly at the inter-monomeric space and pushing the contact residues apart. It also increases the compactness and the rigidity of the monomer. These two factors might play a role for the experimentally observed incremented stability of the monomeric form with increased content of cholesterol. Comparison with simulations on bacterial proteins suggests that the effect of cholesterol is much less pronounced for the latter than for the mouse protein.

Recent studies suggest that Tyr-39 might play a critical role for both the normal function and the pathological dysfunction of α-synuclein (αS), an intrinsically disordered protein involved in Parkinson’s disease. We perform here a comparative analysis between the structural features of human αS and its Y39A, Y39F, and Y39L variants. By the combined application of site-directed mutagenesis, biophysical techniques, and enhanced sampling molecular simulations, we show that removing aromatic functionality at position 39 of monomeric αS leads to protein variants populating more compact conformations, conserving its disordered nature and secondary structure propensities. Contrasting with the subtle changes induced by mutations on the protein structure, removing aromaticity at position 39 impacts strongly on the interaction of αS with the potent amyloid inhibitor phthalocyanine tetrasulfonate (PcTS). Our findings further support the role of Tyr-39 in forming essential inter and intramolecular contacts that might have important repercussions for the function and the dysfunction of αS.

G protein-coupled receptors (GPCRs) are the largest superfamily of signaling proteins. Their activation process is accompanied by conformational changes that have not yet been fully uncovered. Here, we carry out a novel comparative analysis of internal structural fluctuations across a variety of receptors from class A GPCRs, which currently has the richest structural coverage. We infer the local mechanical couplings underpinning the receptors' functional dynamics and finally identify those amino acids whose virtual deletion causes a significant softening of the mechanical network. The relevance of these amino acids is demonstrated by their overlap with those known to be crucial for GPCR function, based on static structural criteria. The differences with the latter set allow us to identify those sites whose functional role is more clearly detected by considering dynamical and mechanical properties. Of these sites with a genuine mechanical/dynamical character, the top ranking is amino acid 7x52, a previously unexplored, and experimentally verifiable key site for GPCR conformational response to ligand binding.

Copper is one of the metals described to bind the Parkinson disease-related protein α-synuclein (aSyn), and to promote its aggregation. Although histidine at position 50 in the aSyn sequence is one of the most studied copper-anchoring sites, its precise role in copper binding and aSyn aggregation is still unclear. Previous studies suggested that this residue does not significantly affect copper-mediated aSyn aggregation. However, our findings showed that the aggregation of the pathological H50Q aSyn mutant is enhanced by copper hints otherwise. Despite the inexistence of a model for aSyn H50Q-copper complexation, we discuss possible mechanisms by which this metal contributes to the misfolding and self-assembly of this particular aSyn mutant. Considering the genetic association of the H50Q mutation with familial forms of Parkinson disease, and the fact that copper homeostasis is deregulated in this disorder, understanding the interplay between both factors will shed light into the molecular and cellular mechanisms triggering the development and spreading of the aSyn pathology.

The N-terminus of the prion protein is a large intrinsically disordered region encompassing approximately 125 amino acids. In this paper, we review its structural and functional properties, with a particular emphasis on its binding to copper ions. The latter is exploited by the region’s conformational flexibility to yield a variety of biological functions. Disease-linked mutations and proteolytic processing of the protein can impact its copper-binding properties, with important structural and functional implications, both in health and disease progression.

The binding between the HIV-1 trans-activator of transcription (Tat) and p300/(CREB-binding protein)-associated factor (PCAF) bromodomain is a crucial step in the HIV-1 life cycle. However, the structure of the full length acetylated Tat bound to PCAF has not been yet determined experimentally. Acetylation of Tat residues can play a critical role in enhancing HIV-1 transcriptional activation. Here, we have combined a fully flexible protein-protein docking approach with molecular dynamics simulations to predict the structural determinants of the complex for the common HIV-1BRU variant. This model reproduces all the crucial contacts between the Tat peptide 46SYGR(AcK)KRRQRC56 and the PCAF bromodomain previously reported by NMR spectroscopy. Additionally, inclusion of the entire Tat protein results in additional contact points at the protein-protein interface. The model is consistent with the available experimental data reported and adds novel information to our previous structural predictions of the PCAF bromodomain in complex with the rare HIVZ2 variant, which was obtained with a less accurate computational method. This improved characterization of Tat.PCAF bromodomain binding may help in defining the structural determinants of other protein interactions involving lysine acetylation.

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.

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.

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.

Protein electrospray ionization (ESI) mass spectrometry (MS)-based techniques are widely used to provide insight into structural proteomics under the assumption that non-covalent protein complexes being transferred into the gas phase preserve basically the same intermolecular interactions as in solution. Here we investigate the applicability of this assumption by extending our previous structural prediction protocol for single proteins in ESI-MS to protein complexes. We apply our protocol to the human insulin dimer (hIns2) as a test case. Our calculations reproduce the main charge and the collision cross section (CCS) measured in ESI-MS experiments. Molecular dynamics simulations for 0.075 ms show that the complex maximizes intermolecular non-bonded interactions relative to the structure in water, without affecting the cross section. The overall gas-phase structure of hIns2 does exhibit differences with the one in aqueous solution, not inferable from a comparison with calculated CCS. Hence, care should be exerted when interpreting ESI-MS proteomics data based solely on NMR and/or X-ray structural information.

We report a girl with a de novo distal deletion of 9p affected by idiopathic central precocious puberty and intellectual disability. Genome-wide array-CGH revealed a terminal deletion of about 11 Mb, allowing to define her karyotype as 46; XX, del(9)(p23-pter). To our knowledge, this is the second reported case of precocious puberty associated with 9p distal deletion. A third case associates precocious puberty with a more proximal 9p deletion del(9)(p12p13,3). In our case, more than 40 genes were encompassed in the deleted region, among which, DMRT1 which is gonad-specific and has a sexually dimorphic expression pattern and ERMP1 which is required in rats for the organization of somatic cells and oocytes into discrete follicular structures. Although we cannot exclude that precocious puberty in our del(9p) patient is a coincidental finding, the report of the other two patients with 9p deletions and precocious puberty indeed suggests a causative relationship.