Search

Background and Purpose: Castration-resistant prostate cancer (CRPC) is a common male malignancy that requires new therapeutic strategies due to acquired resistance to its first-line treatment, docetaxel. The benefits of vitamin D on prostate cancer (PCa) progression have been previously reported. This study aimed to investigate the effects of vitamin D on chemoresistance in CRPC., Experimental Approach: Structure function relationships of potent vitamin D analogues were determined. The combination of the most potent analogue and docetaxel was explored in chemoresistant primary PCa spheroids and in a xenograft mouse model derived from a patient with a chemoresistant CRPC., Key Results: Here, we show that Xe4MeCF3 is more potent than the natural ligand to induce vitamin D receptor (VDR) transcriptional activities and that it has a larger therapeutic window. Moreover, we demonstrate that VDR agonists restore docetaxel sensitivity in PCa spheroids. Importantly, Xe4MeCF3 reduces tumour growth in a chemoresistant CRPC patient-derived xenograft. In addition, this treatment targets signalling pathways associated with cancer progression in the remaining cells., Conclusion and Implications: Taken together, these results unravel the potency of VDR agonists to overcome chemoresistance in CRPC and open new avenues for the clinical management of PCa., (© 2024 The Author(s). British Journal of Pharmacology published by John Wiley & Sons Ltd on behalf of British Pharmacological Society.)

The Vitamin D receptor (VDR) plays a key role in calcium homeostasis, as well as in cell proliferation and differentiation. Among the large number of VDR ligands that have been developed, we have previously shown that BXL-62 and Gemini-72, two C-20-modified vitamin D analogs are highly potent VDR agonists. In this study, we show that both VDR ligands restore the transcriptional activities of VDR variants unresponsive to the natural ligand and identified in patients with rickets. The elucidated mechanisms of action underlying the activities of these C-20-modified analogs emphasize the mutual adaptation of the ligand and the VDR ligand-binding pocket.

Persistent infection with some mucosal α-genus human papillomaviruses (HPVs; the most prevalent one being HPV16) can induce cervical carcinoma, anogenital cancers, and a subset of head and neck squamous cell carcinoma (HNSCC). Cutaneous β-genus HPVs (such as HPV5 and HPV8) associate with skin lesions that can progress into squamous cell carcinoma with sun exposure in Epidermodysplasia verruciformis patients and immunosuppressed patients. Here, we analyzed mechanisms used by E6 proteins from the α- and β-genus to inhibit the interferon-β (IFNB1) response. HPV16 E6 mediates this effect by a strong direct interaction with interferon regulatory factor 3 (IRF3). The binding site of E6 was localized within a flexible linker between the DNA-binding domain and the IRF-activation domain of IRF3 containing an LxxLL motif. The crystallographic structure of the complex between HPV16 E6 and the LxxLL motif of IRF3 was solved and compared with the structure of HPV16 E6 interacting with the LxxLL motif of the ubiquitin ligase E6AP. In contrast, cutaneous HPV5 and HPV8 E6 proteins bind to the IRF3-binding domain (IBiD) of the CREB-binding protein (CBP), a key transcriptional coactivator in IRF3-mediated IFN-β expression. IMPORTANCE Persistent HPV infections can be associated with the development of several cancers. The ability to persist depends on the ability of the virus to escape the host immune system. The type I interferon (IFN) system is the first-line antiviral defense strategy. HPVs carry early proteins that can block the activation of IFN-I. Among mucosal α-genus HPV types, the HPV16 E6 protein has a remarkable property to strongly interact with the transcription factor IRF3. Instead, cutaneous HPV5 and HPV8 E6 proteins bind to the IRF3 cofactor CBP. These results highlight the versatility of E6 proteins to interact with different cellular targets. The interaction between the HPV16 E6 protein and IRF3 might contribute to the higher prevalence of HPV16 than that of other high-risk mucosal HPV types in HPV-associated cancers.

Nuclear receptors are ligand-activated transcription factors that modulate gene regulatory networks from embryonic development to adult physiology and thus represent major targets for clinical interventions in many diseases. Most nuclear receptors function either as homodimers or as heterodimers. The dimerization is crucial for gene regulation by nuclear receptors, by extending the repertoire of binding sites in the promoters or the enhancers of target genes via combinatorial interactions. Here, we focused our attention on an unusual structural variation of the α-helix, called π-turn that is present in helix H7 of the ligand-binding domain of RXR and HNF4. By tracing back the complex evolutionary history of the π-turn, we demonstrate that it was present ancestrally and then independently lost in several nuclear receptor lineages. Importantly, the evolutionary history of the π-turn motif is parallel to the evolutionary diversification of the nuclear receptor dimerization ability from ancestral homodimers to derived heterodimers. We then carried out structural and biophysical analyses, in particular through point mutation studies of key RXR signature residues and showed that this motif plays a critical role in the network of interactions stabilizing homodimers. We further showed that the π-turn was instrumental in allowing a flexible heterodimeric interface of RXR in order to accommodate multiple interfaces with numerous partners and critical for the emergence of high affinity receptors. Altogether, our work allows to identify a functional role for the π-turn in oligomerization of nuclear receptors and reveals how this motif is linked to the emergence of a critical biological function. We conclude that the π-turn can be viewed as a structural exaptation that has contributed to enlarging the functional repertoire of nuclear receptors., Competing Interests: The authors declare that they have no conflict of interest.

Non-canonical residues, caps, crosslinks, and nicks are important to many functions of DNAs, RNAs, proteins, and complexes. However, we do not fully understand how networks of such non-canonical macromolecules generate behavior. One barrier is our limited formats for describing macromolecules. To overcome this barrier, we develop BpForms and BcForms, a toolkit for representing the primary structure of macromolecules as combinations of residues, caps, crosslinks, and nicks. The toolkit can help omics researchers perform quality control and exchange information about macromolecules, help systems biologists assemble global models of cells that encompass processes such as post-translational modification, and help bioengineers design cells.

Opportunistic pathogens such as Candida species can use carboxylic acids, like acetate and lactate, to survive and successfully thrive in different environmental niches. These nonfermentable substrates are frequently the major carbon sources present in certain human body sites, and their efficient uptake by regulated plasma membrane transporters plays a critical role in such nutrient-limited conditions. Here, we cover the physiology and regulation of these proteins and their potential role in Candida virulence. This review also presents an evolutionary analysis of orthologues of the Saccharomyces cerevisiae Jen1 lactate and Ady2 acetate transporters, including a phylogenetic analysis of 101 putative carboxylate transporters in twelve medically relevant Candida species. These proteins are assigned to distinct clades according to their amino acid sequence homology and represent the major carboxylic acid uptake systems in yeast. While Jen transporters belong to the sialate:H + symporter (SHS) family, the Ady2 homologue members are assigned to the acetate uptake transporter (AceTr) family. Here, we reclassify the later members as ATO (acetate transporter ortholog). The new nomenclature will facilitate the study of these transporters, as well as the analysis of their relevance for Candida pathogenesis., (Copyright © 2020 Alves et al.)

Depending on the amino acid sequence, as well as the local environment, some peptides have the capability to fold into multiple secondary structures. Conformational switching between such structures is a key element of protein folding and aggregation. Specifically, understanding the molecular mechanism underlying the transition from an α-helix to a β-hairpin is critical because it is thought to be a harbinger of amyloid assembly. In this study, we explore the energy landscape for an 18-residue peptide (DP5), designed by Araki and Tamura to exhibit equal propensities for the α-helical and β-hairpin forms. We find that the degeneracy is encoded in the multifunnel nature of the underlying free energy landscape. In agreement with experiment, we also observe that mutation of tyrosine at position 12 to serine shifts the equilibrium in favor of the α-helix conformation, by altering the landscape topography. The transition from the α-helix to the β-hairpin is a complex stepwise process, and occurs via collapsed coil-like intermediates. Our findings suggest that even a single mutation can tune the emergent features of the landscape, providing an efficient route to protein design. Interestingly, the transition pathways for the conformational switch seem to be minimally perturbed upon mutation, suggesting that there could be universal microscopic features that are conserved among different switch-competent protein sequences.

Background: In the last two decades, a new phenotype termed Sarcopenic Obesity (SO), in which sarcopenia and obesity coexist, has emerged., Objective: The aim of this systematic review and meta-analysis was first to assess the prevalence of Metabolic syndrome (Mets) among individuals with and without SO, and second, to determine if SO may increase the relative risk of Mets., Methods: This study was conducted in adherence to the Preferred Reporting Items for Systematic Review and Meta-Analyses (PRISMA) guidelines and the data were collated by means of metaanalysis and narrative synthesis., Results: Twelve studies including a total of 11,308 adults with overweight or obesity of both genders met the inclusion criteria and were reviewed, revealing two main findings. First, a similar overall prevalence of Mets in individuals with SO (61.49%; 95% CI: 52.19-70.40) when compared to those without SO (56.74%; 95% CI: 47.32-65.93) was identified. Second, the presence of SO appears not to increase the risk of Mets with respect to those without SO (RR = 1.08, 95% CI: 0.99- 1.17, p = 0.07)., Conclusion: No higher prevalence of Mets among individuals with SO when compared to those with obesity only, nor a significant association between SO and a higher risk of Mets was found., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

Coumermycin A1 is a natural aminocoumarin that inhibits bacterial DNA gyrase, a member of the GHKL proteins superfamily. We report here the first cocrystal structures of gyrase B bound to coumermycin A1, revealing that one coumermycin A1 molecule traps simultaneously two ATP-binding sites. The inhibited dimers from different species adopt distinct sequence-dependent conformations, alternative to the ATP-bound form. These structures provide a basis for the rational development of coumermycin A1 derivatives for antibiotherapy and biotechnology 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., (Copyright © 2019. Published by Elsevier Ltd.)

Retinoid X Receptors (RXRs) act as dimer partners for several nuclear receptors including itself, binding to genomic DNA response elements and regulating gene transcription with cell and gene specificity. As homodimers, RXRs bind direct repeats of the half-site (A/G)G(G/T)TCA separated by 1 nucleotide (DR1) and little variability of this consensus site is observed for natural DR1s. However, these variations are responsible of the modulation of RXR receptors function through differential binding affinity and conformational changes. To further our understanding of the molecular mechanisms underlying RXR-DNA interactions, we examined how RXR DBDs bind to different DR1s using thermodynamics, X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy. We show that the half-site sequences modulate the binding cooperativity that results from the protein-protein contacts between the two DBDs. Chemical shifts perturbation NMR experiments revealed that sequence variations in half-sites induce changes that propagate from the protein-DNA interface to the dimerization interface throughout the DBD fold., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Steroid hormone receptors are important regulators of development and physiology in bilaterian animals, but the role of steroid signaling in cnidarians has been contentious. Cnidarians produce steroids, including A-ring aromatic steroids with a side-chain, but these are probably made through pathways different than the one used by vertebrates to make their A-ring aromatic steroids. Here we present comparative genomic analyses indicating the presence of a previously undescribed nuclear receptor family within medusozoan cnidarians, that we propose to call NR3E. This family predates the diversification of ERR/ER/SR in bilaterians, indicating that the first NR3 evolved in the common ancestor of the placozoan and cnidarian-bilaterian with lineage-specific loss in the anthozoans, even though multiple species in this lineage have been shown to produce aromatic steroids, whose function remain unclear. We discovered serendipitously that a cytoplasmic factor within epidermal cells of transgenic Hydra vulgaris can trigger the nuclear translocation of heterologously expressed human ERα. This led us to hypothesize that aromatic steroids may also be present in the medusozoan cnidarian lineage, which includes Hydra, and may explain the translocation of human ERα. Docking experiments with paraestrol A, a cnidarian A-ring aromatic steroid, into the ligand-binding pocket of Hydra NR3E indicates that, if an aromatic steroid is indeed the true ligand, which remains to be demonstrated, it would bind to the pocket through a partially distinct mechanism from the manner in which estradiol binds to vertebrate ER., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Candida albicans is an opportunistic fungal pathogen that is highly resistant to different oxidative stresses. How reactive sulfur species (RSS) such as sulfite regulate gene expression and the role of the transcription factor Zcf2 and the sulfite exporter Ssu1 in such responses are not known. Here, we show that C. albicans specifically adapts to sulfite stress and that Zcf2 is required for that response as well as induction of genes predicted to remove sulfite from cells and to increase the intracellular amount of a subset of nitrogen metabolites. Analysis of mutants in the sulfate assimilation pathway show that sulfite conversion to sulfide accounts for part of sulfite toxicity and that Zcf2-dependent expression of the SSU1 sulfite exporter is induced by both sulfite and sulfide. Mutations in the SSU1 promoter that selectively inhibit induction by the reactive nitrogen species (RNS) nitrite, a previously reported activator of SSU1 , support a model for C. albicans in which Cta4-dependent RNS induction and Zcf2-dependent RSS induction are mediated by parallel pathways, different from S. cerevisiae in which the transcription factor Fzf1 mediates responses to both RNS and RSS. Lastly, we found that endogenous sulfite production leads to an increase in resistance to exogenously added sulfite. These results demonstrate that C. albicans has a unique response to sulfite that differs from the general oxidative stress response, and that adaptation to internal and external sulfite is largely mediated by one transcription factor and one effector gene., (Copyright © 2017 by the Genetics Society of America.)

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., Competing Interests: The authors have declared that no competing interests exist.

The E6 oncoproteins of high-risk mucosal (hrm) human papillomaviruses (HPVs) contain a pocket that captures LxxLL motifs and a C-terminal motif that recruits PDZ domains, with both functions being crucial for HPV-induced oncogenesis. A chimeric protein was built by fusing a PDZ domain and an LxxLL motif, both known to bind E6. NMR spectroscopy, calorimetry and a mammalian protein complementation assay converged to show that the resulting PDZ-LxxLL chimera is a bivalent nanomolar ligand of E6, while its separated PDZ and LxxLL components are only micromolar binders. The chimera binds to all of the hrm-HPV E6 proteins tested but not to low-risk mucosal or cutaneous HPV E6. Adenovirus-mediated expression of the chimera specifically induces the death of HPV-positive cells, concomitant with increased levels of the tumour suppressor P53, its transcriptional target p21, and the apoptosis marker cleaved caspase 3. The bifunctional PDZ-LxxLL chimera opens new perspectives for the diagnosis and treatment of HPV-induced cancers., (© 2015 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA. This is an open access article under the terms of the Creative Commons Attribution Non-Commercial NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.)

Analysis of an intrinsically disordered protein (IDP) reveals an underlying multifunnel structure for the energy landscape. We suggest that such 'intrinsically disordered' landscapes, with a number of very different competing low-energy structures, are likely to characterise IDPs, and provide a useful way to address their properties. In particular, IDPs are present in many cellular protein interaction networks, and several questions arise regarding how they bind to partners. Are conformations resembling the bound structure selected for binding, or does further folding occur on binding the partner in a induced-fit fashion? We focus on the p53 upregulated modulator of apoptosis (PUMA) protein, which adopts an α-helical conformation when bound to its partner, and is involved in the activation of apoptosis. Recent experimental evidence shows that folding is not necessary for binding, and supports an induced-fit mechanism. Using a variety of computational approaches we deduce the molecular mechanism behind the instability of the PUMA peptide as a helix in isolation. We find significant barriers between partially folded states and the helix. Our results show that the favoured conformations are molten-globule like, stabilised by charged and hydrophobic contacts, with structures resembling the bound state relatively unpopulated in equilibrium.

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

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.

Background: PGC-1α is a crucial regulator of cellular metabolism and energy homeostasis that functionally acts together with the estrogen-related receptors (ERRα and ERRγ) in the regulation of mitochondrial and metabolic gene networks. Dimerization of the ERRs is a pre-requisite for interactions with PGC-1α and other coactivators, eventually leading to transactivation. It was suggested recently (Devarakonda et al) that PGC-1α binds in a strikingly different manner to ERRγ ligand-binding domains (LBDs) compared to its mode of binding to ERRα and other nuclear receptors (NRs), where it interacts directly with the two ERRγ homodimer subunits., Methods/principal Findings: Here, we show that PGC-1α receptor interacting domain (RID) binds in an almost identical manner to ERRα and ERRγ homodimers. Microscale thermophoresis demonstrated that the interactions between PGC-1α RID and ERR LBDs involve a single receptor subunit through high-affinity, ERR-specific L3 and low-affinity L2 interactions. NMR studies further defined the limits of PGC-1α RID that interacts with ERRs. Consistent with these findings, the solution structures of PGC-1α/ERRα LBDs and PGC-1α/ERRγ LBDs complexes share an identical architecture with an asymmetric binding of PGC-1α to homodimeric ERR., Conclusions/significance: These studies provide the molecular determinants for the specificity of interactions between PGC-1α and the ERRs, whereby negative cooperativity prevails in the binding of the coactivators to these receptors. Our work indicates that allosteric regulation may be a general mechanism controlling the binding of the coactivators to homodimers.

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.

Human prion diseases are a heterogeneous group of fatal neurodegenerative disorders, characterized by the deposition of the partially protease-resistant prion protein (PrP(res)), astrocytosis, neuronal loss and spongiform change in the brain. Among inherited forms that represent 15% of patients, different phenotypes have been described depending on the variations detected at different positions within the prion protein gene. Here, we report a new mechanism governing the phenotypic variability of inherited prion diseases. First, we observed that the substitution at residue 211 with either Gln or Asp leads to distinct disorders at the clinical, neuropathological and biochemical levels (Creutzfeldt-Jakob disease or Gerstmann-Sträussler-Scheinker syndrome with abundant amyloid plaques and tau neurofibrillar pathology). Then, using molecular dynamics simulations and biophysical characterization of mutant proteins and an in vitro model of PrP conversion, we found evidence that each substitution impacts differently the stability of PrP and its propensity to produce different protease resistant fragments that may contribute to the phenotypical switch. Thus, subtle differences in the PrP primary structure and stability are sufficient to control amyloid plaques formation and tau abnormal phosphorylation and fibrillation. This mechanism is unique among neurodegenerative disorders and is consistent with the prion hypothesis that proposes a conformational change as the key pathological event in prion disorders.

Actual use of the active form of vitamin D (calcitriol or 1α,25-dihydroxyvitamin D(3)) to treat hyperproliferative disorders is hampered by calcemic effects, hence the continuous development of chemically modified analogues with dissociated profiles. Structurally distinct nonsecosteroidal analogues have been developed to mimic calcitriol activity profiles with low calcium serum levels. Here, we report the crystallographic study of vitamin D nuclear receptor (VDR) ligand binding domain in complexes with six nonsecosteroidal analogues harboring two or three phenyl rings. These compounds induce a stimulated transcription in the nanomolar range, similar to calcitriol. Examination of the protein-ligand interactions reveals the mode of binding of these nonsecosteroidal compounds and highlights the role of the various chemical modifications of the ligands to VDR binding and activity, notably (de)solvation effects. The structures with the tris-aromatic ligands exhibit a rearrangement of a novel region of the VDR ligand binding pocket, helix H6.

Coarse-grained protein models with various levels of granularity and degrees of freedom offer the possibility to explore many phenomena including folding, assembly, and recognition in terms of dynamics and thermodynamics that are inaccessible to all-atom representations in explicit aqueous solution. Here, we present a refined version of the coarse-grained optimized potential for efficient protein structure prediction (OPEP) based on a six-bead representation. The OPEP version 4.0 parameter set, which uses a new analytical formulation for the nonbonded interactions and adds specific side-chain-side-chain interactions for α-helix, is subjected to three tests. First, we show that molecular dynamics simulations at 300 K preserve the experimental rigid conformations of 17 proteins with 37-152 amino acids within a root-mean-square deviation (RMSD) of 3.1 Å after 30 ns. Extending the simulation time to 100 ns for five proteins does not change the RMSDs. Second, replica exchange molecular dynamics (REMD) simulations recover the NMR structures of three prototypical β-hairpin and α-helix peptides and the NMR three-stranded β-sheet topology of a 37-residue WW domain, starting from randomly chosen states. Third, REMD simulations on the ccβ peptide show a temperature transition from a three-stranded coiled coil to amyloid-like aggregates consistent with experiments, while simulations on low molecular weight aggregates of the prion protein helix 1 do not. Overall, these studies indicate the effectiveness of our OPEP4 coarse-grained model for protein folding and aggregation, and report two future directions for improvement.

The amyloid-β protein (Aβ) oligomers are believed to be the main culprits in the cytoxicity of Alzheimer's disease (AD) and p3 peptides (Aβ17-42 fragments) are present in AD amyloid plaques. Many small-molecule or peptide-based inhibitors are known to slow down Aβ aggregation and reduce the toxicity in vitro, but their exact modes of action remain to be determined since there has been no atomic level of Aβ(p3)-drug oligomers. In this study, we have determined the structure of Aβ17-42 trimers both in aqueous solution and in the presence of five small-molecule inhibitors using a multiscale computational study. These inhibitors include 2002-H20, curcumin, EGCG, Nqtrp, and resveratrol. First, we used replica exchange molecular dynamics simulations coupled to the coarse-grained (CG) OPEP force field. These CG simulations reveal that the conformational ensemble of Aβ17-42 trimer can be described by 14 clusters with each peptide essentially adopting turn/random coil configurations, although the most populated cluster is characterized by one peptide with a β-hairpin at Phe19-Leu31. Second, these 14 dominant clusters and the less-frequent fibril-like state with parallel register of the peptides were subjected to atomistic Autodock simulations. Our analysis reveals that the drugs have multiple binding modes with different binding affinities for trimeric Aβ17-42 although they interact preferentially with the CHC region (residues 17-21). The compounds 2002-H20 and Nqtrp are found to be the worst and best binders, respectively, suggesting that the drugs may interfere at different stages of Aβ oligomerization. Finally, explicit solvent molecular dynamics of two predicted Nqtrp-Aβ17-42 conformations describe at atomic level some possible modes of action for Nqtrp.

Many proteins display a strand-loop-strand motif in their amyloid fibrillar states. For instance, the amyloid beta-protein, Abeta1-40, associated with Alzheimer's disease, displays a loop at positions 22-28 in its amyloid fibril state. It has been suggested that this loop could appear early in the aggregation process, but quantitative information regarding its presence in small oligomers remains scant. Because residues 1-15 are disordered in Abeta1-42 fibrils and Abeta10-35 forms fibrils in vitro, we select the peptide Abeta16-35, centered on residues 22-28 and determine the structures and thermodynamics of the monomer and dimer using coarse-grained implicit solvent replica exchange molecular dynamics simulations. Our simulations totalling 5 mus for the monomer and 12 micros for the dimer show no sign of strong secondary structure signals in both instances and the significant impact of dimerization on the global structure of Abeta16-35. They reveal however that the loop 22-28 acts as a quasi-independent unit in both species. The loop structure ensemble we report in Abeta16-35 monomer and dimer has high similarity to the loop formed by the Abeta21-30 peptide in solution and, to a lesser extent, to the loop found in Abeta1-40 fibrils. We discuss the implications of our findings on the assembly of full-length Abeta.

The conversion of the prion protein (PrP) from its cellular form, PrPC, to its pathogenic scrapie form, PrPSc, is a key event in neurodegenerative transmissible spongiform encephalopathies such as Creutzfeldt-Jakob disease (CJD). PrPC is characterized by three helices (H1-H3) and a small antiparallel beta-sheet. One working hypothesis for TSE causation is that oligomeric forms of PrP are the proximate neurotoxic agents. Because these states are transient in character, current experimental studies have failed to provide atomic structures. To gain insights into these intermediates, we have studied PrP125-228 and its CJD-causing T183A variant in their monomer and dimer forms by means of coarse-grained protein molecular dynamics simulations. Our 1.5 microsecond simulations show that the decrease in the thermodynamic stability of PrP monomer upon T183A, consistent with experimental studies, results from a destabilization of the H2H3 subdomain. Comparison of the monomer and dimer properties from wild-type and T183A PrP reveals that helix H1 is robust and the H2H3 subdomain displays a much higher propensity for intra- and inter-beta-sheets in T183A than in the wild-type sequence under denaturing conditions. However, both species display negligible beta-sheet structure. Implications of our simulations on prion propagation are discussed.

Aggregation of the Abeta1-40/Abeta1-42 peptides is a key factor in Alzheimer's disease. Though the inhibitory effect of N-methylated Abeta16-22 (mAbeta16-22) peptides is well characterized in vitro, there is little information on how they disassemble full-length Abeta fibrils or block fibril formation. Here, we report coarse-grained implicit solvent molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulations on Abeta16-22 and mAbeta16-22 monomers, and then a preformed six-chain Abeta16-22 bilayer with either four copies of Abeta16-22 or four copies of mAbeta16-22. Our simulations show that the effect of N-methylation on mAbeta16-22 monomer is to reduce the density of compact forms. While 100 ns MD trajectories do not reveal any significant differences between the two ten-chain systems, the REMD simulations totaling 1 micros help understand the first steps of Abeta16-22 protofibril disassembly by N-methylated inhibitors. Notably, we find that mAbeta16-22 preferentially interacts with Abeta16-22 by blocking both beta-sheet extension and lateral association of layers, but also by intercalation of the inhibitors allowing sequestration of Abeta16-22 peptides. This third binding mode is particularly appealing for blocking Abeta fibrillogenesis.

Current approaches aimed at determining the free energy surface of all-atom medium-size proteins in explicit solvent are slow and are not sufficient to converge to equilibrium properties. To ensure a proper sampling of the configurational space, it is preferable to use reduced representations such as implicit solvent and/or coarse-grained protein models, which are much lighter computationally. Each model must be verified, however, to ensure that it can recover experimental structures and thermodynamics. Here we test the coarse-grained implicit solvent OPEP model with replica exchange molecular dynamics (REMD) on six peptides ranging in length from 10 to 28 residues: two alanine-based peptides, the second beta-hairpin from protein G, the Trp-cage and zinc-finger motif, and a dimer of a coiled coil peptide. We show that REMD-OPEP recovers the proper thermodynamics of the systems studied, with accurate structural description of the beta-hairpin and Trp-cage peptides (within 1-2 A from experiments). The light computational burden of REMD-OPEP, which enables us to generate many hundred nanoseconds at each temperature and fully assess convergence to equilibrium ensemble, opens the door to the determination of the free energy surface of larger proteins and assemblies.