RNA binding induces an allosteric switch in Cyp33 to repress MLL1-mediated transcription.

Mixed-lineage leukemia 1 (MLL1) is a transcription activator of the HOX family, which binds to specific epigenetic marks on histone H3 through its third plant homeodomain (PHD3) domain. Through an unknown mechanism, MLL1 activity is repressed by cyclophilin 33 (Cyp33), which binds to MLL1 PHD3. We determined solution structures of Cyp33 RNA recognition motif (RRM) free, bound to RNA, to MLL1 PHD3, and to both MLL1 and the histone H3 lysine N6-trimethylated. We found that a conserved α helix, amino-terminal to the RRM domain, adopts three different positions facilitating a cascade of binding events. These conformational changes are triggered by Cyp33 RNA binding and ultimately lead to MLL1 release from the histone mark. Together, our mechanistic findings rationalize how Cyp33 binding to MLL1 can switch chromatin to a transcriptional repressive state triggered by RNA binding as a negative feedback loop.

Paradoxical Increase of Permeability and Lipophilicity with the Increasing Topological Polar Surface Area within a Series of PRMT5 Inhibitors.

An imidazolone → triazolone replacement addressed the limited passive permeability of a series of protein arginine methyl transferase 5 (PRMT5) inhibitors. This increase in passive permeability was unexpected given the increase in the hydrogen bond acceptor (HBA) count and topological polar surface area (TPSA), two descriptors that are typically inversely correlated with permeability. Quantum mechanics (QM) calculations revealed that this unusual effect was due to an electronically driven disconnect between TPSA and 3D-PSA, which manifests in a reduction in overall HBA strength as indicated by the HBA moment descriptor from COSMO-RS (conductor-like screening model for real solvation). HBA moment was subsequently deployed as a design parameter leading to the discovery of inhibitors with not only improved passive permeability but also reduced P-glycoprotein (P-gp) transport. Our case study suggests that hidden polarity as quantified by TPSA-3DPSA can be rationally designed through QM calculations.

Optimization of the In Vivo Potency of Pyrazolopyrimidine MALT1 Protease Inhibitors by Reducing Metabolism and Increasing Potency in Whole Blood.

The paracaspase MALT1 has gained increasing interest as a target for the treatment of subsets of lymphomas as well as autoimmune diseases, and there is a need for suitable compounds to explore the therapeutic potential of this target. Here, we report the optimization of the in vivo potency of pyrazolopyrimidines, a class of highly selective allosteric MALT1 inhibitors. High doses of the initial lead compound led to tumor stasis in an activated B-cell-like (ABC) diffuse large B-cell lymphoma (DLBCL) xenograft model, but this compound suffered from a short in vivo half-life and suboptimal potency in whole blood. Guided by metabolism studies, we identified compounds with reduced metabolic clearance and increased in vivo half-life. In the second optimization step, masking one of the hydrogen-bond donors of the central urea moiety through an intramolecular interaction led to improved potency in whole blood. This was associated with improved in vivo potency in a mechanistic model of B cell activation. The optimized compound led to tumor regression in a CARD11 mutant ABC-DLBCL lymphoma xenograft model.

Morphing of Amphipathic Helices to Explore the Activity and Selectivity of Membranolytic Antimicrobial Peptides.

Naturally occurring membranolytic antimicrobial peptides (AMPs) are rarely cell-type selective and highly potent at the same time. Template-based peptide design can be used to generate AMPs with improved properties de novo. Following this approach, 18 linear peptides were obtained by computationally morphing the natural AMP Aurein 2.2d2 GLFDIVKKVVGALG into the synthetic model AMP KLLKLLKKLLKLLK. Eleven of the 18 chimeric designs inhibited the growth of Staphylococcus aureus, and six peptides were tested and found to be active against one resistant pathogenic strain or more. One of the peptides was broadly active against bacterial and fungal pathogens without exhibiting toxicity to certain human cell lines. Solution nuclear magnetic resonance and molecular dynamics simulation suggested an oblique-oriented membrane insertion mechanism of this helical de novo peptide. Temperature-resolved circular dichroism spectroscopy pointed to conformational flexibility as an essential feature of cell-type selective AMPs.

Rationalization of the Membrane Permeability Differences in a Series of Analogue Cyclic Decapeptides.

Cyclization and selected backbone N-methylations are found to be often necessary but not sufficient conditions for peptidic drugs to have a good bioavailability. Thus, the design of cyclic peptides with good passive membrane permeability and good solubility remains a challenge. The backbone scaffold of a recently published series of cyclic decapeptides with six selected backbone N-methylations was designed to favor the adoption of a closed conformation with ß-turns and four transannular hydrogen bonds. Although this conformation was indeed adopted by the peptides as determined by NMR measurements, substantial differences in the membrane permeability were observed. In this work, we aim to rationalize the impact of discrete side chain modifications on membrane permeability for six of these cyclic decapeptides. The thermodynamic and kinetic properties were investigated using molecular dynamics simulations and Markov state modeling in water and chloroform. The study highlights the influence that side-chain modifications can have on the backbone conformation. Peptides with a d-proline in the ß-turns were more likely to adopt, even in water, the closed conformation with transannular hydrogen bonds, which facilitates transition through the membrane. The population of the closed conformation in water was found to correlate positively with PAMPA log Pe.

rbFOX1/MBNL1 competition for CCUG RNA repeats binding contributes to myotonic dystrophy type 1/type 2 differences.

Myotonic dystrophy type 1 and type 2 (DM1, DM2) are caused by expansions of CTG and CCTG repeats, respectively. RNAs containing expanded CUG or CCUG repeats interfere with the metabolism of other RNAs through titration of the Muscleblind-like (MBNL) RNA binding proteins. DM2 follows a more favorable clinical course than DM1, suggesting that specific modifiers may modulate DM severity. Here, we report that the rbFOX1 RNA binding protein binds to expanded CCUG RNA repeats, but not to expanded CUG RNA repeats. Interestingly, rbFOX1 competes with MBNL1 for binding to CCUG expanded repeats and overexpression of rbFOX1 partly releases MBNL1 from sequestration within CCUG RNA foci in DM2 muscle cells. Furthermore, expression of rbFOX1 corrects alternative splicing alterations and rescues muscle atrophy, climbing and flying defects caused by expression of expanded CCUG repeats in a Drosophila model of DM2.

Interconversion Rates between Conformational States as Rationale for the Membrane Permeability of Cyclosporines.

Cyclic peptides have regained interest as potential inhibitors of challenging targets but have often a low bioavailability. The natural product cyclosporine A (CsA) is the textbook exception. Despite its size and polar backbone, it is able to passively cross membranes. This ability is hypothesized to be due to a conformational change from the low-energy conformation in water to a "congruent" conformation that is populated both in water and inside the membrane. Here, we use a combination of NMR measurements and kinetic models based on molecular dynamics simulations to rationalize the difference in the membrane permeability of cyclosporine E (CsE) and CsA. The structure of CsE differs only in a backbone methylation, but its membrane permeability is one order of magnitude lower. The most striking difference is found in the interconversion rates between the conformational states favored in water and in chloroform, which are up to one order of magnitude slower for CsE compared to CsA.

Rational Design of Membrane-Pore-Forming Peptides.

Specific interactions of peptides with lipid membranes are essential for cellular communication and constitute a central aspect of the innate host defense against pathogens. A computational method for generating innovative membrane-pore-forming peptides inspired by natural templates is presented. Peptide representation in terms of sequence- and topology-dependent hydrophobic moments is introduced. This design concept proves to be appropriate for the de novo generation of first-in-class membrane-active peptides with the anticipated mode of action. The designed peptides outperform the natural template in terms of their antibacterial activity. They form a kinked helical structure and self-assemble in the membrane by an entropy-driven mechanism to form dynamically growing pores that are dependent on the lipid composition. The results of this study demonstrate the unique potential of natural template-based peptide design for chemical biology and medicinal chemistry.

Peptide-Membrane Interaction between Targeting and Lysis.

Certain cationic peptides interact with biological membranes. These often-complex interactions can result in peptide targeting to the membrane, or in membrane permeation, rupture, and cell lysis. We investigated the relationship between the structural features of membrane-active peptides and these effects, to better understand these processes. To this end, we employed a computational method for morphing a membranolytic antimicrobial peptide into a nonmembranolytic mitochondrial targeting peptide by "directed simulated evolution." The results obtained demonstrate that superficially subtle sequence modifications can strongly affect the peptides' membranolytic and membrane-targeting abilities. Spectroscopic and computational analyses suggest that N- and C-terminal structural flexibility plays a crucial role in determining the mode of peptide-membrane interaction.

Structural study of the Fox-1 RRM protein hydration reveals a role for key water molecules in RRM-RNA recognition.

The Fox-1 RNA recognition motif (RRM) domain is an important member of the RRM protein family. We report a 1.8 Å X-ray structure of the free Fox-1 containing six distinct monomers. We use this and the nuclear magnetic resonance (NMR) structure of the Fox-1 protein/RNA complex for molecular dynamics (MD) analyses of the structured hydration. The individual monomers of the X-ray structure show diverse hydration patterns, however, MD excellently reproduces the most occupied hydration sites. Simulations of the protein/RNA complex show hydration consistent with the isolated protein complemented by hydration sites specific to the protein/RNA interface. MD predicts intricate hydration sites with water-binding times extending up to hundreds of nanoseconds. We characterize two of them using NMR spectroscopy, RNA binding with switchSENSE and free-energy calculations of mutant proteins. Both hydration sites are experimentally confirmed and their abolishment reduces the binding free-energy. A quantitative agreement between theory and experiment is achieved for the S155A substitution but not for the S122A mutant. The S155 hydration site is evolutionarily conserved within the RRM domains. In conclusion, MD is an effective tool for predicting and interpreting the hydration patterns of protein/RNA complexes. Hydration is not easily detectable in NMR experiments but can affect stability of protein/RNA complexes.

Kinetic Models of Cyclosporin A in Polar and Apolar Environments Reveal Multiple Congruent Conformational States.

The membrane permeability of cyclic peptides and peptidomimetics, which are generally larger and more complex than typical drug molecules, is likely strongly influenced by the conformational behavior of these compounds in polar and apolar environments. The size and complexity of peptides often limit their bioavailability, but there are known examples of peptide natural products such as cyclosporin A (CsA) that can cross cell membranes by passive diffusion. CsA is an undecapeptide with seven methylated backbone amides. Its crystal structure shows a "closed" twisted ß-pleated sheet conformation with four intramolecular hydrogen bonds that is also observed in NMR measurements of CsA in chloroform. When binding to its target cyclophilin, on the other hand, CsA adopts an "open" conformation without intramolecular hydrogen bonds. In this study, we attempted to sample the complete conformational space of CsA in chloroform and in water by molecular dynamics simulations in order to better understand its conformational behavior in these two environments and to rationalize the good membrane permeability of CsA observed experimentally. From 10 µs molecular dynamics simulations in each solvent, Markov state models were constructed to characterize the metastable conformational states. The model in chloroform is compared to nuclear Overhauser effect NMR spectroscopy data reported in this study and taken from the literature. The conformational landscapes in the two solvents show significant overlap but also clearly distinct features.

Synergy between NMR measurements and MD simulations of protein/RNA complexes: application to the RRMs, the most common RNA recognition motifs.

RNA recognition motif (RRM) proteins represent an abundant class of proteins playing key roles in RNA biology. We present a joint atomistic molecular dynamics (MD) and experimental study of two RRM-containing proteins bound with their single-stranded target RNAs, namely the Fox-1 and SRSF1 complexes. The simulations are used in conjunction with NMR spectroscopy to interpret and expand the available structural data. We accumulate more than 50 µs of simulations and show that the MD method is robust enough to reliably describe the structural dynamics of the RRM-RNA complexes. The simulations predict unanticipated specific participation of Arg142 at the protein-RNA interface of the SRFS1 complex, which is subsequently confirmed by NMR and ITC measurements. Several segments of the protein-RNA interface may involve competition between dynamical local substates rather than firmly formed interactions, which is indirectly consistent with the primary NMR data. We demonstrate that the simulations can be used to interpret the NMR atomistic models and can provide qualified predictions. Finally, we propose a protocol for 'MD-adapted structure ensemble' as a way to integrate the simulation predictions and expand upon the deposited NMR structures. Unbiased µs-scale atomistic MD could become a technique routinely complementing the NMR measurements of protein-RNA complexes.

Design and Development of a Cyclic Decapeptide Scaffold with Suitable Properties for Bioavailability and Oral Exposure.

Permeability and oral bioavailability of macrocyclic peptides still represent difficult challenges in drug discovery. Despite the recognized potential of macrocyclic peptides as therapeutics, their use is still restricted to extracellular targets and intravenous administration. Indeed, macrocyclic peptides generally suffer from limited proteolytic stability, high clearance, and poor membrane permeability, and this leads to the absence of systemic exposure after oral administration. To overcome these limitations, we started to investigate the development of a general cyclic decapeptide scaffold that possesses ideal features for cell permeability and oral exposure. On the basis of a rigid hairpin structure, the scaffold design aimed to decrease the overall polarity of the compound, thereby limiting the energetic cost of NH desolvation and the entropy penalty during cell penetration. The results of this study also demonstrate the importance of rigidity for the ß-turn design regarding clearance. To stabilize the scaffold in the desired ß-hairpin conformation, the introduction of d-proline at the i+1 turn position proved to be beneficial for both permeability and clearance. As a result, cyclopeptide decamers with unprecedented high values for oral bioavailability and exposure are reported herein. NMR spectroscopy conformation and dynamic analysis confirmed, for selected examples, the rigidity of the scaffold and the presence of transannular hydrogen bonds in polar and apolar environments. Furthermore, we showed, for one compound, that its transition from a polar environment to an apolar one was accompanied by an increased molecular motion, revealing an entropy contribution to membrane permeation.

Pharmacokinetic Studies around the Mono- and Difunctionalization of a Bioavailable Cyclic Decapeptide Scaffold.

We previously reported the design of several cyclic decapeptides based on a generic scaffold that achieved favorable oral bioavailability and exposure. With the goal to further investigate the potential of this approach, we describe herein the effect of mono- and difunctionalization of this scaffold. A series of cyclic decapeptides were therefore subjected to a range of in vitro assays and pharmacokinetic (PK) studies to investigate whether the introduction of polar or charged groups could be tolerated by the "engineered" scaffold while maintaining good PK profiles. Whereas the introduction of charged amino acids proved-besides maintaining low clearance-to conceal the inherent PK properties of the scaffold, the introduction of polar amino acids (i.e., threonine and pyridyl alanine) led to several cyclic decapeptides exhibiting excellent PK profiles together with a solubility that was significantly improved relative to that of previously reported cyclic decapeptides.

The Signature of the Five-Stranded vRRM Fold Defined by Functional, Structural and Computational Analysis of the hnRNP L Protein.

The RNA recognition motif (RRM) is the far most abundant RNA binding domain. In addition to the typical ß1α1ß2ß3α2ß4 fold, various sub-structural elements have been described and reportedly contribute to the high functional versatility of RRMs. The heterogeneous nuclear ribonucleoprotein L (hnRNP L) is a highly abundant protein of 64 kDa comprising four RRM domains. Involved in many aspects of RNA metabolism, hnRNP L specifically binds to RNAs containing CA repeats or CA-rich clusters. However, a comprehensive structural description of hnRNP L including its sub-structural elements is missing. Here, we present the structural characterization of the RRM domains of hnRNP L and demonstrate their function in repressing exon 4 of SLC2A2. By comparison of the sub-structural elements between the two highly similar paralog families of hnRNP L and PTB, we defined signatures underlying interacting C-terminal coils (ICCs), the RRM34 domain interaction and RRMs with a C-terminal fifth ß-strand, a variation we denoted vRRMs. Furthermore, computational analysis revealed new putative ICC-containing RRM families and allowed us to propose an evolutionary scenario explaining the origins of the ICC and fifth ß-strand sub-structural extensions. Our studies provide insights of domain requirements in alternative splicing mediated by hnRNP L and molecular descriptions for the sub-structural elements. In addition, the analysis presented may help to classify other abundant RRM extensions and to predict structure-function relationships.

Solution structure of the YTH domain in complex with N6-methyladenosine RNA: a reader of methylated RNA.

N(6)A methylation is the most abundant RNA modification occurring within messenger RNA. Impairment of methylase or demethylase functions are associated with severe phenotypes and diseases in several organisms. Beside writer and eraser enzymes of this dynamic RNA epigenetic modification, reader proteins that recognize this modification are involved in numerous cellular processes. Although the precise characterization of these reader proteins remains unknown, preliminary data showed that most potential reader proteins contained a conserved YT521-B homology (YTH) domain. Here we define the YTH domain of rat YT521-B as a N(6)-methylated adenosine reader domain and report its solution structure in complex with a N(6)-methylated RNA. The structure reveals a binding preference for NGANNN RNA hexamer and a deep hydrophobic cleft for m(6)A recognition. These findings establish a molecular function for YTH domains as m(6)A reader domains and should guide further studies into the biological functions of YTH-containing proteins in m(6)A recognition.

Solution and crystal structures of a C-terminal fragment of the neuronal isoform of the polypyrimidine tract binding protein (nPTB).

The eukaryotic polypyrimidine tract binding protein (PTB) serves primarily as a regulator of alternative splicing of messenger RNA, but is also co-opted to other roles such as RNA localisation and translation initiation from internal ribosome entry sites. The neuronal paralogue of PTB (nPTB) is 75% identical in amino acid sequence with PTB. Although the two proteins have broadly similar RNA binding specificities and effects on RNA splicing, differential expression of PTB and nPTB can lead to the generation of alternatively spliced mRNAs. RNA binding by PTB and nPTB is mediated by four RNA recognition motifs (RRMs). We present here the crystal and solution structures of the C-terminal domain of nPTB (nPTB34) which contains RRMs 3 and 4. As expected the structures are similar to each other and to the solution structure of the equivalent fragment from PTB (PTB34). The result confirms that, as found for PTB, RRMs 3 and 4 of nPTB interact with one another to form a stable unit that presents the RNA-binding surfaces of the component RRMs on opposite sides that face away from each other. The major differences between PTB34 and nPTB34 arise from amino acid side chain substitutions on the exposed ß-sheet surfaces and adjoining loops of each RRM, which are likely to modulate interactions with RNA.

Automated and assisted RNA resonance assignment using NMR chemical shift statistics.

The three-dimensional structure determination of RNAs by NMR spectroscopy relies on chemical shift assignment, which still constitutes a bottleneck. In order to develop more efficient assignment strategies, we analysed relationships between sequence and (1)H and (13)C chemical shifts. Statistics of resonances from regularly Watson-Crick base-paired RNA revealed highly characteristic chemical shift clusters. We developed two approaches using these statistics for chemical shift assignment of double-stranded RNA (dsRNA): a manual approach that yields starting points for resonance assignment and simplifies decision trees and an automated approach based on the recently introduced automated resonance assignment algorithm FLYA. Both strategies require only unlabeled RNAs and three 2D spectra for assigning the H2/C2, H5/C5, H6/C6, H8/C8 and H1'/C1' chemical shifts. The manual approach proved to be efficient and robust when applied to the experimental data of RNAs with a size between 20 nt and 42 nt. The more advanced automated assignment approach was successfully applied to four stem-loop RNAs and a 42 nt siRNA, assigning 92-100% of the resonances from dsRNA regions correctly. This is the first automated approach for chemical shift assignment of non-exchangeable protons of RNA and their corresponding (13)C resonances, which provides an important step toward automated structure determination of RNAs.

RNA recognition motifs: boring? Not quite.

The RNA recognition motif (RRM) is one of the most abundant protein domains in eukaryotes. While the structure of this domain is well characterized by the packing of two alpha-helices on a four-stranded beta-sheet, the mode of protein and RNA recognition by RRMs is not clear owing to the high variability of these interactions. Here we report recent structural data on RRM-RNA and RRM-protein interactions showing the ability of this domain to modulate its binding affinity and specificity using each of its constitutive elements (beta-strands, loops, alpha-helices). The extreme structural versatility of the RRM interactions explains why RRM-containing proteins have so diverse biological functions.