Computational analysis of amino acids and their sidechain analogs in crowded solutions of RNA nucleobases with implications for the mRNA–protein complementarity hypothesis

Many critical processes in the cell involve direct binding between RNAs and proteins, making it imperative to fully understand the physicochemical principles behind such interactions at the atomistic level. Here, we use molecular dynamics simulations and 15 μs of sampling to study the behavior of amino acids and amino acid sidechain analogs in high-concentration aqueous solutions of standard RNA nucleobases. Structural and energetic analysis of simulated systems allows us to derive interaction propensity scales for different amino acid/nucleobase combinations. The derived scales closely match and greatly extend the available experimental data, providing a comprehensive foundation for studying RNA–protein interactions in different contexts. By using these scales, we demonstrate a statistically significant connection between nucleobase composition of human mRNA coding sequences and nucleobase interaction propensities of their cognate protein sequences. For example, pyrimidine density profiles of mRNAs match uracil-propensity profiles of their cognate proteins with a median Pearson correlation coefficient of R = −0.70. Our results provide support for the recently proposed hypotheses that mRNAs and their cognate proteins may be physicochemically complementary to each other and bind, especially if unstructured, with the complementarity level being negatively influenced by mRNA adenine content. Finally, we utilize the derived scales to refine the complementarity hypothesis and closely examine its physicochemical underpinnings.     

Nucleic acid fragmentation on the millisecond timescale using a conventional X-ray rotating anode source: application to protein–DNA footprinting

Nucleic acid fragmentation (footprinting) by ·OH radicals is used often as a tool to probe nucleic acid structure and nucleic acid–protein interactions. This method has proven valuable because it provides structural information with single base pair resolution. Recent developments in the field introduced the ‘synchrotron X-ray footprinting’ method, which uses a high-flux X-ray source to produce single base pair fragmentation of nucleic acid in tens of milliseconds. We developed a complementary method that utilizes X-rays generated from a conventional rotating anode machine in which nucleic acid footprints can be generated by X-ray exposures as short as 100–300 ms. Our theoretical and experimental studies indicate that efficient cleavage of nucleic acids by X-rays depends upon sample preparation, energy of the X-ray source and the beam intensity. In addition, using this experimental set up, we demonstrated the feasibility of conducting X-ray footprinting to produce protein–DNA protection portraits at sub-second timescales.     

Nucleic acid binding properties of the nucleic acid chaperone domain of hepatitis delta antigen

The N terminal region of hepatitis delta antigen (HDAg), referred to here as NdAg, has a nucleic acid chaperone activity that modulates the ribozyme activity of hepatitis delta virus (HDV) RNA and stimulates hammerhead ribozyme catalysis. We characterized the nucleic acid binding properties of NdAg, identified the structural and sequence domains important for nucleic acid binding, and studied the correlation between the nucleic acid binding ability and the nucleic acid chaperone activity. NdAg does not recognize the catalytic core of HDV ribozyme specifically. Instead, NdAg interacts with a variety of nucleic acids and has higher affinities to longer nucleic acids. The studies with RNA homopolymers reveal that the binding site size of NdAg is around nine nucleotides long. The extreme N terminal portion of NdAg, the following coiled-coil domain and the basic amino acid clusters in these regions are important for nucleic acid binding. The nucleic acid–NdAg complex is stabilized largely by electrostatic interactions. The formation of RNA–protein complex appears to be a prerequisite for facilitating hammerhead ribozyme catalysis of NdAg and its derivatives. Mutations that reduce the RNA binding activity or high ionic strength that destabilizes the RNA–protein complex, reduce the nucleic acid chaperone activity of NdAg.     

A review of the effects of ATP and hydroxychloroquine on the phase separation of the SARS-CoV-2 nucleocapsid protein

SARS-CoV-2 is the coronavirus causing the ongoing pandemic with > 460 millions of infections and > 6 millions of deaths. SARS-CoV-2 nucleocapsid (N) is the only structural protein which plays essential roles in almost all key steps of the viral life cycle with its diverse functions depending on liquid–liquid phase separation (LLPS) driven by interacting with various nucleic acids. The 419-residue N protein is highly conserved in all variants including delta and omicron, and composed of both folded N-/C-terminal domains (NTD/CTD) as well as three long intrinsically disordered regions (IDRs). Recent results have suggested that its CTD and IDRs are also cryptic nucleic acid–binding domains. In this context, any small molecules capable of interfering in its interaction with nucleic acids are anticipated to modulate its LLPS and associated functions. Indeed, ATP, the energy currency existing at very high concentrations (2–12 mM) in all living cells but absent in viruses, modulates LLPS of N protein, and consequently appears to be evolutionarily hijacked by SARS-CoV-2 to promote its life cycle. Hydroxychloroquine (HCQ) has been also shown to specifically bind NTD and CTD to inhibit their interactions with nucleic acids, as well as to disrupt LLPS. Particularly, the unique structure of the HCQ-CTD complex offers a promising strategy for further design of anti-SARS-CoV-2 drugs with better affinity and specificity. The finding may indicate that LLPS is indeed druggable by small molecules, thus opening up a promising direction for drug discovery/design by targeting LLPS in general.     

Amino acid–base interactions: a three-dimensional analysis of protein–DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid–base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein–DNA complex structures. DNA–backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base–amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein–DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein–DNA interface. In conclusion, the majority of amino acid–base interactions observed follow general principles that apply across all protein–DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are ‘universal’ and ‘context-dependent’. An interactive Web-based atlas of side chain–base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.     

A unified computational view of DNA duplex,triplex, quadruplex and their donor–acceptor interactions

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.     

Synthesis of alanyl nucleobase amino acids and their incorporation into proteins

Proteins which bind to nucleic acids and regulate their structure and functions are numerous and exceptionally important. Such proteins employ a variety of strategies for recognition of the relevant structural elements in their nucleic acid substrates, some of which have been shown to involve rather subtle interactions which might have been difficult to design from first principles. In the present study, we have explored the preparation of proteins containing unnatural amino acids having nucleobase side chains. In principle, the introduction of multiple nucleobase amino acids into the nucleic acid binding domain of a protein should enable these modified proteins to interact with their nucleic acid substrates using Watson-Crick and other base pairing interactions. We describe the synthesis of five alanyl nucleobase amino acids protected in a fashion which enabled their attachment to a suppressor tRNA, and their incorporation into each of two proteins with acceptable efficiencies. The nucleobases studied included cytosine, uracil, thymine, adenine and guanine, i.e. the major nucleobase constituents of DNA and RNA. Dihydrofolate reductase was chosen as one model protein to enable direct comparison of the facility of incorporation of the nucleobase amino acids with numerous other unnatural amino acids studied previously. The Klenow fragment of DNA polymerase I was chosen as a representative DNA binding protein whose mode of action has been studied in detail.     

Direct visualization of Bcl-2 family protein interactions using live cell fluorescent protein redistribution assays

Bcl-2 family proteins have important roles in tumor initiation, progression and resistance to therapy. Pro-survival Bcl-2 proteins are regulated by their interactions with pro-death BH3-only proteins making these protein–protein interactions attractive therapeutic targets. Although these interactions have been extensively characterized biochemically, there is a paucity of tools to assess these interactions in cells. Here, we address this limitation by developing quantitative, high throughput microscopy assays to characterize Bcl-2 and BH3-only protein interactions in live cells. We use fluorescent proteins to label the interacting proteins of interest, enabling visualization and quantification of their mitochondria-localized interactions. Using tool compounds, we demonstrate the suitability of our assays to characterize the cellular activity of putative therapeutic molecules that target the interaction between pro-survival Bcl-2 and pro-death BH3-only proteins. In addition to the relevance of our assays for drug discovery, we anticipate that our work will contribute to an improved understanding of the mechanisms that regulate these important protein–protein interactions within the cell.     

Mass spectrometry‐based protein–protein interaction networks for the study of human diseases

A better understanding of the molecular mechanisms underlying disease is key for expediting the development of novel therapeutic interventions. Disease mechanisms are often mediated by interactions between proteins. Insights into the physical rewiring of protein–protein interactions in response to mutations, pathological conditions, or pathogen infection can advance our understanding of disease etiology, progression, and pathogenesis and can lead to the identification of potential druggable targets. Advances in quantitative mass spectrometry (MS)‐based approaches have allowed unbiased mapping of these disease‐mediated changes in protein–protein interactions on a global scale. Here, we review MS techniques that have been instrumental for the identification of protein–protein interactions at a system‐level, and we discuss the challenges associated with these methodologies as well as novel MS advancements that aim to address these challenges. An overview of examples from diverse disease contexts illustrates the potential of MS‐based protein–protein interaction mapping approaches for revealing disease mechanisms, pinpointing new therapeutic targets, and eventually moving toward personalized applications.     

Photoactivatable Reagents Based on Aryl(trifluoromethyl)diazirines: Synthesis and Application for Studying Nucleic Acid–Protein Interactions

Recent data on the synthesis of photoactivatable derivatives of nucleic acids and proteins on the basis of aryl(trifluoromethyl)diazirines—analogs of nucleosides, nucleotides, oligonucleotides, as well as amino acids and peptides—are reviewed. The synthesis of bi- and polyfunctional photoactivatable reagents, including those containing a cleavable function, designed for postsynthetic modification of biopolymers is described. Data are given on the use of the photoactivatable derivatives for studying nucleic acid–protein interactions by the method of photoaffinity labeling. Special consideration is paid to the results obtained by the authors" team in cooperation with other researchers as well as graduate students of the Chemistry of Natural Products Chair, Chemical Faculty, Moscow State University and destined to solve various scientific tasks in the domain of nucleic acid–protein recognition with the use of photoaffinity crosslinking.     

Switching DNA-binding specificity by unnatural amino acid substitution

The specificity of protein–nucleic acid recognition is believed to originate largely from hydrogen bonding between protein polar atoms, primarily side-chain and polar atoms of nucleic acid bases. One way to design new nucleic acid binding proteins of novel specificity is by structure-guided alterations of the hydrogen bonding patterns of a nucleic acid–protein complex. We have used cI repressor of bacteriophage λ as a model system. In the λ-repressor–DNA complex, the -NH₂ group (hydrogen bond donor) of lysine-4 of λ-repressor forms hydrogen bonds with the amide carbonyl atom of asparagine-55 (acceptor) and the O6 (acceptor) of CG6 of operator site O_L1. Substitution of lysine-4 (two donors) by iso-steric S-(2-hydroxyethyl)-cysteine (one donor and one acceptor), by site-directed mutagenesis and chemical modification, leads to switch of binding specificity of λ-repressor from C:G to T:A at position 6 of O_L1. This suggests that unnatural amino acid substitutions could be a simple way of generating nucleic acid binding proteins of altered specificity.     

An RNA aptamer possessing a novel monovalent cation-mediated fold inhibits lysozyme catalysis by inhibiting the binding of long natural substrates

RNA aptamers are being developed as inhibitors of macromolecular and cellular function, diagnostic tools, and potential therapeutics. Our understanding of the physical nature of this emerging class of nucleic acid–protein complexes is limited; few atomic resolution structures have been reported for aptamers bound to their protein target. Guided by chemical mapping, we systematically minimized an RNA aptamer (Lys1) selected against hen egg white lysozyme. The resultant 59-nucleotide compact aptamer (Lys1.2minE) retains nanomolar binding affinity and the ability to inhibit lysozyme''s catalytic activity. Our 2.0-Å crystal structure of the aptamer–protein complex reveals a helical stem stabilizing two loops to form a protein binding platform that binds lysozyme distal to the catalytic cleft. This structure along with complementary solution analyses illuminate a novel protein–nucleic acid interface; (1) only 410 Ų of solvent accessible surface are buried by aptamer binding; (2) an unusually small fraction (∼18%) of the RNA-protein interaction is electrostatic, consistent with the limited protein phosphate backbone contacts observed in the structure; (3) a single Na⁺ stabilizes the loops that constitute the protein-binding platform, and consistent with this observation, Lys1.2minE–lysozyme complex formation takes up rather than displaces cations at low ionic strength; (4) Lys1.2minE inhibits catalysis of large cell wall substrates but not catalysis of small model substrates; and (5) the helical stem of Lys1.2minE can be shortened to four base pairs (Lys1.2minF) without compromising binding affinity, yielding a 45-nucleotide aptamer whose structure may be an adaptable protein binding platform.     

Compositional complementarity between genomic RNA and coat proteins in positive-sense single-stranded RNA viruses

During packaging in positive-sense single-stranded RNA (+ssRNA) viruses, coat proteins (CPs) interact directly with multiple regions in genomic RNA (gRNA), but the underlying physicochemical principles remain unclear. Here we analyze the high-resolution cryo-EM structure of bacteriophage MS2 and show that the gRNA/CP binding sites, including the known packaging signal, overlap significantly with regions where gRNA nucleobase-density profiles match the corresponding CP nucleobase-affinity profiles. Moreover, we show that the MS2 packaging signal corresponds to the global minimum in gRNA/CP interaction energy in the unstructured state as derived using a linearly additive model and knowledge-based nucleobase/amino-acid affinities. Motivated by this, we predict gRNA/CP interaction sites for a comprehensive set of 1082 +ssRNA viruses. We validate our predictions by comparing them with site-resolved information on gRNA/CP interactions derived in SELEX and CLIP experiments for 10 different viruses. Finally, we show that in experimentally studied systems CPs frequently interact with autologous coding regions in gRNA, in agreement with both predicted interaction energies and a recent proposal that proteins in general tend to interact with own mRNAs, if unstructured. Our results define a self-consistent framework for understanding packaging in +ssRNA viruses and implicate interactions between unstructured gRNA and CPs in the process.     

Discarding functional residues from the substitution table improves predictions of active sites within three-dimensional structures

Substitutions of individual amino acids in proteins may be under very different evolutionary restraints depending on their structural and functional roles. The Environment Specific Substitution Table (ESST) describes the pattern of substitutions in terms of amino acid location within elements of secondary structure, solvent accessibility, and the existence of hydrogen bonds between side chains and neighbouring amino acid residues. Clearly amino acids that have very different local environments in their functional state compared to those in the protein analysed will give rise to inconsistencies in the calculation of amino acid substitution tables. Here, we describe how the calculation of ESSTs can be improved by discarding the functional residues from the calculation of substitution tables. Four categories of functions are examined in this study: protein–protein interactions, protein–nucleic acid interactions, protein–ligand interactions, and catalytic activity of enzymes. Their contributions to residue conservation are measured and investigated. We test our new ESSTs using the program CRESCENDO, designed to predict functional residues by exploiting knowledge of amino acid substitutions, and compare the benchmark results with proteins whose functions have been defined experimentally. The new methodology increases the Z-score by 98% at the active site residues and finds 16% more active sites compared with the old ESST. We also find that discarding amino acids responsible for protein–protein interactions helps in the prediction of those residues although they are not as conserved as the residues of active sites. Our methodology can make the substitution tables better reflect and describe the substitution patterns of amino acids that are under structural restraints only.     

Interactions of Escherichia coli RNA with bacteriophage MS2 coat protein: genomic SELEX

Genomic SELEX is a method for studying the network of nucleic acid–protein interactions within any organism. Here we report the discovery of several interesting and potentially biologically important interactions using genomic SELEX. We have found that bacteriophage MS2 coat protein binds several Escherichia coli mRNA fragments more tightly than it binds the natural, well-studied, phage mRNA site. MS2 coat protein binds mRNA fragments from rffG (involved in formation of lipopolysaccharide in the bacterial outer membrane), ebgR (lactose utilization repressor), as well as from several other genes. Genomic SELEX may yield experimentally induced artifacts, such as molecules in which the fixed sequences participate in binding. We describe several methods (annealing of oligonucleotides complementary to fixed sequences or switching fixed sequences) to eliminate some, or almost all, of these artifacts. Such methods may be useful tools for both randomized sequence SELEX and genomic SELEX.     

Inhibitory effects of archetypical nucleic acid ligands on the interactions of HIV-1 nucleocapsid protein with elements of Psi-RNA

Disrupting the interactions between human immunodeficiency virus type 1 (HIV-1) nucleocapsid (NC) protein and structural elements of the packaging signal (Ψ-RNA) could constitute an ideal strategy to inhibit the functions of this region of the genome leader in the virus life cycle. We have employed electrospray ionization (ESI) Fourier transform mass spectrometry (FTMS) to assess the ability of a series of nucleic acid ligands to bind selected structures of Ψ-RNA and inhibit their specific interactions with NC in vitro. We found that the majority of the ligands included in the study were able to form stable non-covalent complexes with stem–loop 2, 3 and 4 (SL2–4), consistent with their characteristic nucleic acid binding modes. However, only aminoglycosidic antibiotics were capable of dissociating preformed NC•SL3 and NC•SL4 complexes, but not NC•SL2. The apparent specificity of these inhibitory effects is closely dependent on distinctive structural features of the different NC•RNA complexes. The trends observed for the IC₅₀ values correlate very well with those provided by the ligand binding affinities and the dissociation constants of target NC•RNA complexes. This systematic investigation of archetypical nucleic acid ligands provides a valid framework to support the design of novel ligand inhibitors for HIV-1 treatment.