3D-QSAR study for DNA cleavage proteins with a potential anti-tumor ATCUN-like motif

Genomics projects have elucidated several genes that encode protein sequences. Subsequently, the advent of the proteomics age has enabled the synthesis and 3D structure determination for these protein sequences. Some of these proteins incorporate metal atoms but it is often not known whether they are metal-binding proteins and the nature of the biological activity is not understood. Consequently, the development of methods to predict metal-mediated biological activity of proteins from the 3D structure of metal-unbound proteins is a goal of major importance. More specifically, the amino terminal Cu(II)- and Ni(II)-binding (ATCUN) motif is a small metal-binding site found in the N-terminus of many naturally occurring proteins. The ATCUN motif participates in DNA cleavage and has anti-tumor activity. In this study, we calculated average 3D electrostatic potentials (xi(k)) for 265 different proteins including 133 potential ATCUN anti-tumor proteins. We also calculated xi(k) values for the total protein or for the following specific protein regions: the core, inner, middle, and outer orbits. A linear discriminant analysis model was subsequently developed to assign proteins into two groups called ATCUN DNA-cleavage proteins and non-active proteins. The best model found was: ATCUN=1.15.xi(1)(inner)+2.18.xi(5)(middle)+27.57.xi(0)(outer)-27.57.xi(0)(total)+0.09. The model correctly classified 182 out of 197 (91.4%) and 61 out of 66 (92.4%) proteins in training and external predicting series', respectively. Finally, desirability analysis was used to predict the values for the electrostatic potential in one single region and the combined values in two regions that are desirable for ATCUN-like proteins. To the best of our knowledge, the present work is the first study in which desirability analysis has been used in protein quantitative-structure-activity-relationship (QSAR).     

Effect of mutations involving charged residues on the stability of staphylococcal nuclease: a continuum electrostatics study

A continuum electrostatics model is used to calculate the relative stabilities of 117 mutants of staphylococcal nuclease (SNase) involving the mutation of a charged residue to an uncharged residue. The calculations are based on the crystallographic structure of the wild-type protein and attempt to take implicitly into account the effect of the mutations in the denatured state by assuming a linear relationship between the free energy changes caused by the mutation in the native and denatured states. A good correlation (linear correlation coefficient of approximately 0.8) is found with published experimental relative stabilities of these mutants. The results suggest that in the case of SNase (i) charged residues contribute to the stability of the native state mainly through electrostatic interactions, and (ii) native-like electrostatic interactions may persist in the denatured state. The continuum electrostatics method is only moderately sensitive to model parameters and leads to quasi-predictive results for the relative mutant stabilities (error of 2-3 kJ mol(-1) or of the order of k(B)T), except for mutants in which a charged residue is mutated to glycine.     

Mutation of the pseudo-EF-hand of calbindin D9k into a normal EF-hand. Biophysical studies.

The two Ca(2+)-binding sites in calbindin D9k, a protein belonging to the calmodulin superfamily of intracellular proteins, have slightly different structure. The C-terminal site (amino acids 54-65) is a normal EF-hand as in the other proteins of the calmodulin superfamily, while the N-terminal site (amino acids 14-27) contains two additional amino acids, one of which is a proline. We have constructed and studied five mutants of calbindin D9k modified in the N-terminal site. In normal EF-hand structures the first amino acid to coordinate calcium is invariantly an Asp. For this reason Ala15, is exchanged by an Asp in all mutants and the mutants also contain various other changes in this site. The mutants have been characterized by 43Ca, 113Cd and 1H NMR and by the determination of the calcium binding constants using absorption chelators. In two of the mutants (one where Ala14 is deleted, Ala15 is replaced by Asp and Pro20 is replaced by Gly, the other where, in addition, Asn21 is deleted), we find that the structure has changed considerably compared to the wild-type calbindin. The NMR results indicate that the calcium coordination has changed to mainly side-chain carboxyls, from being octahedrally coordinated by mainly back-bone carbonyls, and/or that the coordination number has decreased. The N-terminal site has thus been turned into a normal EF-hand, in which the calcium ion is coordinated by side-chain carboxyls. Furthermore, the calcium binding constants of these two mutant proteins are almost as high as in the wild-type calbindin D9k. That is, the extensive alterations in the N-terminal site have not disrupted the calcium binding ability of the proteins.     

Effect of amino acid substitutions and deletions on the thermal stability, the pH stability and unfolding by urea of bovine calbindin D9k

The influence of amino acid substitutions and deletions on the stability of bovine calbindin D9k, the smallest protein known with a pair of EF-hand calcium-binding sites, has been studied using circular dichroism and ultraviolet absorption spectroscopy. The five modifications are confined to one of the two Ca2+ -binding sites. The Ca2+-loaded forms of the wild-type and mutant calbindins are too stable to be significantly denatured by heating at 90 degrees C or by adding 8 M urea. For the Ca2+-free (apo) forms thermal unfolding appears to be only half complete at 90 degrees C, while denaturation is complete in 7-8 M urea. Four of the mutant proteins show reduced resistance towards unfolding by urea, but one of the modified proteins (Glu-17----Gln) shows an increased stability, presumably because of a reduced electrostatic repulsion in the native state. According to X-ray crystallographic data the OH group of the single tyrosine of calbindin (Tyr-13) is hydrogen-bonded to the carboxyl group of Glu-35, thus linking the two alpha helices flanking the N-terminal Ca2+ site. The pK of ionization of the Tyr-13 hydroxyl group was over 13 for calcium forms of the wild-type protein, between 12.3 and 12.8 for the calcium form of three mutants and between 11.5 and 11.7 for the apoproteins. Significant differences in pH stability between wild type and mutants were observed in the calcium forms, but were not apparent in the apo forms.     

Conferring thermostability to mesophilic proteins through optimized electrostatic surfaces

Recently, there have been several experimental reports of proteins displaying appreciable stability gains through mutation of one or two amino acid residues. Here, we employ a simple theoretical model to quickly screen mutant structures for increased thermostability through optimization of the protein's electrostatic surface. Our results are able to reproduce the experimental observation that elimination of like-charge repulsions and creation of opposite-charge attractions on the protein surface is an efficient method to confer thermostability to a mesophilic protein. Using Poisson-Boltzmann electrostatics, we calculate relative protein stabilities for the exhaustive surface mutagenesis of the cold shock, RNase T1, and CheY proteins. Comparison with 25 experimentally characterized cold shock protein mutants reveals an average correlation of 0.86. The model is also quantitatively accurate when reproducing the experimental D49A and D49H mutant stabilities of RNase T1. This work represents the first comprehensive in silico screening of mutant candidates likely to confer thermostability to mesophilic proteins through optimization of surface electrostatics. Systematic single mutant, followed by double mutant, screening yields a limited number of mutant structures displaying significant stability gains suitable for subsequent experimental characterization.     

Automated computational framework for the analysis of electrostatic similarities of proteins

Charge plays an important role in protein-protein interactions. In the case of excessively charged proteins, their electrostatic potentials contribute to the processes of recognition and binding with other proteins or ligands. We present an automated computational framework for determining the contribution of each charged amino acid to the electrostatic properties of proteins, at atomic resolution level. This framework involves computational alanine scans, calculation of Poisson-Boltzmann electrostatic potentials, calculation of electrostatic similarity distances (ESDs), hierarchical clustering analysis of ESDs, calculation of solvation free energies of association, and visualization of the spatial distributions of electrostatic potentials. The framework is useful to classify families of mutants with similar electrostatic properties and to compare them with the parent proteins in the complex. The alanine scan mutants introduce perturbations in the local electrostatic properties of the proteins and aim in delineating the contribution of each mutated amino acid in the spatial distribution of electrostatic potential, and in biological function when electrostatics is a dominant contributing factor in protein-protein interactions. The framework can be used to design new proteins with tailored electrostatic properties, such as immune system regulators, inhibitors, and vaccines, and in guiding experimental studies. We present an example for the interaction of the immune system protein C3d (the d-fragment of complement protein C3) with its receptor CR2, and we discuss our data in view of a binding site controversy.     

Association between flower stalk elongation,an Arabidopsis developmental trait,and the subcellular location and movement dynamics of the nonstructural protein P3 of Turnip mosaic virus

Virus infections affect plant developmental traits but this aspect of the interaction has not been extensively studied so far. Two strains of Turnip mosaic virus differentially affect Arabidopsis development, especially flower stalk elongation, which allowed phenotypical, cellular, and molecular characterization of the viral determinant, the P3 protein. Transiently expressed wild-type green fluorescent protein-tagged P3 proteins of both strains and selected mutants of them revealed important differences in their behaviour as endoplasmic reticulum (ER)-associated peripheral proteins flowing along the reticulum, forming punctate accumulations. Three-dimensional (3D) model structures of all expressed P3 proteins were computationally constructed through I-TASSER protein structure predictions, which were used to compute protein surfaces and map electrostatic potentials to characterize the effect of amino acid changes on features related to protein interactions and to phenotypical and subcellular results. The amino acid at position 279 was the main determinant affecting stalk development. It also determined the speed of ER-flow of the expressed proteins and their final location. A marked change in the protein surface electrostatic potential correlated with changes in subcellular location. One single amino acid in the P3 viral protein determines all the analysed differential characteristics between strains differentially affecting flower stalk development. A model proposing a role of the protein in the intracellular movement of the viral replication complex, in association with the viral 6K2 protein, is proposed. The type of association between both viral proteins could differ between the strains.     

Effect of Cleavage Mutants on Syncytium Formation Directed by the Wild-Type Fusion Protein of Newcastle Disease Virus

The effects of Newcastle disease virus (NDV) fusion (F) glycoprotein cleavage mutants on the cleavage and syncytium-forming activity of the wild-type F protein were examined. F protein cleavage mutants were made by altering amino acids in the furin recognition region (amino acids 112 to 116) in the F protein of a virulent strain of NDV. Four mutants were made: Q114P replaced the glutamine residue with proline; K115G replaced lysine with glycine; double mutant K115G, R113G replaced both a lysine and an arginine with glycine residues; and a triple mutant, R112G, K115G, F117L, replaced three amino acids to mimic the sequence found in avirulent strains of NDV. All mutants except Q114P were cleavage negative and fusion negative. However, addition of exogenous trypsin cleaved all mutant F proteins and activated fusion. As expected for an oligomeric protein, the fusion-negative mutants had a dominant negative phenotype: cotransfection of wild-type and mutant F protein cDNAs resulted in an inhibition of syncytium formation. The presence of the mutant F protein did not inhibit cleavage of the wild-type protein. Furthermore, evidence is presented that suggests that the mutant protein and the wild-type protein formed heterooligomers. By measuring the syncytium-forming activity of the wild-type protein at various ratios of expression of mutant and wild-type protein, results were obtained that are most consistent with the notion that the size of the functionally active NDV F protein in these assays is a single oligomer, likely a trimer. That a larger oligomer, containing a mix of both wild-type and mutant F proteins, has partial activity cannot, however, be ruled out.     

Local stability identification and the role of key acidic amino acid residues in staphylococcal nuclease unfolding

Staphylococcal nuclease is a single domain protein with 149 amino acids. It has no disulfide bonds, which makes it a simple model for the study of protein folding. In this study, 20 mutants of this protein were generated each with a single base substitution of glycine for negatively charged glutamic acid or aspartic acid. Using differential scanning microcalorimetry in thermal denaturation experiments, we identified two mutants, E75G and E129G, having approximately 43% and 44%, respectively, lower DeltaH(cal) values than the wild-type protein. Furthermore, two mutants, E75Q and E129Q, were created and the results imply that substitution of the Gly residue has little influence on destabilization of the secondary structure that leads to the large perturbation of the tertiary protein structure stability. Two local stable areas formed by the charge-charge interactions around E75 and E129 with particular positively charged amino acids are thus identified as being significant in maintenance of the three-dimensional structure of the protein.     

The central hydrophobic domain of the bovine papillomavirus E5 transforming protein can be functionally replaced by many hydrophobic amino acid sequences containing a glutamine.

The 44-amino-acid E5 transforming protein of bovine papillomavirus can induce growth transformation of cultured rodent fibroblast cell lines. Previous studies revealed that efficient transformation of mouse C127 cells by the E5 protein required a central core of hydrophobic amino acids and several specific carboxyl-terminal amino acids. Although a randomly derived sequence of hydrophobic amino acids could functionally replace the wild-type hydrophobic core, most such sequences could not. We show here that the conserved glutamine at position 17 in the hydrophobic domain is also important for transformation and that insertion of the glutamine can rescue the transforming activity of many but not all otherwise defective mutants containing random hydrophobic sequences. However, a class of mutants was identified that transform efficiently even in the absence of glutamine, demonstrating that the presence of this amino acid is not absolutely required for efficient transformation. E5 proteins containing the glutamine appear to display increased homodimer formation compared with mutant proteins lacking the glutamine, but this amino acid has no apparent effect on protein stability.     

Electrostatics significantly affect the stability of designed homeodomain variants

The role of electrostatic interactions in determining the stability of designed proteins was studied by constructing and analyzing a set of designed variants of the Drosophila engrailed homeodomain. Computational redesign of 29 surface positions results in a 25-fold mutant with moderate stability, similar to the wild-type protein. Incorporating helix dipole and N-capping considerations into the design algorithm by restricting amino acid composition at the helix termini and N-capping positions yields a ninefold mutant of the initial design (a 23-fold mutant of wild-type) that is over 3 kcal mol(-1) more stable than the protein resulting from the unbiased design. Four additional proteins were constructed and analyzed to isolate the effects of helix dipole and N-capping interactions in each helix. Based on the results of urea-denaturation experiments and calculations using the finite difference Poisson-Boltzmann method, both classes of interaction are found to increase the stability of the designed proteins significantly. The simple electrostatic model used in the optimization of rotamers by iterative techniques (ORBIT) force-field, which is similar to the electrostatic models used in other protein design force-fields, is unable to predict the experimentally determined stabilities of the designed variants. The helix dipole and N-capping restrictions provide a simple but effective method to incorporate two types of electrostatic interactions that impact protein stability significantly.     

An evaluation of Poisson-Boltzmann electrostatic free energy calculations through comparison with experimental mutagenesis data

For systems involving highly and oppositely charged proteins, electrostatic forces dominate association and contribute to biomolecular complex stability. Using experimental or theoretical alanine-scanning mutagenesis, it is possible to elucidate the contribution of individual ionizable amino acids to protein association. We evaluated our electrostatic free energy calculations by comparing calculated and experimental data for alanine mutants of five protein complexes. We calculated Poisson-Boltzmann electrostatic free energies based on a thermodynamic cycle, which incorporates association in a reference (Coulombic) and solvated (solution) state, as well as solvation effects. We observe that Coulombic and solvation free energy values correlate with experimental data in highly and oppositely charged systems, but not in systems comprised of similarly charged proteins. We also observe that correlation between solution and experimental free energies is dependent on dielectric coefficient selection for the protein interior. Free energy correlations improve as protein dielectric coefficient increases, suggesting that the protein interior experiences moderate dielectric screening, despite being shielded from solvent. We propose that higher dielectric coefficients may be necessary to more accurately predict protein-protein association. Additionally, our data suggest that Coulombic potential calculations alone may be sufficient to predict relative binding of protein mutants.