Conformations of Macrocyclic Peptides Sampled by Nuclear Magnetic Resonance: Models for Cell-Permeability.

The biological activities and pharmacological properties of peptides and peptide mimetics are determined by their conformational states. Therefore, a detailed understanding of the conformational landscape is crucial for rational drug design. Nuclear magnetic resonance (NMR) is the only method for structure determination in solution. However, it remains challenging to determine the structures of peptides using NMR because of very weak nuclear Overhauser effects (NOEs), the semiquantitative nature of the rotating frame Overhauser effect (ROE), and the low number of NOEs/ROEs in N-methylated peptides. In this study, we introduce a new approach to investigating the structures of modified macrocyclic peptides. We utilize exact NOEs (eNOEs) in viscous solvent mixtures to replicate various cellular environments. eNOEs provide detailed structural information for highly dynamic modified peptides. Structures of high precision were obtained for cyclosporin A, with a backbone atom rmsd of 0.10 Å. Distinct conformational states in different environments were identified for omphalotin A (OmphA), a fungal nematotoxic and multiple backbone N-methylated macrocyclic peptides. A model for cell-permeation is presented for OmphA, based on its structures in polar, apolar, and mixed polarity solvents. During the transition from a polar to an apolar environment, OmphA undergoes a rearrangement of its H-bonding network, accompanied by a cis to trans isomerization of the ω torsion angle within a type VIa ß-turn. We hypothesize that the kinetics of these conformational transitions play a crucial role in determining the membrane-permeation capabilities of OmphA.

An Analysis of Nucleotide-Amyloid Interactions Reveals Selective Binding to Codon-Sized RNA.

Interactions between RNA and proteins are the cornerstone of many important biological processes from transcription and translation to gene regulation, yet little is known about the ancient origin of said interactions. We hypothesized that peptide amyloids played a role in the origin of life and that their repetitive structure lends itself to building interfaces with other polymers through avidity. Here, we report that short RNA with a minimum length of three nucleotides binds in a sequence-dependent manner to peptide amyloids. The 3'-5' linked RNA backbone appears to be well-suited to support these interactions, with the phosphodiester backbone and nucleobases both contributing to the affinity. Sequence-specific RNA-peptide interactions of the kind identified here may provide a path to understanding one of the great mysteries rooted in the origin of life: the origin of the genetic code.

Fibrils Emerging from Droplets: Molecular Guiding Principles behind Phase Transitions of a Short Peptide-Based Condensate Studied by Solid-State NMR.

Biochemical reactions occurring in highly crowded cellular environments require different means of control to ensure productivity and specificity. Compartmentalization of reagents by liquid-liquid phase separation is one of these means. However, extremely high local protein concentrations of up to 400 mg/ml can result in pathological aggregation into fibrillar amyloid structures, a phenomenon that has been linked to various neurodegenerative diseases. Despite its relevance, the process of liquid-to-solid transition inside condensates is still not well understood at the molecular level. We thus herein use small peptide derivatives that can undergo both liquid-liquid and subsequent liquid-to-solid phase transition as model systems to study both processes. Using solid-state nuclear magnetic resonance (NMR) and transmission electron microscopy (TEM), we compare the structure of condensed states of leucine, tryptophan and phenylalanine containing derivatives, distinguishing between liquid-like condensates, amorphous aggregates and fibrils, respectively. A structural model for the fibrils formed by the phenylalanine derivative was obtained by an NMR-based structure calculation. The fibrils are stabilised by hydrogen bonds and side-chain π-π interactions, which are likely much less pronounced or absent in the liquid and amorphous state. Such noncovalent interactions are equally important for the liquid-to-solid transition of proteins, particularly those related to neurodegenerative diseases.

Atomic resolution protein allostery from the multi-state structure of a PDZ domain.

Recent methodological advances in solution NMR allow the determination of multi-state protein structures and provide insights into structurally and dynamically correlated protein sites at atomic resolution. This is demonstrated in the present work for the well-studied PDZ2 domain of protein human tyrosine phosphatase 1E for which protein allostery had been predicted. Two-state protein structures were calculated for both the free form and in complex with the RA-GEF2 peptide using the exact nuclear Overhauser effect (eNOE) method. In the apo protein, an allosteric conformational selection step comprising almost 60% of the domain was detected with an "open" ligand welcoming state and a "closed" state that obstructs the binding site by changing the distance between the ß-sheet 2, α-helix 2, and sidechains of residues Lys38 and Lys72. The observed induced fit-type apo-holo structural rearrangements are in line with the previously published evolution-based analysis covering ~25% of the domain with only a partial overlap with the protein allostery of the open form. These presented structural studies highlight the presence of a dedicated highly optimized and complex dynamic interplay of the PDZ2 domain owed by the structure-dynamics landscape.

A B-factor for NOEs?

Nuclear Overhauser effects (NOEs) are influenced by motion. Here, we derive exact, analytical results for a model of isotropic, harmonic fluctuations of atom positions that corresponds to the one underlying crystallographic B-factors. The model includes steric repulsion and yields closed-form expressions for the expected value of general invertible functions of the distance between two atoms, with the special case r-6 for NOEs. We discuss the implications for the definition of an NOE-based B-factor in solution NMR.

On the Entropy of a One-Dimensional Gas with and without Mixing Using Sinai Billiard.

A one-dimensional gas comprising N point particles undergoing elastic collisions within a finite space described by a Sinai billiard generating identical dynamical trajectories are calculated and analyzed with regard to strict extensivity of the entropy definitions of Boltzmann-Gibbs. Due to the collisions, trajectories of gas particles are strongly correlated and exhibit both chaotic and periodic properties. Probability distributions for the position of each particle in the one-dimensional gas can be obtained analytically, elucidating that the entropy in this special case is extensive at any given number N. Furthermore, the entropy obtained can be interpreted as a measure of the extent of interactions between molecules. The results obtained for the non-mixable one-dimensional system are generalized to mixable one- and two-dimensional systems, the latter by a simple example only providing similar findings.

1H, 13C and 15N resonance assignment of the YTH domain of YTHDC2.

In humans, YTH (YT521-B homology) domain containing protein 2 (YTHDC2) plays a crucial role in the phase-shift from mitosis to meiosis. YTH domains bind to methylated adenosine nucleotides such as m6A. In a phylogenic tree, the YTH domain of YTHDC2 (YTH2) and that of the YTH containing protein YTHDC1 (YTH1) belong to the same sub-group. However, the binding affinity of m6A differs between these proteins. Here, we report 1H, 13C and 15N resonance assignment of YTH2 and its solution structure to examine the difference of the structural architecture and the dynamic properties of YTH1 and YTH2. YTH2 adopts a ß1-α1-ß2-α2-ß3-ß4-ß5-α3-ß6-α4 topology, which was also observed in YTH1. However, the ß4-ß5 loops of YTH1 and YTH2 are distinct in length and amino acid composition. Our data revealed that, unlike in YTH1, the structure of m6A-binding pocket of YTH2 formed by the ß4-ß5 loop is stabilized by electrostatic interaction. This assignment and the structural information for YTH2 will provide the insight on the further functional research of YTHDC2.

An atypical LIR motif within UBA5 (ubiquitin like modifier activating enzyme 5) interacts with GABARAP proteins and mediates membrane localization of UBA5.

Short linear motifs, known as LC3-interacting regions (LIRs), interact with mactoautophagy/autophagy modifiers (Atg8/LC3/GABARAP proteins) via a conserved universal mechanism. Typically, this includes the occupancy of 2 hydrophobic pockets on the surface of Atg8-family proteins by 2 specific aromatic and hydrophobic residues within the LIR motifs. Here, we describe an alternative mechanism of Atg8-family protein interaction with the non-canonical UBA5 LIR, an E1-like enzyme of the ufmylation pathway that preferentially interacts with GABARAP but not LC3 proteins. By solving the structures of both GABARAP and GABARAPL2 in complex with the UBA5 LIR, we show that in addition to the binding to the 2 canonical hydrophobic pockets (HP1 and HP2), a conserved tryptophan residue N-terminal of the LIR core sequence binds into a novel hydrophobic pocket on the surface of GABARAP proteins, which we term HP0. This mode of action is unique for UBA5 and accompanied by large rearrangements of key residues including the side chains of the gate-keeping K46 and the adjacent K/R47 in GABARAP proteins. Swapping mutations in LC3B and GABARAPL2 revealed that K/R47 is the key residue in the specific binding of GABARAP proteins to UBA5, with synergetic contributions of the composition and dynamics of the loop L3. Finally, we elucidate the physiological relevance of the interaction and show that GABARAP proteins regulate the localization and function of UBA5 on the endoplasmic reticulum membrane in a lipidation-independent manner.Abbreviations: ATG: AuTophaGy-related; EGFP: enhanced green fluorescent protein; GABARAP: GABA-type A receptor-associated protein; ITC: isothermal titration calorimetry; KO: knockout; LIR: LC3-interacting region; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; NMR: nuclear magnetic resonance; RMSD: root-mean-square deviation of atomic positions; TKO: triple knockout; UBA5: ubiquitin like modifier activating enzyme 5.

The Solution Structure and Dynamics of Cd-Metallothionein from Helix pomatia Reveal Optimization for Binding Cd over Zn.

Metallothioneins (MTs) are cysteine-rich polypeptides that are naturally found coordinated to monovalent and/or divalent transition metal ions. Three metallothionein isoforms from the Roman snail Helix pomatia are known. They differ in their physiological metal load and in their specificity for transition metal ions such as Cd2+ (HpCdMT isoform) and Cu+ (HpCuMT isoform) or in the absence of a defined metal specificity (HpCd/CuMT isoform). We have determined the solution structure of the Cd-specific isoform (HpCdMT) by nuclear magnetic resonance spectroscopy using recombinant isotopically labeled protein loaded with Zn2+ or Cd2+. Both structures display two-domain architectures, where each domain comprises a characteristic three-metal cluster similar to that observed in the ß-domains of vertebrate MTs. The polypeptide backbone is well-structured over the entire sequence, including the interdomain linker. Interestingly, the two domains display mutual contacts, as observed before for the metallothionein of the snail Littorina littorea, to which both N- and C-terminal domains are highly similar. Increasing the length of the linker motionally decouples both domains and removes mutual contacts between them without having a strong effect on the stability of the individual domains. The structures of Cd6- and Zn6-HpCdMT are nearly identical. However, 15N relaxation, in particular 15N R2 rates, is accelerated for many residues of Zn6-HpCdMT but not for Cd6-HpCdMT, revealing the presence of conformational exchange effects. We suggest that this snail MT isoform is evolutionarily optimized for binding Cd rather than Zn.

High-Resolution Protein 3D Structure Determination in Living Eukaryotic Cells.

Proteins in living cells interact specifically or nonspecifically with an enormous number of biomolecules. To understand the behavior of proteins under intracellular crowding conditions, it is indispensable to observe their three-dimensional (3D) structures at the atomic level in a physiologically natural environment. We demonstrate the first de novo protein structure determinations in eukaryotes with the sf9 cell/baculovirus system using NMR data from living cells exclusively. The method was applied to five proteins, rat calmodulin, human HRas, human ubiquitin, T. thermophilus HB8 TTHA1718, and Streptococcus protein G B1 domain. In all cases, we could obtain structural information from well-resolved in-cell 3D nuclear Overhauser effect spectroscopy (NOESY) data, suggesting that our method can be a standard tool for protein structure determinations in living eukaryotic cells. For three proteins, we achieved well-converged 3D structures. Among these, the in-cell structure of protein G B1 domain was most accurately determined, demonstrating that a helix-loop region is tilted away from a ß-sheet compared to the conformation in diluted solution.

Proteome-wide analysis of phospho-regulated PDZ domain interactions.

A key function of reversible protein phosphorylation is to regulate protein-protein interactions, many of which involve short linear motifs (3-12 amino acids). Motif-based interactions are difficult to capture because of their often low-to-moderate affinities. Here, we describe phosphomimetic proteomic peptide-phage display, a powerful method for simultaneously finding motif-based interaction and pinpointing phosphorylation switches. We computationally designed an oligonucleotide library encoding human C-terminal peptides containing known or predicted Ser/Thr phosphosites and phosphomimetic variants thereof. We incorporated these oligonucleotides into a phage library and screened the PDZ (PSD-95/Dlg/ZO-1) domains of Scribble and DLG1 for interactions potentially enabled or disabled by ligand phosphorylation. We identified known and novel binders and characterized selected interactions through microscale thermophoresis, isothermal titration calorimetry, and NMR We uncover site-specific phospho-regulation of PDZ domain interactions, provide a structural framework for how PDZ domains accomplish phosphopeptide binding, and discuss ligand phosphorylation as a switching mechanism of PDZ domain interactions. The approach is readily scalable and can be used to explore the potential phospho-regulation of motif-based interactions on a large scale.

Regulation of the Activity in the p53 Family Depends on the Organization of the Transactivation Domain.

Despite high sequence homology among the p53 family members, the regulation of their transactivation potential is based on strikingly different mechanisms. Previous studies revealed that the activity of TAp63α is regulated via an autoinhibitory mechanism that keeps inactive TAp63α in a dimeric conformation. While all p73 isoforms are constitutive tetramers, their basal activity is much lower compared with tetrameric TAp63. We show that the dimeric state of TAp63α not only reduces DNA binding affinity, but also suppresses interaction with the acetyltransferase p300. Exchange of the transactivation domains is sufficient to transfer the regulatory characteristics between p63 and p73. Structure determination of the transactivation domains of p63 and p73 in complex with the p300 Taz2 domain further revealed that, in contrast to p53 and p73, p63 has a single transactivation domain. Sequences essential for stabilizing the closed dimer of TAp63α have evolved into a second transactivation domain in p73 and p53.

Protein aggregation of the p63 transcription factor underlies severe skin fragility in AEC syndrome.

The p63 gene encodes a master regulator of epidermal commitment, development, and differentiation. Heterozygous mutations in the C-terminal domain of the p63 gene can cause ankyloblepharon-ectodermal defects-cleft lip/palate (AEC) syndrome, a life-threatening disorder characterized by skin fragility and severe, long-lasting skin erosions. Despite deep knowledge of p63 functions, little is known about mechanisms underlying disease pathology and possible treatments. Here, we show that multiple AEC-associated p63 mutations, but not those causative of other diseases, lead to thermodynamic protein destabilization, misfolding, and aggregation, similar to the known p53 gain-of-function mutants found in cancer. AEC mutant proteins exhibit impaired DNA binding and transcriptional activity, leading to dominant negative effects due to coaggregation with wild-type p63 and p73. Importantly, p63 aggregation occurs also in a conditional knock-in mouse model for the disorder, in which the misfolded p63 mutant protein leads to severe epidermal defects. Variants of p63 that abolish aggregation of the mutant proteins are able to rescue p63's transcriptional function in reporter assays as well as in a human fibroblast-to-keratinocyte conversion assay. Our studies reveal that AEC syndrome is a protein aggregation disorder and opens avenues for therapeutic intervention.

Structural and functional dissection of the DH and PH domains of oncogenic Bcr-Abl tyrosine kinase.

The two isoforms of the Bcr-Abl tyrosine kinase, p210 and p190, are associated with different leukemias and have a dramatically different signaling network, despite similar kinase activity. To provide a molecular rationale for these observations, we study the Dbl-homology (DH) and Pleckstrin-homology (PH) domains of Bcr-Abl p210, which constitute the only structural differences to p190. Here we report high-resolution structures of the DH and PH domains and characterize conformations of the DH-PH unit in solution. Our structural and functional analyses show no evidence that the DH domain acts as a guanine nucleotide exchange factor, whereas the PH domain binds to various phosphatidylinositol-phosphates. PH-domain mutants alter subcellular localization and result in decreased interactions with p210-selective interaction partners. Hence, the PH domain, but not the DH domain, plays an important role in the formation of the differential p210 and p190 Bcr-Abl signaling networks.

NMR Investigation of Structures of G-protein Coupled Receptor Folding Intermediates.

Folding of G-protein coupled receptors (GPCRs) according to the two-stage model (Popot, J. L., and Engelman, D. M. (1990) Biochemistry 29, 4031-4037) is postulated to proceed in 2 steps: partitioning of the polypeptide into the membrane followed by diffusion until native contacts are formed. Herein we investigate conformational preferences of fragments of the yeast Ste2p receptor using NMR. Constructs comprising the first, the first two, and the first three transmembrane (TM) segments, as well as a construct comprising TM1-TM2 covalently linked to TM7 were examined. We observed that the isolated TM1 does not form a stable helix nor does it integrate well into the micelle. TM1 is significantly stabilized upon interaction with TM2, forming a helical hairpin reported previously (Neumoin, A., Cohen, L. S., Arshava, B., Tantry, S., Becker, J. M., Zerbe, O., and Naider, F. (2009) Biophys. J. 96, 3187-3196), and in this case the protein integrates into the hydrophobic interior of the micelle. TM123 displays a strong tendency to oligomerize, but hydrogen exchange data reveal that the center of TM3 is solvent exposed. In all GPCRs so-far structurally characterized TM7 forms many contacts with TM1 and TM2. In our study TM127 integrates well into the hydrophobic environment, but TM7 does not stably pack against the remaining helices. Topology mapping in microsomal membranes also indicates that TM1 does not integrate in a membrane-spanning fashion, but that TM12, TM123, and TM127 adopt predominantly native-like topologies. The data from our study would be consistent with the retention of individual helices of incompletely synthesized GPCRs in the vicinity of the translocon until the complete receptor is released into the membrane interior.

Mechanism of TAp73 inhibition by ΔNp63 and structural basis of p63/p73 hetero-tetramerization.

Members of the p53 tumor-suppressor family are expressed as multiple isoforms. Isoforms with an N-terminal transactivation domain are transcriptionally active, while those ones lacking this domain often inhibit the transcriptional activity of other family members. In squamous cell carcinomas, the high expression level of ΔNp63α inhibits the tumor-suppressor function of TAp73ß. This can in principle be due to blocking of the promoter or by direct interaction between both proteins. p63 and p73 can hetero-oligomerize through their tetramerization domains and a hetero-tetramer consisting of two p63 and two p73 molecules is thermodynamically more stable than both homo-tetramers. Here we show that cells expressing both p63 and p73 exist in mouse epidermis and hair follicle and that hetero-tetramer complexes can be detected by immunoprecipitation in differentiating keratinocytes. Through structure determination of the hetero-tetramer, we reveal why this hetero-tetramer is the thermodynamically preferred species. We have created mutants that exclusively form either hetero-tetramers or homo-tetramers, allowing to investigate the function of these p63/p73 hetero-tetramers. Using these tools, we show that inhibition of TAp73ß in squamous cell carcinomas is due to promoter squelching and not direct interaction.

NMR structure calculation for all small molecule ligands and non-standard residues from the PDB Chemical Component Dictionary.

An algorithm, CYLIB, is presented for converting molecular topology descriptions from the PDB Chemical Component Dictionary into CYANA residue library entries. The CYANA structure calculation algorithm uses torsion angle molecular dynamics for the efficient computation of three-dimensional structures from NMR-derived restraints. For this, the molecules have to be represented in torsion angle space with rotations around covalent single bonds as the only degrees of freedom. The molecule must be given a tree structure of torsion angles connecting rigid units composed of one or several atoms with fixed relative positions. Setting up CYANA residue library entries therefore involves, besides straightforward format conversion, the non-trivial step of defining a suitable tree structure of torsion angles, and to re-order the atoms in a way that is compatible with this tree structure. This can be done manually for small numbers of ligands but the process is time-consuming and error-prone. An automated method is necessary in order to handle the large number of different potential ligand molecules to be studied in drug design projects. Here, we present an algorithm for this purpose, and show that CYANA structure calculations can be performed with almost all small molecule ligands and non-standard amino acid residues in the PDB Chemical Component Dictionary.

Structural basis for phosphorylation-triggered autophagic clearance of Salmonella.

Selective autophagy is mediated by the interaction of autophagy modifiers and autophagy receptors that also bind to ubiquitinated cargo. Optineurin is an autophagy receptor that plays a role in the clearance of cytosolic Salmonella. The interaction between receptors and modifiers is often relatively weak, with typical values for the dissociation constant in the low micromolar range. The interaction of optineurin with autophagy modifiers is even weaker, but can be significantly enhanced through phosphorylation by the TBK1 {TANK [TRAF (tumour-necrosis-factor-receptor-associated factor)-associated nuclear factor κB activator]-binding kinase 1}. In the present study we describe the NMR and crystal structures of the autophagy modifier LC3B (microtubule-associated protein light chain 3 beta) in complex with the LC3 interaction region of optineurin either phosphorylated or bearing phospho-mimicking mutations. The structures show that the negative charge induced by phosphorylation is recognized by the side chains of Arg¹¹ and Lys5¹ in LC3B. Further mutational analysis suggests that the replacement of the canonical tryptophan residue side chain of autophagy receptors with the smaller phenylalanine side chain in optineurin significantly weakens its interaction with the autophagy modifier LC3B. Through phosphorylation of serine residues directly N-terminally located to the phenylalanine residue, the affinity is increased to the level normally seen for receptor-modifier interactions. Phosphorylation, therefore, acts as a switch for optineurin-based selective autophagy.

Spatial elucidation of motion in proteins by ensemble-based structure calculation using exact NOEs.

Proteins are inherently dynamic systems whose motions cover large ranges in both magnitude and timescale. Because of the omnipresence of motion, it is likely that dynamics have important roles in the function of biomolecules. For detailed understanding of a protein's function, the three-dimensional structure and description of its dynamics are therefore required. Structure determination methods are well established, and NMR-relaxation phenomena provide insights into local molecular dynamics; moreover, recently several attempts have been made to detect concerted motion. Here, we present an ensemble-based structure-determination protocol using ensemble-averaged distance restraints obtained from exact NOE rates. Application to the model protein GB3 establishes an ensemble of structures that reveals correlated motion across the ß-sheet, concerted motion between the backbone and side chains localized in the structure core, and a lack of concerted conformational exchange between the ß-sheet and the α-helix.