DNA-protein π-interactions in nature: abundance, structure, composition and strength of contacts between aromatic amino acids and DNA nucleobases or deoxyribose sugar

Targeting and Manipulating Tryptophan Interactions on Proteins

Tryptophan, commonly regarded as buried within the interior cores of proteins to maintain secondary structures, is now being recognized for its significant contributions to protein functionality. However, investigating functional tryptophan-involved interactions across the proteome and manipulating these interactions in live cells are considerable challenges. In this In Focus article, we summarize emerging advances in the field, describing innovative chemistries that leverage distinctive biochemical properties of the indole moiety for targeting and functionally manipulating tryptophan interactions.

Symmetry-breaking malachite green as a near-infrared light-activated fluorogenic photosensitizer for RNA proximity labeling

Cellular RNA is asymmetrically distributed in cells and the regulation of RNA localization is crucial for proper cellular functions. However, limited chemical tools are available to capture dynamic RNA localization in complex biological systems with high spatiotemporal resolution. Here, we developed a new method for RNA proximity labeling activated by near-infrared (NIR) light, which holds the potential for deep penetration. Our method, termed FAP-seq, utilizes a genetically encoded fluorogen activating protein (FAP) that selectively binds to a set of substrates known as malachite green (MG). FAP binding restricts the rotation of MG and rapidly activates its fluorescence in a wash-free manner. By introducing a monoiodo modification to MG, we created a photosensitizer (MG-HI) with the highest singlet oxygen generation ability among various MG derivatives, enabling both protein and RNA proximity labeling in live cells. New insights are provided in the transcriptome analysis with FAP-seq, while a deeper understanding of the symmetry-breaking structural arrangement of FAP-MG-HI was obtained through molecular dynamics simulations. Overall, our wash-free and NIR light-inducible RNA proximity labeling method (FAP-seq) offers a powerful and versatile approach for investigating complex mechanisms underlying RNA-related biological processes.

Genetically encoded bioorthogonal tryptophan decaging in living cells

Tryptophan (Trp) plays a critical role in the regulation of protein structure, interactions and functions through its π system and indole N-H group. A generalizable method for blocking and rescuing Trp interactions would enable the gain-of-function manipulation of various Trp-containing proteins in vivo, but generating such a platform remains challenging. Here we develop a genetically encoded N1-vinyl-caged Trp capable of rapid and bioorthogonal decaging through an optimized inverse electron-demand Diels-Alder reaction, allowing site-specific activation of Trp on a protein of interest in living cells. This chemical activation of a genetically encoded caged-tryptophan (Trp-CAGE) strategy enables precise activation of the Trp of interest underlying diverse important molecular interactions. We demonstrate the utility of Trp-CAGE across various protein families, such as catalase-peroxidases and kinases, as translation initiators and posttranslational modification readers, allowing the modulation of epigenetic signalling in a temporally controlled manner. Coupled with computer-aided prediction, our strategy paves the way for bioorthogonal Trp activation on more than 28,000 candidate proteins within their native cellular settings.

The RNA-Binding Domain of hnRNP U Extends beyond the RGG/RG Motifs

Heterogeneous nuclear ribonucleoprotein U (hnRNP U) is a ubiquitously expressed protein that regulates chromatin architecture through its interactions with numerous DNA, protein, and RNA partners. The RNA-binding domain (RBD) of hnRNP U was previously mapped to an RGG/RG motif within its disordered C-terminal region, but little is understood about its binding mode and potential for selective RNA recognition. Analysis of publicly available hnRNP U enhanced UV cross-linking and immunoprecipitation (eCLIP) data identified high-confidence binding sites within human RNAs. We synthesized a set of diverse RNAs encompassing 11 of these identified cross-link sites for biochemical characterization using a combination of fluorescence anisotropy and electrophoretic mobility shift assays. These in vitro binding experiments with a rationally designed set of RNAs and hnRNP U domains revealed that the RGG/RG motif is a small part of a more expansive RBD that encompasses most of the disordered C-terminal region. This RBD contains a second, previously experimentally uncharacterized RGG/RG motif with RNA-binding properties comparable to those of the canonical RGG/RG motif. These RGG/RG motifs serve redundant functions, with neither serving as the primary RBD. While in isolation, each RGG/RG motif has modest affinity for RNA, together they significantly enhance the association of hnRNP U with RNA, enabling the binding of most of the designed RNA set with low to midnanomolar binding affinities. Identification and characterization of the complete hnRNP U RBD highlight the perils of a reductionist approach to defining biochemical activities in this system and pave the way for a detailed investigation of its RNA-binding specificity.

Single molecule technique unveils the role of electrostatic interactions in ssDNA-gp32 molecular complex stability

The exploration of single-strand DNA-binding protein (SSB)-ssDNA interactions and their crucial roles in essential biological processes lagged behind other types of protein-nucleic acid interactions, such as protein-dsDNA and protein-RNA interactions. The ssDNA binding protein gene product 32 (gp32) of the T4 bacteriophage is a central integrating component of the replication complex that must continuously bind to and unbind from transiently exposed template strands during the DNA synthesis. To gain deeper insights into the electrostatic conditions influencing the stability of the ssDNA-gp32 molecular complex, like the salt concentration or some metal ions proven to specifically bind to gp32, we employed a method that performs rapid measurements of the DNA-protein stability using an α-Hemolysin (α-HL) protein nanopore. We indirectly probed the stability of a protein-nucleic acid complex by monitoring the dissociation process between the gp32 protein and the ssDNA molecular complex in single-molecular electrophysiology experiments, but also through fluorescence spectroscopy techniques. We have shown that the complex is more stable in 0.5 M KCl solution than in 2 M KCl solution and that the presence of Zn2+ ions further increases this stability for any salt used in the present study. This method can be applied to other nucleic acid-protein molecular complexes, as well as for an accurate determination of the drug-protein carrier stability.

An Atlas of the base inter-RNA stacks involved in bacterial translation

Nucleobase-specific noncovalent interactions play a crucial role in translation. Herein, we provide a comprehensive analysis of the stacks between different RNA components in the crystal structures of the bacterial ribosome caught at different translation stages. Analysis of tRNA||rRNA stacks reveals distinct behaviour; both the A-and E-site tRNAs exhibit unique stacking patterns with 23S rRNA bases, while P-site tRNAs stack with 16S rRNA bases. Furthermore, E-site stacks exhibit diverse face orientations and ring topologies-rare for inter-chain RNA interactions-with higher average interaction energies than A or P-site stacks. This suggests that stacking may be essential for stabilizing tRNA progression through the E-site. Additionally, mRNA||rRNA stacks reveal other geometries, which depend on the tRNA binding site, whereas 16S rRNA||23S rRNA stacks highlight the importance of specific bases in maintaining the integrity of the translational complex by linking the two rRNAs. Furthermore, tRNA||mRNA stacks exhibit distinct geometries and energetics at the E-site, indicating their significance during tRNA translocation and elimination. Overall, both A and E-sites display a more diverse distribution of inter-RNA stacks compared to the P-site. Stacking interactions in the active ribosome are not simply accidental byproducts of biochemistry but are likely invoked to compensate and support the integrity and dynamics of translation.

PRA-Pred: Structure-based prediction of protein-RNA binding affinity

Understanding crucial factors that affect the binding affinity of protein-RNA complexes is vital for comprehending their recognition mechanisms. This study involved compiling experimentally measured binding affinity (ΔG) values of 217 protein-RNA complexes and extracting numerous structure-based features, considering RNA, protein, and interactions between protein and RNA. Our findings indicate the significance of RNA base-step parameters, interaction energies, number of atomic contacts in the complex, hydrogen bonds, and contact potentials in understanding the binding affinity. Further, we observed that these factors are influenced by the type of RNA strand and the function of the protein in a protein-RNA complex. Multiple regression equations were developed for different classes of complexes to perform the prediction of the binding affinity between the protein and RNA. We evaluated the models using the jack-knife test and achieved an overall correlation 0.77 between the experimental and predicted binding affinities with a mean absolute error of 1.02 kcal/mol. Furthermore, we introduced a web server, PRA-Pred, intended for the prediction of protein-RNA binding affinity, and it is freely accessible through https://web.iitm.ac.in/bioinfo2/prapred/. We propose that our approach could function as a potential resource for investigating protein-RNA recognitions and developing therapeutic strategies.

Analyzing the aromatic-aromatic interactions in proteins: A2ID 2.0

The Aromatic-Aromatic Interactions Database (A2ID) is a comprehensive repository dedicated to documenting aromatic-aromatic (π-π) networks observed in experimentally determined protein structures. The first version of A2ID was reported in 2011 [Int J Biol Macromol, 2011, 48, 540]. It has undergone a series of significant updates, leading to its current version, which focuses on the identification and analysis of 3,444,619 π-π networks from proteins. The geometrical parameters such as centroid-centroid distances (r) and interplanar angles (ϕ) were used to identify and characterize π-π networks. It was observed that among the 84,500 proteins with at least one aromatic π-π network, about 92.50 % of the instances are found to be either 2π (77.34 %) or 3π (15.23 %) networks. The analysis of interacting amino acid pairs in 2π networks indicated a dominance of PHE residues followed by TYR. The updated version of A2ID incorporates analysis of π-π networks based on SCOP2 and ECOD classifiers, in addition to the existing SCOP, CATH, and EC classifications. This expanded scope allows researchers to explore the characteristics and functional implications of π-π networks in protein structures from multiple perspectives. The current version of A2ID along with its extensive dataset and detailed geometric information is publicly accessible using https://acds.neist.res.in/a2idv2.

Neutrophil Extracellular Traps have DNAzyme activity that drives bactericidal potential

The mechanisms of bacterial killing by neutrophil extracellular traps (NETs) are unclear. DNA, the largest component of NETs is believed to merely be a scaffold with minimal antimicrobial activity through the charge of the backbone. Here, we report that NETs DNA is beyond a scaffold and produces hydroxyl free radicals through the spatially concentrated G-quadruplex/hemin DNAzyme complexes, driving bactericidal effects. Immunofluorescence staining showed colocalization of G-quadruplex and hemin in extruded NETs DNA, and Amplex UltraRed assay portrayed its peroxidase activity. Proximity labeling of bacteria revealed localized concentration of radicals resulting from NETs bacterial trapping. Ex vivo bactericidal assays revealed that G-quadruplex/hemin DNAzyme is the primary driver of bactericidal activity in NETs. NETs are DNAzymes that may have important biological consequences.

PAM-flexible genome editing with an engineered chimeric Cas9

CRISPR enzymes require a defined protospacer adjacent motif (PAM) flanking a guide RNA-programmed target site, limiting their sequence accessibility for robust genome editing applications. In this study, we recombine the PAM-interacting domain of SpRY, a broad-targeting Cas9 possessing an NRN > NYN (R = A or G, Y = C or T) PAM preference, with the N-terminus of Sc + +, a Cas9 with simultaneously broad, efficient, and accurate NNG editing capabilities, to generate a chimeric enzyme with highly flexible PAM preference: SpRYc. We demonstrate that SpRYc leverages properties of both enzymes to specifically edit diverse PAMs and disease-related loci for potential therapeutic applications. In total, the approaches to generate SpRYc, coupled with its robust flexibility, highlight the power of integrative protein design for Cas9 engineering and motivate downstream editing applications that require precise genomic positioning.

The RNA-binding domain of hnRNP U extends beyond the RGG/RG motifs

Heterogeneous nuclear ribonucleoprotein U (hnRNP U) is a ubiquitously expressed protein that regulates chromatin architecture through its interactions with numerous DNA, protein, and RNA partners. The RNA-binding domain (RBD) of hnRNP U was previously mapped to an RGG/RG element within its disordered C-terminal region, but little is understood about its binding mode and potential for selective RNA recognition. Analysis of publicly available hnRNP U enhanced UV crosslinking and immunoprecipitation (eCLIP) data identified high-confidence binding sites within human RNAs. We synthesized a set of diverse RNAs encompassing eleven of these identified crosslink sites for biochemical characterization using a combination of fluorescence anisotropy and electrophoretic mobility shift assays. These in vitro binding experiments with a rationally designed set of RNAs and hnRNP U domains revealed that the RGG/RG element is a small part of a more expansive RBD that encompasses most of the disordered C-terminal region. This RBD contains a second, previously experimentally uncharacterized RGG/RG element with RNA-binding properties comparable to the canonical RGG/RG element. These RGG/RG elements serve redundant functions, with neither serving as the primary RBD. While in isolation each RGG/RG element has modest affinity for RNA, together they significantly enhance the association of hnRNP U with RNA, enabling binding of most of the designed RNA set with low to mid-nanomolar binding affinities. Identification and characterization of the complete hnRNP U RBD highlights the perils of a reductionist approach to defining biochemical activities in this system and paves the way for a detailed investigation of its RNA-binding specificity.

Aromatic Residues Dictate the Transcriptional Repressor and Single-Stranded DNA Binding Activities of Purine-Rich Element Binding Protein B

Purine-rich element binding protein B (Purβ) is a single-stranded DNA (ssDNA) and RNA-binding protein that functions as a transcriptional repressor of genes encoding certain muscle-restricted contractile proteins in the setting of cellular stress or tissue injury. A prior report from our laboratory implicated specific basic amino acid residues in the physical and functional interaction of Purβ with the smooth muscle-α actin gene (Acta2) promoter. Independent structural analysis of fruit fly Purα uncovered a role for several aromatic residues in the binding of this related protein to ssDNA. Herein, we examine the functional importance of a comparable set of hydrophobic residues that are positionally conserved in the repeat I (Y59), II (F155), and III (F256) domains of murine Purβ. Site-directed Y/F to alanine substitutions were engineered, and the resultant Purβ point mutants were tested in various biochemical and cell-based assays. None of the mutations affected the cellular expression, structural stability, or dimerization capacity of Purβ. However, the Y59A and F155A mutants demonstrated weaker Acta2 repressor activity in transfected fibroblasts and reduced binding affinity for the purine-rich strand of an Acta2 cis-regulatory element in vitro. Mutation of Y59 and F155 also altered the multisite binding properties of Purβ for ssDNA and diminished the interaction of Purβ with Y-box binding protein 1, a co-repressor of Acta2. Collectively, these findings suggest that some of the same aromatic residues, which govern the specific and high-affinity binding of Purβ to ssDNA, also mediate certain heterotypic protein interactions underlying the Acta2 repressor function of Purβ.

Occurrence and classification of T-shaped interactions between nucleobases in RNA structures

Understanding the frequency and structural context of discrete noncovalent interactions between nucleotides is of pivotal significance in establishing the rules that govern RNA structure and dynamics. Although T-shaped contacts (i.e., perpendicular stacking contacts) between aromatic amino acids and nucleobases at the nucleic acid-protein interface have recently garnered attention, the analogous contacts within the nucleic acid structures have not been discussed. In this work, we have developed an automated method for identifying and unambiguously classifying T-shaped interactions between nucleobases. Using this method, we identified a total of 3261 instances of T-shaped (perpendicular stacking) contacts between two nucleobases in an array of RNA structures from an up-to-date data set of ≤3.5 Å resolution crystal structures deposited in the Protein Data Bank.

Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA

Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Seven-membered ring nucleobases as inhibitors of human cytidine deaminase and APOBEC3A

The APOBEC3 (APOBEC3A-H) enzyme family as a part of the human innate immune system deaminates cytosine to uracil in single-stranded DNA (ssDNA) and thereby prevents the spread of pathogenic genetic information. However, APOBEC3-induced mutagenesis promotes viral and cancer evolution, thus enabling the progression of diseases and development of drug resistance. Therefore, APOBEC3 inhibition offers a possibility to complement existing antiviral and anticancer therapies and prevent the emergence of drug resistance, thus making such therapies effective for longer periods of time. Here, we synthesised nucleosides containing seven-membered nucleobases based on azepinone and compared their inhibitory potential against human cytidine deaminase (hCDA) and APOBEC3A with previously described 2'-deoxyzebularine (dZ) and 5-fluoro-2'-deoxyzebularine (FdZ). The nanomolar inhibitor of wild-type APOBEC3A was obtained by the incorporation of 1,3,4,7-tetrahydro-2H-1,3-diazepin-2-one in the TTC loop of a DNA hairpin instead of the target 2'-deoxycytidine providing a Ki of 290 ± 40 nM, which is only slightly weaker than the Ki of the FdZ-containing inhibitor (117 ± 15 nM). A less potent but notably different inhibition of human cytidine deaminase (CDA) and engineered C-terminal domain of APOBEC3B was observed for 2'-deoxyribosides of the S and R isomers of hexahydro-5-hydroxy-azepin-2-one: the S-isomer was more active than the R-isomer. The S-isomer shows resemblance in the position of the OH-group observed recently for the hydrated dZ and FdZ in the crystal structures with APOBEC3G and APOBEC3A, respectively. This shows that 7-membered ring analogues of pyrimidine nucleosides can serve as a platform for further development of modified ssDNAs as powerful A3 inhibitors.

PDA-Pred: Predicting the binding affinity of protein-DNA complexes using machine learning techniques and structural features

Protein-DNA interactions play an important role in various biological processes such as gene expression, replication, and transcription. Understanding the important features that dictate the binding affinity of protein-DNA complexes and predicting their affinities is important for elucidating their recognition mechanisms. In this work, we have collected the experimental binding free energy (ΔG) for a set of 391 Protein-DNA complexes and derived several structure-based features such as interaction energy, contact potentials, volume and surface area of binding site residues, base step parameters of the DNA and contacts between different types of atoms. Our analysis on relationship between binding affinity and structural features revealed that the important factors mainly depend on the number of DNA strands as well as functional and structural classes of proteins. Specifically, binding site properties such as number of atom contacts between the DNA and protein, volume of protein binding sites and interaction-based features such as interaction energies and contact potentials are important to understand the binding affinity. Further, we developed multiple regression equations for predicting the binding affinity of protein-DNA complexes belonging to different structural and functional classes. Our method showed an average correlation and mean absolute error of 0.78 and 0.98 kcal/mol, respectively, between the experimental and predicted binding affinities on a jack-knife test. We have developed a webserver, PDA-PreD (Protein-DNA Binding affinity predictor), for predicting the affinity of protein-DNA complexes and it is freely available at https://web.iitm.ac.in/bioinfo2/pdapred/.

PAM-Flexible Genome Editing with an Engineered Chimeric Cas9

CRISPR enzymes require a defined protospacer adjacent motif (PAM) flanking a guide RNA-programmed target site, limiting their sequence accessibility for robust genome editing applications. In this study, we recombine the PAM-interacting domain of SpRY, a broad-targeting Cas9 possessing an NRN > NYN PAM preference, with the N-terminus of Sc++, a Cas9 with simultaneously broad, efficient, and accurate NNG editing capabilities, to generate a chimeric enzyme with highly flexible PAM preference: SpRYc. We demonstrate that SpRYc leverages properties of both enzymes to specifically edit diverse NNN PAMs and disease-related loci for potential therapeutic applications. In total, the unique approaches to generate SpRYc, coupled with its robust flexibility, highlight the power of integrative protein design for Cas9 engineering and motivate downstream editing applications that require precise genomic positioning.

Structural and Energetic Features of Base-Base Stacking Contacts in RNA

Nucleobase π-π stacking is one of the crucial organizing interactions within three-dimensional (3D) RNA architectures. Characterizing the structural variability of these contacts in RNA crystal structures will help delineate their subtleties and their role in determining function. This analysis of different stacking geometries found in RNA X-ray crystal structures is the largest such survey to date; coupled with quantum-mechanical calculations on typical representatives of each possible stacking arrangement, we determined the distribution of stacking interaction energies. A total of 1,735,481 stacking contacts, spanning 359 of the 384 theoretically possible distinct stacking geometries, were identified. Our analysis reveals preferential occurrences of specific consecutive stacking arrangements in certain regions of RNA architectures. Quantum chemical calculations suggest that 88 of the 359 contacts possess intrinsically stable stacking geometries, whereas the remaining stacks require the RNA backbone or surrounding macromolecular environment to force their formation and maintain their stability. Our systematic analysis of π-π stacks in RNA highlights trends in the occurrence and localization of these noncovalent interactions and may help better understand the structural intricacies of functional RNA-based molecular architectures.

qNABpredict: Quick, accurate, and taxonomy-aware sequence-based prediction of content of nucleic acid binding amino acids

Protein sequence-based predictors of nucleic acid (NA)-binding include methods that predict NA-binding proteins and NA-binding residues. The residue-level tools produce more details but suffer high computational cost since they must predict every amino acid in the input sequence and rely on multiple sequence alignments. We propose an alternative approach that predicts content (fraction) of the NA-binding residues, offering more information than the protein-level prediction and much shorter runtime than the residue-level tools. Our first-of-its-kind content predictor, qNABpredict, relies on a small, rationally designed and fast-to-compute feature set that represents relevant characteristics extracted from the input sequence and a well-parametrized support vector regression model. We provide two versions of qNABpredict, a taxonomy-agnostic model that can be used for proteins of unknown taxonomic origin and more accurate taxonomy-aware models that are tailored to specific taxonomic kingdoms: archaea, bacteria, eukaryota, and viruses. Empirical tests on a low-similarity test dataset show that qNABpredict is 100 times faster and generates statistically more accurate content predictions when compared to the content extracted from results produced by the residue-level predictors. We also show that qNABpredict's content predictions can be used to improve results generated by the residue-level predictors. We release qNABpredict as a convenient webserver and source code at http://biomine.cs.vcu.edu/servers/qNABpredict/. This new tool should be particularly useful to predict details of protein-NA interactions for large protein families and proteomes.

DNA and RNA Binding Proteins: From Motifs to Roles in Cancer

DNA and RNA binding proteins (DRBPs) are a broad class of molecules that regulate numerous cellular processes across all living organisms, creating intricate dynamic multilevel networks to control nucleotide metabolism and gene expression. These interactions are highly regulated, and dysregulation contributes to the development of a variety of diseases, including cancer. An increasing number of proteins with DNA and/or RNA binding activities have been identified in recent years, and it is important to understand how their activities are related to the molecular mechanisms of cancer. In addition, many of these proteins have overlapping functions, and it is therefore essential to analyze not only the loss of function of individual factors, but also to group abnormalities into specific types of activities in regard to particular cancer types. In this review, we summarize the classes of DNA-binding, RNA-binding, and DRBPs, drawing particular attention to the similarities and differences between these protein classes. We also perform a cross-search analysis of relevant protein databases, together with our own pipeline, to identify DRBPs involved in cancer. We discuss the most common DRBPs and how they are related to specific cancers, reviewing their biochemical, molecular biological, and cellular properties to highlight their functions and potential as targets for treatment.

The 5′-phosphate enhances the DNA-binding and exonuclease activities of human mitochondrial genome maintenance nuclease 1 (MGME1)

In higher eukaryotes, mitochondria play multiple roles in energy production, signaling, and biosynthesis. Mitochondria possess multiple copies of mitochondrial DNA (mtDNA), which encodes 37 genes that are essential for mitochondrial and cellular function. When mtDNA is challenged by endogenous and exogenous factors, mtDNA undergoes repair, degradation, and compensatory synthesis. mtDNA degradation is an emerging pathway in mtDNA damage response and maintenance. A key factor involved is the human mitochondrial genome maintenance exonuclease 1 (MGME1). Despite previous biochemical and functional studies, controversies exist regarding the polarity of MGME1-mediated DNA cleavage. Also, how DNA sequence may affect the activities of MGME1 remains elusive. Such information is not only fundamental to the understanding of MGME1 but critical for deciphering the mechanism of mtDNA degradation. Herein, we use quantitative assays to examine the effects of substrate structure and sequence on the DNA-binding and enzymatic activities of MGME1. We demonstrate that MGME1 binds to and cleaves from the 5′-end of single-stranded DNA substrates, especially in the presence of 5′-phosphate, which plays an important role in DNA binding and optimal cleavage by MGME1. In addition, MGME1 tolerates certain modifications at the terminal end, such as a 5′-deoxyribosephosphate intermediate formed in base excision repair. We show that MGME1 processes different sequences with varying efficiencies, with dT and dC sequences being the most and least efficiently digested, respectively. Our results provide insights into the enzymatic properties of MGME1 and a rationale for the coordination of MGME1 with the 3′–5′ exonuclease activity of DNA polymerase γ in mtDNA degradation.

Decoding the Mechanism of Specific RNA Targeting by Ribosomal Methyltransferases

Methylation of specific nucleotides is integral for ribosomal biogenesis and also serves as a common mechanism to confer antibiotic resistance by pathogenic bacteria. Here, by determining the high-resolution structure of the 30S-KsgA complex by cryo-electron microscopy, a state was captured, where KsgA juxtaposes between helices h44 and h45 of the 30S ribosome, separating them, thereby enabling remodeling of the surrounded rRNA and allowing the cognate site to enter the methylation pocket. With the structure as a guide, several mutant versions of the ribosomes, where interacting bases in the catalytic helix h45 and surrounding helices h44, h24, and h27, were mutated and evaluated for their methylation efficiency revealing factors that direct the enzyme to its cognate site with high fidelity. The biochemical studies show that the three-dimensional environment of the ribosome enables the interaction of select loop regions in KsgA with the ribosome helices paramount to maintain selectivity.

Transcriptomic and Drug Discovery Analyses Reveal Natural Compounds Targeting the KDM4 Subfamily as Promising Adjuvant Treatments in Cancer

KDM4 proteins are a subfamily of histone demethylases that target the trimethylation of lysines 9 and 36 of histone H3, which are associated with transcriptional repression and elongation respectively. Their deregulation in cancer may lead to chromatin structure alteration and transcriptional defects that could promote malignancy. Despite that KDM4 proteins are promising drug targets in cancer therapy, only a few drugs have been described as inhibitors of these enzymes, while studies on natural compounds as possible inhibitors are still needed. Natural compounds are a major source of biologically active substances and many are known to target epigenetic processes such as DNA methylation and histone deacetylation, making them a rich source for the discovery of new histone demethylase inhibitors. Here, using transcriptomic analyses we determined that the KDM4 family is deregulated and associated with a poor prognosis in multiple neoplastic tissues. Also, by molecular docking and molecular dynamics approaches, we screened the COCONUT database to search for inhibitors of natural origin compared to FDA-approved drugs and DrugBank databases. We found that molecules from natural products presented the best scores in the FRED docking analysis. Molecules with sugars, aromatic rings, and the presence of OH or O- groups favor the interaction with the active site of KDM4 subfamily proteins. Finally, we integrated a protein-protein interaction network to correlate data from transcriptomic analysis and docking screenings to propose FDA-approved drugs that could be used as multitarget therapies or in combination with the potential natural inhibitors of KDM4 enzymes. This study highlights the relevance of the KDM4 family in cancer and proposes natural compounds that could be used as potential therapies.

Structural basis of R-loop recognition by the S9.6 monoclonal antibody

R-loops are ubiquitous, dynamic nucleic-acid structures that play fundamental roles in DNA replication and repair, chromatin and transcription regulation, as well as telomere maintenance. The DNA-RNA hybrid–specific S9.6 monoclonal antibody is widely used to map R-loops. Here, we report crystal structures of a S9.6 antigen-binding fragment (Fab) free and bound to a 13-bp hybrid duplex. We demonstrate that S9.6 exhibits robust selectivity in binding hybrids over double-stranded (ds) RNA and in categorically rejecting dsDNA. S9.6 asymmetrically recognizes a compact epitope of two consecutive RNA nucleotides via their 2′-hydroxyl groups and six consecutive DNA nucleotides via their backbone phosphate and deoxyribose groups. Recognition is mediated principally by aromatic and basic residues of the S9.6 heavy chain, which closely track the curvature of the hybrid minor groove. These findings reveal the molecular basis for S9.6 recognition of R-loops, detail its binding specificity, identify a new hybrid-recognition strategy, and provide a framework for S9.6 protein engineering.

The S9.6 monoclonal antibody is widely used to map R-loops genome wide. Here, Bou-Nader et al., define the nucleic acid-binding specificity of S9.6 and report its crystal structures free and bound to a hybrid, which reveal the asymmetric recognition of the RNA and DNA strands and its A-form conformation.

Quantum biology. π-π entanglement signatures in Protein-DNA interactions

DNA biological functions are carried out by individual proteins that interact with specific sequences along DNA to prime molecular processes required by cellular metabolism. Protein-DNA interactions include DNA replication, gene expression and its regulation, DNA repair, DNA restriction and modification by endonucleases, generally classified as enzymatic functions, or transcription factors functions. To find specific binding target sequences and finalize their activities, proteins must operate in symbiosis with cellular crowded environment identifying extremely small cognate sequences along the DNA chain, ranging from 15-20 bps for repressors to 4-6 bps for restriction enzymes in less than one second.

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.

NMR-Based Studies on Odorant-Melanoidin Interactions in Coffee Beverages

A quantitative ¹H NMR-based approach was established, which allowed the direct and noninvasive analysis of molecular interactions between key coffee odorants and high-molecular-weight (HMW) melanoidin polymers. A clear distinction between covalent and noncovalent interactions was achieved by monitoring the time dependency of odorant-polymer interactions, resulting in four scenarios: covalent, π-π, covalent and π-π-, as well as no interactions. Evaluation of temperature influence on e.g. 2-furfurylthiol (FFT), revealed an altered behavior with increased π-π stacking at lower temperatures and accelerated covalent interactions at higher temperatures. Human sensory experiments with HMW material and a coffee aroma reconstitution model showed a drastic reduction of "roasty/sulfury" aroma notes, as well as an increased "sweetish/caramel-like" flavor. The lack of interactions between the "sweetish/caramel" smelling 4-hydroxy-2,5-dimethyl-3(2H)-furanone with the HMW melanoidins in combination with the high binding affinity of coffee thiols explains the sensory evaluation and is obviously the reason for the fast disappearance of the typical "roasty/sulfury" aroma impressions of a freshly prepared coffee brew.

Exploring the Nature of Hydrogen Bonding between RNA and Proteins: A Comprehensive Analysis of RNA : Protein Complexes

A nonredundant dataset of ∼300 high (up to 2.5 Å) resolution X-ray structures of RNA:protein complexes were analyzed for hydrogen bonds between amino-acid residues and canonical ribonucleotides (rNs). The identified 17100 contacts were classified based on the identity (rA, rC, rG or rU) and interacting fragment (base, sugar, or ribose) of the rN, the nature (polar or nonpolar) and interacting moiety (main chain or side chain) of the amino-acid residue, as well as the rN and amino-acid atoms participating in the hydrogen bonding. 80 possible hydrogen-bonding combinations (4 (rNs) X 20 (amino acids)) involve a wide variety of RNA and protein types and are present in multiple occurrences in almost all PDB files. Comparison with the analogously-selected DNA:protein complexes reveals that the absence of 2'-OH group in DNA mainly accounts for the differences in DNA:protein and RNA:protein hydrogen bonding. Search for intrinsically-stable base:amino acid pairs containing single or multiple hydrogen bonds reveals 37 unique pairs, which may act as well-defined RNA:protein interaction motifs. Overall, our work collectively analyzes the largest set of nucleic acid-protein hydrogen bonds to date, and therefore highlights several trends that may help frame structural rules governing the physiochemical characteristics of RNA:protein recognition.

Amyloid-like Self-Assembly of a Hydrophobic Cell-Penetrating Peptide and Its Use as a Carrier for Nucleic Acids

Cell-penetrating peptides (CPPs) are a topical subject potentially exploitable for creating nanotherapeutics for the delivery of bioactive loads. These compounds are often classified into three major categories according to their physicochemical characteristics: cationic, amphiphilic, and hydrophobic. Among them, the group of hydrophobic CPPs has received increasing attention in recent years due to toxicity concerns posed by highly cationic CPPs. The hexapeptide PFVYLI (P, proline; F, phenylalanine; V, valine; Y, tyrosine; L, leucine; and I, isoleucine), a fragment derived from the C-terminal portion of α1-antitrypsin, is a prototypal example of hydrophobic CPP. This sequence shows reduced cytotoxicity and a capacity of nuclear localization, and its small size readily hints at its suitability as a building block to construct nanostructured materials. In this study, we examine the self-assembling properties of PFVYLI and investigate its ability to form noncovalent complexes with nucleic acids. By using a combination of biophysical tools including synchrotron small-angle X-ray scattering and atomic force microscopy-based infrared spectroscopy, we discovered that this CPP self-assembles into discrete nanofibrils with remarkable amyloidogenic features. Over the course of days, these fibrils coalesce into rodlike crystals that easily reach the micrometer range. Despite lacking cationic residues in the composition, PFVYLI forms noncovalent complexes with nucleic acids that retain β-sheet pairing found in amyloid aggregates. In vitro vectorization experiments performed with double-stranded DNA fragments indicate that complexes promote the internalization of nucleic acids, revealing that tropism toward cell membranes is preserved upon complexation. On the other hand, transfection assays with splice-correction oligonucleotides (SCOs) for luciferase expression show limited bioactivity across a narrow concentration window, suggesting that the propensity to form amyloidogenic aggregates may trigger endosomal entrapment. We anticipate that the findings presented here open perspectives for using this archetypical hydrophobic CPP in the fabrication of nanostructured scaffolds, which potentially integrate properties of amyloids and translocation capabilities of CPPs.

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.

A comparative study of protein-ssDNA interactions

Single-stranded DNA-binding proteins (SSBs) play crucial roles in DNA replication, recombination and repair, and serve as key players in the maintenance of genomic stability. While a number of SSBs bind single-stranded DNA (ssDNA) non-specifically, the others recognize and bind specific ssDNA sequences. The mechanisms underlying this binding discrepancy, however, are largely unknown. Here, we present a comparative study of protein-ssDNA interactions by annotating specific and non-specific SSBs and comparing structural features such as DNA-binding propensities and secondary structure types of residues in SSB-ssDNA interactions, protein-ssDNA hydrogen bonding and π-π interactions between specific and non-specific SSBs. Our results suggest that protein side chain-DNA base hydrogen bonds are the major contributors to protein-ssDNA binding specificity, while π-π interactions may mainly contribute to binding affinity. We also found the enrichment of aspartate in the specific SSBs, a key feature in specific protein-double-stranded DNA (dsDNA) interactions as reported in our previous study. In addition, no significant differences between specific and non-specific groups with respect of conformational changes upon ssDNA binding were found, suggesting that the flexibility of SSBs plays a lesser role than that of dsDNA-binding proteins in conferring binding specificity.

A more efficient CRISPR-Cas12a variant derived from Lachnospiraceae bacterium MA2020

CRISPR effector proteins introduce double-stranded breaks into the mammalian genome, facilitating gene editing by non-homologous end-joining or homology-directed repair. Unlike the more commonly studied Cas9, the CRISPR effector protein Cas12a/Cpf1 recognizes a T-rich protospacer adjacent motif (PAM) and can process its own CRISPR RNA (crRNA) array, simplifying the use of multiple guide RNAs. We observed that the Cas12a ortholog of Lachnospiraceae bacterium MA2020 (Lb2Cas12a) edited mammalian genes with efficiencies comparable to those of AsCas12a and LbCas12a. Compared to these well-characterized Cas12a orthologs, Lb2Cas12a is smaller and recognizes a narrow set of PAM TTTV. We introduced two mutations into Lb2Cas12a, Q571K and C1003Y, that increased its cleavage efficiency for a range of target sequences beyond those of the commonly used Cas12a orthologs AsCas12a and LbCas12a. In addition to the canonical TTTV PAM, this variant, Lb2-KY, also efficiently cleaved target regions with CTTN PAMs. Finally, we demonstrated that Lb2-KY ribonucleoprotein (RNP) complexes edited two hemoglobin target regions useful for correcting common forms of sickle-cell anemia more efficiently than commercial AsCas12a RNP complexes. Thus, Lb2-KY has distinctive properties useful for modifying a range of clinically relevant targets in the human genome.

Investigating the lipid fingerprint of SLC6 neurotransmitter transporters: a comparison of dDAT, hDAT, hSERT, and GlyT2

The local lipid annulus, or "fingerprint", of four SLC6 transporters (dDAT, hDAT, hSERT, and GlyT2) embedded in a complex neuronal membrane were compared and characterised using molecular dynamics. Our analysis included the development of new tools to improve membrane leaflet detection and the analysis of leaflet-dependent properties. Overall, the lipid fingerprints of the four transporters are comprised of similar lipids when grouped by headgroup or tail saturation. The enrichment and depletion of specific lipids, including sites of cholesterol contacts, varies between transporters. The subtle differences in lipid fingerprints results in varying membrane biophysical properties near the transporter. Our results highlight that the lipid-fingerprint of SLC6 transporters in complex membranes is highly dependent on membrane composition. Our results further characterize how the presence and identity of membrane proteins affects the complex interplay of lipid-protein interactions, influencing the local lipid environment and membrane biophysical properties.

Protein-DNA target search relies on quantum walk

Protein-DNA interactions play a fundamental role in all life systems. A critical issue of such interactions is given by the strategy of protein search for specific targets on DNA. The mechanisms by which the protein are able to find relatively small cognate sequences, typically 15-20 base pairs (bps) for repressors, and 4-6 bps for restriction enzymes among the millions of bp of non-specific chromosomal DNA have hardly engaged researchers for decades. Recent experimental studies have generated new insights on the basic processes of protein-DNA interactions evidencing the underlying complex dynamic phenomena involved, which combine three-dimensional and one-dimensional motion along the DNA chain. It has been demonstrated that protein molecules have an extraordinary ability to find the target very quickly on the DNA chain, in some cases, with two orders of magnitude faster than the diffusion limit. This unique property of protein-DNA search mechanism is known as facilitated diffusion. Several theoretical mechanisms have been suggested to describe the origin of facilitated diffusion. However, none of such models currently has the ability to fully describe the protein search strategy. In this paper, we suggest that the ability of proteins to identify consensus sequences on DNA is based on the entanglement of π-π electrons between DNA nucleotides and protein amino acids. The π-π entanglement is based on Quantum Walk (QW), through Coin-position entanglement (CPE). First, the protein identifies a dimer belonging to the consensus sequence, and localize a π on such dimer, hence, the other π electron scans the DNA chain until the sequence is identified. Focusing on the example of recognition of consensus sequences of EcoRV or EcoRI, we will describe the quantum features of QW on protein-DNA complexes during the search strategy, such as walker quadratic spreading on a coherent superposition of different vertices and environment-supported long-time survival probability of the walker. We will employ both discrete- or continuous-time versions of QW. Biased and unbiased classical Random Walk (CRW) have been used for a long time to describe the Protein-DNA search strategy. QW, the quantum version of CRW, has been widely studied for its applications in quantum information applications. In our biological application, the walker (the protein) resides at a vertex in a graph (the DNA structural topology). Differently to CRW, where the walker moves randomly, the quantum walker can hop along the edges in the graph to reach other vertices entering coherently a superposition across different vertices spreading quadratically faster than CRW analogous evidencing the typical speed up features of the QW. When applied to a protein-DNA target search problem, QW gives the possibility to achieve the experimental diffusional motion of proteins over diffusion classical limits experienced along DNA chains exploiting quantum features such as CPE and long-time survival probability supported by the environment. In turn, we come to the conclusion that, under quantum picture, the protein search strategy does not distinguish between one-dimensional (1D) and three-dimensional (3D) cases.

Protein Side-Chain-DNA Contacts Probed by Fast Magic-Angle Spinning NMR

Protein–nucleic acid interactions are essential in a variety of biological events ranging from the replication of genomic DNA to the synthesis of proteins. Noncovalent interactions guide such molecular recognition events, and protons are often at the center of them, particularly due to their capability of forming hydrogen bonds to the nucleic acid phosphate groups. Fast magic-angle spinning experiments (100 kHz) reduce the proton NMR line width in solid-state NMR of fully protonated protein–DNA complexes to such an extent that resolved proton signals from side-chains coordinating the DNA can be detected. We describe a set of NMR experiments focusing on the detection of protein side-chains from lysine, arginine, and aromatic amino acids and discuss the conclusions that can be obtained on their role in DNA coordination. We studied the 39 kDa enzyme of the archaeal pRN1 primase complexed with DNA and characterize protein–DNA contacts in the presence and absence of bound ATP molecules.

Evolution of DNA replication origin specification and gene silencing mechanisms

DNA replication in eukaryotic cells initiates from replication origins that bind the Origin Recognition Complex (ORC). Origin establishment requires well-defined DNA sequence motifs in Saccharomyces cerevisiae and some other budding yeasts, but most eukaryotes lack sequence-specific origins. A 3.9 Å structure of S. cerevisiae ORC-Cdc6-Cdt1-Mcm2-7 (OCCM) bound to origin DNA revealed that a loop within Orc2 inserts into a DNA minor groove and an α-helix within Orc4 inserts into a DNA major groove. Using a massively parallel origin selection assay coupled with a custom mutual-information-based modeling approach, and a separate analysis of whole-genome replication profiling, here we show that the Orc4 α-helix contributes to the DNA sequence-specificity of origins in S. cerevisiae and Orc4 α-helix mutations change genome-wide origin firing patterns. The DNA sequence specificity of replication origins, mediated by the Orc4 α-helix, has co-evolved with the gain of ORC-Sir4-mediated gene silencing and the loss of RNA interference.

Contrary to most eukaryotes that lack sequence-specific origins of replication, S. cerevisiae origins are defined by specific DNA sequence motifs. Here the authors reveal that multiple subunits of ORC, including Orc2 and Orc4, contribute to the sequence-specificity of origins in S. cerevisiae.

Insights into the Direct Oxidative Repair of Etheno Lesions: MD and QM/MM Study on the Substrate Scope of ALKBH2 and AlkB

E. coli AlkB and human ALKBH2 belong to the AlkB family enzymes, which contain several α-ketoglutarate (α-KG)/Fe(II)-dependent dioxygenases that repair alkylated DNA. Specifically, the AlkB enzymes catalyze decarboxylation of α-KG to generate a high-valent Fe(IV)-oxo species that oxidizes alkyl groups on DNA adducts. AlkB and ALKBH2 have been reported to differentially repair select etheno adducts, with preferences for 1,N⁶-ethenoadenine (1,N⁶-εA) and 3,N⁴-ethenocytosine (3,N⁴-εC) over 1,N²-ethenoguanine (1,N²-εG). However, N²,3-ethenoguanine (N²,3-εG), the most common etheno adduct, is not repaired by the AlkB enzymes. Unfortunately, a structural understanding of the differential activity of E. coli AlkB and human ALKBH2 is lacking due to challenges acquiring atomistic details for a range of substrates using experiments. This study uses both molecular dynamics (MD) simulations and ONIOM(QM:MM) calculations to determine how the active site changes upon binding each etheno adduct and characterizes the corresponding catalytic impacts. Our data reveal that the preferred etheno substrates (1,N⁶-εA and 3,N⁴-εC) form favorable interactions with catalytic residues that situate the lesion near the Fe(IV)-oxo species and permit efficient oxidation. In contrast, although the damage remains correctly aligned with respect to the Fe(IV)-oxo moiety, repair of 1,N²-εG is mitigated by increased solvation of the active site and a larger distance between Fe(IV)-oxo and the aberrant carbons. Binding of non-substrate N²,3-εG in the active site disrupts key DNA-enzyme interactions, and positions the aberrant carbon atoms even further from the Fe(IV)-oxo species, leading to prohibitively high barriers for oxidative catalysis. Overall, our calculations provide the first structural insight required to rationalize the experimentally-reported substrate specificities of AlkB and ALKBH2 and thereby highlight the roles of several active site residues in the repair of etheno adducts that directly correlates with available experimental data. These proposed catalytic strategies can likely be generalized to other α-KG/Fe(II)-dependent dioxygenases that play similar critical biological roles, including epigenetic and post-translational regulation.

Ratiometric Detection of ATP by Fluorescent Cyclophanes with Bellows-Type Sensing Mechanism

Pyrene-based cyclophanes have been synthesized with the aim to realize a bellows-type sensing mechanism for the ratiometric detection of nucleotide concentrations in a buffered aqueous solution. The sensing mechanism involves the encapsulation of a nucleobase between two pyrene rings, which affects the monomer-excimer equilibrium of the receptor in the excited state. The nature of the spacer and its connection pattern to pyrene rings have been varied to achieve high selectivity for ATP. The 1,8-substituted pyrene-based cyclophane with the 2,2'-diaminodiethylamine spacer demonstrates the best selectivity for ATP showing a 50-fold increase in the monomer-excimer emission ratio upon saturation with the nucleotide. The receptor can detect ATP within the biological concentrations range over a wide pH range. NMR and spectroscopic studies have revealed the importance of hydrogen bonding and stacking interactions for achieving a required receptor selectivity. The probe has been successfully applied for the real-time monitoring of creatine kinase activity.

PlantPAN3.0: a new and updated resource for reconstructing transcriptional regulatory networks from ChIP-seq experiments in plants

The Plant Promoter Analysis Navigator (PlantPAN; http://PlantPAN.itps.ncku.edu.tw/) is an effective resource for predicting regulatory elements and reconstructing transcriptional regulatory networks for plant genes. In this release (PlantPAN 3.0), 17 230 TFs were collected from 78 plant species. To explore regulatory landscapes, genomic locations of TFBSs have been captured from 662 public ChIP-seq samples using standard data processing. A total of 1 233 999 regulatory linkages were identified from 99 regulatory factors (TFs, histones and other DNA-binding proteins) and their target genes across seven species. Additionally, this new version added 2449 matrices extracted from ChIP-seq peaks for cis-regulatory element prediction. In addition to integrated ChIP-seq data, four major improvements were provided for more comprehensive information of TF binding events, including (i) 1107 experimentally verified TF matrices from the literature, (ii) gene regulation network comparison between two species, (iii) 3D structures of TFs and TF-DNA complexes and (iv) condition-specific co-expression networks of TFs and their target genes extended to four species. The PlantPAN 3.0 can not only be efficiently used to investigate critical cis- and trans-regulatory elements in plant promoters, but also to reconstruct high-confidence relationships among TF–targets under specific conditions.

New insights into protein-DNA binding specificity from hydrogen bond based comparative study

Knowledge of protein-DNA binding specificity has important implications in understanding DNA metabolism, transcriptional regulation and developing therapeutic drugs. Previous studies demonstrated hydrogen bonds between amino acid side chains and DNA bases play major roles in specific protein-DNA interactions. In this paper, we investigated the roles of individual DNA strands and protein secondary structure types in specific protein-DNA recognition based on side chain-base hydrogen bonds. By comparing the contribution of each DNA strand to the overall binding specificity between DNA-binding proteins with different degrees of binding specificity, we found that highly specific DNA-binding proteins show balanced hydrogen bonding with each of the two DNA strands while multi-specific DNA binding proteins are generally biased towards one strand. Protein-base pair hydrogen bonds, in which both bases of a base pair are involved in forming hydrogen bonds with amino acid side chains, are more prevalent in the highly specific protein-DNA complexes than those in the multi-specific group. Amino acids involved in side chain-base hydrogen bonds favor strand and coil secondary structure types in highly specific DNA-binding proteins while multi-specific DNA-binding proteins prefer helices.

Interaction of enzymes with lignocellulosic materials: causes, mechanism and influencing factors

For the production of biofuel (bioethanol), enzymatic adsorption onto a lignocellulosic biomass surface is a prior condition for the enzymatic hydrolysis process to occur. Lignocellulosic substances are mainly composed of cellulose, hemicellulose and lignin. The polysaccharide matrix (cellulose and hemicellulose) is capable of producing bioethanol. Therefore, lignin is removed or its concentration is reduced from the adsorption substrates by pretreatments. Selected enzymes are used for the production of reducing sugars from cellulosic materials, which in turn are converted to bioethanol. Adsorption of enzymes onto the substrate surface is a complicated process. A large number of research have been performed on the adsorption process, but little has been done to understand the mechanism of adsorption process. This article reviews the mechanisms of adsorption of enzymes onto the biomass surfaces. A conceptual adsorption mechanism is presented which will fill the gaps in literature and help researchers and industry to use adsorption more efficiently. The process of enzymatic adsorption starts with the reciprocal interplay of enzymes and substrates and ends with the establishment of molecular and cellular binding. The kinetics of an enzymatic reaction is almost the same as that of a characteristic chemical catalytic reaction. The influencing factors discussed in detail are: surface characteristics of the participating materials, the environmental factors, such as the associated flow conditions, temperature, concentration, etc. Pretreatment of lignocellulosic materials and optimum range of shear force and temperature for getting better results of adsorption are recommended.

The physical basis and practical consequences of biological promiscuity

Proteins interact with metabolites, nucleic acids, and other proteins to orchestrate the myriad catalytic, structural and regulatory functions that support life from the simplest microbes to the most complex multicellular organisms. These molecular interactions are often exquisitely specific, but never perfectly so. Adventitious "promiscuous" interactions are ubiquitous due to the thousands of macromolecules and small molecules crowded together in cells. Such interactions may perturb protein function at the molecular level, but as long as they do not compromise organismal fitness, they will not be removed by natural selection. Although promiscuous interactions are physiologically irrelevant, they are important because they can provide a vast reservoir of potential functions that can provide the starting point for evolution of new functions, both in nature and in the laboratory.

Structural Patterns and Stabilities of Hydrogen-Bonded Pairs Involving Ribonucleotide Bases and Arginine, Glutamic Acid, or Glutamine Residues of Proteins from Quantum Mechanical Calculations

Ribonucleotide:protein interactions play crucial roles in a number of biological processes. Unlike the RNA:protein interface where van der Waals contacts are prevalent, the recognition of a single ribonucleotide such as ATP by a protein occurs predominantly through hydrogen-bonding interactions. As a first step toward understanding the role of hydrogen bonding in ribonucleotide:protein recognition, the present work employs density functional theory to provide a detailed quantum-mechanical analysis of the structural and energetic characteristics of 18 unique hydrogen-bonded pairs involving the nucleobase/nucleoside moiety of four canonical ribonucleotides and the side chains of three polar amino-acid residues (arginine, glutamine, and glutamic acid) of proteins. In addition, we model five new pairs that are till now not observed in crystallographically identified ribonucleotide:protein complexes but may be identified in complexes crystallized in the future. We critically examine the characteristics of each pair in its ribonucleotide:protein crystal structure occurrence and (gas phase and water phase) optimized intrinsic structure. We further evaluated the interaction energy of each pair and characterized the associated hydrogen bonds using a number of quantum mechanics-based relationships including natural bond orbital analysis, quantum theory atoms in molecules analysis, Iogansen relationships, Nikolaienko-Bulavin-Hovorun relationships, and noncovalent interaction-reduced density gradient analysis. Our analyses reveal rich variability in hydrogen bonds in the crystallographic as well as intrinsic structure of each pair, which includes conventional O/N-H···N/O and C-H···O hydrogen bonds as well as donor/acceptor-bifurcated hydrogen bonds. Further, we identify five combinations of nucleobase and amino acid moieties; each of which exhibits at least two alternate (i.e., multimodal) structures that interact through the same nucleobase edge. In fact, one such pair exhibits four multimodal structures; one of which possesses unconventional "amino-acceptor" hydrogen bonding with comparable (-9.4 kcal mol^-1) strength to the corresponding conventional (i.e., amino:donor) structure (-9.2 kcal mol^-1). This points to the importance of amino-acceptor hydrogen bonds in RNA:protein interactions and suggests that such interactions must be considered in the future while studying the dynamics in the context of molecular recognition. Overall, our study provides preliminary insights into the intrinsic features of ribonucleotide:amino acid interactions, which may help frame a clearer picture of the molecular basis of RNA:protein recognition and further appreciate the role of such contacts in biology.

Nucleotide Binding Modes in a Motor Protein Revealed by 31 P- and 1 H-Detected MAS Solid-State NMR Spectroscopy

Protein-nucleic acid interactions play important roles not only in energy-providing reactions, such as ATP hydrolysis, but also in reading, extending, packaging, or repairing genomes. Although they can often be analyzed in detail with X-ray crystallography, complementary methods are needed to visualize them in complexes, which are not crystalline. Here, we show how solid-state NMR spectroscopy can detect and classify protein-nucleic interactions through site-specific 1 H- and 31 P-detected spectroscopic methods. The sensitivity of 1 H chemical-shift values on noncovalent interactions involved in these molecular recognition processes is exploited allowing us to probe directly the chemical bonding state, an information, which is not directly accessible from an X-ray structure. We show that these methods can characterize interactions in easy-to-prepare sediments of the 708 kDa dodecameric DnaB helicase in complex with ADP:AlF4- :DNA, and this despite the very challenging size of the complex.

Revisiting the Potential Energy Surface of the Stacked Cytosine Dimer: FNO-CCSD(T) Interaction Energies, SAPT Decompositions, and Benchmarking

Nucleobase stacking interactions are crucial for the stability of nucleic acids. This study investigates base stacking energies of the cytosine homodimer in different configurations, including intermolecular separation plots, detailed twist dependence, and displaced structures. Highly accurate ab initio quantum chemical single point energies using an energy function based on MP2 complete basis set extrapolation ([6 → 7]ZaPa-NR) and a CCSD(T)/cc-pVTZ-F12 high-level correction are presented as new reference data, providing the most accurate stacking energies of nucleobase dimers currently available. Accurate SAPT2+(3)δMP2 energy decomposition is used to obtain detailed insights into the nature of base stacking interactions at varying vertical distances and twist values. The ab initio symmetry adapted perturbation theory (SAPT) energy decomposition suggests that the base stacking originates from an intricate interplay between dispersion attraction, short-range exchange-repulsion, and Coulomb interaction. The interpretation of the SAPT data is a complex issue as key energy terms vary substantially in the region of optimal (low energy) base stacking geometries. Thus, attempts to highlight one leading stabilizing SAPT base stacking term may be misleading and the outcome strongly depends on the used geometries within the range of geometries sampled in nucleic acids upon thermal fluctuations. Modern dispersion-corrected density functional theory (among them DSD-BLYP-D3, ωB97M-V, and ωB97M-D3BJ) is benchmarked and often reaches up to spectroscopic accuracy (below 1 kJ/mol). The classical AMBER force field is benchmarked with multiple different sets of point-charges (e.g. HF, DFT, and MP2-based) and is found to produce reasonable agreement with the benchmark data.