Experimental and numerical study of the poplar plastocyanin isoforms using Tyr as a probe for electrostatic similarity and dissimilarity

Modeling pKas of unfolded proteins to probe structural models of unfolded state

Modeling unfolded states of proteins has implications for protein folding and stability. Since in unfolded state proteins adopt multiple conformations, any experimentally measured quantity is ensemble averaged, therefore the computed quantity should be ensemble averaged as well. Here, we investigate the possibility that one can model an unfolded state ensemble with the coil model approach, algorithm such as “flexible-meccano” [Ozenne V et al., Flexible-meccano: A tool for the generation of explicit ensemle descriptions of intrinsically disordered proteins and their associated experimental observables, Bioinformatics 28:1463–1470, 2012], developed to generate structures for intrinsically disordered proteins. We probe such a possibility by using generated structures to calculate pKas of titratable groups and compare with experimental data. It is demonstrated that even with a small number of representative structures of unfolded state, the average calculated pKas are in very good agreement with experimentally measured pKas. Also, predictions are made for titratable groups for which there is no experimental data available. This suggests that the coil model approach is suitable for generating 3D structures of unfolded state of proteins. To make the approach suitable for large-scale modeling, which requires limited number of structures, we ranked the structures according to their solvent accessible surface area (SASA). It is shown that in the majority of cases, the top structures with smallest SASA are enough to represent unfolded state.

Structural comparison of the poplar plastocyanin isoforms PCa and PCb sheds new light on the role of the copper site geometry in interactions with redox partners in oxygenic photosynthesis

Plastocyanin (PC) from poplar leaves is present in two isoforms, PCa and PCb, which differ in sequence by amino acid replacements at locations remote from the copper center and simultaneously act in the photosynthetic electron-transport chain. We describe ultra-high resolution structures of PCa and high-resolution structures of PCb, both under oxidizing and reducing conditions at pH 4, 6 and 8. The docking on cytochrome f and photosystem I, respectively, has been modeled for both isoforms. PCa and PCb exhibit closely similar overall and active-site structures, except for a difference in the relative orientation of the acidic patches. The isoforms exhibit substantial differences in the dependence of the reduced (Cu(I)) geometry on pH. In PCa, the decrease in pH causes a gradual dissociation of His87 from Cu(I) at low pH, probably adopting a neutral tautomeric state. In PCb, the histidine remains covalently bound to Cu(I) and may adopt a doubly protonated state at low pH. The fact that both isoforms have similar although not identical functions in photosynthetic electron flows suggests that the His87 imidazole does not play a crucial role for the pathway of electron transport from cytochrome f to oxidized PC.

A fast and accurate method for predicting pKa of residues in proteins

Predicting the pH-activities of residues in proteins is an important problem in enzyme engineering and protein design. A novel predictor called 'Pred-pK(a)' was developed based on the physicochemical properties of amino acids and protein 3D structure. The Pred-pK(a) approach considers the influence of all other residues of the protein to predict the pK(a) value of an ionizable residue. An empirical equation was formulated, in which the pK(a) value was a distance-dependent function of physicochemical parameters of 20 amino acid types, describing their electrostatic and van der Waals interaction, as well as the effects of hydrogen bonds and solvation. Two sets of coefficients, {a(alpha)} and {b(l)}, were used in the predictor: {a(alpha)} is the weight factors of 20 amino acid types and {b(l)} is the weight factors of physicochemical properties of amino acids. An iterative double least square procedure was proposed to solve the two sets of weight factors alternately and iteratively in a training set. The two coefficient sets {a(alpha)} and {b(l)} thus obtained were used to predict the pK(a) values of residues in a protein. The average predictive error is +/-0.6 pH in less than a minute in common personal computer.

Mutants, overexpressors, and interactors of Arabidopsis plastocyanin isoforms: revised roles of plastocyanin in photosynthetic electron flow and thylakoid redox state

Two homologous plastocyanin isoforms are encoded by the genes PETE1 and PETE2 in the nuclear genome of Arabidopsis thaliana. The PETE2 transcript is expressed at considerably higher levels and the PETE2 protein is the more abundant isoform. Null mutations in the PETE genes resulted in plants, designated pete1 and pete2, with decreased plastocyanin contents. However, despite reducing plastocyanin levels by over approximately 90%, a pete2 null mutation on its own affects rates of photosynthesis and growth only slightly, whereas pete1 knockout plants, with about 60-80% of the wild-type plastocyanin level, did not show any alteration. Hence, plastocyanin concentration is not limiting for photosynthetic electron flow under optimal growth conditions, perhaps implying other possible physiological roles for the protein. Indeed, plastocyanin has been proposed previously to cooperate with cytochrome c(6A) (Cyt c(6A)) in thylakoid redox reactions, but we find no evidence for a physical interaction between the two proteins, using interaction assays in yeast. We observed homodimerization of Cyt c(6A) in yeast interaction assays, but also Cyt c(6A) homodimers failed to interact with plastocyanin. Moreover, phenotypic analysis of atc6-1 pete1 and atc6-1 pete2 double mutants, each lacking Cyt c(6A) and one of the two plastocyanin-encoding genes, failed to reveal any genetic interaction. Overexpression of either PETE1 or PETE2 in the pete1 pete2 double knockout mutant background results in essentially wild-type photosynthetic performance, excluding the possibility that the two plastocyanin isoforms could have distinct functions in thylakoid electron flow.

Comparative study of the stability of poplar plastocyanin isoforms

The stability of the two isoforms of poplar plastocyanin (PCa and PCb) was studied with differential scanning calorimetry (DSC) technique. It was shown that the thermal unfolding of both isoforms is an irreversible process with two endothermic and one exothermic peaks. The melting temperature of PCb was found to be 1.3+/-0.2 K degrees higher than of PCa, which indicates that PCb is more stable. The enthalpy of unfolding was estimated from the heat capacity curves and was found to be significantly higher for PCb at salt concentration I=0.1 M. In addition, PCb unfolding enthalpy and melting temperature are much more sensitive to the changes in the salt concentration as found in the experiments done at different ionic strength. The experiments were complemented with numerical calculations. The salt effect on the stability was modeled using the X-ray structure of PCa and a homology modeled structure of PCb. It was found, in agreement with the experimental data, that the stability of PCb changes by 4.7 kJ more than PCa, as the salt concentration increases from zero to 0.1 M. Thus, the differences in only 12 amino acid positions between "a" and "b" isoforms result in a measurable difference in the folding enthalpy and a significant difference in the salt dependence. The optimization of the electrostatic energies of PCa and PCb were studied and it was shown that PCb is better electrostatically optimized.