Occurrence and stability of anion–π interactions between phosphate and nucleobases in functional RNA molecules

We present a systematic structural and energetic characterization of phosphate(OP)–nucleobase anion…π stacking interactions in RNAs. We observed OP–nucleobase stacking contacts in a variety of structural motifs other than regular helices and spanning broadly diverse sequence distances. Apart from the stacking between a phosphate and a guanine or a uracil two-residue upstream in specific U-turns, such interactions in RNA have been scarcely characterized to date. Our QM calculations showed an energy minimum at a distance between the OP atom and the nucleobase plane centroid slightly below 3 Å for all the nucleobases. By sliding the OP atom over the nucleobase plane we localized the optimal mutual positioning of the stacked moieties, corresponding to an energy minimum below -6 kcal•mol⁻¹, for all the nucleobases, consistently with the projections of the OP atoms over the different π-rings we observed in experimental occurrences. We also found that the strength of the interaction clearly correlates with its electrostatic component, pointing to it as the most relevant contribution. Finally, as OP–uracil and OP–guanine interactions represent together 86% of the instances we detected, we also proved their stability under dynamic conditions in model systems simulated by state-of-the art DFT-MD calculations.     

Anatomy of noncovalent interactions between the nucleobases or ribose and π-containing amino acids in RNA–protein complexes

A set of >300 nonredundant high-resolution RNA–protein complexes were rigorously searched for π-contacts between an amino acid side chain (W, H, F, Y, R, E and D) and an RNA nucleobase (denoted π–π interaction) or ribose moiety (denoted sugar–π). The resulting dataset of >1500 RNA–protein π-contacts were visually inspected and classified based on the interaction type, and amino acids and RNA components involved. More than 80% of structures searched contained at least one RNA–protein π-interaction, with π–π contacts making up 59% of the identified interactions. RNA–protein π–π and sugar–π contacts exhibit a range in the RNA and protein components involved, relative monomer orientations and quantum mechanically predicted binding energies. Interestingly, π–π and sugar–π interactions occur more frequently with RNA (4.8 contacts/structure) than DNA (2.6). Moreover, the maximum stability is greater for RNA–protein contacts than DNA–protein interactions. In addition to highlighting distinct differences between RNA and DNA–protein binding, this work has generated the largest dataset of RNA–protein π-interactions to date, thereby underscoring that RNA–protein π-contacts are ubiquitous in nature, and key to the stability and function of RNA–protein complexes.     

Insights into the DNA stabilizing contributions of a bicyclic cytosine analogue: crystal structures of DNA duplexes containing 7,8-dihydropyrido [2,3-d]pyrimidin-2-one

The incorporation of the bicyclic cytosine analogue 7,8-dihydropyrido[2,3-d]pyrimidin-2-one (X) into DNA duplexes results in a significant enhancement of their stability (3–4 K per modification). To establish the effects of X on the local hydrogen-bonding and base stacking interactions and the overall DNA conformation, and to obtain insights into the correlation between the structure and stability of X-containing DNA duplexes, the crystal structures of [d(CGCGAATT-X-GCG)]₂ and [d(CGCGAAT-X-CGCG)]₂ have been determined at 1.9–2.9 Å resolutions. In all of the structures, the analogue X base pairs with the purine bases on the opposite strands through Watson–Crick and/or wobble type hydrogen bonds. The additional ring of the X base is stacked on the thymine bases at the 5′-side and overall exhibits greatly enhanced stacking interactions suggesting that this is a major contribution to duplex stabilization.     

New 3-O-substituted xanthone derivatives as promising acetylcholinesterase inhibitors

A new series of 3-O-substituted xanthone derivatives were synthesised and evaluated for their anti-cholinergic activities against acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). The results indicated that the xanthone derivatives possessed good AChE inhibitory activity with eleven of them (5, 8, 11, 17, 19, 21-23, 26-28) exhibited significant effects with the IC₅₀ values ranged 0.88 to 1.28 µM. The AChE enzyme kinetic study of 3-(4-phenylbutoxy)-9H-xanthen-9-one (23) and ethyl 2-((9-oxo-9H-xanthen-3-yl)oxy)acetate (28) showed a mixed inhibition mechanism. Molecular docking study showed that 23 binds to the active site of AChE and interacts via extensive π–π stacking with the indole and phenol side chains of Trp86 and Tyr337, besides the hydrogen bonding with the hydration site and π–π interaction with the phenol side chain of Y72. This study revealed that 3-O-alkoxyl substituted xanthone derivatives are potential lead structures, especially 23 and 28 which can be further developed into potent AChE inhibitors.     

Residual Structure of Unfolded Ubiquitin as Revealed by Hydrogen/Deuterium-Exchange 2D NMR

The characterization of residual structures persistent in unfolded proteins in concentrated denaturant solution is currently an important issue in studies of protein folding because the residual structure present, if any, in the unfolded state may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the hydrogen/deuterium (H/D)-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride. We employed a dimethylsulfoxide (DMSO)-quenched H/D-exchange NMR technique with the use of spin desalting columns, which allowed us to perform a quick medium exchange from 6 M guanidinium chloride to a quenching DMSO solution. Based on the backbone resonance assignment of ubiquitin in the DMSO solution, we successfully investigated the H/D-exchange kinetics of 60 identified peptide amide groups in the ubiquitin sequence. Although a majority of these amide groups were not protected, certain amide groups involved in a middle helix (residues 23–34) and an N-terminal β-hairpin (residues 2–16) were significantly protected with a protection factor of 2.1–4.2, indicating that there were residual structures in unfolded ubiquitin and that these amide groups were more than 52% hydrogen bonded in the residual structures. We show that the hydrogen-bonded residual structures in the α-helix and the β-hairpin are formed even in 6 M guanidinium chloride, suggesting that these residual structures may function as a folding initiation site to guide the subsequent folding reactions of ubiquitin.     

Exploring the Influence of Carbon Nanoparticles on the Formation of β-Sheet-Rich Oligomers of IAPP22–28 Peptide by Molecular Dynamics Simulation

Recent advances in nanotechnologies have led to wide use of nanomaterials in biomedical field. However, nanoparticles are found to interfere with protein misfolding and aggregation associated with many human diseases. It is still a controversial issue whether nanoparticles inhibit or promote protein aggregation. In this study, we used molecular dynamics simulations to explore the effects of three kinds of carbon nanomaterials including graphene, carbon nanotube and C₆₀ on the aggregation behavior of islet amyloid polypeptide fragment 22–28 (IAPP_22–28). The diverse behaviors of IAPP_22–28 peptides on the surfaces of carbon nanomaterials were studied. The results suggest these nanomaterials can prevent β-sheet formation in differing degrees and further affect the aggregation of IAPP_22–28. The π–π stacking and hydrophobic interactions are different in the interactions between peptides and different nanoparticles. The subtle differences in the interaction are due to the difference in surface curvature and area. The results demonstrate the adsorption interaction has competitive advantages over the interactions between peptides. Therefore, the fibrillation of IAPP_22–28 may be inhibited at its early stage by graphene or SWCNT. Our study can not only enhance the understanding about potential effects of nanomaterials to amyloid formation, but also provide valuable information to develop potential β-sheet formation inhibitors against type II diabetes.     

Structure,Dynamics, and Interaction of Mycobacterium tuberculosis (Mtb) DprE1 and DprE2 Examined by Molecular Modeling,Simulation, and Electrostatic Studies

The enzymes decaprenylphosphoryl-β-D-ribose oxidase (DprE1) and decaprenylphosphoryl-β-D-ribose-2-epimerase (DprE2) catalyze epimerization of decaprenylphosporyl ribose (DPR) todecaprenylphosporyl arabinose (DPA) and are critical for the survival of Mtb. Crystal structures of DprE1 so far reported display significant disordered regions and no structural information is known for DprE2. We used homology modeling, protein threading, molecular docking and dynamics studies to investigate the structural and dynamic features of Mtb DprE1 and DprE2 and DprE1-DprE2 complex. A three-dimensional model for DprE2 was generated using the threading approach coupled with ab initio modeling. A 50 ns simulation of DprE1 and DprE2 revealed the overall stability of the structures. Principal Component Analysis (PCA) demonstrated the convergence of sampling in both DprE1 and DprE2. In DprE1, residues in the 269–330 area showed considerable fluctuation in agreement with the regions of disorder observed in the reported crystal structures. In DprE2, large fluctuations were detected in residues 95–113, 146–157, and 197–226. The study combined docking and MD simulation studies to map and characterize the key residues involved in DprE1-DprE2 interaction. A 60 ns MD simulation for DprE1-DprE2 complex was also performed. Analysis of data revealed that the docked complex is stabilized by H-bonding, hydrophobic and ionic interactions. The key residues of DprE1 involved in DprE1-DprE2 interactions belong to the disordered region. We also examined the docked complex of DprE1-BTZ043 to investigate the binding pocket of DprE1 and its interactions with the inhibitor BTZ043. In summary, we hypothesize that DprE1-DprE2 interaction is crucial for the synthesis of DPA and DprE1-DprE2 complex may be a new therapeutic target amenable to pharmacological validation. The findings have important implications in tuberculosis (TB) drug discovery and will facilitate drug development efforts against TB.     

Solution and Crystallographic Structures of the Central Region of the Phosphoprotein from Human Metapneumovirus

Human metapneumovirus (HMPV) of the family Paramyxoviridae is a major cause of respiratory illness worldwide. Phosphoproteins (P) from Paramyxoviridae are essential co-factors of the viral RNA polymerase that form tetramers and possess long intrinsically disordered regions (IDRs). We located the central region of HMPV P (P_ced) which is involved in tetramerization using disorder analysis and modeled its 3D structure ab initio using Rosetta fold-and-dock. We characterized the solution-structure of P_ced using small angle X-ray scattering (SAXS) and carried out direct fitting to the scattering data to filter out incorrect models. Molecular dynamics simulations (MDS) and ensemble optimization were employed to select correct models and capture the dynamic character of P_ced. Our analysis revealed that oligomerization involves a compact central core located between residues 169-194 (P_core), that is surrounded by flexible regions with α-helical propensity. We crystallized this fragment and solved its structure at 3.1 Å resolution by molecular replacement, using the folded core from our SAXS-validated ab initio model. The RMSD between modeled and experimental tetramers is as low as 0.9 Å, demonstrating the accuracy of the approach. A comparison of the structure of HMPV P to existing mononegavirales P_ced structures suggests that P_ced evolved under weak selective pressure. Finally, we discuss the advantages of using SAXS in combination with ab initio modeling and MDS to solve the structure of small, homo-oligomeric protein complexes.     

Bacillus subtilis tau subunit of DNA polymerase III interacts with bacteriophage SPP1 replicative DNA helicase G40P

Genetic evidence suggests that the Bacillus subtilis dnaX gene only encodes for the τ subunit of both DNA polymerases III (Pol IIIs). The B.subtilis full-length protein and their mutant derivatives τ(373– 563) (lacking the N-terminal, domains I–III or amino acid residues 1–372) and τ(1–372) (lacking the C-terminal region or amino acids 373–563) have been purified. The τ protein forms tetramers, τ(373– 563) forms dimers, whereas τ(1–372), depending on the ionic strength, forms trimers or tetramers in solution. In the absence of single-stranded (ss) DNA and a nucleotide cofactor, τ interacts with the SPP1 hexameric replicative G40P DNA helicase in solution or with G40P-ATP bound to ssDNA, with a 1:1 stoichiometry. G40P(109–442), lacking the N-terminal amino acid residues 1–108, interacts with the C-terminal moiety of τ. The data indicate that the interaction of G40P with the τ subunit of Pol III, is relevant for the loading of the Pol IIIs into the SPP1 G38P-promoted open complex.     

Interactions between the 2.4 and 4.2 regions of sigmaS, the stress-specific sigma factor of Escherichia coli, and the -10 and -35 promoter elements

The σ^s subunit of Escherichia coli RNA polymerase holoenzyme (Eσ^S) is a key factor of gene expression upon entry into stationary phase and in stressful conditions. The selectivity of promoter recognition by Eσ^S and the housekeeping Eσ⁷⁰ is as yet not clearly understood. We used a genetic approach to investigate the interaction of σ^S with its target promoters. Starting with down-promoter variants of a σ^S promoter target, osmEp, altered in the –10 or –35 elements, we isolated mutant forms of σ^S suppressing the promoter defects. The activity of these suppressors on variants of osmEp and ficp, another target of σ^S, indicated that σ^S is able to interact with the same key features within a promoter sequence as σ⁷⁰. Indeed, (i) σ^S can recognize the –35 element of some but not all its target promoters, through interactions with its 4.2 region; and (ii) amino acids within the 2.4 region participate in the recognition of the –10 element. More specifically, residues Q152 and E155 contribute to the strong preference of σ^S for a C in position –13 and residue R299 can interact with the –31 nucleotide in the –35 element of the target promoters.     

Aromatic residues in RNase T stack with nucleobases to guide the sequence‐specific recognition and cleavage of nucleic acids

RNase T is a classical member of the DEDDh family of exonucleases with a unique sequence preference in that its 3′‐to‐5′ exonuclease activity is blocked by a 3′‐terminal dinucleotide CC in digesting both single‐stranded RNA and DNA. Our previous crystal structure analysis of RNase T‐DNA complexes show that four phenylalanine residues, F29, F77, F124, and F146, stack with the two 3′‐terminal nucleobases. To elucidate if the π–π stacking interactions between aromatic residues and nucleobases play a critical role in sequence‐specific protein–nucleic acid recognition, here we mutated two to four of the phenylalanine residues in RNase T to tryptophan (W mutants) and tyrosine (Y mutants). The Escherichia coli strains expressing either the W mutants or the Y mutants had slow growth phenotypes, suggesting that all of these mutants could not fully substitute the function of the wild‐type RNase T in vivo. DNA digestion assays revealed W mutants shared similar sequence specificity with wild‐type RNase T. However, the Y mutants exhibited altered sequence‐dependent activity, digesting ssDNA with both 3′‐end CC and GG sequences. Moreover, the W and Y mutants had reduced DNA‐binding activity and lower thermal stability as compared to wild‐type RNase T. Taken together, our results suggest that the four phenylalanine residues in RNase T not only play critical roles in sequence‐specific recognition, but also in overall protein stability. Our results provide the first evidence showing that the π−π stacking interactions between nucleobases and protein aromatic residues may guide the sequence‐specific activity for DNA and RNA enzymes.     

N(omega)-arginine dimethylation modulates the interaction between a Gly/Arg-rich peptide from human nucleolin and nucleic acids

We studied the interaction between a synthetic peptide (sequence Ac-GXGGFGGXGGFXGGXGG-NH₂, where X = arginine, N^ω,N^ω-dimethylarginine, DMA, or lysine) corresponding to residues 676–692 of human nucleolin and several DNA and RNA substrates using double filter binding, melting curve analysis and circular dichroism spectroscopy. We found that despite the reduced capability of DMA in forming hydrogen bonds, N^ω,N^ω-dimethylation does not affect the strength of the binding to nucleic acids nor does it have any effect on stabilization of a double-stranded DNA substrate. However, circular dichroism studies show that unmethylated peptide can perturb the helical structure, especially in RNA, to a much larger extent than the DMA peptide.     

Helical coherence of DNA in crystals and solution

The twist, rise, slide, shift, tilt and roll between adjoining base pairs in DNA depend on the identity of the bases. The resulting dependence of the double helix conformation on the nucleotide sequence is important for DNA recognition by proteins, packaging and maintenance of genetic material, and other interactions involving DNA. This dependence, however, is obscured by poorly understood variations in the stacking geometry of the same adjoining base pairs within different sequence contexts. In this article, we approach the problem of sequence-dependent DNA conformation by statistical analysis of X-ray and NMR structures of DNA oligomers. We evaluate the corresponding helical coherence length—a cumulative parameter quantifying sequence-dependent deviations from the ideal double helix geometry. We find, e.g. that the solution structure of synthetic oligomers is characterized by 100–200 Å coherence length, which is similar to ~150 Å coherence length of natural, salmon-sperm DNA. Packing of oligomers in crystals dramatically alters their helical coherence. The coherence length increases to 800–1200 Å, consistent with its theoretically predicted role in interactions between DNA at close separations.     

Effects of Hydroxylated Carbon Nanotubes on the Aggregation of Aβ16–22 Peptides: A Combined Simulation and Experimental Study

The pathogenesis of Alzheimer’s disease (AD) is associated with the aggregation of amyloid-β (Aβ) peptides into toxic aggregates with β-sheet character. In a previous computational study, we showed that pristine single-walled carbon nanotubes (SWCNTs) can inhibit the formation of β-sheet-rich oligomers in the central hydrophobic core fragment of Aβ (Aβ_16–22). However, the poor solubility of SWCNTs in water hinders their use in biomedical applications and nanomedicine. Here, we investigate the influence of hydroxylated SWCNT, a water-soluble SWCNT derivative, on the aggregation of Aβ_16–22 peptides using all-atom explicit-water replica exchange molecular dynamics simulations. Our results show that hydroxylated SWCNTs can significantly inhibit β-sheet formation and shift the conformations of Aβ_16–22 oligomers from ordered β-sheet-rich structures toward disordered coil aggregates. Detailed analyses of the SWCNT-Aβ interaction reveal that the inhibition of β-sheet formation by hydroxylated SWCNTs mainly results from strong electrostatic interactions between the hydroxyl groups of SWCNTs and the positively charged residue K16 of Aβ_16–22 and hydrophobic and aromatic stacking interactions between SWCNTs and F19 and F20. In addition, our atomic force microscopy and thioflavin T fluorescence experiments confirm the inhibitory effect of both pristine and hydroxylated SWCNTs on Aβ_16–22 fibrillization, in support of our previous and present replica exchange molecular dynamics simulation results. These results demonstrate that hydroxylated SWCNTs efficiently inhibit the aggregation of Aβ_16–22; in addition, they offer molecular insight into the inhibition mechanism, thus providing new clues for the design of therapeutic drugs against amyloidosis.     

Predicting future community-level ocular Chlamydia trachomatis infection prevalence using serological,clinical, molecular,and geospatial data

Trachoma is an infectious disease characterized by repeated exposures to Chlamydia trachomatis (Ct) that may ultimately lead to blindness. Efficient identification of communities with high infection burden could help target more intensive control efforts. We hypothesized that IgG seroprevalence in combination with geospatial layers, machine learning, and model-based geostatistics would be able to accurately predict future community-level ocular Ct infections detected by PCR. We used measurements from 40 communities in the hyperendemic Amhara region of Ethiopia to assess this hypothesis. Median Ct infection prevalence among children 0–5 years old increased from 6% at enrollment, in the context of recent mass drug administration (MDA), to 29% by month 36, following three years without MDA. At baseline, correlation between seroprevalence and Ct infection was stronger among children 0–5 years old (ρ = 0.77) than children 6–9 years old (ρ = 0.48), and stronger than the correlation between active trachoma and Ct infection (0-5y ρ = 0.56; 6-9y ρ = 0.40). Seroprevalence was the strongest concurrent predictor of infection prevalence at month 36 among children 0–5 years old (cross-validated R² = 0.75, 95% CI: 0.58–0.85), though predictive performance declined substantially with increasing temporal lag between predictor and outcome measurements. Geospatial variables, a spatial Gaussian process, and stacked ensemble machine learning did not meaningfully improve predictions. Serological markers among children 0–5 years old may be an objective tool for identifying communities with high levels of ocular Ct infections, but accurate, future prediction in the context of changing transmission remains an open challenge.     

Structure–activity relationships at a nucleobase-stacking tryptophan required for chemomechanical coupling in the DNA resecting motor-nuclease AdnAB

Mycobacterial AdnAB is a heterodimeric helicase-nuclease that initiates homologous recombination by resecting DNA double-strand breaks. The AdnB subunit hydrolyzes ATP to drive single-nucleotide steps of 3′-to-5′ translocation of AdnAB on the tracking DNA strand via a ratchet-like mechanism. Trp325 in AdnB motif III, which intercalates into the tracking strand and makes a π stack on a nucleobase 5′ of a flipped-out nucleoside, is the putative ratchet pawl without which ATP hydrolysis is mechanically futile. Here, we report that AdnAB mutants wherein Trp325 was replaced with phenylalanine, tyrosine, histidine, leucine, or alanine retained activity in ssDNA-dependent ATP hydrolysis but displayed a gradient of effects on DSB resection. The resection velocities of Phe325 and Tyr325 mutants were 90% and 85% of the wild-type AdnAB velocity. His325 slowed resection rate to 3% of wild-type and Leu325 and Ala325 abolished DNA resection. A cryo-EM structure of the DNA-bound Ala325 mutant revealed that the AdnB motif III peptide was disordered and the erstwhile flipped out tracking strand nucleobase reverted to a continuous base-stacked arrangement with its neighbors. We conclude that π stacking of Trp325 on a DNA nucleobase triggers and stabilizes the flipped-out conformation of the neighboring nucleoside that underlies formation of a ratchet pawl.     

Discovery and Structure Activity Relationship of Small Molecule Inhibitors of Toxic β-Amyloid-42 Fibril Formation

Increasing evidence implicates Aβ peptides self-assembly and fibril formation as crucial events in the pathogenesis of Alzheimer disease. Thus, inhibiting Aβ aggregation, among others, has emerged as a potential therapeutic intervention for this disorder. Herein, we employed 3-aminopyrazole as a key fragment in our design of non-dye compounds capable of interacting with Aβ42 via a donor-acceptor-donor hydrogen bond pattern complementary to that of the β-sheet conformation of Aβ42. The initial design of the compounds was based on connecting two 3-aminopyrazole moieties via a linker to identify suitable scaffold molecules. Additional aryl substitutions on the two 3-aminopyrazole moieties were also explored to enhance π-π stacking/hydrophobic interactions with amino acids of Aβ42. The efficacy of these compounds on inhibiting Aβ fibril formation and toxicity in vitro was assessed using a combination of biophysical techniques and viability assays. Using structure activity relationship data from the in vitro assays, we identified compounds capable of preventing pathological self-assembly of Aβ42 leading to decreased cell toxicity.     

Insight into F plasmid DNA segregation revealed by structures of SopB and SopB–DNA complexes

Accurate DNA segregation is essential for genome transmission. Segregation of the prototypical F plasmid requires the centromere-binding protein SopB, the NTPase SopA and the sopC centromere. SopB displays an intriguing range of DNA-binding properties essential for partition; it binds sopC to form a partition complex, which recruits SopA, and it also coats DNA to prevent non-specific SopA–DNA interactions, which inhibits SopA polymerization. To understand the myriad functions of SopB, we determined a series of SopB–DNA crystal structures. SopB does not distort its DNA site and our data suggest that SopB–sopC forms an extended rather than wrapped partition complex with the SopA-interacting domains aligned on one face. SopB is a multidomain protein, which like P1 ParB contains an all-helical DNA-binding domain that is flexibly attached to a compact (β₃–α)₂ dimer-domain. Unlike P1 ParB, the SopB dimer-domain does not bind DNA. Moreover, SopB contains a unique secondary dimerization motif that bridges between DNA duplexes. Both specific and non-specific SopB–DNA bridging structures were observed. This DNA-linking function suggests a novel mechanism for in trans DNA spreading by SopB, explaining how it might mask DNA to prevent DNA-mediated inhibition of SopA polymerization.