Mapping noncovalent ligand binding to stemloop domains of the HIV-1 packaging signal by tandem mass spectrometry

The binding modes and structural determinants of the noncovalent complexes formed by aminoglycoside antibiotics with conserved domains of the HIV-1 packaging signal (Psi-RNA) were investigated using electrospray ionization (ESI) Fourier transform mass spectrometry (FTMS). The location of the aminoglycoside binding sites on the different stemloop structures was revealed by characteristic coverage gaps in the ion series obtained by sustained off-resonance irradiation collision induced dissociation (SORI-CID) of the antibiotic-RNA assemblies. The site positions were confirmed using mutants that eliminated salient structural features of the Psi-RNA domains. The effects of the mutations on the binding properties of the different substrates served to validate the position of the aminoglycoside site on the wild-type structures. Additional information was provided by docking experiments performed on the different aminoglycoside-stemloop complexes. The results have shown that, in the absence of features disrupting the regular A-helix of the double-stranded stem, aminoglycosides tend to bind in an area situated between the upper stem and the loop regions, as demonstrated for stemloop SL3. The presence of a tandem wobbles motif in SL4 modifies the regular geometry of the upper stem, which does not affect the general site location, but greatly increases its solution binding affinity compared with SL3. The platform motif in SL2 locates the binding site in the stem midsection and confers upon this stemloop an intermediate affinity toward aminoglycosides. In SL3 and SL4, the extensive overlap of the antibiotic site with the region used to bind the nucleocapsid (NC) protein provides the basis for a competition mechanism that could explain the aminoglycoside inhibition of the NC.SL3 and NC.SL4 assemblies. In contrast, the minimal overlap between the aminoglycoside and the NC sites in SL2 accounts for the absence of inhibition of the NC.SL2 complex.     

Screening of Minimalist Noncanonical Sites in Duplex DNA and RNA Reveals Context and Motif-Selective Binding by Fluorogenic Base Probes

We hypothesize that programmable hybridization to noncanonical nucleic acid motifs may be achieved by macromolecular display of binders to individual noncanonical pairs (NCPs). As each recognition element may individually have weak binding to an NCP, we developed a semi-rational approach to detect low affinity interactions between selected nitrogenous bases and noncanonical sites in duplex DNA and RNA. A set of fluorogenic probes was synthesized by coupling abiotic (triazines, pyrimidines) and native RNA bases to thiazole orange (TO) dye. This probe library was screened against duplex nucleic acid substrates bearing single abasic, single NCP, and tandem NCP sites. Probe engagement with NCP sites was reported by 100–1000× fluorescence enhancement over background. Binding is strongly context-dependent, reflective of both molecular recognition and stability: less stable motifs are more likely to bind a synthetic probe. Further, DNA and RNA substrates exhibit entirely different abasic and single NCP binding profiles. While probe binding in the abasic and single NCP screens was monotonous, much richer binding profiles were observed with the screen of tandem NCP sites in RNA, in part due to increased steric accessibility. In addition to known binding interactions between the triazine melamine (M) and T/U sites, the NCP screens identified new targeting elements for pyrimidine-rich motifs in single NCPs and 2×2 internal bulges. We anticipate that semi-rational approaches of this type will lead to programmable noncanonical hybridization strategies at the macromolecular level.     

Mass spectrometric investigation of protein alkylation by the RNA footprinting probe kethoxal

The reactivity of the RNA footprinting reagent kethoxal (KT) toward proteins was investigated by electrospray ionization-Fourier transform mass spectrometry. Using standard peptides, KT was shown to selectively modify the guanidino group of arginine side chains at neutral pH, while primary amino groups of lysine and N-terminus were found to be unreactive under these conditions. Gas-phase fragmentation of KT adducts provided evidence for a cyclic 1,2-diol structure. Esterification of the 1,2-diol product was obtained in borate buffer, and its structure was also investigated by tandem mass spectrometry. When model proteins were probed with this RNA footprinting reagent, the adducts proved to be sufficiently stable to allow for the application of different peptide-mapping procedures to identify the location of modified arginines. Probing of proteins under native folding conditions provided modification patterns that very closely matched the structural context of arginines in the global protein structure. A strong correlation was demonstrated between the susceptibility to modification and residue accessibility calculated from the known 3D structure. When the complexes formed by HIV-1 nucleocapsid (NC) protein and RNA stemloops SL2 and SL3 were investigated, KT footprinting provided accurate information regarding the involvement of individual arginines in binding RNA and showed different reactivity according to their mode of interaction.     

Femtomole SHAPE reveals regulatory structures in the authentic XMRV RNA genome

Higher-order structure influences critical functions in nearly all noncoding and coding RNAs. Most single-nucleotide resolution RNA structure determination technologies cannot be used to analyze RNA from scarce biological samples, like viral genomes. To make quantitative RNA structure analysis applicable to a much wider array of RNA structure-function problems, we developed and applied high-sensitivity selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) to structural analysis of authentic genomic RNA of the xenotropic murine leukemia virus-related virus (XMRV). For analysis of fluorescently labeled cDNAs generated in high-sensitivity SHAPE experiments, we developed a two-color capillary electrophoresis approach with zeptomole molecular detection limits and subfemtomole sensitivity for complete SHAPE experiments involving hundreds of individual RNA structure measurements. High-sensitivity SHAPE data correlated closely (R = 0.89) with data obtained by conventional capillary electrophoresis. Using high-sensitivity SHAPE, we determined the dimeric structure of the XMRV packaging domain, examined dynamic interactions between the packaging domain RNA and viral nucleocapsid protein inside virion particles, and identified the packaging signal for this virus. Despite extensive sequence differences between XMRV and the intensively studied Moloney murine leukemia virus, architectures of the regulatory domains are similar and reveal common principles of gammaretrovirus RNA genome packaging.     

Discovering protein-binding RNA motifs with a generative model of RNA sequences

Recent advances in high-throughput experimental technologies have generated a huge amount of data on interactions between proteins and nucleic acids. Motivated by the big experimental data, several computational methods have been developed either to predict binding sites in a sequence or to determine if an interaction exists between protein and nucleic acid sequences. However, most of the methods cannot be used to discover new nucleic acid sequences that bind to a target protein because they are classifiers rather than generators. In this paper we propose a generative model for constructing protein-binding RNA sequences and motifs using a long short-term memory (LSTM) neural network. Testing the model for several target proteins showed that RNA sequences generated by the model have high binding affinity and specificity for their target proteins and that the protein-binding motifs derived from the generated RNA sequences are comparable to the motifs from experimentally validated protein-binding RNA sequences. The results are promising and we believe this approach will help design more efficient in vitro or in vivo experiments by suggesting potential RNA aptamers for a target protein.     

Segmental,Domain‐Selective Perdeuteration and Small‐Angle Neutron Scattering for Structural Analysis of Multi‐Domain Proteins

Multi‐domain proteins play critical roles in fine‐tuning essential processes in cellular signaling and gene regulation. Typically, multiple globular domains that are connected by flexible linkers undergo dynamic rearrangements upon binding to protein, DNA or RNA ligands. RNA binding proteins (RBPs) represent an important class of multi‐domain proteins, which regulate gene expression by recognizing linear or structured RNA sequence motifs. Here, we employ segmental perdeuteration of the three RNA recognition motif (RRM) domains in the RBP TIA‐1 using Sortase A mediated protein ligation. We show that domain‐selective perdeuteration combined with contrast‐matched small‐angle neutron scattering (SANS), SAXS and computational modeling provides valuable information to precisely define relative domain arrangements. The approach is generally applicable to study conformational arrangements of individual domains in multi‐domain proteins and changes induced by ligand binding.     

Segmental,Domain‐Selective Perdeuteration and Small‐Angle Neutron Scattering for Structural Analysis of Multi‐Domain Proteins

Multi‐domain proteins play critical roles in fine‐tuning essential processes in cellular signaling and gene regulation. Typically, multiple globular domains that are connected by flexible linkers undergo dynamic rearrangements upon binding to protein, DNA or RNA ligands. RNA binding proteins (RBPs) represent an important class of multi‐domain proteins, which regulate gene expression by recognizing linear or structured RNA sequence motifs. Here, we employ segmental perdeuteration of the three RNA recognition motif (RRM) domains in the RBP TIA‐1 using Sortase A mediated protein ligation. We show that domain‐selective perdeuteration combined with contrast‐matched small‐angle neutron scattering (SANS), SAXS and computational modeling provides valuable information to precisely define relative domain arrangements. The approach is generally applicable to study conformational arrangements of individual domains in multi‐domain proteins and changes induced by ligand binding.     

Noncovalent interactions: defining cooperativity. Ligand binding aided by reduced dynamic behavior of receptors. Binding of bacterial cell wall analogues to ristocetin A

Changes in the relative populations of the monomer and asymmetric dimer forms of ristocetin A, upon binding of two molecules of ligand, suggest that ligand binding is negatively cooperative with respect to dimerization. However, strong hydrogen bonds formed in the binding sites of the ligands are reinforced in the dimer relative to the monomer, and the barrier to dissociation of the dimer is increased upon binding of the ligands. It is concluded that the interactions which are common in the binding of both ligands are made with positive cooperativity with respect to those involved in dimerization. The conclusions are relevant to the binding of ligands to proteins, where ligand binding energy can be derived from stabilization of the protein in its ligand-bound form.     

Fabrication of Pascal‐triangle Lattice of Proteins by Inducing Ligand Strategy

A protein Pascal triangle has been constructed as new type of supramolecular architecture by using the inducing ligand strategy that we previously developed for protein assemblies. Although mathematical studies on this famous geometry have a long history, no work on such Pascal triangles fabricated from native proteins has been reported so far due to their structural complexity. In this work, by carefully tuning the specific interactions between the native protein building block WGA and the inducing ligand R‐SL , a 2D Pascal‐triangle lattice with three types of triangular voids has been assembled. Moreover, a 3D crystal structure was obtained based on the 2D Pascal triangles. The distinctive carbohydrate binding sites of WGA and the intralayer as well as interlayer dimerization of RhB was the key to facilitate nanofabrication in solution. This strategy may be applied to prepare and explore various sophisticated assemblies based on native proteins.     

Immunomodulatory peptides from IgSF proteins: a review

The immunoglobulin (Ig) domain is a highly conserved domain predominantly observed in cell surface proteins due to its ability to resist proteolysis. By mutation and selection the Ig domain has evolved to serve diverse biological functions including growth and development, signaling, adhesion and protein-carbohydrate interactions. Collectively, proteins with Ig-like domain constitute the immunoglobulin superfamily (IgSF). The IgSF proteins make up over 2% of human genes constituting the largest gene family in the human genome. Analogous to the complementarity determining regions (CDR)s that form the antigen combining sites of the antibody, the high specificity of the IgSF receptor-ligand interaction is attributed to the sequence and structure of the CDR-like regions unique to each IgSF protein. Hence, CDR-like regions provide ideal templates for the design of mimetics that can potentially perturb specific IgSF receptor/ligand interactions. The determinants of binding are localized near the CDR-like regions, conformation is determined locally and is unique for each loop. In structure based drug design one of the approaches to identify lead agents is to map the receptor/ligand binding epitope onto a small peptide. Data from theoretical, structural and functional studies have been adopted in the design of novel peptide antagonists of the IgSF protein-protein interactions. Many peptide antagonists have shown significant therapeutic potential in multiple animal models. The design of the IgSF peptide analogs, rationale as therapeutic targets, functional efficacy and the clinical benefits are reviewed here.     

Why not consider a spherical protein? Implications of backbone hydrogen bonding for protein structure and function

The intrinsic ability of protein structures to exhibit the geometric features required for molecular function in the absence of evolution is examined in the context of three systems: the reference set of real, single domain protein structures, a library of computationally generated, compact homopolypeptides, artificial structures with protein-like secondary structural elements, and quasi-spherical random proteins packed at the same density as proteins but lacking backbone secondary structure and hydrogen bonding. Without any evolutionary selection, the library of artificial structures has similar backbone hydrogen bonding, global shape, surface to volume ratio and statistically significant structural matches to real protein global structures. Moreover, these artificial structures have native like ligand binding cavities, and a tiny subset has interfacial geometries consistent with native-like protein-protein interactions and DNA binding. In contrast, the quasi-spherical random proteins, being devoid of secondary structure, have a lower surface to volume ratio and lack ligand binding pockets and intermolecular interaction interfaces. Surprisingly, these quasi-spherical random proteins exhibit protein like distributions of virtual bond angles and almost all have a statistically significant structural match to real protein structures. This implies that it is local chain stiffness, even without backbone hydrogen bonding, and compactness that give rise to the likely completeness of the library solved single domain protein structures. These studies also suggest that the packing of secondary structural elements generates the requisite geometry for intermolecular binding. Thus, backbone hydrogen bonding plays an important role not only in protein structure but also in protein function. Such ability to bind biological molecules is an inherent feature of protein structure; if combined with appropriate protein sequences, it could provide the non-zero background probability for low-level function that evolution requires for selection to occur.     

Validation of an empirical RNA-ligand scoring function for fast flexible docking using RiboDock®

We report the design and validation of a fast empirical function for scoring RNA-ligand interactions, and describe its implementation within RiboDock, a virtual screening system for automated flexible docking. Building on well-known protein-ligand scoring function foundations, features were added to describe the interactions of common RNA-binding functional groups that were not handled adequately by conventional terms, to disfavour non-complementary polar contacts, and to control non-specific charged interactions. The results of validation experiments against known structures of RNA-ligand complexes compare favourably with previously reported methods. Binding modes were well predicted in most cases and good discrimination was achieved between native and non-native ligands for each binding site, and between native and non-native binding sites for each ligand. Further evidence of the ability of the method to identify true RNA binders is provided by compound selection ('enrichment factor') experiments based around a series of HIV-1 TAR RNA-binding ligands. Significant enrichment in true binders was achieved amongst high scoring docking hits, even when selection was from a library of structurally related, positively charged molecules. Coupled with a semi-automated cavity detection algorithm for identification of putative ligand binding sites, also described here, the method is suitable for the screening of very large databases of molecules against RNA and RNA-protein interfaces, such as those presented by the bacterial ribosome.     

Dynamic properties of water around a protein-DNA complex from molecular dynamics simulations

Formation of protein-DNA complex is an important step in regulation of genes in living organisms. One important issue in this problem is the role played by water in mediating the protein-DNA interactions. In this work, we have carried out atomistic molecular dynamics simulations to explore the heterogeneous dynamics of water molecules present in different regions around a complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. It is demonstrated that such heterogeneous water motions around the complex are correlated with the relaxation time scales of hydrogen bonds formed by those water molecules with the protein and DNA. The calculations reveal the existence of a fraction of extraordinarily restricted water molecules forming a highly rigid thin layer in between the binding motifs of the protein and DNA. It is further proved that higher rigidity of water layers around the complex originates from more frequent reformations of broken water-water hydrogen bonds. Importantly, it is found that the formation of the complex affects the transverse and longitudinal degrees of freedom of surrounding water molecules in a nonuniform manner.     

Control of the stability of a protein-RNA complex by the position of fluorine in a base analogue

The effects of modifying the electronic characteristics of nonpolar base analogues substituted at positions involved in stacking interactions between SL2 RNA and the U1A protein are described. A surprisingly large difference in the stability between complexes formed with base analogues that differ only in the position of substitution of a single fluorine atom is observed. The results of high-level ab initio calculations of the interactions between the nonpolar base analogue and the amino acid side chain correlate with the experimentally observed trends in complex stability, which suggests that changes in stacking interactions that result from varying the position and degree of fluorine substitution contribute to the effects of fluorine substitution on the stability of the U1A-SL2 RNA complex.     

Site mapping and small molecule blind docking reveal a possible target site on the SARS-CoV-2 main protease dimer interface

The SARS-CoV-2 virus is causing COVID-19 resulting in an ongoing pandemic with serious health, social, and economic implications. Much research is focused in repurposing or identifying new small molecules which may interact with viral or host-cell molecular targets. An important SARS-CoV-2 target is the main protease (M^pro), and the peptidomimetic α-ketoamides represent prototypical experimental inhibitors. The protease is characterised by the dimerization of two monomers each which contains the catalytic dyad defined by Cys¹⁴⁵ and His⁴¹ residues (active site). Dimerization yields the functional homodimer. Here, our aim was to investigate small molecules, including lopinavir and ritonavir, α-ketoamide 13b, and ebselen, for their ability to interact with the M^pro. The sirtuin 1 agonist SRT1720 was also used in our analyses. Blind docking to each monomer individually indicated preferential binding of the ligands in the active site. Site-mapping of the dimeric protease indicated a highly reactive pocket in the dimerization region at the domain III apex. Blind docking consistently indicated a strong preference of ligand binding in domain III, away from the active site. Molecular dynamics simulations indicated that ligands docked both to the active site and in the dimerization region at the apex, formed relatively stable interactions. Overall, our findings do not obviate the superior potency with respect to inhibition of protease activity of covalently-linked inhibitors such as α-ketoamide 13b in the M^pro active site. Nevertheless, along with those from others, our findings highlight the importance of further characterisation of the M^pro active site and any potential allosteric sites.     

Cooperativity and anti-cooperativity between ligand binding and the dimerization of ristocetin A: asymmetry of a homodimer complex and implications for signal transduction

Background: Recent work has indicated that dimerization is important in the mode of action of the vancomycin group of glycopeptide antibiotics. NMR studies have shown that one member of this group, ristocetin A₁ forms an asymmetric dimer with two physically different binding sites for cell wall peptides. Ligand binding by ristocetin A and dimerization are slightly anti-cooperative. In contrast, for the other glycopeptide antibiotics of the vancomycin group that have been examined so far, binding of cell wall peptides and dimerization are cooperative.Results: Here we show that the two halves of the asymmetric homodimer formed by ristocetin A have different affinities for ligand binding. One of these sites is preferentially filled before the other, and binding to this site is cooperative with dimerization. Ligand binding to the other, less favored half of the dimer, is anti-cooperative with dimerization.Conclusions: In dinner complexes, anti-cooperativity of dimerization upon ligand binding can be a result of asymmetry, in which two binding sites have different affinities for ligands. Such a system, in which one binding site is filled preferentially, may be a mechanism by which the cooperativity between ligand binding and dimerization is fine tuned and may thus have relevance to the control of signal transduction in biological systems.     

Novel features for identifying A-minors in three-dimensional RNA molecules

RNA tertiary interactions or tertiary motifs are conserved structural patterns formed by pairwise interactions between nucleotides. They include base-pairing, base-stacking, and base-phosphate interactions. A-minor motifs are the most common tertiary interactions in the large ribosomal subunit. The A-minor motif is a nucleotide triple in which minor groove edges of an adenine base are inserted into the minor groove of neighboring helices, leading to interaction with a stabilizing base pair. We propose here novel features for identifying and predicting A-minor motifs in a given three-dimensional RNA molecule. By utilizing the features together with machine learning algorithms including random forests and support vector machines, we show experimentally that our approach is capable of predicting A-minor motifs in the given RNA molecule effectively, demonstrating the usefulness of the proposed approach. The techniques developed from this work will be useful for molecular biologists and biochemists to analyze RNA tertiary motifs, specifically A-minor interactions.     

Phosphate recognition in structural biology

Drug-discovery research in the past decade has seen an increased selection of targets with phosphate recognition sites, such as protein kinases and phosphatases, in the past decade. This review attempts, with the help of database-mining tools, to give an overview of the most important principles in molecular recognition of phosphate groups by enzymes. A total of 3003 X-ray crystal structures from the RCSB Protein Data Bank with bound organophosphates has been analyzed individually, in particular for H-bonding interactions between proteins and ligands. The various known binding motifs for phosphate binding are reviewed, and similarities to phosphate complexation by synthetic receptors are highlighted. An analysis of the propensities of amino acids in various classes of phosphate-binding enzymes showed characteristic distributions of amino acids used for phosphate binding. This review demonstrates that structure-based lead development and optimization should carefully address the phosphate-binding-site environment and also proposes new alternatives for filling such sites.