Role of aromatic amino acids in protein-nucleic acid recognition

Isomorphic Fluorescent Nucleosides

In 1960, Weber prophesied that "There are many ways in which the properties of the excited state can be utilized to study points of ignorance of the structure and function of proteins". This has been realized, illustrating that an intrinsic and highly responsive fluorophore such as tryptophan can alter the course of an entire scientific discipline. But what about RNA and DNA? Adapting Weber's protein photophysics prophecy to nucleic acids requires the development of intrinsically emissive nucleoside surrogates as, unlike Trp, the canonical nucleobases display unusually low emission quantum yields, which render nucleosides, nucleotides, and oligonucleotides practically dark for most fluorescence-based applications.Over the past decades, we have developed emissive nucleoside surrogates that facilitate the monitoring of nucleoside-, nucleotide-, and nucleic acid-based transformations at a nucleobase resolution in real time. The premise underlying our approach is the identification of minimal atomic/structural perturbations that endow the synthetic analogs with favorable photophysical features while maintaining native conformations and pairing. As illuminating probes, the photophysical parameters of such isomorphic nucleosides display sensitivity to microenvironmental factors. Responsive isomorphic analogs that function similarly to their native counterparts in biochemical contexts are defined as isofunctional.Early analogs included pyrimidines substituted with five-membered aromatic heterocycles at their 5 position and have been used to assess the polarity of the major groove in duplexes. Polarized quinazolines have proven useful in assembling FRET pairs with established fluorophores and have been used to study RNA-protein and RNA-small-molecule binding. Completing a fluorescent ribonucleoside alphabet, composed of visibly emissive purine (thA, thG) and pyrimidine (thU, thC) analogs, all derived from thieno[3,4-d]pyrimidine as the heterocyclic nucleus, was a major breakthrough. To further augment functionality, a second-generation emissive RNA alphabet based on an isothiazolo[4,3-d]pyrimidine core (thA, tzG, tzU, and tzC) was fabricated. This single-atom "mutagenesis" restored the basic/coordinating nitrogen corresponding to N7 in the purine skeleton and elevated biological recognition.The isomorphic emissive nucleosides and nucleotides, particularly the purine analogs, serve as substrates for diverse enzymes. Beyond polymerases, we have challenged the emissive analogs with metabolic and catabolic enzymes, opening optical windows into the biochemistry of nucleosides and nucleotides as metabolites as well as coenzymes and second messengers. Real-time fluorescence-based assays for adenosine deaminase, guanine deaminase, and cytidine deaminase have been fabricated and used for inhibitor discovery. Emissive cofactors (e.g., SthAM), coenzymes (e.g., NtzAD+), and second messengers (e.g., c-di-tzGMP) have been enzymatically synthesized, using xyNTPs and native enzymes. Both their biosynthesis and their transformations can be fluorescently monitored in real time.Highly isomorphic and isofunctional emissive surrogates can therefore be fabricated and judiciously implemented. Beyond their utility, side-by-side comparison to established analogs, particularly to 2-aminopurine, the workhorse of nucleic acid biophysics over 5 decades, has proven prudent as they refined the scope and limitations of both the new analogs and their predecessors. Challenges, however, remain. Associated with such small heterocycles are relatively short emission wavelengths and limited brightness. Recent advances in multiphoton spectroscopy and further structural modifications have shown promise for overcoming such barriers.

DNA-Binding Activity of CAMTA3 Is Essential for Its Function: Identification of Critical Amino Acids for Its Transcriptional Activity

Calmodulin-binding transcription activators (CAMTAs), a small family of highly conserved transcription factors, function in calcium-mediated signaling pathways. Of the six CAMTAs in Arabidopsis, CAMTA3 regulates diverse biotic and abiotic stress responses. A recent study has shown that CAMTA3 is a guardee of NLRs (Nucleotide-binding, Leucine-rich repeat Receptors) in modulating plant immunity, raising the possibility that CAMTA3 transcriptional activity is dispensable for its function. Here, we show that the DNA-binding activity of CAMTA3 is essential for its role in mediating plant immune responses. Analysis of the DNA-binding (CG-1) domain of CAMTAs in plants and animals showed strong conservation of several amino acids. We mutated six conserved amino acids in the CG-1 domain to investigate their role in CAMTA3 function. Electrophoretic mobility shift assays using these mutants with a promoter of its target gene identified critical amino acid residues necessary for DNA-binding activity. In addition, transient assays showed that these residues are essential for the CAMTA3 function in activating the Rapid Stress Response Element (RSRE)-driven reporter gene expression. In line with this, transgenic lines expressing the CG-1 mutants of CAMTA3 in the camta3 mutant failed to rescue the mutant phenotype and restore the expression of CAMTA3 downstream target genes. Collectively, our results provide biochemical and genetic evidence that the transcriptional activity of CAMTA3 is indispensable for its function.

Structural Protein Effects Underpinning Cognitive Developmental Delay of the PURA p.Phe233del Mutation Modelled by Artificial Intelligence and the Hybrid Quantum Mechanics–Molecular Mechanics Framework

A whole-exome capture and next-generation sequencing was applied to an 11 y/o patient with a clinical history of congenital hypotonia, generalized motor and cognitive neurodevelopmental delay, and severe cognitive deficit, and without any identifiable Syndromic pattern, and to her parents, we disclosed a de novo heterozygous pathogenic mutation, c.697_699del p.Phe233del (rs786204835)(ACMG classification PS2, PM1, PM2, PP5), harbored in the PURA gene (MIM*600473) (5q31.3), associated with Autosomal Dominant Mental Retardation 31 (MIM # 616158). We used the significant improvement in the accuracy of protein structure prediction recently implemented in AlphaFold that incorporates novel neural network architectures and training procedures based on the evolutionary, physical, and geometric constraints of protein structures. The wild-type (WT) sequence and the mutated sequence, missing the Phe233, were reconstructed. The predicted local Distance Difference Test (lDDT) for the PURAwt and the PURA–Phe233del showed that the occurrence of the Phe233del affects between 220–320 amino acids. The distortion in the PURA structural conformation in the ~5 Å surrounding area after the p.Phe233del produces a conspicuous disruption of the repeat III, where the DNA and RNA helix unwinding capability occurs. PURA Protein–DNA docking corroborated these results in an in silico analysis that showed a loss of the contact of the PURA–Phe233del III repeat domain model with the DNA. Together, (i) the energetic and stereochemical, (ii) the hydropathic indexes and polarity surfaces, and (iii) the hybrid Quantum Mechanics–Molecular Mechanics (QM–MM) analyses of the PURA molecular models demarcate, at the atomic resolution, the specific surrounding region affected by these mutations and pave the way for future cell-based functional analysis. To the best of our knowledge, this is the first report of a de novo mutation underpinning a PURA syndrome in a Latin American patient and highlights the importance of predicting the molecular effects in protein structure using artificial intelligence algorithms and molecular and atomic resolution stereochemical analyses.

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.

Crystal structure and self-assembly on graphite of a pyrazolo[1,5-c]pyrimidine derivative

The crystal structure of the heterocyclic compound 2-(4-methoxyphenyl)-7-phenylpyrazolo[1,5-c]pyrimidine, C₁₉H₁₅N₃O, has been determined and its self-assembly on the surface of graphite has been examined using atomic force microscopy (AFM). The title compound crystallized in the monoclinic space group P2₁/c, with two independent molecules in the asymmetric unit. The packing of the L-shaped molecules in the crystal is governed by arene interactions, in the absence of any conventional hydrogen-bonding interactions. The packing arrangement reveals four types of dimeric motifs stabilized by π-π and C-H...π interactions. At low coverage, molecules assemble into long needle-like islands on the graphite surface. High-resolution AFM images reveal that the molecules interact through weak noncovalent interactions between the aromatic H atoms and the methoxy O atoms.

Correlation in Domain Fluctuations Navigates Target Search of a Viral Peptide along RNA

Biological macromolecules often exhibit correlations in fluctuations involving distinct domains. This study decodes their functional implications in RNA-protein recognition and target-specific binding. The target search of a peptide along RNA in a viral TAR-Tat complex is closely monitored using atomistic simulations, steered molecular dynamics simulations, free energy calculations, and a machine-learning-based clustering technique. An anticorrelated domain fluctuation is identified between the tetraloop and the bulge region in the apo form of TAR RNA that sets a hierarchy in the domain-specific fluctuations at each binding event and that directs the succeeding binding footsteps. Thus, at each binding footstep, the dynamic partner selects an RNA location for binding where it senses a higher fluctuation, which is conventionally reduced upon binding. This event stimulates an alternate domain fluctuation, which then dictates sequential binding footstep/s and thus the search progresses. Our cross-correlation maps show that the fluctuations relay from one domain to another specific domain until the anticorrelation between those interdomain fluctuations sustains. Artificial attenuation of that hierarchical domain fluctuation inhibits specific RNA binding. The binding is completed with the arrival of a few long-lived water molecules that mediate slightly distant RNA-protein sites and finally stabilize the overall complex. The study underscores the functional importance of naturally designed fluctuating RNA motifs (bulge, tetraloop) and their interplay in dictating the directionality of the search in a highly dynamic environment.

3D Interaction Homology: Hydropathic Analyses of the "π-Cation" and "π-π" Interaction Motifs in Phenylalanine, Tyrosine, and Tryptophan Residues

Three-dimensional (3D) maps of the hydropathic environments of protein amino acid residues are information-rich descriptors of preferred conformations, interaction types and energetics, and solvent accessibility. The interactions made by each residue are the primary factor for rotamer selection and the secondary, tertiary, and even quaternary protein structure. Our evolving basis set of environmental data for each residue type can be used to understand the protein structure. This work focuses on the aromatic residues phenylalanine, tyrosine, and tryptophan and their structural roles. We calculated and analyzed side chain-to-environment 3D maps for over 70,000 residues of these three types that reveal, with respect to hydrophobic and polar interactions, the environment around each. After binning with backbone ϕ/ψ and side chain χ₁, we clustered each bin by 3D similarities between map-map pairs. For each of the three residue types, four bins were examined in detail: one in the β-pleat, two in the right-hand α-helix, and one in the left-hand α-helix regions of the Ramachandran plot. For high degrees of side chain burial, encapsulation of the side chain by hydrophobic interactions is ubiquitous. The more solvent-exposed side chains are more likely to be involved in polar interactions with their environments. Evidence for π-π interactions was observed in about half of the residues surveyed [phenylalanine (PHE): 53.3%, tyrosine (TYR): 34.1%, and tryptophan (TRP): 55.7%], but on an energy basis, this contributed to only ∼4% of the total. Evidence for π-cation interactions was observed in 14.1% of PHE, 8.3% of TYR, and 26.8% of TRP residues, but on an energy basis, this contributed to only ∼1%. This recognition of even these subtle interactions in the 3D hydropathic environment maps is key support for our interaction homology paradigm of protein structure elucidation and possibly prediction.

ATP regulates RNA-driven cold inducible RNA binding protein phase separation

Intrinsically disordered proteins and proteins containing intrinsically disordered regions are highly abundant in the proteome of eukaryotes and are extensively involved in essential biological functions. More recently, their role in the organization of biomolecular condensates has become evident and along with their misregulation in several neurologic disorders. Currently, most studies involving these proteins are carried out in vitro and using purified proteins. Given that in cells, condensate‐forming proteins are exposed to high, millimolar concentrations of cellular metabolites, we aimed to reveal the interactions of cellular metabolites and a representative condensate‐forming protein. Here, using the arginine–glycine/arginine–glycine–glycine (RG/RGG)‐rich cold inducible RNA binding protein (CIRBP) as paradigm, we studied binding of the cellular metabolome to CIRBP. We found that most of the highly abundant cellular metabolites, except nucleotides, do not directly bind to CIRBP. ATP, ADP, and AMP as well as NAD⁺, NADH, NADP⁺, and NADPH directly interact with CIRBP, involving both the folded RNA‐recognition motif and the disordered RG/RGG region. ATP binding inhibited RNA‐driven phase separation of CIRBP. Thus, it might be beneficial to include cellular metabolites in in vitro liquid–liquid phase separation studies of RG/RGG and other condensate‐forming proteins in order to better mimic the cellular environment in the future.

Bacterial Growth Inhibition Screen (BGIS) identifies a loss-of-function mutant of the DEK oncogene, indicating DNA modulating activities of DEK in chromatin

The DEK oncoprotein regulates cellular chromatin function via a number of protein-protein interactions. However, the biological relevance of its unique pseudo-SAP/SAP-box domain, which transmits DNA modulating activities in vitro, remains largely speculative. As hypothesis-driven mutations failed to yield DNA-binding null (DBN) mutants, we combined random mutagenesis with the Bacterial Growth Inhibition Screen (BGIS) to overcome this bottleneck. Re-expression of a DEK-DBN mutant in newly established human DEK knockout cells failed to reduce the increase in nuclear size as compared to wild type, indicating roles for DEK-DNA interactions in cellular chromatin organization. Our results extend the functional roles of DEK in metazoan chromatin and highlight the predictive ability of recombinant protein toxicity in E. coli for unbiased studies of eukaryotic DNA modulating protein domains.

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.

Evaluation of aromatic amino acids as potential biomarkers in breast cancer by Raman spectroscopy analysis

The scope of the work undertaken in this paper was to explore the feasibility and reliability of using the Raman signature of aromatic amino acids as a marker in the detection of the presence of breast cancer and perhaps, even the prediction of cancer development in very early stages of cancer onset. To be able to assess this hypothesis, we collected most recent and relevant literature in which Raman spectroscopy was used as an analytical tool in the evaluation of breast cell lines and breast tissue, re-analyzed all the Raman spectra, and extracted all spectral bands from each spectrum that were indicative of aromatic amino acids. The criteria for the consideration of the various papers for this study, and hence, the inclusion of the data that they contained were two-fold: (1) The papers had to focus on the characterization of breast tissue with Raman spectroscopy, and (2) the spectra provided within these papers included the spectral range of 500-1200 cm-1, which constitutes the characteristic region for aromatic amino acid vibrational modes. After all the papers that satisfied these criteria were collected, the relevant spectra from each paper were extracted, processed, normalized. All data were then plotted without bias in order to decide whether there is a pattern that can shed light on a possible diagnostic classification. Remarkably, we have been able to demonstrate that cancerous breast tissues and cells decidedly exhibit overexpression of aromatic amino acids and that the difference between the extent of their presence in cancerous cells and healthy cells is overwhelming. On the basis of this analysis, we conclude that it is possible to use the signature Raman bands of aromatic amino acids as a biomarker for the detection, evaluation and diagnosis of breast cancer.

Interdependent YpsA- and YfhS-Mediated Cell Division and Cell Size Phenotypes in Bacillus subtilis

Although many bacterial cell division factors have been uncovered over the years, evidence from recent studies points to the existence of yet-to-be-discovered factors involved in cell division regulation. Thus, it is important to identify factors and conditions that regulate cell division to obtain a better understanding of this fundamental biological process. We recently reported that in the Gram-positive organisms Bacillus subtilis and Staphylococcus aureus, increased production of YpsA resulted in cell division inhibition. In this study, we isolated spontaneous suppressor mutations to uncover critical residues of YpsA and the pathways through which YpsA may exert its function. Using this technique, we were able to isolate four unique intragenic suppressor mutations in ypsA (E55D, P79L, R111P, and G132E) that rendered the mutated YpsA nontoxic upon overproduction. We also isolated an extragenic suppressor mutation in yfhS, a gene that encodes a protein of unknown function. Subsequent analysis confirmed that cells lacking yfhS were unable to undergo filamentation in response to YpsA overproduction. We also serendipitously discovered that YfhS may play a role in cell size regulation. Finally, we provide evidence showing a mechanistic link between YpsA and YfhS.IMPORTANCEBacillus subtilis is a rod-shaped Gram-positive model organism. The factors fundamental to the maintenance of cell shape and cell division are of major interest. We show that increased expression of ypsA results in cell division inhibition and impairment of colony formation on solid medium. Colonies that do arise possess compensatory suppressor mutations. We have isolated multiple intragenic (within ypsA) mutants and an extragenic suppressor mutant. Further analysis of the extragenic suppressor mutation led to a protein of unknown function, YfhS, which appears to play a role in regulating cell size. In addition to confirming that the cell division phenotype associated with YpsA is disrupted in a yfhS-null strain, we also discovered that the cell size phenotype of the yfhS knockout mutant is abolished in a strain that also lacks ypsA This highlights a potential mechanistic link between these two proteins; however, the underlying molecular mechanism remains to be elucidated.

Genome-Wide Identification of RNA Silencing-Related Genes and Their Expressional Analysis in Response to Heat Stress in Barley (Hordeum vulgare L.)

Barley (Hordeum vulgare L.) is an economically important crop cultivated in temperate climates all over the world. Adverse environmental factors negatively affect its survival and productivity. RNA silencing is a conserved pathway involved in the regulation of growth, development and stress responses. The key components of RNA silencing are the Dicer-like proteins (DCLs), Argonautes (AGOs) and RNA-dependent RNA polymerases (RDRs). Despite its economic importance, there is no available comprehensive report on barley RNA silencing machinery and its regulation. In this study, we in silico identified five DCL (HvDCL), eleven AGO (HvAGO) and seven RDR (HvRDR) genes in the barley genome. Genomic localization, phylogenetic analysis, domain organization and functional/catalytic motif identification were also performed. To understand the regulation of RNA silencing, we experimentally analysed the transcriptional changes in response to moderate, persistent or gradient heat stress treatments: transcriptional accumulation of siRNA- but not miRNA-based silencing factor was consistently detected. These results suggest that RNA silencing is dynamically regulated and may be involved in the coordination of development and environmental adaptation in barley. In summary, our work provides information about barley RNA silencing components and will be a ground for the selection of candidate factors and in-depth functional/mechanistic analyses.

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.

Structural basis of chimpanzee APOBEC3H dimerization stabilized by double-stranded RNA

APOBEC3H (A3H) is a mammal-specific cytidine deaminase that potently restricts the replication of retroviruses. Primate A3Hs are known to exert key selective pressures against the cross-species transmission of primate immunodeficiency viruses from chimpanzees to humans. Despite recent advances, the molecular structures underlying the functional mechanisms of primate A3Hs have not been fully understood. Here, we reveal the 2.20-Å crystal structure of the chimpanzee A3H (cpzA3H) dimer bound to a short double-stranded RNA (dsRNA), which appears to be similar to two recently reported structures of pig-tailed macaque A3H and human A3H. In the structure, the dsRNA-binding interface forms a specialized architecture with unique features. The analysis of the dsRNA nucleotides in the cpzA3H complex revealed the GC-rich palindrome-like sequence preference for dsRNA interaction, which is largely determined by arginine residues in loop 1. In cells, alterations of the cpzA3H residues critical for the dsRNA interaction severely reduce intracellular protein stability due to proteasomal degradation. This suggests that cpzA3H stability is regulated by the dsRNA-mediated dimerization as well as by unknown cellular machinery through proteasomal degradation in cells. Taken together, these findings highlight unique structural features of primate A3Hs that are important to further understand their cellular functions and regulation.

DNA-protein interaction: identification, prediction and data analysis

Life in living organisms is dependent on specific and purposeful interaction between other molecules. Such purposeful interactions make the various processes inside the cells and the bodies of living organisms possible. DNA-protein interactions, among all the types of interactions between different molecules, are of considerable importance. Currently, with the development of numerous experimental techniques, diverse methods are convenient for recognition and investigating such interactions. While the traditional experimental techniques to identify DNA-protein complexes are time-consuming and are unsuitable for genome-scale studies, the current high throughput approaches are more efficient in determining such interaction at a large-scale, but they are clearly too costly to be practice for daily applications. Hence, according to the availability of much information related to different biological sequences and clearing different dimensions of conditions in which such interactions are formed, with the developments related to the computer, mathematics, and statistics motivate scientists to develop bioinformatics tools for prediction the interaction site(s). Until now, there has been much progress in this field. In this review, the factors and conditions governing the interaction and the laboratory techniques for examining such interactions are addressed. In addition, developed bioinformatics tools are introduced and compared for this reason and, in the end, several suggestions are offered for the promotion of such tools in prediction with much more precision.

Aromatic SOFAST-HMBC for proteins at natural 13C abundance

We propose here SOFAST-HMBC as a new complementary NMR tool for aromatic side chain assignment of protein samples at natural ¹³C abundance. The characteristic peak patterns detected in SOFAST-HMBC for each aromatic side chain allow straightforward assignment of all protons and carbons (including quaternary ones) of the aromatic ring, and for tyrosine and phenylalanine, connection to the CB of the aliphatic chain. The performance of SOFAST-HMBC is demonstrated for three small proteins (7-14 kDa) at millimolar sample concentration using modern high-field NMR instruments equipped with cryogenically cooled probes. Despite the low amount of NMR-active ¹³C nuclei in these samples, ¹H-¹³C multiple-bond correlation spectra of good quality were obtained in reasonable experimental times of typically less than 24 h.

Trinucleotide Base Pair Stacking Free Energy for Understanding TF-DNA Recognition and the Functions of SNPs

Single nucleotide polymorphisms (SNPs) affect base pair stacking, which is the primary factor for maintaining the stability of DNA. However, the mechanism of how SNPs lead to phenotype variations is still unclear. In this work, we connected SNPs and base pair stacking by a 3-mer base pair stacking free energy matrix. The SNPs with large base pair stacking free energy differences led to phenotype variations. A molecular dynamics (MD) simulation was then applied. Our results showed that base pair stacking played an important role in the transcription factor (TF)-DNA interaction. Changes in DNA structure mainly originate from TF-DNA interactions, and with the increased base pair stacking free energy, the structure of DNA approaches its free type, although its binding affinity was increased by the SNP. In addition, quantitative models using base pair stacking features revealed that base pair stacking can be used to predict TF binding specificity. As such, our work combined knowledge from bioinformatics and structural biology and provided a new understanding of the relationship between SNPs and phenotype variations. The 3-mer base pair stacking free energy matrix is useful in high-throughput screening of SNPs and predicting TF-DNA binding affinity.

BEST and SOFAST experiments for resonance assignment of histidine and tyrosine side chains in 13C/15N labeled proteins

Aromatic amino-acid side chains are essential components for the structure and function of proteins. We present herein a set of NMR experiments for time-efficient resonance assignment of histidine and tyrosine side chains in uniformly ¹³C/¹⁵N-labeled proteins. The use of band-selective ¹³C pulses allows to deal with linear chains of coupled spins, thus avoiding signal loss that occurs in branched spin systems during coherence transfer. Furthermore, our pulse schemes make use of longitudinal ¹H relaxation enhancement, Ernst-angle excitation, and simultaneous detection of ¹H and ¹³C steady-state polarization to achieve significant signal enhancements.

Structural and functional dissection of differentially expressed tomato WRKY transcripts in host defense response against the vascular wilt pathogen (Fusarium oxysporum f. sp. lycopersici)

The WRKY transcription factors have indispensable role in plant growth, development and defense responses. The differential expression of WRKY genes following the stress conditions has been well demonstrated. We investigated the temporal and tissue-specific (root and leaf tissues) differential expression of plant defense-related WRKY genes, following the infection of Fusarium oxysporum f. sp. lycopersici (Fol) in tomato. The genome-wide computational analysis revealed that during the Fol infection in tomato, 16 different members of WRKY gene superfamily were found to be involved, of which only three WRKYs (SolyWRKY4, SolyWRKY33, and SolyWRKY37) were shown to have clear-cut differential gene expression. The quantitative real time PCR (qRT-PCR) studies revealed different gene expression profile changes in tomato root and leaf tissues. In root tissues, infected with Fol, an increased expression for SolyWRKY33 (2.76 fold) followed by SolyWRKY37 (1.93 fold) gene was found at 24 hrs which further increased at 48 hrs (5.0 fold). In contrast, the leaf tissues, the expression was more pronounced at an earlier stage of infection (24 hrs). However, in both cases, we found repression of SolyWRKY4 gene, which further decreased at an increased time interval. The biochemical defense programming against Fol pathogenesis was characterized by the highest accumulation of H₂O₂ (at 48 hrs) and enhanced lignification. The functional diversity across the characterized WRKYs was explored through motif scanning using MEME suite, and the WRKYs specific gene regulation was assessed through the DNA protein docking studies The functional WRKY domain modeled had β sheets like topology with coil and turns. The DNA-protein interaction results revealed the importance of core residues (Tyr, Arg, and Lys) in a feasible WRKY-W-box DNA interaction. The protein interaction network analysis revealed that the SolyWRKY33 could interact with other proteins, such as mitogen-activated protein kinase 5 (MAPK), sigma factor binding protein1 (SIB1) and with other WRKY members including WRKY70, WRKY1, and WRKY40, to respond various biotic and abiotic stresses. The STRING results were further validated through Predicted Tomato Interactome Resource (PTIR) database. The CELLO2GO web server revealed the functional gene ontology annotation and protein subcellular localization, which predicted that SolyWRKY33 is involved in amelioration of biological stress (39.3%) and other metabolic processes (39.3%). The protein (SolyWRKY33) most probably located inside the nucleus (91.3%) with having transcription factor binding activity. We conclude that the defense response following the Fol challenge was accompanied by differential expression of the SolyWRKY4(↓), SolyWRKY33(↑) and SolyWRKY37(↑) transcripts. The biochemical changes are occupied by elicitation of H₂O₂ generation and accumulation and enhanced lignified tissues.

Comprehensive review and empirical analysis of hallmarks of DNA-, RNA- and protein-binding residues in protein chains

Proteins interact with a variety of molecules including proteins and nucleic acids. We review a comprehensive collection of over 50 studies that analyze and/or predict these interactions. While majority of these studies address either solely protein–DNA or protein–RNA binding, only a few have a wider scope that covers both protein–protein and protein–nucleic acid binding. Our analysis reveals that binding residues are typically characterized with three hallmarks: relative solvent accessibility (RSA), evolutionary conservation and propensity of amino acids (AAs) for binding. Motivated by drawbacks of the prior studies, we perform a large-scale analysis to quantify and contrast the three hallmarks for residues that bind DNA-, RNA-, protein- and (for the first time) multi-ligand-binding residues that interact with DNA and proteins, and with RNA and proteins. Results generated on a well-annotated data set of over 23 000 proteins show that conservation of binding residues is higher for nucleic acid- than protein-binding residues. Multi-ligand-binding residues are more conserved and have higher RSA than single-ligand-binding residues. We empirically show that each hallmark discriminates between binding and nonbinding residues, even predicted RSA, and that combining them improves discriminatory power for each of the five types of interactions. Linear scoring functions that combine these hallmarks offer good predictive performance of residue-level propensity for binding and provide intuitive interpretation of predictions. Better understanding of these residue-level interactions will facilitate development of methods that accurately predict binding in the exponentially growing databases of protein sequences.

Using Separation-of-Function Mutagenesis To Define the Full Spectrum of Activities Performed by the Est1 Telomerase Subunit in Vivo

A leading objective in biology is to identify the complete set of activities performed by each gene. Identification of a comprehensive set of separation...

A leading objective in biology is to identify the complete set of activities that each gene performs in vivo. In this study, we have asked whether a genetic approach can provide an efficient means of achieving this goal, through the identification and analysis of a comprehensive set of separation-of-function (sof⁻) mutations in a gene. Toward this goal, we have subjected the Saccharomyces cerevisiae EST1 gene, which encodes a regulatory subunit of telomerase, to intensive mutagenesis (with an average coverage of one mutation for every 4.5 residues), using strategies that eliminated those mutations that disrupted protein folding/stability. The resulting set of sof⁻ mutations defined four biochemically distinct activities for the Est1 telomerase protein: two temporally separable steps in telomerase holoenzyme assembly, a telomerase recruitment activity, and a fourth newly discovered regulatory function. Although biochemically distinct, impairment of each of these four different activities nevertheless conferred a common phenotype (critically short telomeres) comparable to that of an est1-∆ null strain. This highlights the limitations of gene deletions, even for nonessential genes; we suggest that employing a representative set of sof⁻ mutations for each gene in future high- and low-throughput investigations will provide deeper insights into how proteins interact inside the cell.

The Structure of Ethylbenzene, Styrene and Phenylacetylene Determined by Total Neutron Scattering

Organic solvents such as phenylacetylene, styrene and ethylbenzene are widely used in industrial processes, especially in the production of rubber or thermoplastics. Despite their important applications detailed knowledge about their structure is limited. In this paper the structures of these three aromatic solvents were investigated using neutron diffraction. The results show that many of their structural characteristics are similar, although the structure of phenylacetylene is more ordered and has a smaller solvation sphere than either ethylbenzene or styrene. Two regions within the first coordination sphere, in which the surrounding molecules show different preferable orientations with respect to the central molecule, were found for each liquid. Additionally, the localisation of the aliphatic chains reveals that they tend to favour closer interactions with each other than to the aromatic rings of the adjacent molecules.

Structure-based prediction of transcription factor binding specificity using an integrative energy function

Transcription factors (TFs) regulate gene expression through binding to specific target DNA sites. Accurate annotation of transcription factor binding sites (TFBSs) at genome scale represents an essential step toward our understanding of gene regulation networks. In this article, we present a structure-based method for computational prediction of TFBSs using a novel, integrative energy (IE) function. The new energy function combines a multibody (MB) knowledge-based potential and two atomic energy terms (hydrogen bond and π interaction) that might not be accurately captured by the knowledge-based potential owing to the mean force nature and low count problem. We applied the new energy function to the TFBS prediction using a non-redundant dataset that consists of TFs from 12 different families. Our results show that the new IE function improves the prediction accuracy over the knowledge-based, statistical potentials, especially for homeodomain TFs, the second largest TF family in mammals.

CONTACT

jguo4@uncc.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Direct Comparison of Amino Acid and Salt Interactions with Double-Stranded and Single-Stranded DNA from Explicit-Solvent Molecular Dynamics Simulations

Given the ubiquitous nature of protein-DNA interactions, it is important to understand the interaction thermodynamics of individual amino acid side chains for DNA. One way to assess these preferences is to perform molecular dynamics (MD) simulations. Here we report MD simulations of 20 amino acid side chain analogs interacting simultaneously with both a 70-base-pair double-stranded DNA and with a 70-nucleotide single-stranded DNA. The relative preferences of the amino acid side chains for dsDNA and ssDNA match well with values deduced from crystallographic analyses of protein-DNA complexes. The estimated apparent free energies of interaction for ssDNA, on the other hand, correlate well with previous simulation values reported for interactions with isolated nucleobases, and with experimental values reported for interactions with guanosine. Comparisons of the interactions with dsDNA and ssDNA indicate that, with the exception of the positively charged side chains, all types of amino acid side chain interact more favorably with ssDNA, with intercalation of aromatic and aliphatic side chains being especially notable. Analysis of the data on a base-by-base basis indicates that positively charged side chains, as well as sodium ions, preferentially bind to cytosine in ssDNA, and that negatively charged side chains, and chloride ions, preferentially bind to guanine in ssDNA. These latter observations provide a novel explanation for the lower salt dependence of DNA duplex stability in GC-rich sequences relative to AT-rich sequences.

Identification of Position-Specific Correlations between DNA-Binding Domains and Their Binding Sites. Application to the MerR Family of Transcription Factors

The large and increasing volume of genomic data analyzed by comparative methods provides information about transcription factors and their binding sites that, in turn, enables statistical analysis of correlations between factors and sites, uncovering mechanisms and evolution of specific protein-DNA recognition. Here we present an online tool, Prot-DNA-Korr, designed to identify and analyze crucial protein-DNA pairs of positions in a family of transcription factors. Correlations are identified by analysis of mutual information between columns of protein and DNA alignments. The algorithm reduces the effects of common phylogenetic history and of abundance of closely related proteins and binding sites. We apply it to five closely related subfamilies of the MerR family of bacterial transcription factors that regulate heavy metal resistance systems. We validate the approach using known 3D structures of MerR-family proteins in complexes with their cognate DNA binding sites and demonstrate that a significant fraction of correlated positions indeed form specific side-chain-to-base contacts. The joint distribution of amino acids and nucleotides hence may be used to predict changes of specificity for point mutations in transcription factors.

DNA-HMGB1 interaction: The nuclear aggregates of polyamine mediation

Nuclear aggregates of polyamines (NAPs) are supramolecular compounds generated by the self-assembly of protonated nuclear polyamines (spermine, spermidine and putrescine) and phosphate ions. In the presence of genomic DNA, the hierarchical process of self-structuring ultimately produces nanotube-like polymers that envelop the double helix. Because of their modular nature and their aggregation-disaggregation dynamics, NAPs confer plasticity and flexibility to DNA. Through the disposition of charges, NAPs also enable a bidirectional stream of information between the genome and interacting moieties. High mobility group (HMG) B1 is a non-histone chromosomal protein that binds to DNA and that influences multiple nuclear processes. Because genomic DNA binds to either NAPs or HMGB1 protein, we explored the ability of in vitro self-assembled NAPs (ivNAPs) to mediate the DNA-HMGB1 interaction. To this end, we structured DNA-NAPs-HMGB1 and DNA-HMGB1-NAPs ternary complexes in vitro through opportune sequential incubations. Mobility shift electrophoresis and atomic force microscopy showed that the DNA-ivNAPs-HGMB1 complex had conformational assets supposedly more suitable those of the DNA-HGMB1-ivNAPs to comply with the physiological and functional requirements of DNA. Our findings indicated that ivNAPs act as mediators of the DNA-HMGB1 interaction.

Statistical analysis of structural determinants for protein-DNA-binding specificity

DNA-binding proteins play critical roles in biological processes including gene expression, DNA packaging and DNA repair. They bind to DNA target sequences with different degrees of binding specificity, ranging from highly specific (HS) to nonspecific (NS). Alterations of DNA-binding specificity, due to either genetic variation or somatic mutations, can lead to various diseases. In this study, a comparative analysis of protein-DNA complex structures was carried out to investigate the structural features that contribute to binding specificity. Protein-DNA complexes were grouped into three general classes based on degrees of binding specificity: HS, multispecific (MS), and NS. Our results show a clear trend of structural features among the three classes, including amino acid binding propensities, simple and complex hydrogen bonds, major/minor groove and base contacts, and DNA shape. We found that aspartate is enriched in HS DNA binding proteins and predominately binds to a cytosine through a single hydrogen bond or two consecutive cytosines through bidentate hydrogen bonds. Aromatic residues, histidine and tyrosine, are highly enriched in the HS and MS groups and may contribute to specific binding through different mechanisms. To further investigate the role of protein flexibility in specific protein-DNA recognition, we analyzed the conformational changes between the bound and unbound states of DNA-binding proteins and structural variations. The results indicate that HS and MS DNA-binding domains have larger conformational changes upon DNA-binding and larger degree of flexibility in both bound and unbound states. Proteins 2016; 84:1147-1161. © 2016 Wiley Periodicals, Inc.

Topology of RNA-protein nucleobase-amino acid π-π interactions and comparison to analogous DNA-protein π-π contacts

The present work analyzed 120 high-resolution X-ray crystal structures and identified 335 RNA-protein π-interactions (154 nonredundant) between a nucleobase and aromatic (W, H, F, or Y) or acyclic (R, E, or D) π-containing amino acid. Each contact was critically analyzed (including using a visual inspection protocol) to determine the most prevalent composition, structure, and strength of π-interactions at RNA-protein interfaces. These contacts most commonly involve F and U, with U:F interactions comprising one-fifth of the total number of contacts found. Furthermore, the RNA and protein π-systems adopt many different relative orientations, although there is a preference for more parallel (stacked) arrangements. Due to the variation in structure, the strength of the intermolecular forces between the RNA and protein components (as determined from accurate quantum chemical calculations) exhibits a significant range, with most of the contacts providing significant stability to the associated RNA-protein complex (up to -65 kJ mol(-1)). Comparison to the analogous DNA-protein π-interactions emphasizes differences in RNA- and DNA-protein π-interactions at the molecular level, including the greater abundance of RNA contacts and the involvement of different nucleobase/amino acid residues. Overall, our results provide a clearer picture of the molecular basis of nucleic acid-protein binding and underscore the important role of these contacts in biology, including the significant contribution of π-π interactions to the stability of nucleic acid-protein complexes. Nevertheless, more work is still needed in this area in order to further appreciate the properties and roles of RNA nucleobase-amino acid π-interactions in nature.

Signature motifs of GDP polyribonucleotidyltransferase, a non-segmented negative strand RNA viral mRNA capping enzyme, domain in the L protein are required for covalent enzyme-pRNA intermediate formation

The unconventional mRNA capping enzyme (GDP polyribonucleotidyltransferase, PRNTase; block V) domain in RNA polymerase L proteins of non-segmented negative strand (NNS) RNA viruses (e.g. rabies, measles, Ebola) contains five collinear sequence elements, Rx(3)Wx(3–8)ΦxGxζx(P/A) (motif A; Φ, hydrophobic; ζ, hydrophilic), (Y/W)ΦGSxT (motif B), W (motif C), HR (motif D) and ζxxΦx(F/Y)QxxΦ (motif E). We performed site-directed mutagenesis of the L protein of vesicular stomatitis virus (VSV, a prototypic NNS RNA virus) to examine participation of these motifs in mRNA capping. Similar to the catalytic residues in motif D, G1100 in motif A, T1157 in motif B, W1188 in motif C, and F1269 and Q1270 in motif E were found to be essential or important for the PRNTase activity in the step of the covalent L-pRNA intermediate formation, but not for the GTPase activity that generates GDP (pRNA acceptor). Cap defective mutations in these residues induced termination of mRNA synthesis at position +40 followed by aberrant stop–start transcription, and abolished virus gene expression in host cells. These results suggest that the conserved motifs constitute the active site of the PRNTase domain and the L-pRNA intermediate formation followed by the cap formation is essential for successful synthesis of full-length mRNAs.

Implications of aromatic-aromatic interactions: From protein structures to peptide models

With increasing structural information on proteins, the opportunity to understand physical forces governing protein folding is also expanding. One of the significant non-covalent forces between the protein side chains is aromatic-aromatic interactions. Aromatic interactions have been widely exploited and thoroughly investigated in the context of folding, stability, molecular recognition, and self-assembly processes. Through this review, we discuss the contribution of aromatic interactions to the activity and stability of thermophilic, mesophilic, and psychrophilic proteins. Being hydrophobic, aromatic amino acids tend to reside in the protein hydrophobic interior or transmembrane segments of proteins. In such positions, it can play a diverse role in soluble and membrane proteins, and in α-helix and β-sheet stabilization. We also highlight here some excellent investigations made using peptide models and several approaches involving aryl-aryl interactions, as an increasingly popular strategy in protein and peptide engineering. A recent survey described the existence of aromatic clusters (trimer, tetramer, pentamer, and higher order assemblies), revealing the self-associating property of aryl groups, even in folded protein structures. The application of this self-assembly of aromatics in the generation of modern bionanomaterials is also discussed.

The peptide NCbz-Val-Tyr-OMe and aromatic π-π interactions

The peptide N-benzyloxycarbonyl-L-valyl-L-tyrosine methyl ester or NCbz-Val-Tyr-OMe (where NCbz is N-benzyloxycarbonyl and OMe indicates the methyl ester), C(23)H(28)N(2)O(6), has an extended backbone conformation. The aromatic rings of the Tyr residue and the NCbz group are involved in various attractive intra- and intermolecular aromatic π-π interactions which stabilize the conformation and packing in the crystal structure, in addition to N-H...O and O-H...O hydrogen bonds. The aromatic π-π interactions include parallel-displaced, perpendicular T-shaped, perpendicular L-shaped and inclined orientations.

A thermodynamic insight into the recognition of hydrophilic and hydrophobic amino acids in pure water by aza-scorpiand type receptors

Thermodynamic studies about the interaction of scorpiand aza-macrocycles with amino acids in water show entropy driven stabilisations often associated with solvation/desolvation processes.

Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur

Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.

Distinguishing closely related amyloid precursors using an RNA aptamer

Although amyloid fibrils assembled in vitro commonly involve a single protein, fibrils formed in vivo can contain multiple protein sequences. The amyloidogenic protein human β2-microglobulin (hβ2m) can co-polymerize with its N-terminally truncated variant (ΔN6) in vitro to form hetero-polymeric fibrils that differ from their homo-polymeric counterparts. Discrimination between the different assembly precursors, for example by binding of a biomolecule to one species in a mixture of conformers, offers an opportunity to alter the course of co-assembly and the properties of the fibrils formed. Here, using hβ2m and its amyloidogenic counterpart, ΔΝ6, we describe selection of a 2'F-modified RNA aptamer able to distinguish between these very similar proteins. SELEX with a N30 RNA pool yielded an aptamer (B6) that binds hβ2m with an EC50 of ∼200 nM. NMR spectroscopy was used to assign the (1)H-(15)N HSQC spectrum of the B6-hβ2m complex, revealing that the aptamer binds to the face of hβ2m containing the A, B, E, and D β-strands. In contrast, binding of B6 to ΔN6 is weak and less specific. Kinetic analysis of the effect of B6 on co-polymerization of hβ2m and ΔN6 revealed that the aptamer alters the kinetics of co-polymerization of the two proteins. The results reveal the potential of RNA aptamers as tools for elucidating the mechanisms of co-assembly in amyloid formation and as reagents able to discriminate between very similar protein conformers with different amyloid propensity.

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

Four hundred twenty-eight high-resolution DNA-protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π-contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π-π stacking, π-π T-shaped or sugar-π contacts. Three hundred forty-four nucleobase-amino acid π-π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA-protein sugar-π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C-H···π and/or lone-pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π-π and sugar-π interactions display a range of relative monomer orientations and therefore interaction energies (up to -50 (-70) kJ mol(-1) for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA-protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.

Aromatic interactions modulate the 5'-base selectivity of the DNA-binding autoantibody ED-10

We present detailed computational analyses of the binding of four dinucleotides to a highly sequence-selective single-stranded DNA (ssDNA) binding antibody (ED-10) and selected point mutants. Anti-DNA antibodies are central to the pathogenesis of systemic lupus erythematosus (SLE), and a more complete understanding of the mode of binding of DNA and other ligands will be necessary to elucidate the role of anti-DNA antibodies in the kidney inflammation associated with SLE. Classical molecular mechanics based molecular dynamics simulations and density functional theory (DFT) computations were applied to pinpoint the origin of selectivity for the 5'-nucleotide. In particular, the strength of interactions between each nucleotide and the surrounding residues were computed using MMGBSA as well as DFT applied to a cluster model of the binding site. The results agree qualitatively with experimental binding free energies, and indicate that π-stacking, CH/π, NH/π, and hydrogen-bonding interactions all contribute to 5'-base selectivity in ED-10. Most importantly, the selectivity for dTdC over dAdC arises primarily from differences in the strength of π-stacking and XH/π interactions with the surrounding aromatic residues; hydrogen bonds play little role. These data suggest that a key Tyr residue, which is not present in other anti-DNA antibodies, plays a key role in the 5'-base selectivity, while we predict that the mutation of a single Trp residue can tune the selectivity for dTdC over dAdC.

Examination of tyrosine/adenine stacking interactions in protein complexes

The π-stacking interactions between tyrosine amino acid side chains and adenine-bearing ligands are examined. Crystalline protein structures from the protein data bank (PDB) exhibiting face-to-face tyrosine/adenine arrangements were used to construct 20 unique 4-methylphenol/N9-methyladenine (p-cresol/9MeA) model systems. Full geometry optimization of the 20 crystal structures with the M06-2X density functional theory method identified 11 unique low-energy conformations. CCSD(T) complete basis set (CBS) limit interaction energies were estimated for all of the structures to determine the magnitude of the interaction between the two ring systems. CCSD(T) computations with double-ζ basis sets (e.g., 6-31G*(0.25) and aug-cc-pVDZ) indicate that the MP2 method overbinds by as much as 3.07 kcal mol(-1) for the crystal structures and 3.90 kcal mol(-1) for the optimized structures. In the 20 crystal structures, the estimated CCSD(T) CBS limit interaction energy ranges from -4.00 to -6.83 kcal mol(-1), with an average interaction energy of -5.47 kcal mol(-1), values remarkably similar to the corresponding data for phenylalanine/adenine stacking interactions. Geometry optimization significantly increases the interaction energies of the p-cresol/9MeA model systems. The average estimated CCSD(T) CBS limit interaction energy of the 11 optimized structures is 3.23 kcal mol(-1) larger than that for the 20 crystal structures.

Significant strength of charged DNA-protein π-π interactions: a preliminary study of cytosine

The present work characterized the preferred gas-phase structure and optimum interaction energy of both parallel stacked and perpendicular T-shaped dimers between cytosine (C), as a representative nucleobase, and aspartic/glutamic acid (DE), aspartate/glutamate (DE(-)) or arginine (R(+)), using detailed M06-2X/6-31+G(d,p) potential energy surface scans as a function of the relative monomer orientation. Through comparison to previous literature on the π-π interactions between the DNA nucleobases and the aromatic amino acid residues, this work will allow for comparisons between DNA-protein interactions involving aromatic and acyclic R-side chains, as well as comparisons of the relative geometric dependence and magnitude of π-π (C:DE), πcation-π (C:R(+)), and πanion-π (C:DE(-)) interactions. Our results show that the preferred relative monomer orientation is highly dependent on the monomer composition and charge, and is dictated by electrostatic-driven interactions. More importantly, for the first time, we report that the π-π interactions between cytosine and (neutral) aspartic/glutamic acid are up to approximately -40 kJ mol(-1), while the πcation-π or πanion-π interactions between cytosine and arginine or aspartate/glutamate are up to approximately -90 and -99 kJ mol(-1), respectively. An extensive investigation of the effects of the computational methodology implemented, including comparisons to detailed CCSD(T)/CBS potential energy surfaces and interaction energies, supports the use of M06-2X, as well as ωB97X-D, to study DNA-protein π-π interactions of varying composition and charge. Most importantly, the CCSD(T)/CBS results verify the strong nature of these DNA-protein π-π interactions, as well as the unique nature of the πcation-π and πanion-π counterparts. Therefore, our results emphasize that a wide variety of different types of noncovalent interactions between both cyclic and acyclic π-containing components can significantly contribute to the stability of DNA-protein complexes and likely play a larger role in biology than currently accepted.

Uses of phage display in agriculture: sequence analysis and comparative modeling of late embryogenesis abundant client proteins suggest protein-nucleic acid binding functionality

A group of intrinsically disordered, hydrophilic proteins-Late Embryogenesis Abundant (LEA) proteins-has been linked to survival in plants and animals in periods of stress, putatively through safeguarding enzymatic function and prevention of aggregation in times of dehydration/heat. Yet despite decades of effort, the molecular-level mechanisms defining this protective function remain unknown. A recent effort to understand LEA functionality began with the unique application of phage display, wherein phage display and biopanning over recombinant Seed Maturation Protein homologs from Arabidopsis thaliana and Glycine max were used to retrieve client proteins at two different temperatures, with one intended to represent heat stress. From this previous study, we identified 21 client proteins for which clones were recovered, sometimes repeatedly. Here, we use sequence analysis and homology modeling of the client proteins to ascertain common sequence and structural properties that may contribute to binding affinity with the protective LEA protein. Our methods uncover what appears to be a predilection for protein-nucleic acid interactions among LEA client proteins, which is suggestive of subcellular residence. The results from this initial computational study will guide future efforts to uncover the protein protective mechanisms during heat stress, potentially leading to phage-display-directed evolution of synthetic LEA molecules.

On the importance and origin of aromatic interactions in chemistry and biodisciplines

Aromatic systems contain both σ- and π-electrons, which in turn constitute σ- and π-molecular orbitals (MOs). In discussing the properties of these systems, researchers typically refer to the highest occupied and lowest unoccupied MOs, which are π MOs. The characteristic properties of aromatic systems, such as their low ionization potentials and electron affinities, high polarizabilities and stabilities, and small band gaps (in spectroscopy called the N → V1 space), can easily be explained based on their electronic structure. These one-electron properties point to characteristic features of how aromatic systems interact with each other. Unlike hydrogen bonding systems, which primarily interact through electrostatic forces, complexes containing aromatic systems, especially aromatic stacked pairs, are predominantly stabilized by dispersion attraction. The stabilization energy in the benzene dimer is rather small (~2.5 kcal/mol) but strengthens with heteroatom substitution. The stacked interaction of aromatic nucleic acid bases is greater than 10 kcal/mol, and for the most stable stacked pair, guanine and cytosine, it reaches approximately 17 kcal/mol. Although these values do not equal the planar H-bonded interactions of these bases (~29 kcal/mol), stacking in DNA is more frequent than H-bonding and, unlike H-bonding, is not significantly weakened when passing from the gas phase to a water environment. Consequently, the stacking of aromatic systems represents the leading stabilization energy contribution in biomacromolecules and in related nanosystems. Therefore stacking (dispersion) interactions predominantly determine the double helical structure of DNA, which underlies its storage and transfer of genetic information. Similarly, dispersion is the dominant contributor to attractive interactions involving aromatic amino acids within the hydrophobic core of a protein, which is critical for folding. Therefore, understanding the nature of aromatic interactions, which depend greatly on quantum mechanical (QM) calculations, is of key importance in biomolecular science. This Account shows that accurate binding energies for aromatic complexes should be based on computations made at the (estimated) CCSD(T)/complete basis set limit (CBS) level of theory. This method is the least computationally intensive one that can give accurate stabilization energies for all common classes of noncovalent interactions (aromatic-aromatic, H-bonding, ionic, halogen bonding, charge-transfer, etc.). These results allow for direct comparison of binding energies between different interaction types. Conclusions based on lower-level QM calculations should be considered with care.

Recognition of platinum-DNA adducts by HMGB1a

Cisplatin (CP) and oxaliplatin (OX), platinum-based drugs used widely in chemotherapy, form adducts on intrastrand guanines (5'GG) in genomic DNA. DNA damage recognition proteins, transcription factors, mismatch repair proteins, and DNA polymerases discriminate between CP- and OX-GG DNA adducts, which could partly account for differences in the efficacy, toxicity, and mutagenicity of CP and OX. In addition, differential recognition of CP- and OX-GG adducts is highly dependent on the sequence context of the Pt-GG adduct. In particular, DNA binding protein domain HMGB1a binds to CP-GG DNA adducts with up to 53-fold greater affinity than to OX-GG adducts in the TGGA sequence context but shows much smaller differences in binding in the AGGC or TGGT sequence contexts. Here, simulations of the HMGB1a-Pt-DNA complex in the three sequence contexts revealed a higher number of interface contacts for the CP-DNA complex in the TGGA sequence context than in the OX-DNA complex. However, the number of interface contacts was similar in the TGGT and AGGC sequence contexts. The higher number of interface contacts in the CP-TGGA sequence context corresponded to a larger roll of the Pt-GG base pair step. Furthermore, geometric analysis of stacking of phenylalanine 37 in HMGB1a (Phe37) with the platinated guanines revealed more favorable stacking modes correlated with a larger roll of the Pt-GG base pair step in the TGGA sequence context. These data are consistent with our previous molecular dynamics simulations showing that the CP-TGGA complex was able to sample larger roll angles than the OX-TGGA complex or either CP- or OX-DNA complexes in the AGGC or TGGT sequences. We infer that the high binding affinity of HMGB1a for CP-TGGA is due to the greater flexibility of CP-TGGA compared to OX-TGGA and other Pt-DNA adducts. This increased flexibility is reflected in the ability of CP-TGGA to sample larger roll angles, which allows for a higher number of interface contacts between the Pt-DNA adduct and HMGB1a.