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.

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.

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.

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.

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.

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.

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.

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.

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.

Tex19 and Sectm1 concordant molecular phylogenies support co-evolution of both eutherian-specific genes.

BACKGROUND: Transposable elements (TE) have attracted much attention since they shape the genome and contribute to species evolution. Organisms have evolved mechanisms to control TE activity. Testis expressed 19 (Tex19) represses TE expression in mouse testis and placenta. In the human and mouse genomes, Tex19 and Secreted and transmembrane 1 (Sectm1) are neighbors but are not homologs. Sectm1 is involved in immunity and its molecular phylogeny is unknown. METHODS: Using multiple alignments of complete protein sequences (MACS), we inferred Tex19 and Sectm1 molecular phylogenies. Protein conserved regions were identified and folds were predicted. Finally, expression patterns were studied across tissues and species using RNA-seq public data and RT-PCR. RESULTS: We present 2 high quality alignments of 58 Tex19 and 58 Sectm1 protein sequences from 48 organisms. First, both genes are eutherian-specific, i.e., exclusively present in mammals except monotremes (platypus) and marsupials. Second, Tex19 and Sectm1 have both duplicated in Sciurognathi and Bovidae while they have remained as single copy genes in all further placental mammals. Phylogenetic concordance between both genes was significant (p-value < 0.05) and supported co-evolution and functional relationship. At the protein level, Tex19 exhibits 3 conserved regions and 4 invariant cysteines. In particular, a CXXC motif is present in the N-terminal conserved region. Sectm1 exhibits 2 invariant cysteines and an Ig-like domain. Strikingly, Tex19 C-terminal conserved region was lost in Haplorrhini primates while a Sectm1 C-terminal extra domain was acquired. Finally, we have determined that Tex19 and Sectm1 expression levels anti-correlate across the testis of several primates (ρ = -0.72) which supports anti-regulation. CONCLUSIONS: Tex19 and Sectm1 co-evolution and anti-regulated expressions support a strong functional relationship between both genes. Since Tex19 operates a control on TE and Sectm1 plays a role in immunity, Tex19 might suppress an immune response directed against cells that show TE activity in eutherian reproductive tissues.

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.

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.

Evolution of nuclear retinoic acid receptor alpha (RARα) phosphorylation sites. Serine gain provides fine-tuned regulation.

The human nuclear retinoic acid (RA) receptor alpha (hRARα) is a ligand-dependent transcriptional regulator, which is controlled by a phosphorylation cascade. The cascade starts with the RA-induced phosphorylation of a serine residue located in the ligand-binding domain, S(LBD), allowing the recruitment of the cdk7/cyclin H/MAT1 subcomplex of TFIIH through the docking of cyclin H. It ends by the subsequent phosphorylation by cdk7 of an other serine located in the N-terminal domain, S(NTD). Here, we show that this cascade relies on an increase in the flexibility of the domain involved in cyclin H binding, subsequently to the phosphorylation of S(LBD). Owing to the functional importance of RARα in several vertebrate species, we investigated whether the phosphorylation cascade was conserved in zebrafish (Danio rerio), which expresses two RARα genes: RARα-A and RARα-B. We found that in zebrafish RARαs, S(LBD) is absent, whereas S(NTD) is conserved and phosphorylated. Therefore, we analyzed the pattern of conservation of the phosphorylation sites and traced back their evolution. We found that S(LBD) is most often absent outside mammalian RARα and appears late during vertebrate evolution. In contrast, S(NTD) is conserved, indicating that the phosphorylation of this functional site has been under ancient high selection constraint. This suggests that, during evolution, different regulatory circuits control RARα activity.

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.

Force field parameters for the simulation of modified histone tails.

We describe the development of force field parameters for methylated lysines and arginines, and acetylated lysine for the CHARMM all-atom force field. We also describe a CHARMM united-atom force field for modified sidechains suitable for use with fragment-based docking methods. The development of these parameters is based on results of ab initio quantum mechanics calculations of model compounds with subsequent refinement and validation by molecular mechanics and molecular dynamics simulations. The united-atom parameters are tested by fragment docking to target proteins using the MCSS procedure. The all-atom force field is validated by molecular dynamics simulations of multiple experimental structures. In both sets of calculations, the computational predictions using the force field were compared to the corresponding experimental structures. We show that the parameters yield an accurate reproduction of experimental structures. Together with the existing CHARMM force field, these parameters will enable the general modeling of post-translational modifications of histone tails.

Dynamic correlation networks in human peroxisome proliferator-activated receptor-γ nuclear receptor protein.

Peroxisome proliferator-activated receptor-γ nuclear receptor (PPAR-γ) belongs to the superfamily of nuclear receptor proteins that function as ligand-dependent transcription factors and plays a specific physiological role as a regulator of lipid metabolism. A number of experimental studies have suggested that allostery plays an important role in the functioning of PPAR-γ. Here we use normal-mode analysis of PPAR-γ to characterize a network of dynamically coupled amino acids that link physiologically relevant binding surfaces such as the ligand-dependent activation domain AF-2 with the ligand binding site and the heterodimer interface. Multiple calculations were done in both the presence and absence of the agonist rosiglitazone, and the differences in dynamics were characterized. The global dynamics of the ligand binding domain were affected by the ligand, and in particular, changes to the network of dynamically correlated amino acids were observed with only small changes in conformation. These results suggest that changes in dynamic couplings can be functionally significant with respect to the transmission of allosteric signals.

Far infrared spectra of solid state aliphatic amino acids in different protonation states.

Far infrared spectra of zwitterionic, cationic, and anionic forms of aliphatic amino acids in solid state have been studied experimentally. Measurements were done on glycine, L-alanine, L-valine, L-leucine, and L-isoleucine powder samples and film samples obtained from dried solutions prepared at pH ranging from 1 to 13. Solid state density functional theory calculations were also performed, and detailed potential energy distributions were obtained from normal mode results. A good correspondence between experimental and simulated spectra was achieved and this allowed us to propose an almost complete band assignment for the far infrared spectra of zwitterionic forms. In the 700-50 cm(-1) range, three regions were identified, each corresponding to a characteristic set of normal modes. A first region between 700 and 450 cm(-1) mainly contained the carboxylate bending, rocking, and wagging modes as well as the ammonium torsional mode. The 450-250 cm(-1) region was representative of backbone and sidechain skeletal bending modes. At last, the low wavenumber zone, below 250 cm(-1), was characteristic of carboxylate and skeletal torsional modes and of lattice modes. Assignments are also proposed for glycine cationic and anionic forms, but could not be obtained for all aliphatic amino acids due to the lack of structural data. This work is intended to provide fundamental information for the understanding of peptides vibrational properties.

Dynamics of beta3 integrin I-like and hybrid domains: insight from simulations on the mechanism of transition between open and closed forms.

The conformational dynamics of the I-like and Hybrid domains from the beta3 integrin headpiece were studied by molecular dynamics simulation and normal mode analysis. Crystallographic structures of integrins show that the integrin headpiece can exist in largely different conformations manifested by a significant difference in the angle between the I-like and Hybrid domains. The relative orientation of these two domains is believed to be a crucial element of integrin function, as it may relate local structural modifications induced by ligand binding into large-scale conformational changes. To investigate the detailed mechanisms responsible for this coupling, we carried out molecular dynamics simulations of the I-like/Hybrid system and employed quasi-harmonic and normal mode analyses to characterize the large-scale motions. Our results show that the conformational transition of I-like and Hybrid domains inferred from crystallographic data is contained in the low-frequency dynamics of the system. Using targeted molecular dynamics simulations, we investigated the roles played by two structural elements of the I-like domain, the alpha7 and alpha1 helices, in the interdomain transition. From our results, we propose that these two helices function in tandem to initiate large-scale, interdomain conformational transition apparent in integrin activation and signaling.

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.