Association of protein-DNA recognition complexes: electrostatic and nonelectrostatic effects

PNBACE: an ensemble algorithm to predict the effects of mutations on protein-nucleic acid binding affinity

BACKGROUND

Mutations occurring in nucleic acids or proteins may affect the binding affinities of protein-nucleic acid interactions. Although many efforts have been devoted to the impact of protein mutations, few computational studies have addressed the effect of nucleic acid mutations and explored whether the identical methodology could be applied to the prediction of binding affinity changes caused by these two mutation types.

RESULTS

Here, we developed a generalized algorithm named PNBACE for both DNA and protein mutations. We first demonstrated that DNA mutations could induce varying degrees of changes in binding affinity from multiple perspectives. We then designed a group of energy-based topological features based on different energy networks, which were combined with our previous partition-based energy features to construct individual prediction models through feature selections. Furthermore, we created an ensemble model by integrating the outputs of individual models using a differential evolution algorithm. In addition to predicting the impact of single-point mutations, PNBACE could predict the influence of multiple-point mutations and identify mutations significantly reducing binding affinities. Extensive comparisons indicated that PNBACE largely performed better than existing methods on both regression and classification tasks.

CONCLUSIONS

PNBACE is an effective method for estimating the binding affinity changes of protein-nucleic acid complexes induced by DNA or protein mutations, therefore improving our understanding of the interactions between proteins and DNA/RNA.

Physicochemical property based computational scheme for classifying DNA sequence elements of Saccharomyces cerevisiae

GenerationE of huge "omics" data necessitates the development and application of computational methods to annotate the data in terms of biological features. In the context of DNA sequence, it is important to unravel the hidden physicochemical signatures. For this purpose, we have considered various sequence elements such as promoter, ACS, LTRs, telomere, and retrotransposon of the model organism Saccharomyces cerevisiae. Contributions due to di-nucleotides play a major role in studying the DNA conformation profile. The physicochemical parameters used are hydrogen bonding energy, stacking energy and solvation energy per base pair. Our computational study shows that all sequence elements in this study have distinctive physicochemical signatures and the same can be exploited for prediction experiments. The order that we see in a DNA sequence is dictated by biological regions and hence, there exists role of dependency in the sequence makeup, keeping this in mind we are proposing two computational schemes (a) using a windowing block size procedure and (b) using di-nucleotide transitions. We obtained better discriminating profile when we analyzed the sequence data in windowing manner. In the second novel approach, we introduced the di-nucleotide transition probability matrix (DTPM) to study the hidden layer of information embedded in the sequences. DTPM has been used as weights for scanning and predictions. This proposed computational scheme incorporates the memory property which is more realistic to study the physicochemical properties embedded in DNA sequences. Our analysis shows that the DTPM scheme performs better than the existing method in this applied region. Characterization of these elements will be a key to genome editing applications and advanced machine learning approaches may also require such distinctive profiles as useful input features.

Proteins Recognizing DNA: Structural Uniqueness and Versatility of DNA-Binding Domains in Stem Cell Transcription Factors

Proteins in the form of transcription factors (TFs) bind to specific DNA sites that regulate cell growth, differentiation, and cell development. The interactions between proteins and DNA are important toward maintaining and expressing genetic information. Without knowing TFs structures and DNA-binding properties, it is difficult to completely understand the mechanisms by which genetic information is transferred between DNA and proteins. The increasing availability of structural data on protein-DNA complexes and recognition mechanisms provides deeper insights into the nature of protein-DNA interactions and therefore, allows their manipulation. TFs utilize different mechanisms to recognize their cognate DNA (direct and indirect readouts). In this review, we focus on these recognition mechanisms as well as on the analysis of the DNA-binding domains of stem cell TFs, discussing the relative role of various amino acids toward facilitating such interactions. Unveiling such mechanisms will improve our understanding of the molecular pathways through which TFs are involved in repressing and activating gene expression.

Electrostatic hot spot on DNA-binding domains mediates phosphate desolvation and the pre-organization of specificity determinant side chains

A major obstacle towards elucidating the molecular basis of transcriptional regulation is the lack of a detailed understanding of the interplay between non-specific and specific protein–DNA interactions. Based on molecular dynamics simulations of C₂H₂ zinc fingers (ZFs) and engrailed homeodomain transcription factors (TFs), we show that each of the studied DNA-binding domains has a set of highly constrained side chains in preset configurations ready to form hydrogen bonds with the DNA backbone. Interestingly, those domains that bury their recognition helix into the major groove are found to have an electrostatic hot spot for Cl⁻ ions located on the same binding cavity as the most buried DNA phosphate. The spot is characterized by three protein hydrogen bond donors, often including two basic side chains. If bound, Cl⁻ ions, likely mimicking phosphates, steer side chains that end up forming specific contacts with bases into bound-like conformations. These findings are consistent with a multi-step DNA-binding mechanism in which a pre-organized set of TF side chains assist in the desolvation of phosphates into well defined sites, prompting the re-organization of specificity determining side chains into conformations suitable for the recognition of their cognate sequence.

CHARMM: the biomolecular simulation program

CHARMM (Chemistry at HARvard Molecular Mechanics) is a highly versatile and widely used molecular simulation program. It has been developed over the last three decades with a primary focus on molecules of biological interest, including proteins, peptides, lipids, nucleic acids, carbohydrates, and small molecule ligands, as they occur in solution, crystals, and membrane environments. For the study of such systems, the program provides a large suite of computational tools that include numerous conformational and path sampling methods, free energy estimators, molecular minimization, dynamics, and analysis techniques, and model-building capabilities. The CHARMM program is applicable to problems involving a much broader class of many-particle systems. Calculations with CHARMM can be performed using a number of different energy functions and models, from mixed quantum mechanical-molecular mechanical force fields, to all-atom classical potential energy functions with explicit solvent and various boundary conditions, to implicit solvent and membrane models. The program has been ported to numerous platforms in both serial and parallel architectures. This article provides an overview of the program as it exists today with an emphasis on developments since the publication of the original CHARMM article in 1983.

The human cytomegalovirus UL44 C clamp wraps around DNA

Processivity factors tether the catalytic subunits of DNA polymerases to DNA so that continuous synthesis of long DNA strands is possible. The human cytomegalovirus DNA polymerase subunit UL44 forms a C clamp-shaped dimer intermediate in structure between monomeric herpes simplex virus UL42, which binds DNA directly via a basic surface, and the trimeric sliding clamp PCNA, which encircles DNA. To investigate how UL44 interacts with DNA, calculations were performed in which a 12 bp DNA oligonucleotide was docked to UL44. The calculations suggested that UL44 encircles DNA, which interacts with basic residues both within the cavity of the C clamp and in flexible loops of UL44 that complete the "circle." The results of mutational and crosslinking studies were consistent with this model. Thus, UL44 is a "hybrid" of UL42 and PCNA: its structure is intermediate between the two and its mode of interaction with DNA has elements of both.

GBr6NL: a generalized Born method for accurately reproducing solvation energy of the nonlinear Poisson-Boltzmann equation

The nonlinear Poisson-Boltzmann (NLPB) equation can provide accurate modeling of electrostatic effects for nucleic acids and highly charged proteins. Generalized Born methods have been developed to mimic the linearized Poisson-Boltzmann (LPB) equation at substantially reduced cost. The computer time for solving the NLPB equation is approximately fivefold longer than for the LPB equation, thus presenting an even greater obstacle. Here we present the first generalized Born method, GBr(6)NL, for mimicking the NLPB equation. GBr(6)NL is adapted from GBr(6), a generalized Born method recently developed to reproduce the solvation energy of the LPB equation [Tjong and Zhou, J. Phys. Chem. B 111, 3055 (2007)]. Salt effects predicted by GBr(6)NL on 55 proteins overall deviate from NLPB counterparts by 0.5 kcal/mol from ionic strengths from 10 to 1000 mM, which is approximately 10% of the average magnitudes of the salt effects. GBr(6)NL predictions for the salts effects on the electrostatic interaction energies of two protein:RNA complexes are very promising.

Do electrostatic interactions destabilize protein-nucleic acid binding?

The negatively charged phosphates of nucleic acids are often paired with positively charged residues upon binding proteins. It was thus counter-intuitive when previous Poisson-Boltzmann (PB) calculations gave positive energies from electrostatic interactions, meaning that they destabilize protein-nucleic acid binding. Our own PB calculations on protein-protein binding have shown that the sign and the magnitude of the electrostatic component are sensitive to the specification of the dielectric boundary in PB calculations. A popular choice for the boundary between the solute low dielectric and the solvent high dielectric is the molecular surface; an alternative is the van der Waals (vdW) surface. In line with results for protein-protein binding, in this article, we found that PB calculations with the molecular surface gave positive electrostatic interaction energies for two protein-RNA complexes, but the signs are reversed when the vdW surface was used. Therefore, whether destabilizing or stabilizing effects are predicted depends on the choice of the dielectric boundary. The two calculation protocols, however, yielded similar salt effects on the binding affinity. Effects of charge mutations differentiated the two calculation protocols; PB calculations with the vdW surface had smaller deviations overall from experimental data.

Role of hydrogen bonds in protein-DNA recognition: effect of nonplanar amino groups

Amino groups are one of the various types of hydrogen bond donors, abundantly found in protein main chains, protein side chains, and DNA bases. The polar hydrogen atoms of these groups exhibit short ranged, specific, and directional hydrogen bonds, which play a decisive role in the specificity and stability of protein-DNA complexes. To date, planar amino groups are only considered for the analysis of protein-DNA interfacial hydrogen bonds. This assumption regarding hydrogen atom positions possibly failed to establish the expected role of hydrogen bonds in protein-DNA recognition. We have performed ab initio quantum chemical studies on amino acid side chains and DNA bases containing amino groups as well as on specific hydrogen bonded residue pairs selected from high-resolution protein-DNA complex crystal structures. Our results suggest that occurrences of pyramidal amino groups are more probable in comparison with the usually adopted planar geometry. This increases the quality of the existing hydrogen bonds in almost all cases. Further, detailed analysis of protein-DNA interfacial hydrogen bonds in 107 crystal structures using the in-house program "pyrHBfind" indicates that consideration of energetically more preferred nonplanar amino groups improves the geometry of hydrogen bonds and also gives rise to new contacts amounting to nearly 14.5% of the existing interactions. Large improvements have been observed specifically for the amino groups of guanine, which faces the DNA minor groove and thus helps to resolve the problem of insufficient directional contacts observed in many minor groove binding complexes. Apart from guanine, improvement observed for asparagine, glutamine, adenine, or cytosine also indicates that the consideration of nonplanar amino groups leads to a more realistic scenario of hydrogen bonds occurring between protein and DNA residues.

Structural Studies on Pax-8 Prd Domain/DNA Complex

Pax-8 is a member of the Pax family of transcription factors and is essential in the development of thyroid follicular cells. Pax-8 has two DNA-binding domains: the paired domain and the homeo domain. In this study, a preliminary X-ray diffraction analysis of the mammalian Pax-8 paired domain in complex with the C-site of the thyroglobulin promoter was achieved. The Pax-8 paired domain was crystallized by the hanging-drop vapor-diffusion method in complex with both a blunt-ended 26 bp DNA fragment and with a sticky-ended 24 bp DNA fragment with two additional overhanging bases. Crystallization experiments make clear that the growth of transparent crystals with large dimensions and regular shape is particularly influenced by ionic strength. The crystals of Pax-8 complex with blunt-ended and sticky-ended DNA, diffracted synchrotron radiation to 6.0 and 8.0 A resolution and belongs both to the C centered monoclinic system with cell dimensions: a = 89.88 A, b = 80.05 A, c = 67.73 A, and beta = 124.3 degrees and a = 256.56, b = 69.07, c = 99.32 A, and beta = 98.1 degrees , respectively. Fluorescence experiments suggest that the crystalline disorder, deduced by the poor diffraction, can be attributed to the low homogeneity of the protein-DNA sample. The theoretical comparative model of the Pax-8 paired domain complexed with the C-site of the thyroglobulin promoter shows the probable presence of some specific protein-DNA interactions already observed in other Pax proteins and the important role of the cysteine residues of PAI subdomain in the redox control of the DNA recognition.

Critical jump sizes in DNA-protein interactions

Interaction of a protein molecule with a specific-site on the DNA lattice can be modeled as an unbiased random jump process. Here we show that there exists a critical jump size (kc) beyond which site-specific association of a protein molecule with a DNA lattice cannot be facilitated. The maximum achievable association rate is predicted to be approximately 10(10) mol-1 s-1. This critical jump size scales with the total length of DNA lattice (N) as kc proportional, variantN2/3. Beyond kc the mean first passage time MFPT (denoted as T) required for the protein molecule to target the specific-site follows a linear scaling law as T proportional, variantN rather than the usual T proportional, variantN2 scaling law. On the basis of these results we argue that the evolution of the super coiled structures of the genomic DNA must be a consequence of the existence of this critical jump sizes. We finally show that the random jump method of searching for the specific-site by the protein molecule on the DNA lattice itself introduce an abstract linear type potential favoring the site-specific association rate.

Thermodynamic and electrostatic properties of ternary Oxytricha nova TEBP-DNA complex

Telomeres constitute the nucleoprotein ends of eukaryotic chromosomes which are essential for their proper function. Telomere end binding protein (TEBP) from Oxytricha nova was among the first telomeric proteins, which were well characterized biologically. TEBP consists of two protein subunits (alpha, beta) and forms a ternary complex with single stranded telomeric DNA containing tandem repeats TTTTGGGG. This work presents the characterization of the thermodynamic and electrostatic properties of this complex by computational chemistry methods (continuum Poisson-Boltzmann and solvent accessible surface calculations). Our calculations give a new insight into molecular properties of studied system. Based on the thermodynamic analysis we provide a rationale for the experimental observation that alpha and ssDNA forms a binary complex and the beta subunit joins alpha:ssDNA complex only after the latter is formed. Calculations of distribution of the molecular electrostatic potential for protein subunits alone and for all possible binary complexes revealed the important role of the "guiding funnel" potential generated by alpha:ssDNA complex. This potential may help the beta subunit to dock to the already formed alpha:DNA intermediate in highly steric and electrostatic favorable manner. Our pK(a) calculations of TEBP are able to explain the experimental mobility shifts of the complex in electrophoretic non-denaturating gels.

Solvation change and ion release during aminoacylation by aminoacyl-tRNA synthetases

Discrimination between cognate and non-cognate tRNAs by aminoacyl-tRNA synthetases occurs at several steps of the aminoacylation pathway. We have measured changes of solvation and counter-ion distribution at various steps of the aminoacylation pathway of glutamyl- and glutaminyl-tRNA synthetases. The decrease in the association constant with increasing KCl concentration is relatively small for cognate tRNA binding when compared to known DNA-protein interactions. The electro-neutral nature of the tRNA binding domain may be largely responsible for this low ion release stoichiometry, suggesting that a relatively large electrostatic component of the DNA-protein interaction free energy may have evolved for other purposes, such as, target search. Little change in solvation upon tRNA binding is seen. Non-cognate tRNA binding actually increases with increasing KCl concentration indicating that charge repulsion may be a significant component of binding free energy. Thus, electrostatic interactions may have been used to discriminate between cognate and non-cognate tRNAs in the binding step. The catalytic constant of glutaminyl-tRNA synthetase increases with increasing osmotic pressure indicating a water release of 8.4 +/- 1.4 mol/mol in the transition state, whereas little change is seen in the case of glutamyl-tRNA synthetase. We propose that the significant amount of water release in the transition state, in the case of glutaminyl-tRNA synthetase, is due to additional contact of the protein with the tRNA in the transition state.

DNA-protein interactions under random jump conditions

We model the site-specific association of a protein molecule with DNA as a random walk with random jumps. Results show that the simultaneous occurrence of processes such as sliding, hopping, and intersegmental transfer can facilitate the diffusion-controlled site-specific association rate. We have also shown that sliding would dominate at lower DNA length, whereas at higher lengths hopping and intersegmental transfer would dominate. Apart from this, we predict that the association rate is directly proportional to the size of nonspecific DNA that flanks the specific site. These results are consistent with the experimental observations.