Structure-Guided Peptidomimetic Design Leads to Nanomolar β-Hairpin Inhibitors of the Tat−TAR Interaction of Bovine Immunodeficiency Virus,

The Tat protein of immunodeficiency viruses is the main activator of viral gene expression. By binding specifically to its cognate site, the transactivator response element (TAR), Tat mediates a strong induction of the production of all viral transcripts. In seeking a new chemical solution to inhibiting viral protein-RNA interactions, we recently identified inhibitors of the viral Tat protein from the bovine immunodeficiency virus (BIV) using conformationally constrained beta-hairpin peptidomimetics. We identified a micromolar ligand, called BIV2, and the structure of its complex with BIV TAR was determined by NMR. In this work, we demonstrate that this chemistry can rapidly yield highly potent and selective ligands. On the basis of the structure, we synthesized and assayed libraries of mutant peptidomimetics. Remarkably, we were able in just a few rounds of design and synthesis to discover nanomolar inhibitors of the Tat-TAR interaction in BIV that selectively bind the BIV TAR RNA compared to RNA structures as closely related as the HIV-1 TAR or RRE elements. The molecular recognition principles developed in this study have been exploited in discovering related peptidomimetic inhibitors of the Tat-TAR interaction in HIV-1.

High-resolution structural and thermodynamic analysis of extreme stabilization of human procarboxypeptidase by computational protein design

Recent efforts to design de novo or redesign the sequence and structure of proteins using computational techniques have met with significant success. Most, if not all, of these computational methodologies attempt to model atomic-level interactions, and hence high-resolution structural characterization of the designed proteins is critical for evaluating the atomic-level accuracy of the underlying design force-fields. We previously used our computational protein design protocol RosettaDesign to completely redesign the sequence of the activation domain of human procarboxypeptidase A2. With 68% of the wild-type sequence changed, the designed protein, AYEdesign, is over 10 kcal/mol more stable than the wild-type protein. Here, we describe the high-resolution crystal structure and solution NMR structure of AYEdesign, which show that the experimentally determined backbone and side-chains conformations are effectively superimposable with the computational model at atomic resolution. To isolate the origins of the remarkable stabilization, we have designed and characterized a new series of procarboxypeptidase mutants that gain significant thermodynamic stability with a minimal number of mutations; one mutant gains more than 5 kcal/mol of stability over the wild-type protein with only four amino acid changes. We explore the relationship between force-field smoothing and conformational sampling by comparing the experimentally determined free energies of the overall design and these focused subsets of mutations to those predicted using modified force-fields, and both fixed and flexible backbone sampling protocols.

The C-terminal domains of vertebrate CstF-64 and its yeast orthologue Rna15 form a new structure critical for mRNA 3'-end processing

Yeast Rna15 and its vertebrate orthologue CstF-64 play critical roles in mRNA 3 '-end processing and in transcription termination downstream of poly(A) sites. These proteins contain N-terminal domains that recognize the poly(A) site, but little is known about their highly conserved C-terminal regions. Here we show by NMR that the C-terminal domains of CstF-64 and Rna15 fold into a three-helix bundle with an uncommon topological arrangement. The structure defines a cluster of evolutionary conserved yet exposed residues we show to be essential for the interaction between Pcf11 and Rna15. Furthermore, we demonstrate that this interaction is critical for the function of Rna15 in 3 '-end processing but dispensable for transcription termination. The C-terminal domain of the Rna15 homologue Pti1 contains critical sequence alterations within this region that are predicted to prevent Pcf11 interaction, providing an explanation for the distinct functions of these two closely related proteins in the 3 '-end formation of RNA polymerase II transcripts. These results define the role of the C-terminal half of Rna15 and provide insight into the network of protein/protein interactions responsible for assembly of the 3 '-end processing apparatus.

An all-atom, distance-dependent scoring function for the prediction of protein-DNA interactions from structure

We have developed an all-atom statistical potential function for the prediction of protein-DNA interactions from their structures, and show that this method outperforms similar, lower-resolution statistical potentials in a series of decoy discrimination experiments. The all-atom formalism appears to capture details of atomic interactions that are missed by the lower-resolution methods, with the majority of the discriminatory power arising from its description of short-range atomic contacts. We show that, on average, the method is able to identify 90% of near-native docking decoys within the best-scoring 10% of structures in a given decoy set, and it compares favorably with an optimized physical potential function in a test of structure-based identification of DNA binding-sequences. These results demonstrate that all-atom statistical functions specific to protein-DNA interactions can achieve great discriminatory power despite the limited size of the structural database. They also suggest that the statistical scores may soon be able to achieve accuracy on par with more complex, physical potential functions.

Mis-translation of a computationally designed protein yields an exceptionally stable homodimer: implications for protein engineering and evolution

We recently used computational protein design to create an extremely stable, globular protein, Top7, with a sequence and fold not observed previously in nature. Since Top7 was created in the absence of genetic selection, it provides a rare opportunity to investigate aspects of the cellular protein production and surveillance machinery that are subject to natural selection. Here we show that a portion of the Top7 protein corresponding to the final 49 C-terminal residues is efficiently mis-translated and accumulates at high levels in Escherichia coli. We used circular dichroism, size-exclusion chromatography, small-angle X-ray scattering, analytical ultra-centrifugation, and NMR spectroscopy to show that the resulting C-terminal fragment (CFr) protein adopts a compact, extremely stable, homo-dimeric structure. Based on the solution structure, we engineered an even more stable variant of CFr by disulfide-induced covalent circularisation that should be an excellent platform for design of novel functions. The accumulation of high levels of CFr exposes the high error rate of the protein translation machinery. The rarity of correspondingly stable fragments in natural proteins coupled with the observation that high quality ribosome binding sites are found to occur within E. coli protein-coding regions significantly less often than expected by random chance implies a stringent evolutionary pressure against protein sub-fragments that can independently fold into stable structures. The symmetric self-association between two identical mis-translated CFr sub-domains to generate an extremely stable structure parallels a mechanism for natural protein-fold evolution by modular recombination of protein sub-structures.

Decoding RNA motional codes

When proteins and small molecules bind to RNA, they often alter its conformation. These structural changes are an essential aspect of the ability of RNA to sense signaling molecules and modulate gene expression. Thus far, few studies have been dedicated to understanding how RNA moves at a residue level and how these motions change upon complex formation. A recent report highlights how intrinsic motions in RNA correlate with its ability to bind to cognate ligands.

High-resolution structural validation of the computational redesign of human U1A protein

Achieving atomic-level resolution in the computational design of a protein structure remains a challenging problem despite recent progress. Rigorous experimental tests are needed to improve protein design algorithms, yet studies of the structure and dynamics of computationally designed proteins are very few. The NMR structure and backbone dynamics of a redesigned protein of 96 amino acids are compared here with the design target, human U1A protein. We demonstrate that the redesigned protein reproduces the target structure to within the uncertainty of the NMR coordinates, even as 65 out of 96 amino acids were simultaneously changed by purely computational methods. The dynamics of the backbone of the redesigned protein also mirror those of human U1A, suggesting that the protein design algorithm captures the shape of the potential energy landscape in addition to the local energy minimum.

Structure of stem-loop IV of Tetrahymena telomerase RNA

Conserved domains within the RNA component of telomerase provide the template for reverse transcription, recruit protein components to the holoenzyme and are required for enzymatic activity. Among the functionally essential domains in ciliate telomerase RNA is stem-loop IV, which strongly stimulates telomerase activity and processivity even when provided in trans. The NMR structure of Tetrahymena thermophila stem-loop IV shows a highly structured distal stem-loop linked to a conformationally flexible template-proximal region by a bulge that severely kinks the entire RNA. Through extensive structure-function studies, we identify residues that contribute to both these structural features and to enzymatic activity, with no apparent effect on the binding of TERT protein. We propose that the bending induced by the GA bulge and the flexibility of the template-proximal region allow positioning of the prestructured apical loop during the catalytic cycle.

A model for C74 addition by CCA-adding enzymes: C74 addition, like C75 and A76 addition, does not involve tRNA translocation

The CCA-adding enzyme adds CCA to the 3'-end of tRNA one nucleotide at a time, using CTP and ATP as substrates. We found previously that tRNA does not rotate or translocate on the enzyme during the addition of C75 and A76. We therefore predicted that the growing 3'-end of tRNA must, upon addition of each nucleotide, refold to reposition the new 3'-hydroxyl equivalently relative to the solitary nucleotidyltransferase motif. Cocrystal structures of the class I archaeal Archaeoglobus fulgidus enzyme, poised for addition of C75 and A76, confirmed this prediction. We have also demonstrated that an evolutionarily flexible beta-turn facilitates progressive refolding of the 3'-terminal C74 and C75 residues during C75 and A76 addition. Although useful cocrystals corresponding to C74 addition have not yet been obtained, we now show experimentally that tRNA does not rotate or translocate during C74 addition. We therefore propose, based on the existing A. fulgidus cocrystal structures, that the same flexible beta-turn functions as a wedge between the discriminator base (N73) and the terminal base pair of the acceptor stem, unstacking and repositioning N73 to attack the incoming CTP. Thus a single flexible beta-turn would orchestrate consecutive addition of all three nucleotides without significant movement of the tRNA on the enzyme surface.

The Cbf5-Nop10 complex is a molecular bracket that organizes box H/ACA RNPs

Box H/ACA ribonucleoprotein particles (RNPs) catalyze RNA pseudouridylation and direct processing of ribosomal RNA, and are essential architectural components of vertebrate telomerases. H/ACA RNPs comprise four proteins and a multihelical RNA. Two proteins, Cbf5 and Nop10, suffice for basal enzymatic activity in an archaeal in vitro system. We now report their cocrystal structure at 1.95-A resolution. We find that archaeal Cbf5 can assemble with yeast Nop10 and with human telomerase RNA, consistent with the high sequence identity of the RNP components between archaea and eukarya. Thus, the Cbf5-Nop10 architecture is phylogenetically conserved. The structure shows how Nop10 buttresses the active site of Cbf5, and it reveals two basic troughs that bidirectionally extend the active site cleft. Mutagenesis results implicate an adjacent basic patch in RNA binding. This tripartite RNA-binding surface may function as a molecular bracket that organizes the multihelical H/ACA and telomerase RNAs.

TAR RNA recognition by a cyclic peptidomimetic of Tat protein

The search for new antiviral drugs that repress HIV viral replication by blocking transactivation of viral RNA transcription has long been advocated as an approach to novel antiviral therapy. However, research in this area has so far failed to yield attractive lead compounds because of the insufficient development of RNA-based medicinal chemistry. One difficulty in efforts to inhibit protein-RNA interactions using small druglike molecules is the large surface areas typically found at these interfaces. To overcome this problem, we sought to identify constrained peptidomimetic inhibitors that would provide potential new drug leads. We previously reported the discovery of a cyclic peptide mimic of the RNA-binding domain of BIV Tat protein based on a designed beta-hairpin scaffold. We demonstrated that the cyclic peptide bound BIV TAR RNA with an affinity comparable to that of the RNA-binding domain of the Tat protein and inhibited protein binding to the RNA. In this study, we report the structure of the complex of the cyclic peptide bound to BIV TAR RNA determined using heteronuclear NMR methods. The structure reveals a beta-hairpin conformation in the bound peptide, which adopts an unexpected orientation in the major groove of the RNA opposite those observed for peptides derived from the Tat protein. This structure suggests many ways in which to optimize the compound and enhance its activity and pharmacological potential and represents a further step in the rational design of a new class of HIV-1 virus replication inhibitors based on peptidomimetic chemistry.

Identification of the optimal third generation antifolate against P. falciparum and P. vivax

Inhibitors of dihydrofolate reductase (DHFR) have been mainstays in the treatment of falciparum malaria. Resistance to one of these antifolates, pyrimethamine, is now common in Plasmodium falciparum populations. Antifolates have not traditionally been recommended for treatment of vivax malaria. However, recent studies have suggested that a third-generation antifolate, WR99210, is remarkably effective even against highly pyrimethamine-resistant parasites from both species. Two methods were used to identify a compound that is effective against quadruple mutant alleles from P. falciparum (N51I/C59R/S108N/I164L) and from Plasmodium vivax (57L/111L/117T/173F). The first was simple yeast system used to screen a panel of WR99210 analogs. The biguanide prodrug, JPC-2056, of the 2-chloro-4-trifluoromethoxy analog of WR99210 was effective against both the P. falciparum and P. vivax enzymes, and has been selected for further development. The second method compared the analogs in silico by docking them in the known structure of the P. falciparum DHFR-thymidylate synthase. The program reproduced well the position of the triazine ring, but the calculated energies of ligand binding were very similar for different compounds and therefore did not reproduce the observed trends in biological activity. The WR99210 family of molecules is flexible due to a long bridge between the triazine ring and the substituted benzene. During docking, multiple conformations were observed for the benzene ring part of the molecules in the DHFR active site, making computer-based predictions of binding energy less informative than for more rigid ligands. This flexibility is a key factor in their effectiveness against the highly mutant forms of DHFR.

Protein families and RNA recognition

This minireview series examines the structural principles underlying the biological function of RNA-binding proteins. The structural work of the last decade has elucidated the structures of essentially all the major RNA-binding protein families; it has also demonstrated how RNA recognition takes place. The ribosome structures have further integrated this knowledge into principles for the assembly of complex ribonucleoproteins. Structural and biochemical work has revealed unexpectedly that several RNA-binding proteins bind to other proteins in addition to RNA or instead of RNA. This tremendous increase in the structural knowledge has expanded not only our understanding of the RNA recognition principle, but has also provided new insight into the biological function of these proteins and has helped to design better experiments to understand their biological roles.

Monitoring tat peptide binding to TAR RNA by solid-state 31P-19F REDOR NMR

Complexes of the HIV transactivation response element (TAR) RNA with the viral regulatory protein tat are of special interest due in particular to the plasticity of the RNA at this binding site and to the potential for therapeutic targeting of the interaction. We performed REDOR solid-state NMR experiments on lyophilized samples of a 29 nt HIV-1 TAR construct to measure conformational changes in the tat-binding site concomitant with binding of a short peptide comprising the residues of the tat basic binding domain. Peptide binding was observed to produce a nearly 4 A decrease in the separation between phosphorothioate and 2'F labels incorporated at A27 in the upper helix and U23 in the bulge, respectively, consistent with distance changes observed in previous solution NMR studies, and with models showing significant rearrangement in position of bulge residue U23 in the bound-form RNA. In addition to providing long-range constraints on free TAR and the TAR-tat complex, these results suggest that in RNAs known to undergo large deformations upon ligand binding, 31P-19F REDOR measurements can also serve as an assay for complex formation in solid-state samples. To our knowledge, these experiments provide the first example of a solid-state NMR distance measurement in an RNA-peptide complex.

13C NMR relaxation studies of RNA base and ribose nuclei reveal a complex pattern of motions in the RNA binding site for human U1A protein

The widespread importance of induced fit and order-disorder transition in RNA recognition by proteins and small molecules makes it imperative that RNA motional properties are characterized quantitatively. Until now, however, very few studies have been dedicated to the systematic characterization of RNA motion and to their changes upon protein or small-molecule binding. The U1A protein-RNA complexes provide some of the best-studied examples of the role of RNA motional changes upon protein binding. Here, we report (13)C NMR relaxation studies of base and ribose dynamics for the RNA internal loop target of human U1A protein located within the 3'-untranslated region (3'-UTR) of the mRNA coding for U1A itself. We also report the semi-quantitative analysis of both fast (nano- to picosecond) and intermediate (micro- to millisecond) motions for this paradigmatic RNA system. We measure (13)C T(1), T(1rho) and heteronuclear nuclear Overhauser effects (NOEs) for sugar and base nuclei, as well as the power dependence of T(1rho) at 500 MHz and 750 MHz, and analyze these results using the model-free formalism. The results provide a much clearer picture of the type of motions experienced by this RNA in the absence of the protein than was provided by the analysis of the structure based solely on NOEs and scalar couplings. They define a model where the RNA internal loop region "breathes" on a micro- to millisecond timescale with respect to the double-helical regions. Superimposed on this slower motion, the residues at the very tip of the loop undergo faster (nano- to picosecond) motions. We hypothesize that these motions allow the RNA to sample multiple conformations so that the protein can select a structure within the ensemble that optimizes intermolecular contacts.

The structure of an enzyme-activating fragment of human telomerase RNA

The ribonucleoprotein enzyme telomerase ensures the stability and fidelity of linear chromosome ends by elongating the telomeric DNA that is lost during each round of DNA replication. All telomerases contain a catalytic protein component homologous to viral reverse transcriptases (TERT) and an RNA (TR) that provides the template sequence, acts as the scaffold for ribonucleoprotein assembly, and activates the enzyme for catalysis. Vertebrate telomerase RNAs contain three highly conserved structural and functional domains: the template domain, the "CR4-CR5" or "activation" domain essential for activation of the enzymatic activity, and a 3'-terminal "box H/ACA"-homology domain responsible for ribonucleprotein assembly and maturation. Here we report the NMR structure of a functionally essential RNA structural element derived from the human telomerase RNA CR4-CR5 domain. This RNA, referred to as hTR J6, forms a stable hairpin interrupted by a single nucleotide bulge and an asymmetric internal loop. Previous work on telomerase has shown that deletion of the hTR J6 asymmetric internal loop results in an RNA incapable of binding the enzymatic protein component of the RNP and therefore an inactive RNP without telomerase activity. We demonstrate here that the J6 internal loop introduces a twist in the RNA structure that may position the entire domain into the catalytic site of the enzyme.

Protein and RNA Dynamics Play Key Roles in Determining the Specific Recognition of GU-rich Polyadenylation Regulatory Elements by Human Cstf-64 Protein

The N-terminal domain of the 64 kDa subunit of the cleavage stimulation factor (CstF-64) recognizes GU-rich elements within the 3'-untranslated region of eukaryotic mRNAs. This interaction is essential for mRNA 3' end processing and transcription termination, and its strength affects the efficiency of utilization of different polyadenylation sites. The structure of the RNA-binding N-terminal domain of CstF-64 showed how the N-terminal RNA recognition motif of CstF-64 recognizes GU-rich RNAs. However, it is still perplexing how this protein can bind selectively to RNAs that are rich in G and U residues regardless of their detailed sequence composition, yet discriminate effectively against non-GU-RNAs. We investigated by NMR the dynamics of the CstF-64 RNA-binding domain, both free and bound to two GU-rich RNA sequences that represent polyadenylation regulatory elements. While the free protein displays the motional properties typical of a well-folded protein domain and is uniformly rigid, the protein-RNA interface acquires significant mobility on the micro- to millisecond time-scale once GU-rich RNAs binds to it. These motional features, we propose, are intrinsic to the functional requirement to bind all GU-rich sequences and yet to discriminate against non-GU-rich RNAs. This behavior may be a general mechanism by which some RNA-binding proteins are able to bind to classes of sequences, as opposed to a well-defined sequence or consensus.

NMR studies of protein-nucleic acid interactions

Protein-DNA and protein-RNA complexes play key functional roles in every living organism. Therefore, the elucidation of their structure and dynamics is an important goal of structural and molecular biology. Nuclear magnetic resonance (NMR) studies of protein and nucleic acid complexes have common features with studies of protein-protein complexes: the interaction surfaces between the molecules must be carefully delineated, the relative orientation of the two species needs to be accurately and precisely determined, and close intermolecular contacts defined by nuclear Overhauser effects (NOEs) must be obtained. However, differences in NMR properties (e.g., chemical shifts) and biosynthetic pathways for sample productions generate important differences. Chemical shift differences between the protein and nucleic acid resonances can aid the NMR structure determination process; however, the relatively limited dispersion of the RNA ribose resonances makes the process of assigning intermolecular NOEs more difficult. The analysis of the resulting structures requires computational tools unique to nucleic acid interactions. This chapter summarizes the most important elements of the structure determination by NMR of protein-nucleic acid complexes and their analysis. The main emphasis is on recent developments (e.g., residual dipolar couplings and new Web-based analysis tools) that have facilitated NMR studies of these complexes and expanded the type of biological problems to which NMR techniques of structural elucidation can now be applied.

A new hydrogen-bonding potential for the design of protein-RNA interactions predicts specific contacts and discriminates decoys

RNA-binding proteins play many essential roles in the regulation of gene expression in the cell. Despite the significant increase in the number of structures for RNA-protein complexes in the last few years, the molecular basis of specificity remains unclear even for the best-studied protein families. We have developed a distance and orientation-dependent hydrogen-bonding potential based on the statistical analysis of hydrogen-bonding geometries that are observed in high-resolution crystal structures of protein-DNA and protein-RNA complexes. We observe very strong geometrical preferences that reflect significant energetic constraints on the relative placement of hydrogen-bonding atom pairs at protein-nucleic acid interfaces. A scoring function based on the hydrogen-bonding potential discriminates native protein-RNA structures from incorrectly docked decoys with remarkable predictive power. By incorporating the new hydrogen-bonding potential into a physical model of protein-RNA interfaces with full atom representation, we were able to recover native amino acids at protein-RNA interfaces.