Computational investigation of proton transfer,pKa shifts and pH‐optimum of protein–DNA and protein–RNA complexes

Protein–nucleic acid interactions play a crucial role in many biological processes. This work investigates the changes of pKa values and protonation states of ionizable groups (including nucleic acid bases) that may occur at protein–nucleic acid binding. Taking advantage of the recently developed pKa calculation tool DelphiPka, we utilize the large protein–nucleic acid interaction database (NPIDB database) to model pKa shifts caused by binding. It has been found that the protein's interfacial basic residues experience favorable electrostatic interactions while the protein acidic residues undergo proton uptake to reduce the energy cost upon the binding. This is in contrast with observations made for protein–protein complexes. In terms of DNA/RNA, both base groups and phosphate groups of nucleotides are found to participate in binding. Some DNA/RNA bases undergo pKa shifts at complex formation, with the binding process tending to suppress charged states of nucleic acid bases. In addition, a weak correlation is found between the pH‐optimum of protein–DNA/RNA binding free energy and the pH‐optimum of protein folding free energy. Overall, the pH‐dependence of protein–nucleic acid binding is not predicted to be as significant as that of protein–protein association. Proteins 2017; 85:282–295. © 2016 Wiley Periodicals, Inc.     

Poisson-Boltzmann analysis of the lambda repressor-operator interaction.

A theoretical study of the ion atmosphere contribution to the binding free energy of the lambda repressor-operator complex is presented. The finite-difference form of the Poisson-Boltzmann equation was solved to calculate the electrostatic interaction energy of the amino-terminal domain of the lambda repressor with a 9 or 45 base pair oligonucleotide. Calculations were performed at various distances between repressor and operator as well as at different salt concentrations to determine ion atmosphere contributions to the total electrostatic interaction. Details in the distribution of charges on DNA and protein atoms had a strong influence on the calculated total interaction energies. In contrast, the calculated salt contributions are relatively insensitive to changes in the details of the charge distribution. The results indicate that the ion atmosphere contribution favors association at all protein-DNA distances studied. The theoretical number of ions released upon repressor-operator binding appears to be in reasonable agreement with experimental data.     

EF loop conformational change triggers ligand binding in beta-lactoglobulins

Beta-lactoglobulins, belonging to the lipocalin family, are a widely studied group of proteins, characterized by the ability to solubilize and transport hydrophobic ligands, especially fatty acids. Despite many reports, the mechanism of ligand binding and the functional role of these proteins is still unclear, and many contradicting concepts are often encountered in the literature. In the present paper the comparative analysis of the binding properties of beta-lactoglobulins has been performed using sequence-derived information, structure-based electrostatic calculations, docking simulations, and NMR experiments. Our results reveal for the first time the mechanism of beta-lactoglobulin ligand binding, which is completely determined by the opening-closing of EF loop, triggered by Glu89 protonation. The alkaline shift observed for Glu89 pKa in porcine beta-lactoglobulin (pKa 9.7) with respect to the bovine species (pKa 5.5) depends upon the interplay of electrostatic effects of few nearby key residues. Porcine protein is therefore able to bind fatty acids provided that the appropriate pH solution conditions are met (pH > 8.6), where the EF loop conformational change can take place. The unusually high pH of binding detected for porcine beta-lactoglobulin seems to be functional to lipases activity. Theoretical pKa calculations extended to representative beta-lactoglobulins allowed the identification of key residues involved in structurally and functionally important electrostatic interactions. The results presented here provide a strong indication that the described conformational change is a common feature of all beta-lactoglobulins.     

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.     

Predicting DNA-binding amino acid residues from electrostatic stabilization upon mutation to Asp/Glu and evolutionary conservation

Binding of polyanionic DNA depends on the cluster of electropositive atoms in the binding site of a DNA-binding protein. Such a cluster of electropositive protein atoms would be electrostatically unfavorable without stabilizing interactions from the respective electronegative DNA atoms and would likely be evolutionary conserved due to its critical biological role. Consequently, our strategy for predicting DNA-binding residues is based on detecting a cluster of evolutionary conserved surface residues that are electrostatically stabilized upon mutation to negatively charged Asp/Glu residues. The method requires as input the protein structure and sufficient sequence homologs to define each residue's relative conservation, and it yields as output experimentally testable residues that are predicted to bind DNA. By incorporating characteristic DNA-binding site features (i.e., electrostatic strain and amino acid conservation), the new method yields a prediction accuracy of 83%, which is much higher than methods based on only electrostatic strain (57%) or conservation alone (50%). It is also less sensitive to protein conformational changes upon DNA binding than methods that mainly depend on the 3D protein structure.     

Density functional theory for the nonspecific binding of salt to polyelectrolytes: thermodynamic properties

The thermodynamics of the nonspecific binding of salt to a polyelectrolyte molecule is studied using a density functional approach. The polyelectrolyte molecule is modeled as an infinite, inflexible, and impenetrable charged cylinder and the counterions and co-ions are modeled as charged hard spheres of equal diameter. The density functional theory is based on a hybrid approach where the hard-sphere contribution to the one-particle correlation function is evaluated nonperturbatively and the ionic contribution to the one-particle correlation function is evaluated perturbatively. The advantage of the approach is that analytical expressions are available for all the correlation functions. The calculated single ion preferential interaction coefficients, excess free energy, and activity coefficients show a nonmonotonic variation as a function of polyion charge in the presence of divalent ions. These properties display considerable departure from the predictions of the nonlinear Poisson-Boltzmann (NLPB) equation, with qualitative differences in some cases, which may be attributed to correlation effects neglected in the NLPB theory.     

Interaction of nucleic acids with carbon nanotubes and dendrimers

Nucleic acid interaction with nanoscale objects like carbon nanotubes (CNTs) and dendrimers is of fundamental interest because of their potential application in CNT separation, gene therapy and antisense therapy. Combining nucleic acids with CNTs and dendrimers also opens the door towards controllable self-assembly to generate various supra-molecular and nano-structures with desired morphologies. The interaction between these nanoscale objects also serve as a model system for studying DNA compaction, which is a fundamental process in chromatin organization. By using fully atomistic simulations, here we report various aspects of the interactions and binding modes of DNA and small interfering RNA (siRNA) with CNTs, graphene and dendrimers. Our results give a microscopic picture and mechanism of the adsorption of single- and double-strand DNA (ssDNA and dsDNA) on CNT and graphene. The nucleic acid-CNT interaction is dominated by the dispersive van der Waals (vdW) interaction. In contrast, the complexation of DNA (both ssDNA and dsDNA) and siRNA with various generations of poly-amido-amine (PAMAM) dendrimers is governed by electrostatic interactions. Our results reveal that both the DNA and siRNA form stable complex with the PAMAM dendrimer at a physiological pH when the dendrimer is positively charged due to the protonation of the primary amines. The size and binding energy of the complex increase with increase in dendrimer generation. We also give a summary of the current status in these fields and discuss future prospects.     

Calculations of free-energy contributions to protein-RNA complex stabilization

The problem of calculating binding affinities of protein-RNA complexes is addressed by analyzing a computational strategy of modeling electrostatic free energies based on a nonlinear Poisson-Boltzmann (NLPB) model and linear response approximation (LRA). The underlying idea is to treat binding as a two-step process. Solutions to the NLPB equation calculate free energies arising from electronic polarizability and the LRA is constructed from molecular dynamics simulations to model reorganization free energies due to conformational transitions. By implementing a consistency condition of requiring the NLPB model to reproduce the solute-solvent free-energy transitions determined by the LRA, a "macromolecule dielectric constant" (epsilon(m)) for treating reorganization is obtained. The applicability of this hybrid approach was evaluated by calculating the absolute free energy of binding and free-energy changes for amino acid substitutions in the complex between the U1A spliceosomal protein and its cognate RNA hairpin. Depending on the residue substitution, epsilon(m) varied from 3 to 18, and reflected dipolar reorientation not included in the polarization modeled by epsilon(m) = 2. Although the changes in binding affinities from substitutions modeled strictly at the implicit level by the NLPB equation with epsilon(m) = 4 reproduced the experimental values with good overall agreement, substitutions problematic to this simple treatment showed significant improvement when solved by the NLPB-LRA approach.     

Molecular Interactions and Residues Involved in Force Generation in the T4 Viral DNA Packaging Motor

Many viruses utilize molecular motors to package their genomes into preformed capsids. A striking feature of these motors is their ability to generate large forces to drive DNA translocation against entropic, electrostatic, and bending forces resisting DNA confinement. A model based on recently resolved structures of the bacteriophage T4 motor protein gp17 suggests that this motor generates large forces by undergoing a conformational change from an extended to a compact state. This transition is proposed to be driven by electrostatic interactions between complementarily charged residues across the interface between the N- and C-terminal domains of gp17. Here we use atomistic molecular dynamics simulations to investigate in detail the molecular interactions and residues involved in such a compaction transition of gp17. We find that although electrostatic interactions between charged residues contribute significantly to the overall free energy change of compaction, interactions mediated by the uncharged residues are equally if not more important. We identify five charged residues and six uncharged residues at the interface that play a dominant role in the compaction transition and also reveal salt bridging, van der Waals, and solvent hydrogen-bonding interactions mediated by these residues in stabilizing the compact form of gp17. The formation of a salt bridge between Glu309 and Arg494 is found to be particularly crucial, consistent with experiments showing complete abrogation in packaging upon Glu309Lys mutation. The computed contributions of several other residues are also found to correlate well with single-molecule measurements of impairments in DNA translocation activity caused by site-directed mutations.     

Electrostatic contribution to the energy of binding of aromatic ligands with DNA

The method of solution of the nonlinear Poisson-Boltzmann equation was used to calculate electrostatic energy of binding of various aromatic ligands with DNA oligomers of different length. Analysis of the electrostatic contribution was made in terms of a two-step DNA binding process: formation of the intercalation cavity and insertion of the ligand. The total electrostatic energy was also partitioned into components: the energy of atom-atom coulombic interactions and the energy of interaction with surrounding water. The results indicate that electrostatic interactions are, as a whole, unfavorable to the intercalation process and that a correct analysis of structure-energy interrelation for Ligand-DNA interactions should only be accomplished at the level of the components rather than at the level of total electrostatic energy.     

Bindings of RNAse to denatured and native deoxyribonucleic acids

The binding of pancreatic ribonuclease-A by denatured DNA, native DNA, poly-dA, and poly-dT, has been studied by a gel filtration method. With denatured DNA at pH 7.5, ionic strength 0.053M, there is one binding site per 12 nucleotides and the equilibrium binding constant per site is 9.7 × 10⁴ l./mole. The binding constant increases by a factor of 8 as the pH is decreased from 8 to 7. The strength of the binding of denatured DNA increases with decreasing ionic strength. At pH 7.5, native DNA binds about ? as strongly as does denatured DNA. The binding affinity increases in the order poly-dA, denatured DNA, and poly-dT. These results support the view that the binding of denatured DNA involves both electrostatic interactions between the negatively charged polynucleotide and the positively charged protein, and an interaction of the protein with a pyrimidine residue of the denatured DNA, and thus that the binding is basically similar to that between RNAse and its substrate RNA.     

Investigations of Sso7d catalytic residues by NMR titration shifts and electrostatic calculations

Sso7d is a small basic protein consisting of 62 amino acids isolated from the thermoacidophilic archeobacterium Sulfolobus solfataricus. The protein is endowed with DNA binding properties, RNase activity, and the capability of rescuing aggregated proteins in the presence of ATP. In this study, the electrostatic properties of Sso7d are investigated by using the Poisson-Boltzmann calculation of the surface potential distribution and following by NMR spectroscopy the proton chemical shift pH titration of acidic residues. Although the details of the catalytic mechanism still have to be defined, the results from NMR experiments confirm the possible involvement of Glu35 as the proton acceptor in the catalytic reaction, as seen by its abnormally high pK(a) value. Poisson-Boltzmann calculations and NMR titration shifts suggest the presence of a possible hydrogen bond between Glu35 and Tyr33, with a consequent rather rigid arrangement at these positions. Comparison with RNase T1 suggests that Tyr7 may be a good candidate for acting as a proton donor in the active site of Sso7d as shown by its low phenolic pK(a) of approximately 9.3. Titration experiments performed with the UpA, a RNA dinucleotide model, showed that the protein residues affected by the interaction are mainly located in a different region with respect to the surface affected by DNA recognition, in good agreement with the surface potential distribution found with electrostatic calculations.     

Peptide nucleic acid (PNA) amphiphiles: synthesis, self-assembly, and duplex stability

Peptide amphiphiles comprising a class of conjugates of peptide nucleic acid (PNA), natural amino acids, and n-alkanes were synthesized and studied. These PNA amphiphiles (PNAA) self-assemble at concentrations between 10 and 50 muM and exhibit water solubilities above 500 muM. The highly specific, stable DNA binding properties of PNAs are preserved by these modifications, with no significant differences between the thermodynamics of DNA binding of the PNA peptide and the PNA amphiphile. Proper solubilization of the PNAA required the attachment of (Lys)(2) and (Glu)(4) peptides to PNAs, which affected the PNAA-DNA duplex stability by electrostatic interactions between these charged amino acids and the negatively charged DNA backbone. These electrostatic effects did not affect the specificity of DNA binding, however. Electrostatic effects are screened with added salt, in a manner consistent with previous studies of PNA-DNA duplex stability and predictions from a charged-cylinder model for the duplex.     

Comparison of calculation and experiment implicates significant electrostatic contributions to the binding stability of barnase and barstar

The contributions of electrostatic interactions to the binding stability of barnase and barstar were studied by the Poisson-Boltzmann model with three different protocols: a), the dielectric boundary specified as the van der Waals (vdW) surface of the protein along with a protein dielectric constant (epsilon (p)) of 4; b), the dielectric boundary specified as the molecular (i.e., solvent-exclusion (SE)) surface along with epsilon (p) = 4; and c), "SE + epsilon (p) = 20." The "vdW + epsilon (p) = 4" and "SE + epsilon (p) = 20" protocols predicted an overall electrostatic stabilization whereas the "SE + epsilon (p) = 4" protocol predicted an overall electrostatic destabilization. The "vdW + epsilon (p) = 4" protocol was most consistent with experiment. It quantitatively reproduced the observed effects of 17 mutations neutralizing charged residues lining the binding interface and the measured coupling energies of six charge pairs across the interface and reasonably rationalized the experimental ionic strength and pH dependences of the binding constant. In contrast, the "SE + epsilon (p) = 4" protocol predicted significantly larger coupling energies of charge pairs whereas the "SE + epsilon (p) = 20" protocol did not predict any pH dependence. This study calls for further scrutiny of the different Poisson-Boltzmann protocols and demonstrates potential danger in drawing conclusions on electrostatic contributions based on a particular calculation protocol.     

Electrostatic interactions in the association of proteins: an analysis of the thrombin-hirudin complex.

The role of electrostatic interactions in stabilization of the thrombin-hirudin complex has been investigated by means of two macroscopic approaches: the modified Tanford-Kirkwood model and the finite-difference method for numerical solution of the Poisson-Boltzmann equations. The electrostatic potentials around the thrombin and hirudin molecules were asymmetric and complementary, and it is suggested that these fields influence the initial orientation in the process of the complex formation. The change of the electrostatic binding energy due to mutation of acidic residues in hirudin has been calculated and compared with experimentally determined changes in binding energy. In general, the change in electrostatic binding energy for a particular mutation calculated by the modified Tanford-Kirkwood approach agreed well with the experimentally observed change. The finite-difference approach tended to overestimate changes in binding energy when the mutated residues were involved in short-range electrostatic interactions. Decreases in binding energy caused by mutations of amino acids that do not make any direct ionic interactions (e.g., Glu 61 and Glu 62 of hirudin) can be explained in terms of the interaction of these charges with the positive electrostatic potential of thrombin. Differences between the calculated and observed changes in binding energy are discussed in terms of the crystal structure of the thrombin-hirudin complex.     

pH-Dependent conformational stability of the ribotoxin alpha-sarcin and four active site charge substitution variants

Alpha-sarcin is an exquisitely specific ribonuclease that binds and cleaves a single phosphodiester bond in the large rRNA of the eukaryotic ribosome, inactivating it. To better understand this remarkable activity, the contributions of the active site residues (His 50, Glu 96, and His 137) to the conformational stability have been determined as a function of pH using variant proteins containing uncharged substitutes. Wild-type alpha-sarcin and the variants are maximally stable near pH 5.5, coinciding with the pH of optimal activity. A comparison of the stability vs pH profiles determined by thermal denaturation experiments to those calculated on the basis of pKa values shows that the charged forms of Glu 96 and His 137 compromise the enzyme's stability, lowering it. In contrast to barnase, there is little evidence for significant electrostatic interactions in the denatured states of alpha-sarcin or its active site variants between pH 3.5 and pH 8.5. Alpha-sarcin contains a long beta-hairpin and surface loops which are highly positively charged and which play key roles in membrane translocation and in ribosome binding. These positive charges decrease the stability of alpha-sarcin, particularly below pH 5. Hydrogen exchange measurements have been performed at pH 5.5 and reveal that the catalytic residues are firmly anchored in highly stable elements of secondary structure. Significant, though lower, levels of protection are observed for many amide protons in the positively charged beta-hairpin and long loops.