A structural signature motif enlightens the origin and diversification of nuclear receptors.

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.

High-Risk Mucosal Human Papillomavirus 16 (HPV16) E6 Protein and Cutaneous HPV5 and HPV8 E6 Proteins Employ Distinct Strategies To Interfere with Interferon Regulatory Factor 3-Mediated Beta Interferon Expression.

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.

Vitamin D Analogs Bearing C-20 Modifications Stabilize the Agonistic Conformation of Non-Responsive Vitamin D Receptor Variants.

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.

A multifunnel energy landscape encodes the competing α-helix and ß-hairpin conformations for a designed peptide.

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.

Crucial role of nonspecific interactions in amyloid nucleation.

Protein oligomers have been implicated as toxic agents in a wide range of amyloid-related diseases. However, it has remained unsolved whether the oligomers are a necessary step in the formation of amyloid fibrils or just a dangerous byproduct. Analogously, it has not been resolved if the amyloid nucleation process is a classical one-step nucleation process or a two-step process involving prenucleation clusters. We use coarse-grained computer simulations to study the effect of nonspecific attractions between peptides on the primary nucleation process underlying amyloid fibrillization. We find that, for peptides that do not attract, the classical one-step nucleation mechanism is possible but only at nonphysiologically high peptide concentrations. At low peptide concentrations, which mimic the physiologically relevant regime, attractive interpeptide interactions are essential for fibril formation. Nucleation then inevitably takes place through a two-step mechanism involving prefibrillar oligomers. We show that oligomers not only help peptides meet each other but also, create an environment that facilitates the conversion of monomers into the ß-sheet-rich form characteristic of fibrils. Nucleation typically does not proceed through the most prevalent oligomers but through an oligomer size that is only observed in rare fluctuations, which is why such aggregates might be hard to capture experimentally. Finally, we find that the nucleation of amyloid fibrils cannot be described by classical nucleation theory: in the two-step mechanism, the critical nucleus size increases with increases in both concentration and interpeptide interactions, which is in direct contrast with predictions from classical nucleation theory.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems.

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.

Targeting the Two Oncogenic Functional Sites of the HPV E6 Oncoprotein with a High-Affinity Bivalent Ligand.

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.

Substitutions at residue 211 in the prion protein drive a switch between CJD and GSS syndrome, a new mechanism governing inherited neurodegenerative disorders.

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.

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.

BpForms and BcForms: a toolkit for concretely describing non-canonical polymers and complexes to facilitate global biochemical networks.

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.

Targeting the early steps of Abeta16-22 protofibril disassembly by N-methylated inhibitors: a numerical study.

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.

The conversion of helix H2 to beta-sheet is accelerated in the monomer and dimer of the prion protein upon T183A mutation.

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.

Structures and thermodynamics of Alzheimer's amyloid-beta Abeta(16-35) monomer and dimer by replica exchange molecular dynamics simulations: implication for full-length Abeta fibrillation.

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.

Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent.

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.

Structural Basis for DNA Gyrase Interaction with Coumermycin A1.

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.

Modulation of RXR-DNA complex assembly by DNA context.

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.

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.

NR3E receptors in cnidarians: A new family of steroid receptor relatives extends the possible mechanisms for ligand binding.

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.

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.