Comparative study of the stability of poplar plastocyanin isoforms

Effects of protein and phosphate buffer concentrations on thermal denaturation of lysozyme analyzed by isoconversional method

Thermal denaturation of lysozymes was studied as a function of protein concentration, phosphate buffer concentration, and scan rate using differential scanning calorimetry (DSC), which was then analyzed by the isoconversional method. The results showed that lysozyme thermal denaturation was only slightly affected by the protein concentration and scan rate. When the protein concentration and scan rate increased, the denaturation temperature (Tm) also increased accordingly. On the contrary, the Tm decreased with the increase of phosphate buffer concentration. The denaturation process of lysozymes was accelatated and the thermal stability was reduced with the increase of phosphate concentration. One part of degeneration process was not reversible where the aggregation occurred. The other part was reversible. The apparent activation energy (Ea) was computed by the isoconversional method. It decreased with the increase of the conversion ratio (α). The observed denaturation process could not be described by a simple reaction mechanism. It was not a process involving 2 standard reversible states, but a multi-step process. The new opportunities for investigating the kinetics process of protein denaturation can be supplied by this novel isoconversional method.

Isolation and functional characterization of DNA damage repair protein (DRT) from Lepidium latifolium L

We have isolated and in silico characterized a cold regulated plastocyanin encoding gene from Lepidium latifolium L designated as LlaDRT. Its cDNA sequence (JN214346) consists of a 504 bp ORF, 48 and 205 bp of 5' and 3' UTR regions, respectively encoding a protein of 17.07 KDa and pI 4.95. In silico and phylogenetic analysis of LlaDRT suggested that the protein has features of a typical plastocyanin family member and of a nearest relative of the predominant isoform of Arabidopsis (PETE2) plastocyanin. Validation of stress response of LlaDRT by qPCR under different abiotic stress regulators viz salicylic acid, jasmonic acid, calcium chloride, ethylene and abscisic acid revealed its possible regulation and crosstalk amongst different pathways.

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.

Cupricyclins, novel redox-active metallopeptides based on conotoxins scaffold

Highly stable natural scaffolds which tolerate multiple amino acid substitutions represent the ideal starting point for the application of rational redesign strategies to develop new catalysts of potential biomedical and biotechnological interest. The knottins family of disulphide-constrained peptides display the desired characteristics, being highly stable and characterized by hypervariability of the inter-cysteine loops. The potential of knottins as scaffolds for the design of novel copper-based biocatalysts has been tested by engineering a metal binding site on two different variants of an ω-conotoxin, a neurotoxic peptide belonging to the knottins family. The binding site has been designed by computational modelling and the redesigned peptides have been synthesized and characterized by optical, fluorescence, electron spin resonance and nuclear magnetic resonance spectroscopy. The novel peptides, named Cupricyclin-1 and -2, bind one Cu²⁺ ion per molecule with nanomolar affinity. Cupricyclins display redox activity and catalyze the dismutation of superoxide anions with an activity comparable to that of non-peptidic superoxide dismutase mimics. We thus propose knottins as a novel scaffold for the design of catalytically-active mini metalloproteins.

Spatiotemporal analysis of copper homeostasis in Populus trichocarpa reveals an integrated molecular remodeling for a preferential allocation of copper to plastocyanin in the chloroplasts of developing leaves

Plastocyanin, which requires copper (Cu) as a cofactor, is an electron carrier in the thylakoid lumen and essential for photoautotrophic growth of plants. The Cu microRNAs, which are expressed during Cu deprivation, down-regulate several transcripts that encode for Cu proteins. Since plastocyanin is not targeted by the Cu microRNAs, a cofactor economy model has been proposed in which plants prioritize Cu for use in photosynthetic electron transport. However, defects in photosynthesis are classic symptoms of Cu deprivation, and priorities in Cu cofactor delivery have not been determined experimentally. Using hydroponically grown Populus trichocarpa (clone Nisqually-1), we have established a physiological and molecular baseline for the response to Cu deficiency. An integrated analysis showed that Cu depletion strongly reduces the activity of several Cu proteins including plastocyanin, and consequently, photosynthesis and growth are decreased. Whereas plastocyanin mRNA levels were only mildly affected by Cu depletion, this treatment strongly affected the expression of other Cu proteins via Cu microRNA-mediated transcript down-regulation. Polyphenol oxidase was newly identified as Cu regulated and targeted by a novel Cu microRNA, miR1444. Importantly, a spatiotemporal analysis after Cu resupply to previously depleted plants revealed that this micronutrient is preferentially allocated to developing photosynthetic tissues. Plastocyanin and photosynthetic electron transport efficiency were the first to recover after Cu addition, whereas recovery of the other Cu-dependent activities was delayed. Our findings lend new support to the hypothesis that the Cu microRNAs serve to mediate a prioritization of Cu cofactor use. These studies also highlight poplar as an alternative sequenced model for spatiotemporal analyses of nutritional homeostasis.

Contribution of plastocyanin isoforms to photosynthesis and copper homeostasis in Arabidopsis thaliana grown at different copper regimes

In land plants plastocyanin is indispensable and therefore copper (Cu) availability is a prerequisite for growth. When Cu supply is limited, higher plants prioritize the Cu delivery to plastocyanin by down-regulation of other Cu proteins. Arabidopsis has two plastocyanin genes (PETE1 and PETE2). PETE2 is the predominant isoform in soil-grown plants and in hydroponic cultures it is accumulated in response to Cu addition. It functions as a Cu sink when more Cu is available, in addition to its role as an electron carrier. PETE1 is not affected by Cu feeding and it is the isoform that drives electron transport under Cu-deficiency. Cu feeding rescued the defect in photosystem II electron flux (Phi(PSII)) in the pete1 mutant whereas Phi(PSII) was not changed in the pete2 mutant as Cu was added. Plants with mutations in the plastocyanin genes had altered Cu homeostasis. The pete2 mutant accumulated more Cu/Zn superoxide dismutase (CSD2 and CSD1) and Cu chaperone (CCS) whereas the pete1 mutant accumulated less. On the other hand, less iron superoxide dismutase (FeSOD) and microRNA398b were observed in the pete2 mutant, whereas more were accumulated in the pete1 mutant. Our data suggest that plastocyanin isoforms are different in their response to Cu and the absence of either one changes the Cu homeostasis. Also a small amount of plastocyanin is enough to support efficient electron transport and more PETE2 is accumulated as more Cu is added, presumably, to buffer the excess Cu.

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.

Assessing the quality of the homology-modeled 3D structures from electrostatic standpoint: test on bacterial nucleoside monophosphate kinase families

In this study, we address the issue of performing meaningful pK(a) calculations using homology modeled three-dimensional (3D) structures and analyze the possibility of using the calculated pK(a) values to detect structural defects in the models. For this purpose, the 3D structure of each member of five large protein families of a bacterial nucleoside monophosphate kinases (NMPK) have been modeled by means of homology-based approach. Further, we performed pK(a) calculations for the each model and for the template X-ray structures. Each bacterial NMPK family used in the study comprised on average 100 members providing a pool of sequences and 3D models large enough for reliable statistical analysis. It was shown that pK(a) values of titratable groups, which are highly conserved within a family, tend to be conserved among the models too. We demonstrated that homology modeled structures with sequence identity larger than 35% and gap percentile smaller than 10% can be used for meaningful pK(a) calculations. In addition, it was found that some highly conserved titratable groups either exhibit large pK(a) fluctuations among the models or have pK(a) values shifted by several pH units with respect to the pK(a) calculated for the X-ray structure. We demonstrated that such case usually indicates structural errors associated with the model. Thus, we argue that pK(a) calculations can be used for assessing the quality of the 3D models by monitoring fluctuations of the pK(a) values for highly conserved titratable residues within large sets of homologous proteins.

Optimization of electrostatic interactions in protein-protein complexes

In this article, we present a statistical analysis of the electrostatic properties of 298 protein-protein complexes and 356 domain-domain structures extracted from the previously developed database of protein complexes (ProtCom, http://www.ces.clemson.edu/compbio/protcom). For each structure in the dataset we calculated the total electrostatic energy of the binding and its two components, Coulombic and reaction field energy. It was found that in a vast majority of the cases (>90%), the total electrostatic component of the binding energy was unfavorable. At the same time, the Coulombic component of the binding energy was found to favor the complex formation while the reaction field component of the binding energy opposed the binding. It was also demonstrated that the components in a wild-type (WT) structure are optimized/anti-optimized with respect to the corresponding distributions, arising from random shuffling of the charged side chains. The degree of this optimization was assessed through the Z-score of WT energy in respect to the random distribution. It was found that the Z-scores of Coulombic interactions peak at a considerably negative value for all 654 cases considered while the Z-score of the reaction field energy varied among different types of complexes. All these findings indicate that the Coulombic interactions within WT protein-protein complexes are optimized to favor the complex formation while the total electrostatic energy predominantly opposes the binding. This observation was used to discriminate WT structures among sets of structural decoys and showed that the electrostatic component of the binding energy is not a good discriminator of the WT; while, Coulombic or reaction field energies perform better depending upon the decoy set used.

Poisson-Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies

The salt dependence of the binding free energy of five protein-protein hetero-dimers and two homo-dimers/tetramers was calculated from numerical solutions to the Poisson-Boltzmann equation. Overall, the agreement with experimental values is very good. In all cases except one involving the highly charged lactoglobulin homo-dimer, increasing the salt concentration is found both experimentally and theoretically to decrease the binding affinity. To clarify the source of salt effects, the salt-dependent free energy of binding is partitioned into screening terms and to self-energy terms that involve the interaction of the charge distribution of a monomer with its own ion atmosphere. In six of the seven complexes studied, screening makes the largest contribution but self-energy effects can also be significant. The calculated salt effects are found to be insensitive to force-field parameters and to the internal dielectric constant assigned to the monomers. Nonlinearities due to high charge densities, which are extremely important in the binding of proteins to negatively charged membrane surfaces and to nucleic acids, make much smaller contributions to the protein-protein complexes studied here, with the exception of highly charged lactoglobulin dimers. Our results indicate that the Poisson-Boltzmann equation captures much of the physical basis of the nonspecific salt dependence of protein-protein complexation.