1-Hydroxy-2(1H)-pyridinone-Based Chelators with Potential Catechol O-Methyl Transferase Inhibition and Neurorescue Dual Action against Parkinson's Disease.

Two analogues of tolcapone where the nitrocatechol group has been replaced by a 1-hydroxy-2(1H)-pyridinone have been designed and synthesised. These compounds are expected to have a dual mode of action both beneficial against Parkinson's disease: they are designed to be inhibitors of catechol O-methyl transferase, which contribute to the reduction of dopamine in the brain, and to protect neurons against oxidative damage. To assess whether these compounds are worthy of biological assessment to demonstrate these effects, measurement of their pKa and stability constants for Fe(III), in silico modelling of their potential to inhibit COMT and blood-brain barrier scoring were performed. These results demonstrate that the compounds may indeed have the desired properties, indicating they are indeed promising candidates for further evaluation.

Agonistic and potentiating effects of perfluoroalkyl substances (PFAS) on the Atlantic cod (Gadus morhua) peroxisome proliferator-activated receptors (Ppars).

Toxicity mediated by per- and polyfluoroalkyl substances (PFAS), and especially perfluoroalkyl acids (PFAAs), has been linked to activation of peroxisome proliferator-activated receptors (Ppar) in many vertebrates. Here, we present the primary structures, phylogeny, and tissue-specific distributions of the Atlantic cod (Gadus morhua) gmPpara1, gmPpara2, gmPparb, and gmPparg, and demonstrate that the carboxylic acids PFHxA, PFOA, PFNA, as well as the sulfonic acid PFHxS, activate gmPpara1 in vitro, which was also supported by in silico analyses. Intriguingly, a binary mixture of PFOA and the non-activating PFOS produced a higher activation of gmPpara1 compared to PFOA alone, suggesting that PFOS has a potentiating effect on receptor activation. Supporting the experimental data, docking and molecular dynamics simulations of single and double-ligand complexes led to the identification of a putative allosteric binding site, which upon binding of PFOS stabilizes an active conformation of gmPpara1. Notably, binary exposures of gmPpara1, gmPpara2, and gmPparb to model-agonists and PFAAs produced similar potentiating effects. This study provides novel mechanistic insights into how PFAAs may modulate the Ppar signaling pathway by either binding the canonical ligand-binding pocket or by interacting with an allosteric binding site. Thus, individual PFAAs, or mixtures, could potentially modulate the Ppar-signaling pathway in Atlantic cod by interfering with at least one gmPpar subtype.

A computational study of the molecular basis of antibiotic resistance in a DXR mutant.

Proteins of the independent mevalonate pathway for isoprenoid biosynthesis are important targets for the development of new antibacterial compounds as this pathway is present in most pathogenic organisms such as Mycobacterium tuberculosis, DPlasmodium falciparum and Escherichia coli, but is not present in mammalian species, including humans. Deoxy-D-xylulose 5-phosphate reductoisomerase (DXR) is an important target in this pathway and the most effective DXR inhibitor to date is fosmidomycin, which is used to treat malaria and, more recently, tuberculosis. Recently, Armstrong C. M. et al. showed that a mutant of DXR, S222T, induces a loss of the fosmidomycin inhibition efficiency, even though the bacteria culture is still viable and able to produce isoprenoids. As this represents a potential fosmidomycin-resistant mutation, it is important to understand the mechanism of this apparent mutation-induced resistance to fosmidomycin. Here, we used molecular dynamics simulations and Molecular Mechanics/Poisson Boltzmann Surface Area analysis to understand the structural and energetic basis of the resistance. Our results suggest that the point mutation results in changes to the structural dynamics of an active site loop that probably protects the active site and facilitates enzymatic reaction. From the simulation analysis, we also showed that the mutation results in changes in the interaction energy profiles in a way that can explain the observed activity of the mutant protein toward the natural inhibitor deoxy-D-xylulose 5-phosphate. These results should be taken into consideration in future efforts to develop new therapeutic antibiotic compounds that target DXR.

Allostery in Its Many Disguises: From Theory to Applications.

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Recurrent activating mutations of PPARγ associated with luminal bladder tumors.

The upregulation of PPARγ/RXRα transcriptional activity has emerged as a key event in luminal bladder tumors. It renders tumor cell growth PPARγ-dependent and modulates the tumor microenvironment to favor escape from immuno-surveillance. The activation of the pathway has been linked to PPARG gains/amplifications resulting in PPARγ overexpression and to recurrent activating point mutations of RXRα. Here, we report recurrent mutations of PPARγ that also activate the PPARγ/RXRα pathway, conferring PPARγ-dependency and supporting a crucial role of PPARγ in luminal bladder cancer. These mutations are found throughout the protein-including N-terminal, DNA-binding and ligand-binding domains-and most of them enhance protein activity. Structure-function studies of PPARγ variants with mutations in the ligand-binding domain allow identifying structural elements that underpin their gain-of-function. Our study reveals genomic alterations of PPARG that lead to pro-tumorigenic PPARγ/RXRα pathway activation in luminal bladder tumors and may open the way towards alternative options for treatment.

Unrolr: Structural analysis of protein conformations using stochastic proximity embedding.

Molecular dynamics (MD) simulations are widely used to explore the conformational space of biological macromolecules. Advances in hardware, as well as in methods, make the generation of large and complex MD datasets much more common. Although different clustering and dimensionality reduction methods have been applied to MD simulations, there remains a need for improved strategies that handle nonlinear data and/or can be applied to very large datasets. We present an original implementation of the pivot-based version of the stochastic proximity embedding method aimed at large MD datasets using the dihedral distance as a metric. The advantages of the algorithm in terms of data storage and computational efficiency are presented, as well as the implementation realized. Application and testing through the analysis of a 200 ns accelerated MD simulation of a 35-residue villin headpiece is discussed. Analysis of the simulation shows the promise of this method to organize large conformational ensembles. © 2018 Wiley Periodicals, Inc.

Publisher Correction: Post-translational modifications in DNA topoisomerase 2α highlight the role of a eukaryote-specific residue in the ATPase domain.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Post-translational modifications in DNA topoisomerase 2α highlight the role of a eukaryote-specific residue in the ATPase domain.

Type 2 DNA topoisomerases (Top2) are critical components of key protein complexes involved in DNA replication, chromosome condensation and segregation, as well as gene transcription. The Top2 were found to be the main targets of anticancer agents, leading to intensive efforts to understand their functional and physiological role as well as their molecular structure. Post-translational modifications have been reported to influence Top2 enzyme activities in particular those of the mammalian Top2α isoform. In this study, we identified phosphorylation, and for the first time, acetylation sites in the human Top2α isoform produced in eukaryotic expression systems. Structural analysis revealed that acetylation sites are clustered on the catalytic domains of the homodimer while phosphorylation sites are located in the C-terminal domain responsible for nuclear localization. Biochemical analysis of the eukaryotic-specific K168 residue in the ATPase domain shows that acetylation affects a key position regulating ATP hydrolysis through the modulation of dimerization. Our findings suggest that acetylation of specific sites involved in the allosteric regulation of human Top2 may provide a mechanism for modulation of its catalytic activity.

Alternative dimerization interfaces in the glucocorticoid receptor-α ligand binding domain.

BACKGROUND: Nuclear hormone receptors (NRs) constitute a large family of multi-domain ligand-activated transcription factors. Dimerization is essential for their regulation, and both DNA binding domain (DBD) and ligand binding domain (LBD) are implicated in dimerization. Intriguingly, the glucocorticoid receptor-α (GRα) presents a DBD dimeric architecture similar to that of the homologous estrogen receptor-α (ERα), but an atypical dimeric architecture for the LBD. The physiological relevance of the proposed GRα LBD dimer is a subject of debate. METHODS: We analyzed all GRα LBD homodimers observed in crystals using an energetic analysis based on the PISA and on the MM/PBSA methods and a sequence conservation analysis, using the ERα LBD dimer as a reference point. RESULTS: Several dimeric assemblies were observed for GRα LBD. The assembly generally taken to be physiologically relevant showed weak binding free energy and no significant residue conservation at the contact interface, while an alternative homodimer mediated by both helix 9 and C-terminal residues showed significant binding free energy and residue conservation. However, none of the GRα LBD assemblies found in crystals are as stable or conserved as the canonical ERα LBD dimer. GRα C-terminal sequence (F-domain) forms a steric obstacle to the canonical dimer assembly in all available structures. CONCLUSIONS: Our analysis calls for a re-examination of the currently accepted GRα homodimer structure and experimental investigations of the alternative architectures. GENERAL SIGNIFICANCE: This work questions the validity of the currently accepted architecture. This has implications for interpreting physiological data and for therapeutic design pertaining to glucocorticoid research.

Peptide-Protein Binding Investigated by Far-IR Spectroscopy and Molecular Dynamics Simulations.

Molecular dynamics (MD) simulations and far-infrared (far-IR) spectroscopy were combined to study peptide binding by the second PDZ domain (PDZ1) of MAGI1, which has been identified as an important target for the Human Papilloma Virus. PDZ1 recognizes and binds to the C-terminal end of the E6 protein from high-risk Human Papilloma Virus. The far-IR spectra of two forms of the protein, an unbound APO form and a HOLO form (where the PDZ1 is bound to an 11-residue peptide derived from the C terminus of HPV16 E6), were obtained. MD simulations were used to determine the most representative structure of each form and these were used to compute their respective IR spectra by normal mode analysis. Far-UV circular dichroism spectroscopy was used to confirm the secondary structure content and the stability through temperature-dependent studies. Both the experimental and calculated far-IR spectra showed a red shift of the low-frequency peaks upon peptide binding. The calculations show that this is coincident with an increased number of hydrogen bonds formed as the peptide augments the protein ß-sheet. We further identified the contribution of surface-bound water molecules to bands in the far-IR and, through the calculations, identified potential pathways for allosteric communication. Together, these results demonstrate the utility of combining far-IR experiments and MD studies to study peptide binding by proteins.

Allosteric Regulation in the Ligand Binding Domain of Retinoic Acid Receptorγ.

Retinoic acid (RA) plays key roles in cell differentiation and growth arrest through nuclear retinoic acid receptors (RARs), which are ligand-dependent transcription factors. While the main trigger of RAR activation is the binding of RA, phosphorylation of the receptors has also emerged as an important regulatory signal. Phosphorylation of the RARγ N-terminal domain (NTD) is known to play a functional role in neuronal differentiation. In this work, we investigated the phosphorylation of RARγ ligand binding domain (LBD), and present evidence that the phosphorylation status of the LBD affects the phosphorylation of the NTD region. We solved the X-ray structure of a phospho-mimetic mutant of the LBD (RARγ S371E), which we used in molecular dynamics simulations to characterize the consequences of the S371E mutation on the RARγ structural dynamics. Combined with simulations of the wild-type LBD, we show that the conformational equilibria of LBD salt bridges (notably R387-D340) are affected by the S371E mutation, which likely affects the recruitment of the kinase complex that phosphorylates the NTD. The molecular dynamics simulations also showed that a conservative mutation in this salt bridge (R387K) affects the dynamics of the LBD without inducing large conformational changes. Finally, cellular assays showed that the phosphorylation of the NTD of RARγ is differentially regulated by retinoic acid in RARγWT and in the S371N, S371E and R387K mutants. This multidisciplinary work highlights an allosteric coupling between phosphorylations of the LBD and the NTD of RARγ and supports the importance of structural dynamics involving electrostatic interactions in the regulation of RARs activity.

Synthesis and biological evaluation of phosphate isosters of fosmidomycin and analogs as inhibitors of Escherichia coli and Mycobacterium smegmatis 1-deoxyxylulose 5-phosphate reductoisomerases.

Hydroxamate analogs of fosfoxacin, the phosphate homolog of fosmidomycin, have been synthesized and their activity tested on Escherichia coli and Mycobacterium smegmatis DXRs. Except for compound 4b, the IC50 values of phosphate derivatives are approximately 10-fold higher than those of the corresponding phosphonates. Although their inhibitory activity on Escherichia coli DXR is less efficient than their phosphonate analogs, we report the ability of phosphate compounds to inhibit the growth of Escherichia coli. This work points out that the uptake of fosfoxacin and its analogs is taking place via the GlpT and UhpT transporters. As expected, these compounds are inefficient to inhibit the growth of M. smegmatis growth inhibition probably due to a lack of uptake.

Tautomers of a Fluorescent G Surrogate and Their Distinct Photophysics Provide Additional Information Channels.

Thienoguanosine ((th) G) is an isomorphic nucleoside analogue acting as a faithful fluorescent substitute of G, with respectable quantum yield in oligonucleotides. Photophysical analysis of (th) G reveals the existence of two ground-state tautomers with significantly shifted absorption and emission wavelengths, and high quantum yield in buffer. Using (TD)-DFT calculations, the tautomers were identified as the H1 and H3 keto-amino tautomers. When incorporated into the loop of (-)PBS, the (-)DNA copy of the HIV-1 primer binding site, both tautomers are observed and show differential sensitivity to protein binding. The red-shifted H1 tautomer is strongly favored in matched (-)/(+)PBS duplexes, while the relative emission of the H3 tautomer can be used to detect single nucleotide polymorphisms. These tautomers and their distinct environmental sensitivity provide unprecedented information channels for analyzing G residues in oligonucleotides and their complexes.

PSSweb: protein structural statistics web server.

With the increasing number of protein structures available, there is a need for tools capable of automating the comparison of ensembles of structures, a common requirement in structural biology and bioinformatics. PSSweb is a web server for protein structural statistics. It takes as input an ensemble of PDB files of protein structures, performs a multiple sequence alignment and computes structural statistics for each position of the alignment. Different optional functionalities are proposed: structure superposition, Cartesian coordinate statistics, dihedral angle calculation and statistics, and a cluster analysis based on dihedral angles. An interactive report is generated, containing a summary of the results, tables, figures and 3D visualization of superposed structures. The server is available at http://pssweb.org.

Quantitative sampling of conformational heterogeneity of a DNA hairpin using molecular dynamics simulations and ultrafast fluorescence spectroscopy.

Molecular dynamics (MD) simulations and time resolved fluorescence (TRF) spectroscopy were combined to quantitatively describe the conformational landscape of the DNA primary binding sequence (PBS) of the HIV-1 genome, a short hairpin targeted by retroviral nucleocapsid proteins implicated in the viral reverse transcription. Three 2-aminopurine (2AP) labeled PBS constructs were studied. For each variant, the complete distribution of fluorescence lifetimes covering 5 orders of magnitude in timescale was measured and the populations of conformers experimentally observed to undergo static quenching were quantified. A binary quantification permitted the comparison of populations from experimental lifetime amplitudes to populations of aromatically stacked 2AP conformers obtained from simulation. Both populations agreed well, supporting the general assumption that quenching of 2AP fluorescence results from pi-stacking interactions with neighboring nucleobases and demonstrating the success of the proposed methodology for the combined analysis of TRF and MD data. Cluster analysis of the latter further identified predominant conformations that were consistent with the fluorescence decay times and amplitudes, providing a structure-based rationalization for the wide range of fluorescence lifetimes. Finally, the simulations provided evidence of local structural perturbations induced by 2AP. The approach presented is a general tool to investigate fine structural heterogeneity in nucleic acid and nucleoprotein assemblies.

Far infrared spectra of solid state L-serine, L-threonine, L-cysteine, and L-methionine in different protonation states.

In this study, experimental far infrared measurements of L-serine, L-threonine, L-cysteine, and L-methionine are presented showing the spectra for the 1.0-13.0 pH range. In parallel, solid state DFT calculations were performed on the amino acid zwitterions in the crystalline form. We focused on the lowest frequency far infrared normal modes, which required the most precision and convergence of the calculations. Analysis of the computational results, which included the potential energy distribution of the vibrational modes, permitted a detailed and almost complete assignment of the experimental spectrum. In addition to characteristic signals of the two main acid-base couples, CO2H/CO2(-) and NH3(+)/NH2, specific side chain contributions for these amino acids, including CCO and CCS vibrational modes were analyzed. This study is in line with the growing application of FIR measurements to biomolecules.

Molecular dynamics for computational proteomics of methylated histone H3.

BACKGROUND: Post-translational modifications of histones, and in particular of their disordered N-terminal tails, play a major role in epigenetic regulation. The identification of proteins and proteic domains that specifically bind modified histones is therefore of paramount importance to understand the molecular mechanisms of epigenetics. METHODS: We performed an energetic analysis using the MM/PBSA method in order to study known complexes between methylated histone H3 and effector domains of the PHD family. We then developed a simple molecular dynamics based predictive model based on our analysis. RESULTS: We present a thorough validation of our procedure, followed by the computational predictions of new PHD domains specific for binding histone H3 methylated on lysine 4 (K4). CONCLUSIONS: PHD domains recognize methylated K4 on histone H3 in the context of a linear interaction motif (LIM) formed by the first four amino acids of histone H3 as opposed to recognition of a single methylated site. PHD domains with different sequences find chemically equivalent solutions for stabilizing the histone LIM and these can be identified from energetic analysis. This analysis, in turn, allows for the identification of new PHD domains that bind methylated H3K4 using information that cannot be retrieved from sequence comparison alone. GENERAL SIGNIFICANCE: Molecular dynamics simulations can be used to devise computational proteomics protocols that are both easy to implement and interpret, and that yield reliable predictions that compare favorably to and complement experimental proteomics methods. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Structural insights into the molecular mechanism of vitamin D receptor activation by lithocholic acid involving a new mode of ligand recognition.

The vitamin D receptor (VDR), an endocrine nuclear receptor for 1α,25-dihydroxyvitamin D3, acts also as a bile acid sensor by binding lithocholic acid (LCA). The crystal structure of the zebrafish VDR ligand binding domain in complex with LCA and the SRC-2 coactivator peptide reveals the binding of two LCA molecules by VDR. One LCA binds to the canonical ligand-binding pocket, and the second one, which is not fully buried, is anchored to a site located on the VDR surface. Despite the low affinity of the alternative site, the binding of the second molecule promotes stabilization of the active receptor conformation. Biological activity assays, structural analysis, and molecular dynamics simulations indicate that the recognition of two ligand molecules is crucial for VDR agonism by LCA. The unique binding mode of LCA provides clues for the development of new chemical compounds that target alternative binding sites for therapeutic applications.

Protein structural statistics with PSS.

Characterizing the variability within an ensemble of protein structures is a common requirement in structural biology and bioinformatics. With the increasing number of protein structures becoming available, there is a need for new tools capable of automating the structural comparison of large ensemble of structures. We present Protein Structural Statistics (PSS), a command-line program written in Perl for Unix-like environments, dedicated to the calculation of structural statistics for a set of proteins. PSS can perform multiple sequence alignments, structure superpositions, calculate Cartesian and dihedral coordinate statistics, and execute cluster analyses. An HTML report that contains a convenient summary of results with figures, tables, and hyperlinks can also be produced. PSS is a new tool providing an automated way to compare multiple structures. It integrates various types of structural analyses through an user-friendly and flexible interface, facilitating the access to powerful but more specialized programs. PSS is easy to modify and extend and is distributed under a free and open source license. The relevance of PSS is illustrated by examples of application to pertinent biological problems.

Phosphorylation of the retinoic acid receptor alpha induces a mechanical allosteric regulation and changes in internal dynamics.

Nuclear receptor proteins constitute a superfamily of proteins that function as ligand dependent transcription factors. They are implicated in the transcriptional cascades underlying many physiological phenomena, such as embryogenesis, cell growth and differentiation, and apoptosis, making them one of the major signal transduction paradigms in metazoans. Regulation of these receptors occurs through the binding of hormones, and in the case of the retinoic acid receptor (RAR), through the binding of retinoic acid (RA). In addition to this canonical scenario of RAR activity, recent discoveries have shown that RAR regulation also occurs as a result of phosphorylation. In fact, RA induces non-genomic effects, such as the activation of kinase signaling pathways, resulting in the phosphorylation of several targets including RARs themselves. In the case of RARα, phosphorylation of Ser369 located in loop L9-10 of the ligand-binding domain leads to an increase in the affinity for the protein cyclin H, which is part of the Cdk-activating kinase complex of the general transcription factor TFIIH. The cyclin H binding site in RARα is situated more than 40 Å from the phosphorylated serine. Using molecular dynamics simulations of the unphosphorylated and phosphorylated forms of the receptor RARα, we analyzed the structural implications of receptor phosphorylation, which led to the identification of a structural mechanism for the allosteric coupling between the two remote sites of interest. The results show that phosphorylation leads to a reorganization of a local salt bridge network, which induces changes in helix extension and orientation that affects the cyclin H binding site. This results in changes in conformation and flexibility of the latter. The high conservation of the residues implicated in this signal transduction suggests a mechanism that could be applied to other nuclear receptor proteins.