Characterization of p38α autophosphorylation inhibitors that target the non-canonical activation pathway.

p38α is a versatile protein kinase that can control numerous processes and plays important roles in the cellular responses to stress. Dysregulation of p38α signaling has been linked to several diseases including inflammation, immune disorders and cancer, suggesting that targeting p38α could be therapeutically beneficial. Over the last two decades, numerous p38α inhibitors have been developed, which showed promising effects in pre-clinical studies but results from clinical trials have been disappointing, fueling the interest in the generation of alternative mechanisms of p38α modulation. Here, we report the in silico identification of compounds that we refer to as non-canonical p38α inhibitors (NC-p38i). By combining biochemical and structural analyses, we show that NC-p38i efficiently inhibit p38α autophosphorylation but weakly affect the activity of the canonical pathway. Our results demonstrate how the structural plasticity of p38α can be leveraged to develop therapeutic opportunities targeting a subset of the functions regulated by this pathway.

The G protein database, GproteinDb.

Two-thirds of signaling substances, several sensory stimuli and over one-third of drugs act via receptors coupling to G proteins. Here, we present an online platform for G protein research with reference data and tools for analysis, visualization and design of scientific studies across disciplines and areas. This platform may help translate new pharmacological, structural and genomic data into insights on G protein signaling vital for human physiology and medicine. The G protein database is accessible at https://gproteindb.org.

GPCRdb in 2021: integrating GPCR sequence, structure and function.

G protein-coupled receptors (GPCRs) form both the largest family of membrane proteins and drug targets, mediating the action of one-third of medicines. The GPCR database, GPCRdb serves >4 000 researchers every month and offers reference data, analysis of own or literature data, experiment design and dissemination of published datasets. Here, we describe new and updated GPCRdb resources with a particular focus on integration of sequence, structure and function. GPCRdb contains all human non-olfactory GPCRs (and >27 000 orthologs), G-proteins and arrestins. It includes over 2 000 drug and in-trial agents and nearly 200 000 ligands with activity and availability data. GPCRdb annotates all published GPCR structures (updated monthly), which are also offered in a refined version (with re-modeled missing/distorted regions and reverted mutations) and provides structure models of all human non-olfactory receptors in inactive, intermediate and active states. Mutagenesis data in the GPCRdb spans natural genetic variants, GPCR-G protein interfaces, ligand sites and thermostabilising mutations. A new sequence signature tool for identification of functional residue determinants has been added and two data driven tools to design ligand site mutations and constructs for structure determination have been updated extending their coverage of receptors and modifications. The GPCRdb is available at https://gpcrdb.org.

QligFEP: an automated workflow for small molecule free energy calculations in Q.

The process of ligand binding to a biological target can be represented as the equilibrium between the relevant solvated and bound states of the ligand. This which is the basis of structure-based, rigorous methods such as the estimation of relative binding affinities by free energy perturbation (FEP). Despite the growing capacity of computing power and the development of more accurate force fields, a high throughput application of FEP is currently hampered due to the need, in the current schemes, of an expert user definition of the "alchemical" transformations between molecules in the series explored. Here, we present QligFEP, a solution to this problem using an automated workflow for FEP calculations based on a dual topology approach. In this scheme, the starting poses of each of the two ligands, for which the relative affinity is to be calculated, are explicitly present in the MD simulations associated with the (dual topology) FEP transformation, making the perturbation pathway between the two ligands univocal. We show that this generalized method can be applied to accurately estimate solvation free energies for amino acid sidechain mimics, as well as the binding affinity shifts due to the chemical changes typical of lead optimization processes. This is illustrated in a number of protein systems extracted from other FEP studies in the literature: inhibitors of CDK2 kinase and a series of A2A adenosine G protein-coupled receptor antagonists, where the results obtained with QligFEP are in excellent agreement with experimental data. In addition, our protocol allows for scaffold hopping perturbations to identify the binding affinities between different core scaffolds, which we illustrate with a series of Chk1 kinase inhibitors. QligFEP is implemented in the open-source MD package Q, and works with the most common family of force fields: OPLS, CHARMM and AMBER.

Free energy calculations of RNA interactions.

This review discusses the use of molecular dynamics free energy calculations for characterizing RNA interactions, with particular emphasis on molecular recognition events involved in mRNA translation on the ribosome. The general methodology for efficient free energy calculations is outlined and our specific implementation for binding free energy changes due to base mutations in mRNA and tRNA is described. We show that there are a number of key problems related to the accuracy of protein synthesis that can be addressed with this type of computational approach and several such examples are discussed in detail. These include the decoding of mRNA during peptide chain elongation, initiation and termination of translation, as well as the energetic effects of base tautomerization and tRNA modifications. It is shown that free energy calculations can be made sufficiently reliable to allow quantitative conclusions to be drawn regarding the energetics of cognate versus non-cognate interactions and its structural origins.

Computer simulations of the catalytic mechanism of wild-type and mutant ß-phosphoglucomutase.

ß-Phosphoglucomutase (ß-PGM) has served as an important model system for understanding biological phosphoryl transfer. This enzyme catalyzes the isomerization of ß-glucose-1-phosphate to ß-glucose-6-phosphate in a two-step process proceeding via a bisphosphate intermediate. The conventionally accepted mechanism is that both steps are concerted processes involving acid-base catalysis from a nearby aspartate (D10) side chain. This argument is supported by the observation that mutation of D10 leaves the enzyme with no detectable activity. However, computational studies have suggested that a substrate-assisted mechanism is viable for many phosphotransferases. Therefore, we carried out empirical valence bond (EVB) simulations to address the plausibility of this mechanistic alternative, including its role in the abolished catalytic activity of the D10S, D10C and D10N point mutants of ß-PGM. In addition, we considered both of these mechanisms when performing EVB calculations of the catalysis of the wild type (WT), H20A, H20Q, T16P, K76A, D170A and E169A/D170A protein variants. Our calculated activation free energies confirm that D10 is likely to serve as the general base/acid for the reaction catalyzed by the WT enzyme and all its variants, in which D10 is not chemically altered. Our calculations also suggest that D10 plays a dual role in structural organization and maintaining electrostatic balance in the active site. The correct positioning of this residue in a catalytically competent conformation is provided by a functionally important conformational change in this enzyme and by the extensive network of H-bonding interactions that appear to be exquisitely preorganized for the transition state stabilization.

Characterization of Ligand Binding to GPCRs Through Computational Methods.

The recent increase in available G protein-coupled receptor structures now contributes decisively to the structure-based ligand design. In this context, computational approaches in combination with medicinal chemistry and pharmacology are extremely helpful. Here, we provide an update on our structure-based computational protocols, used to answer key questions related to GPCR-ligand binding. All combined, these techniques can shed light on ligand binding modes, determine the molecular basis of conformational selection, for agonists and antagonists, as well as of subtype selectivity. To illustrate each of these questions, we will consider examples from existing projects on three families of class A (rhodopsin-like) GPCRs: one small-molecule (nucleotide-like) family, i.e., the adenosine receptors, and two peptide-binding receptors: neuropeptide-Y and angiotensin II receptors. The successful application of the same computational protocols to investigate this diverse group of receptor families gives an idea of the general applicability of our methodology in the characterization of GPCR-ligand binding.

A close-up view of codon selection in eukaryotic initiation.

When given an option to choose among a set of alternatives and only one selection is right, one might stop and reflect over which one is best. However, the ribosome has no time to stop and make such reflections, proteins need to be produced and very fast. Eukaryotic translation initiation is an example of such a conundrum. Here, scanning for the correct codon match must be fast, efficient and accurate. We highlight our recent computational findings, which show how the initiation machinery manages to recognize one specific codon among many possible challengers, by fine-tuning the energetic landscape of base-pairing with the aid of the initiation factors eIF1 and eIF1A. Using a recent 3-dimensional structure of the eukaryotic initiation complex we have performed simulations of codon recognition in atomic detail. These calculations provide an in-depth energetic and structural view of how discrimination against near-cognate codons is achieved by the initiation complex.

GPCR-ModSim: A comprehensive web based solution for modeling G-protein coupled receptors.

GPCR-ModSim (http://open.gpcr-modsim.org) is a centralized and easy to use service dedicated to the structural modeling of G-protein Coupled Receptors (GPCRs). 3D molecular models can be generated from amino acid sequence by homology-modeling techniques, considering different receptor conformations. GPCR-ModSim includes a membrane insertion and molecular dynamics (MD) equilibration protocol, which can be used to refine the generated model or any GPCR structure uploaded to the server, including if desired non-protein elements such as orthosteric or allosteric ligands, structural waters or ions. We herein revise the main characteristics of GPCR-ModSim and present new functionalities. The templates used for homology modeling have been updated considering the latest structural data, with separate profile structural alignments built for inactive, partially-active and active groups of templates. We have also added the possibility to perform multiple-template homology modeling in a unique and flexible way. Finally, our new MD protocol considers a series of distance restraints derived from a recently identified conserved network of helical contacts, allowing for a smoother refinement of the generated models which is particularly advised when there is low homology to the available templates. GPCR- ModSim has been tested on the GPCR Dock 2013 competition with satisfactory results.

Triple helical DNA in a duplex context and base pair opening.

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson-Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson-Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.

Cambios neuropsccológicos asociados a estimulacón cerebral profunda en enfermeedad de parkinson: revisión teórica / Neuropsychologcaal changes associated with dee brain stimulation of parkinsonn disease: theoreticalrreview

La enfermedad de Parkinson es una patología neurodegenerativa causada por la pérdida de células dopaminergicas en la sustancia negra mesencefálica. Esto produce una disfunción de los núcleos basales que se manifiesta con síntomas motores como temblor, rigidez y bradicinecia, entre otros. Con la estimulación cerebral profunda (ECP) ha resurgido la cirugía como opción terapéutica y es el núcleo subtalámico el área diana predilecta. Los estudios muestran mejoras significativas en los déficits motores, pero no hay claridad sobre los cambios neuropsicológicos de los pacientes sometidos a ECP. Se hace una revisión de los diferentes estudios que han investigado los cambios cognitivos, emocionales y comportamentales, concluyendo que la mayoría de habilidades cognitivas se mantienen o mejoran después de la ECP, pero pueden existir cambios emocionales y comportamentales adversos que están relacionadas con el núcleo cerebral donde se implanta el electrodo y con las características premorbidas de personalidad.

Parkinson's disease is a neurodegenerative disorder attributable to midbrain dopaminergic cell loss within the substantia nigra. This causes a dysfunction of the basal ganglia manifested by motor symptoms such as tremor, rigidity, bradykinesia among others. With Deep Brain Stimulation (DBS), neurosurgery has emerged as a therapeutic option, being the subthalamic nucleus its main target area. Studies show significant improvement in motor deficits, but there is no knowledge on the neuropsychological changes in patients after DBS. A review of several studies that have researched the cognitive, emotional and behavioral changes concluded that most cognitive skills are either maintained or improved after DBS, but there may be adverse emotional and behavioral changes that are related to the core brain where the electrode is implanted and with its premorbid personality characteristics.

New information content in RNA base pairing deduced from quantitative analysis of high-resolution structures.

Non-canonical base pairs play important roles in organizing the complex three-dimensional folding of RNA. Here, we outline methodology developed both to analyze the spatial patterns of interacting base pairs in known RNA structures and to reconstruct models from the collective experimental information. We focus attention on the structural context and deformability of the seven pairing patterns found in greatest abundance in the helical segments in a set of well-resolved crystal structures, including (i-ii) the canonical A.U and G.C Watson-Crick base pairs, (iii) the G.U wobble pair, (iv) the sheared G.A pair, (v) the A.U Hoogsteen pair, (vi) the U.U wobble pair, and (vii) the G.A Watson-Crick-like pair. The non-canonical pairs stand out from the canonical associations in terms of apparent deformability, spanning a broader range of conformational states as measured by the six rigid-body parameters used to describe the spatial arrangements of the interacting bases, the root-mean-square deviations of the base-pair atoms, and the fluctuations in hydrogen-bonding geometry. The deformabilties, the modes of base-pair deformation, and the preferred sites of occurrence depend on sequence. We also characterize the positioning and overlap of the base pairs with respect to the base pairs that stack immediately above and below them in double-helical fragments. We incorporate the observed positions of the bases, base pairs, and intervening phosphorus atoms in models to predict the effects of the non-canonical interactions on overall helical structure.