Bioinformatic analysis of the protein/DNA interface

Insights into DNA solvation found in protein-DNA structures

The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.

Events at the Interface: How Do Interfaces Modulate the Dynamics and Functionalities of Guest Molecules?

Chemical and biological interfaces are of various types, which could be between two materials of the same and/or different states, two phases of the same material, biological substrates and the outer environment, surfactant or polymeric membranes and the bulk, and so forth. Small-molecule guests frequently interact with such interfaces that decide their functionalities. The structural and behavioral properties undergo considerable characteristic changes, which control their final course of action in the targeted application. This Perspective will discuss mainly the chemical interfaces constituted by the surfactants, polymers, lipids, and nucleic acids and their impacts on the dynamics of small-molecule guests. Some specific and interesting phenomena and future prospects will be elucidated in this Perspective.

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks

Database-derived water probability densities around structurally and sequentially distinct DNA dinucleotide fragments reproduce the known hydration motifs, which thus can be used as building blocks to predict DNA hydration.

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the ‘spine of hydration’ in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

A deterministic code for transcription factor-DNA recognition through computation of binding interfaces

The recognition code between transcription factor (TF) amino acids and DNA bases remains poorly understood. Here, the determinants of TF amino acid-DNA base binding selectivity were identified through the analysis of crystals of TF-DNA complexes. Selective, high-frequency interactions were identified for the vast majority of amino acid side chains (‘structural code’). DNA binding specificities were then independently assessed by meta-analysis of random-mutagenesis studies of Zn finger-target DNA sequences. Selective, high-frequency interactions were identified for the majority of mutagenized residues (‘mutagenesis code’). The structural code and the mutagenesis code were shown to match to a striking level of accuracy (P = 3.1 × 10−33), suggesting the identification of fundamental rules of TF binding to DNA bases. Additional insight was gained by showing a geometry-dictated choice among DNA-binding TF residues with overlapping specificity. These findings indicate the existence of a DNA recognition mode whereby the physical-chemical characteristics of the interacting residues play a deterministic role. The discovery of this DNA recognition code advances our knowledge on fundamental features of regulation of gene expression and is expected to pave the way for integration with higher-order complexity approaches.

Variability of the DNA Backbone Geometry in DNA-Protein Complexes: Experimental Data Analysis

We have analyzed and compared the available experimental data (PDB) on the backbone geometry of the DNA in solution (NMR), in crystals (X-rays), and in complexes with proteins (X-rays and cryo-electron microscopy). The deoxyribose (pseudorotational angle τ₀) and ε/ζ (BI-BII transition in phosphates) flexibilities are practically the same in the four samples. The α/γ mobility is minimal in crystalline DNA: on the histograms, there is one canonical and one noncanonical t/t peak. The α/γ mobility increases in DNA solutions (three more noncanonical peaks) and is maximal in DNA-protein complexes (another additional peak). On a large amount of data, we have confirmed that the three main degrees of freedom of the sugar-phosphate backbone are "orthogonal": changes in any of the angles τ₀, (ζ-ε), and (γ-α) occur, as a rule, at a constant (usually canonical) value of any other. In the DNA-protein complexes, none of the geometrical parameters commonly used to distinguish the A and B forms of DNA, except for Zp and its simpler analog Zp', show an unambiguous correlation with τ₀. Proteins, binding to DNA, in 59% of cases change the local shape of the helix up to the characteristic of the A-form without switching the deoxyribose conformation from south to north. However, we have found simple local characteristics of one nucleotide that correlate with the angles τ₀ and (ζ-ε). These are the angles C3'C1'N* and C4'C3'P(2), respectively. They are orthogonal in DNA-protein complexes exactly as the pair τ₀ and (ζ-ε). Most characteristics of DNA in complexes with proteins are the same in X-ray and in cryo-EM data, except for the histogram for the angle τ₀. We offer a possible explanation for this difference. We also discuss the artifacts on the ε/ζ histogram for DNA in solutions caused by the currently used NMR refinement protocols.

PDB-2-PBv3.0: An updated protein block database

Our protein block (PB) sequence database PDB-2-PBv1.0 provides PB sequences and dihedral angles for 74,297 protein structures comprising of 103,252 protein chains of Protein Data Bank (PDB) as on 2011. Since there are a lot of practical applications of PB and also as the size of PDB database increases, it becomes necessary to provide the PB sequences for all PDB protein structures. The current updated PDB-2-PBv3.0 contains PB sequences for 147,602 PDB structures comprising of 400,355 protein chains as on October 2019. When compared to our previous version PDB-2-PBv1.0, the current PDB-2-PBv3.0 contains 2- and 4-fold increase in the number of protein structures and chains, respectively. Notably, it provides PB information for any protein chain, regardless of the missing atom records of protein structure data in PDB. It includes protein interaction information with DNA and RNA along with their corresponding functional classes from Nucleic Acid Database (NDB) and PDB. Now, the updated version allows the user to download multiple PB records by parameter search and/or by a given list. This database is freely accessible at http://bioinfo.bdu.ac.in/pb3.

Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA

Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.

Kinetics of the Conformational Transformation between B- and A-Forms in the Drew-Dickerson Dodecamer

Some DNA sequences in crystals and in complexes with proteins can exist in the forms intermediate between the B- and A-DNA. Based on this, it was implied that the B-to-A transition for any DNA molecule should go through these intermediate forms also in kinetics. More precisely, the helix parameter Slide has to change first, and the molecule should take the E-form. After that, the Roll parameter changes. In the present work, we simulated the kinetics of the B-A transition in the Drew-Dickerson dodecamer, a known B-philic DNA oligomer. We used the "sugar" coarse-grained model that reproduces ribose flexibility, preserves sequence specificity, employs implicit water and explicit ions, and offers the possibility to vary friction. As the control parameter of the transition, we chose the volume available for a counterion and considered the change from a large to a small volume. In the described system, the B-to-A conformational transformation proved to correspond to a first-order phase transition. The molecule behaves like a small cluster in the region of such a transition, jumping between the A- and B-forms in a wide range of available volumes. The viscosity of the solvent does not affect the midpoint of the transition but only the overall mobility of the system. All helix parameters change synchronously on average, we have not observed the sequence "Slide first, Roll later" in kinetics, and the E-DNA is not a necessary step for the transition between the B- and A-forms in the studied system. So, the existence of the intermediate DNA forms requires specific conditions, shifting the common balance of interactions: certain nucleotide sequence in specific solution or/and the interaction with some protein.

Motif orientation matters: Structural characterization of TEAD1 recognition of genomic DNA

TEAD transcription factors regulate gene expression through interactions with DNA and other proteins. They are crucial for the development of eukaryotic organisms and to control the expression of genes involved mostly in cell proliferation and differentiation; however, their deregulation can lead to tumorigenesis. To study the interactions of TEAD1 with M-CAT motifs and their inverted versions, the K_D of each complex was determined, and H/D exchange, quantitative chemical cross-linking, molecular docking, and smFRET were utilized for structural characterization. ChIP-qPCR was employed to correlate the results with a cell line model. The results obtained showed that although the inverted motif has 10× higher K_D, the same residues were affected by the presence of M-CAT in both orientations. Molecular docking and smFRET revealed that TEAD1 binds the inverted motif rotated 180°. In addition, the inverted motif was proven to be occupied by TEAD1 in Jurkat cells, suggesting that the low-affinity binding sites present in the human genome may possess biological relevance.

Base editing: advances and therapeutic opportunities

Base editing - the introduction of single-nucleotide variants (SNVs) into DNA or RNA in living cells - is one of the most recent advances in the field of genome editing. As around half of known pathogenic genetic variants are due to SNVs, base editing holds great potential for the treatment of numerous genetic diseases, through either temporary RNA or permanent DNA base alterations. Recent advances in the specificity, efficiency, precision and delivery of DNA and RNA base editors are revealing exciting therapeutic opportunities for these technologies. We expect the correction of single point mutations will be a major focus of future precision medicine.

Structural alphabets for conformational analysis of nucleic acids available at dnatco.datmos.org

A detailed description of the dnatco.datmos.org web server implementing the universal structural alphabet of nucleic acids is presented. It is capable of processing any mmCIF- or PDB-formatted files containing DNA or RNA molecules; these can either be uploaded by the user or supplied as the wwPDB or PDB-REDO structural database access code. The web server performs an assignment of the nucleic acid conformations and presents the results for the intuitive annotation, validation, modeling and refinement of nucleic acids.

The Involvement of the McsB Arginine Kinase in Clp-Dependent Degradation of the MgsR Regulator in Bacillus subtilis

Regulated ATP-dependent proteolysis is a common feature of developmental processes and plays also a crucial role during environmental perturbations such as stress and starvation. The Bacillus subtilis MgsR regulator controls a subregulon within the stress- and stationary phase σ^B regulon. After ethanol exposition and a short time-window of activity, MgsR is ClpXP-dependently degraded with a half-life of approximately 6 min. Surprisingly, a protein interaction analysis with MgsR revealed an association with the McsB arginine kinase and an in vivo degradation assay confirmed a strong impact of McsB on MgsR degradation. In vitro phosphorylation experiments with arginine (R) by lysine (K) substitutions in McsB and its activator McsA unraveled all R residues, which are essentially needed for the arginine kinase reaction. Subsequently, site directed mutagenesis of the MgsR substrate was used to substitute all arginine residues with glutamate (R-E) to mimic arginine phosphorylation and to test their influence on MgsR degradation in vivo. It turned out, that especially the R33E and R94/95E residues (RRPI motif), the latter are adjacently located to the two redox-sensitive cysteines in a 3D model, have the potential to accelerate MgsR degradation. These results imply that selective arginine phosphorylation may have favorable effects for Clp dependent degradation of short-living regulatory proteins. We speculate that in addition to its kinase activity and adaptor function for the ClpC ATPase, McsB might also serve as a proteolytic adaptor for the ClpX ATPase in the degradation mechanism of MgsR.

Accurate Representation of Protein-Ligand Structural Diversity in the Protein Data Bank (PDB)

The number of available protein structures in the Protein Data Bank (PDB) has considerably increased in recent years. Thanks to the growth of structures and complexes, numerous large-scale studies have been done in various research areas, e.g., protein-protein, protein-DNA, or in drug discovery. While protein redundancy was only simply managed using simple protein sequence identity threshold, the similarity of protein-ligand complexes should also be considered from a structural perspective. Hence, the protein-ligand duplicates in the PDB are widely known, but were never quantitatively assessed, as they are quite complex to analyze and compare. Here, we present a specific clustering of protein-ligand structures to avoid bias found in different studies. The methodology is based on binding site superposition, and a combination of weighted Root Mean Square Deviation (RMSD) assessment and hierarchical clustering. Repeated structures of proteins of interest are highlighted and only representative conformations were conserved for a non-biased view of protein distribution. Three types of cases are described based on the number of distinct conformations identified for each complex. Defining these categories decreases by 3.84-fold the number of complexes, and offers more refined results compared to a protein sequence-based method. Widely distinct conformations were analyzed using normalized B-factors. Furthermore, a non-redundant dataset was generated for future molecular interactions analysis or virtual screening studies.

Structural variations within proteins can be as large as variations observed across their homologues

Understanding the structural plasticity of proteins is key to understanding the intricacies of their functions and mechanistic basis. In the current study, we analyzed the available multiple crystal structures of the same protein for the structural differences. For this purpose we used an abstraction of protein structures referred as Protein Blocks (PBs) that was previously established. We also characterized the nature of the structural variations for a few proteins using molecular dynamics simulations. In both the cases, the structural variations were summarized in the form of substitution matrices of PBs. We show that certain conformational states are preferably replaced by other specific conformational states. Interestingly, these structural variations are highly similar to those previously observed across structures of homologous proteins (r² = 0.923) or across the ensemble of conformations from NMR data (r² = 0.919). Thus our study quantitatively shows that overall trends of structural changes in a given protein are nearly identical to the trends of structural differences that occur in the topologically equivalent positions in homologous proteins. Specific case studies are used to illustrate the nature of these structural variations.

Discrete analyses of protein dynamics

Protein structures are highly dynamic macromolecules. This dynamics is often analysed through experimental and/or computational methods only for an isolated or a limited number of proteins. Here, we explore large-scale protein dynamics simulation to observe dynamics of local protein conformations using different perspectives. We analysed molecular dynamics to investigate protein flexibility locally, using classical approaches such as RMSf, solvent accessibility, but also innovative approaches such as local entropy. First, we focussed on classical secondary structures and analysed specifically how β-strand, β-turns, and bends evolve during molecular simulations. We underlined interesting specific bias between β-turns and bends, which are considered as the same category, while their dynamics show differences. Second, we used a structural alphabet that is able to approximate every part of the protein structures conformations, namely protein blocks (PBs) to analyse (i) how each initial local protein conformations evolve during dynamics and (ii) if some exchange can exist among these PBs. Interestingly, the results are largely complex than simple regular/rigid and coil/flexible exchange. AbbreviationsN_eqnumber of equivalentPBProtein BlocksPDBProtein DataBankRMSfroot mean square fluctuationsCommunicated by Ramaswamy H. Sarma.

DNAproDB: an interactive tool for structural analysis of DNA-protein complexes

Many biological processes are mediated by complex interactions between DNA and proteins. Transcription factors, various polymerases, nucleases and histones recognize and bind DNA with different levels of binding specificity. To understand the physical mechanisms that allow proteins to recognize DNA and achieve their biological functions, it is important to analyze structures of DNA–protein complexes in detail. DNAproDB is a web-based interactive tool designed to help researchers study these complexes. DNAproDB provides an automated structure-processing pipeline that extracts structural features from DNA–protein complexes. The extracted features are organized in structured data files, which are easily parsed with any programming language or viewed in a browser. We processed a large number of DNA–protein complexes retrieved from the Protein Data Bank and created the DNAproDB database to store this data. Users can search the database by combining features of the DNA, protein or DNA–protein interactions at the interface. Additionally, users can upload their own structures for processing privately and securely. DNAproDB provides several interactive and customizable tools for creating visualizations of the DNA–protein interface at different levels of abstraction that can be exported as high quality figures. All functionality is documented and freely accessible at http://dnaprodb.usc.edu.

Chemical Biology Interrogates Protein Arginine Phosphorylation

Protein arginine phosphorylation has emerged as an important regulatory mechanism in Gram-positive bacteria, yet a couple of key mysteries remain. In this issue of Cell Chemical Biology, Fuhrmann et al. (2016) use activity-based probes to demonstrate that the protein arginine phosphatase YwlE is inactivated upon oxidative stress to unleash heat-shock gene expression.

Comparative Activity of Adenosine Deaminase Acting on RNA (ADARs) Isoforms for Correction of Genetic Code in Gene Therapy

Introduction:

Members of the adenosine deaminase acting on RNA (ADAR) family of enzymes consist of double-stranded RNA-binding domains (dsRBDs) and a deaminase domain (DD) that converts adenosine (A) into inosine (I), which acts as guanosine (G) during translation. Using the MS2 system, we engineered the DD of ADAR1 to direct it to a specific target. The aim of this work was to compare the deaminase activities of ADAR1-DD and various isoforms of ADAR2-DD.

Materials and Methods:

We measured the binding affinity of the artificial enzyme system on a Biacore ™ X100. ADARs usually target dsRNA, so we designed a guide RNA complementary to the target RNA, and then fused the guide sequence to the MS2 stem-loop. A mutated amber (TAG) stop codon at 58 amino acid (TGG) of EGFP was targeted. After transfection of these three factors into HEK 293 cells, we observed fluorescence signals of various intensities.

Results:

ADAR2-long without the Alu-cassette yielded a much higher fluorescence signal than ADAR2-long with the Alu-cassette. With another isoform, ADAR2-short, which is 81 bp shorter at the C-terminus, the fluorescence signal was undetectable. A single amino acid substitution of ADAR2-long-DD (E488Q) rendered the enzyme more active than the wild type. The results of fluorescence microscopy suggested that ADAR1-DD is more active than ADAR2-long-DD. Western blots and sequencing confirmed that ADAR1-DD was more active than any other DD.

Conclusion:

This study provides information that should facilitate the rational use of ADAR variants for genetic restoration and treatment of genetic diseases.

Natural Products Containing 'Rare' Organophosphorus Functional Groups

Phosphorous-containing molecules are essential constituents of all living cells. While the phosphate functional group is very common in small molecule natural products, nucleic acids, and as chemical modification in protein and peptides, phosphorous can form P⁻N (phosphoramidate), P⁻S (phosphorothioate), and P⁻C (e.g., phosphonate and phosphinate) linkages. While rare, these moieties play critical roles in many processes and in all forms of life. In this review we thoroughly categorize P⁻N, P⁻S, and P⁻C natural organophosphorus compounds. Information on biological source, biological activity, and biosynthesis is included, if known. This review also summarizes the role of phosphorylation on unusual amino acids in proteins (N- and S-phosphorylation) and reviews the natural phosphorothioate (P⁻S) and phosphoramidate (P⁻N) modifications of DNA and nucleotides with an emphasis on their role in the metabolism of the cell. We challenge the commonly held notion that nonphosphate organophosphorus functional groups are an oddity of biochemistry, with no central role in the metabolism of the cell. We postulate that the extent of utilization of some phosphorus groups by life, especially those containing P⁻N bonds, is likely severely underestimated and has been largely overlooked, mainly due to the technological limitations in their detection and analysis.

Base editing: precision chemistry on the genome and transcriptome of living cells

RNA-guided programmable nucleases from CRISPR systems generate precise breaks in DNA or RNA at specified positions. In cells, this activity can lead to changes in DNA sequence or RNA transcript abundance. Base editing is a newer genome-editing approach that uses components from CRISPR systems together with other enzymes to directly install point mutations into cellular DNA or RNA without making double-stranded DNA breaks. DNA base editors comprise a catalytically disabled nuclease fused to a nucleobase deaminase enzyme and, in some cases, a DNA glycosylase inhibitor. RNA base editors achieve analogous changes using components that target RNA. Base editors directly convert one base or base pair into another, enabling the efficient installation of point mutations in non-dividing cells without generating excess undesired editing by-products. In this Review, we summarize base-editing strategies to generate specific and precise point mutations in genomic DNA and RNA, highlight recent developments that expand the scope, specificity, precision and in vivo delivery of base editors and discuss limitations and future directions of base editing for research and therapeutic applications.

Twist and turn: a revised structural view on the unpaired bubble of class II CPD photolyase in complex with damaged DNA

Cyclobutane pyrimidine dimer (CPD) photolyases harness the energy of blue light to repair UV-induced DNA CPDs. Upon binding, CPD photolyases cause the photodamage to flip out of the duplex DNA and into the catalytic site of the enzyme. This process, called base-flipping, induces a kink in the DNA, as well as an unpaired bubble, which are stabilized by a network of protein-nucleic acid interactions. Previously, several co-crystal structures have been reported in which the binding mode of CPD photolyases has been studied in detail. However, in all cases the internucleoside linkage of the photodamage site was a chemically synthesized formacetal analogue and not the natural phosphodiester. Here, the first crystal structure and conformational analysis via molecular-dynamics simulations of a class II CPD photolyase in complex with photodamaged DNA that contains a natural cyclobutane pyrimidine dimer with an intra-lesion phosphodiester linkage are presented. It is concluded that a highly conserved bubble-intruding region (BIR) mediates stabilization of the open form of CPD DNA when complexed with class II CPD photolyases.

Site-directed RNA editing by adenosine deaminase acting on RNA for correction of the genetic code in gene therapy

Site-directed RNA editing is an important technique for correcting gene sequences and ultimately tuning protein function. In this study, we engineered the deaminase domain of adenosine deaminase acting on RNA (ADAR1) and the MS2 system to target-specific adenosines, with the goal of correcting G-to-A mutations at the RNA level. For this purpose, the ADAR1 deaminase domain was fused downstream of the RNA-binding protein MS2, which has affinity for the MS2 RNA. To direct editing to specific targets, we designed guide RNAs complementary to target RNAs. The guide RNAs directed the ADAR1 deaminase to the desired editing site, where it converted adenosine to inosine. To provide proof of principle, we used an allele of enhanced green fluorescent protein (EGFP) bearing a mutation at the 58th amino acid (TGG), encoding Trp, into an amber (TAG) or ochre (TAA) stop codon. In HEK-293 cells, our system could convert stop codons to read-through codons, thereby turning on fluorescence. We confirmed the specificity of editing at the DNA level by restriction fragment length polymorphism analysis and sequencing, and at the protein level by western blotting. The editing efficiency of this enzyme system was ~5%. We believe that this system could be used to treat genetic diseases resulting from G-to-A point mutations.

Knowledge-based prediction of protein backbone conformation using a structural alphabet

Libraries of structural prototypes that abstract protein local structures are known as structural alphabets and have proven to be very useful in various aspects of protein structure analyses and predictions. One such library, Protein Blocks, is composed of 16 standard 5-residues long structural prototypes. This form of analyzing proteins involves drafting its structure as a string of Protein Blocks. Predicting the local structure of a protein in terms of protein blocks is the general objective of this work. A new approach, PB-kPRED is proposed towards this aim. It involves (i) organizing the structural knowledge in the form of a database of pentapeptide fragments extracted from all protein structures in the PDB and (ii) applying a knowledge-based algorithm that does not rely on any secondary structure predictions and/or sequence alignment profiles, to scan this database and predict most probable backbone conformations for the protein local structures. Though PB-kPRED uses the structural information from homologues in preference, if available. The predictions were evaluated rigorously on 15,544 query proteins representing a non-redundant subset of the PDB filtered at 30% sequence identity cut-off. We have shown that the kPRED method was able to achieve mean accuracies ranging from 40.8% to 66.3% depending on the availability of homologues. The impact of the different strategies for scanning the database on the prediction was evaluated and is discussed. Our results highlight the usefulness of the method in the context of proteins without any known structural homologues. A scoring function that gives a good estimate of the accuracy of prediction was further developed. This score estimates very well the accuracy of the algorithm (R2 of 0.82). An online version of the tool is provided freely for non-commercial usage at http://www.bo-protscience.fr/kpred/.

A DNA Structural Alphabet Distinguishes Structural Features of DNA Bound to Regulatory Proteins and in the Nucleosome Core Particle

We analyzed the structural behavior of DNA complexed with regulatory proteins and the nucleosome core particle (NCP). The three-dimensional structures of almost 25 thousand dinucleotide steps from more than 500 sequentially non-redundant crystal structures were classified by using DNA structural alphabet CANA (Conformational Alphabet of Nucleic Acids) and associations between ten CANA letters and sixteen dinucleotide sequences were investigated. The associations showed features discriminating between specific and non-specific binding of DNA to proteins. Important is the specific role of two DNA structural forms, A-DNA, and BII-DNA, represented by the CANA letters AAA and BB2: AAA structures are avoided in non-specific NCP complexes, where the wrapping of the DNA duplex is explained by the periodic occurrence of BB2 every 10.3 steps. In both regulatory and NCP complexes, the extent of bending of the DNA local helical axis does not influence proportional representation of the CANA alphabet letters, namely the relative incidences of AAA and BB2 remain constant in bent and straight duplexes.

A structural bioinformatics investigation on protein–DNA complexes delineates their modes of interaction

A non-redundant dataset of 629 protein–DNA complexes has been used to investigate on amino acid composition of protein-DNA interfaces. Structural proteins, transcription factors and DNA-related enzymes show specific patterns accounting for different modes of their interaction with DNA.

Implication of Terminal Residues at Protein-Protein and Protein-DNA Interfaces

Terminal residues of protein chains are charged and more flexible than other residues since they are constrained only on one side. Do they play a particular role in protein-protein and protein-DNA interfaces? To answer this question, we considered large sets of non-redundant protein-protein and protein-DNA complexes and analyzed the status of terminal residues and their involvement in interfaces. In protein-protein complexes, we found that more than half of terminal residues (62%) are either modified by attachment of a tag peptide (10%) or have missing coordinates in the analyzed structures (52%). Terminal residues are almost exclusively located at the surface of proteins (94%). Contrary to charged residues, they are not over or under-represented in protein-protein interfaces, but strongly prefer the peripheral region of interfaces when present at the interface (83% of terminal residues). The almost exclusive location of terminal residues at the surface of the proteins or in the rim regions of interfaces explains that experimental methods relying on tail hybridization can be successfully applied without disrupting the complexes under study. Concerning conformational rearrangement in protein-protein complexes, despite their expected flexibility, terminal residues adopt similar locations between the free and bound forms of the docking benchmark. In protein-DNA complexes, N-terminal residues are twice more frequent than C-terminal residues at interfaces. Both N-terminal and C-terminal residues are under-represented in interfaces, in contrast to positively charged residues, which are strongly favored. When located in protein-DNA interfaces, terminal residues prefer the periphery. N-terminal and C-terminal residues thus have particular properties with regard to interfaces, which cannot be reduced to their charged nature.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

DNA Backbone BI/BII Distribution and Dynamics in E2 Protein-Bound Environment Determined by Molecular Dynamics Simulations

BI and BII conformational substates in the DNA backbone typify canonical B-form DNA. The BI and BII substates are important for structural variation of DNA and have been implicated in protein-nucleic acid recognition mechanisms. Recent refinements have been made to nucleic acid force fields employed in molecular dynamics simulations that demonstrate a better ability to model the BI and BII states, leading to overall improved modeling of DNA structure and dynamics. These force field improvements have yet to be significantly demonstrated in the context of a protein-DNA system extended to long time scales. Our plan was to run molecular dynamics simulations of a well-studied protein-DNA system (E2-DNA) into the microsecond time scale and determine the ability of the force field to populate BII states in the DNA backbone consistent with dinucleotide steps crystallized in the BII conformation. The results showed that the dinucleotide steps in the E2-DNA complex with the highest BII populations from simulation trajectories corresponded to the dinucleotide steps crystallized in the BII state and that decoy BI and BII states converge to the same results within approximately one microsecond.

The dynamic epitranscriptome: A to I editing modulates genetic information

Adenosine to inosine editing (A to I editing) is a cotranscriptional process that contributes to transcriptome complexity by deamination of adenosines to inosines. Initially, the impact of A to I editing has been described for coding targets in the nervous system. Here, A to I editing leads to recoding and changes of single amino acids since inosine is normally interpreted as guanosine by cellular machines. However, more recently, new roles for A to I editing have emerged: Editing was shown to influence splicing and is found massively in Alu elements. Moreover, A to I editing is required to modulate innate immunity. We summarize the multiple ways in which A to I editing generates transcriptome variability and highlight recent findings in the field.

RPI-Pred: predicting ncRNA-protein interaction using sequence and structural information

RNA-protein complexes are essential in mediating important fundamental cellular processes, such as transport and localization. In particular, ncRNA-protein interactions play an important role in post-transcriptional gene regulation like mRNA localization, mRNA stabilization, poly-adenylation, splicing and translation. The experimental methods to solve RNA-protein interaction prediction problem remain expensive and time-consuming. Here, we present the RPI-Pred (RNA-protein interaction predictor), a new support-vector machine-based method, to predict protein-RNA interaction pairs, based on both the sequences and structures. The results show that RPI-Pred can correctly predict RNA-protein interaction pairs with ∼94% prediction accuracy when using sequence and experimentally determined protein and RNA structures, and with ∼83% when using sequences and predicted protein and RNA structures. Further, our proposed method RPI-Pred was superior to other existing ones by predicting more experimentally validated ncRNA-protein interaction pairs from different organisms. Motivated by the improved performance of RPI-Pred, we further applied our method for reliable construction of ncRNA-protein interaction networks. The RPI-Pred is publicly available at: http://ctsb.is.wfubmc.edu/projects/rpi-pred.

Local dynamics of proteins and DNA evaluated from crystallographic B factors

The dynamics of protein and nucleic acid structures is as important as their average static picture. The local molecular dynamics concealed in diffraction images is expressed as so-called B factors. To find out how the crystal-derived B factors represent the dynamic behaviour of atoms and residues of proteins and DNA in their complexes, the distributions of scaled B factors from a carefully curated data set of over 700 protein-DNA crystal structures were analyzed [Schneider et al. (2014), Nucleic Acids Res. 42, 3381-3394]. Amino acids and nucleotides were categorized based on their molecular neighbourhood as solvent-accessible, solvent-inaccessible (i.e. forming the protein core) or lying at protein-protein or protein-DNA interfaces; the backbone and side-chain atoms were analyzed separately. The B factors of two types of crystal-ordered water molecules were also analyzed. The analysis confirmed several expected features of protein and DNA dynamics, but also revealed surprising facts. Solvent-accessible amino acids have B factors that are larger than those of residues at the biomolecular interfaces, and core-forming amino acids are the most restricted in their movement. A unique feature of the latter group is that their side-chain and backbone atoms are restricted in their movement to the same extent; in all other amino-acid groups the side chains are more floppy than the backbone. The low values of the B factors of water molecules bridging proteins with DNA and the very large fluctuations of DNA phosphates are surprising. The features discriminating different types of residues are less pronounced in structures with lower crystallographic resolution. Some of the observed trends are likely to be the consequence of improper refinement protocols that may need to be rectified.