The linked conservation of structure and function in a family of high diversity: the monomeric cupredoxins

Molecular evolution of Phytocyanin gene and analysis of expression at different coloring periods in apple (Malus domestica)

BACKGROUND

PC (phytocyanin) is a class of copper-containing electron transfer proteins closely related to plant photosynthesis, abiotic stress responses growth and development in plants, and regulation of the expression of some flavonoids and phenylpropanoids, etc., however, compared with other plants, the PC gene family has not been systematically characterized in apple.

RESULTS

A total of 59 MdPC gene members unevenly distributed across 12 chromosomes were identified at the genome-wide level. The proteins of the MdPC family were classified into four subfamilies based on differences in copper binding sites and glycosylation sites: Apple Early nodulin-like proteins (MdENODLs), Apple Uclacyanin-like proteins (MdUCLs), Apple Stellacyanin-like proteins (MdSCLs), and Apple Plantacyanin-like proteins (MdPLCLs). Some MdPC members with similar gene structures and conserved motifs belong to the same group or subfamily. The internal collinearity analysis revealed 14 collinearity gene pairs among members of the apple MdPC gene. Interspecific collinearity analysis showed that apple had 31 and 35 homologous gene pairs with strawberry and grape, respectively. Selection pressure analysis indicated that the MdPC gene was under purifying selection. Prediction of protein interactions showed that MdPC family members interacted strongly with the Nad3 protein. GO annotation results indicated that the MdPC gene also regulated the biosynthesis of phenylpropanoids. Chip data analysis showed that (MdSCL3, MdSCL7 and MdENODL27) were highly expressed in mature fruits and peels. Many cis-regulatory elements related to light response, phytohormones, abiotic stresses and flavonoid biosynthetic genes regulation were identified 2000 bp upstream of the promoter of the MdPC gene, and qRT-PCR results showed that gene members in Group IV (MdSCL1/3, MdENODL27) were up-regulated at all five stages of apple coloring, but the highest expression was observed at the DAF13 (day after fruit bag removal) stage. The gene members in Group II (MdUCL9, MdPLCL3) showed down-regulated or lower expression in the first four stages of apple coloring but up-regulated and highest expression in the DAF 21 stage.

CONCLUSION

Herein, one objective of these findings is to provide valuable information for understanding the structure, molecular evolution, and expression pattern of the MdPC gene, another major objective in this study was designed to lay the groundwork for further research on the molecular mechanism of PC gene regulation of apple fruit coloration.

Beyond the coupled distortion model: structural analysis of the single domain cupredoxin AcoP, a green mononuclear copper centre with original features

Cupredoxins are widely occurring copper-binding proteins with a typical Greek-key beta barrel fold. They are generally described as electron carriers that rely on a T1 copper centre coordinated by four ligands provided by the folded polypeptide. The discovery of novel cupredoxins demonstrates the high diversity of this family, with variations in terms of copper-binding ligands, copper centre geometry, redox potential, as well as biological function. AcoP is a periplasmic cupredoxin belonging to the iron respiratory chain of the acidophilic bacterium Acidithiobacillus ferrooxidans. AcoP presents original features, including high resistance to acidic pH and a constrained green-type copper centre of high redox potential. To understand the unique properties of AcoP, we undertook structural and biophysical characterization of wild-type AcoP and of two Cu-ligand mutants (H166A and M171A). The crystallographic structures, including native reduced AcoP at 1.65 Å resolution, unveil a typical cupredoxin fold. The presence of extended loops, never observed in previously characterized cupredoxins, might account for the interaction of AcoP with physiological partners. The Cu-ligand distances, determined by both X-ray diffraction and EXAFS, show that the AcoP metal centre seems to present both T1 and T1.5 features, in turn suggesting that AcoP might not fit well to the coupled distortion model. The crystal structures of two AcoP mutants confirm that the active centre of AcoP is highly constrained. Comparative analysis with other cupredoxins of known structures, suggests that in AcoP the second coordination sphere might be an important determinant of active centre rigidity due to the presence of an extensive hydrogen bond network. Finally, we show that other cupredoxins do not perfectly follow the coupled distortion model as well, raising the suspicion that further alternative models to describe copper centre geometries need to be developed, while the importance of rack-induced contributions should not be underestimated.

pLM-BLAST: distant homology detection based on direct comparison of sequence representations from protein language models

MOTIVATION

The detection of homology through sequence comparison is a typical first step in the study of protein function and evolution. In this work, we explore the applicability of protein language models to this task.

RESULTS

We introduce pLM-BLAST, a tool inspired by BLAST, that detects distant homology by comparing single-sequence representations (embeddings) derived from a protein language model, ProtT5. Our benchmarks reveal that pLM-BLAST maintains a level of accuracy on par with HHsearch for both highly similar sequences (with >50% identity) and markedly divergent sequences (with <30% identity), while being significantly faster. Additionally, pLM-BLAST stands out among other embedding-based tools due to its ability to compute local alignments. We show that these local alignments, produced by pLM-BLAST, often connect highly divergent proteins, thereby highlighting its potential to uncover previously undiscovered homologous relationships and improve protein annotation.

AVAILABILITY AND IMPLEMENTATION

pLM-BLAST is accessible via the MPI Bioinformatics Toolkit as a web server for searching precomputed databases (https://toolkit.tuebingen.mpg.de/tools/plmblast). It is also available as a standalone tool for building custom databases and performing batch searches (https://github.com/labstructbioinf/pLM-BLAST).

Overexpression of LpCPC from Lilium pumilum confers saline-alkali stress (NaHCO3) resistance

Lilium Pumilum with wide distribution is highly tolerant to salinity. The blue copper protein LpCPC (Lilium pumilum Cucumber Peeling Cupredoxin) gene was cloned from Lilium pumilum, which has the conserved regions of type I copper protein. Moreover, LpCPC has the closest relation to CPC from Actinidia chinensis using DNAMAN software and MEGA7 software. qRT-PCR indicated that LpCPC expression was higher in root and bulb of Lilium pumilum, and the expression of the LpCPC gene increased and reached the highest level at 12 h in bulbs under 20 mM NaHCO₃. The transgenic yeast was more tolerant compared with the control under NaHCO₃ stress. Compared with the wild type, overexpressing plants indicated a relatively lower degree of wilting. In addition, the chlorophyll content, soluble phenol content, and lignin content of overexpressing lines were higher than that of wild-type, whereas the relative conductivity of overexpressing plants was significantly lower than that of wild-type plants. Expression of essential genes including NHX1 and SOS1 in salt stress response pathways are steadily higher in overexpression tobacco than that in wild-types. Transgenic lines had much higher levels of CCR1 and CAD, which are involved in lignin production, compared with wild-type lines. The yeast two-hybrid technique was applied to screen probable interacting proteins interacting with LpCPC. Eight proteins interacted with LpCPC were screened, and five of which were demonstrated to be associated with plant salinity resistance. Overall, the role of gene LpCPC is mediating molecule responses in increasing saline-alkali stress resistance, indicating that it is an essential gene to enhance salt tolerance in Lilium pumilum.

The gene LpBCP increased NaHCO3 resistance by enhancing lignin or ROS scavenging in the Nicotiana benthamiana

Lilium pumilum is an excellent wildflower germplasm resource with high resistance to salinity stress. The gene LpBCP plays an important role in salinity tolerance of L. pumilum. Studying the molecular mechanism of salinity resistance in L. pumilum will provide insights into multiple aspects, including breeding better varieties, environmental protection, improving soil conditions, etc. Conventional methods were used to determine different physiological indicators of Nicotiana benthamiana after NaHCO₃ treatment, i.e. chlorophyll content, soluble phenol content and lignin content. RT-qPCR was carried out to find expression of LpBCP in different organs and under abiotic stresses. DAB was used to detect H₂ O₂ in leaves in situ. A yeast two-hybrid system was used to screen for LpBCP interacting proteins. LpBCP was cloned from bulbs of L. pumilum. The highest expression of LpBCP was in roots and bulbs of transgenic plants. LpBCP-overexpressed plants showed less wilting, compared to WT plants. LpBCP transgenic plants have higher chlorophyll, soluble phenol and lignin content, and lower relative conductivity under 500 mM NaHCO₃ stress. In addition, H₂ O₂ scavenging in transgenic plants was much improved, indicating increased resistance to NaHCO₃ stress. Thirteen LpBCP-interacting proteins were screened using the yeast two-hybrid method and five were associated with salt stress. Based on our findings, LPBCP could be a key gene that can be used to improve L. pumilum salt tolerance.

Effect of the Met148Leu mutation on the structure and dynamics of the rusticyanin protein from Acidithiobacillus sp. FJ2

The rusticyanin protein, a blue monomeric copper protein type-1, is one of the main components in the iron-electron transfer chain of the Acidithiobacillus ferrooxidans, and is the product of the rus gene expression. Herein, first the bacterial DNA of Acidithiobacillus sp. FJ2 was extracted. Then, the rus gene sequence and the sequence amino acid rusticyanin protein were determined. The Met148Leu mutation increased the oxidase activity of the rusticyanin protein, thereby enhancing the efficiency of the bioleaching process by bacteria Acidithiobacillus ferroxidans. Met148Leu mutation was created in the rusticyanin protein, then molecular dynamics (MD) simulations and structural analysis were performed. The MD analysis of the wild-type and mutant protein demonstrated a slight instability in the mutant protein and significant instability in the active site of the mutant protein. The usefulness of this study is the genetic manipulation of the native Acidithiobacillus sp. FJ2 bacterium, which can boost the bioleaching efficiency of the bacterium to some extent, and investigating its effects on the structure of a mutant protein using computational methods.

Identification of Phytocyanin Gene Family in Legume Plants and their Involvement in Nodulation of Medicago truncatula

The establishment of symbiosis between legume and rhizobium results in the formation of nodule. Phytocyanins (PCs) are a class of plant-specific blue copper proteins, playing critical roles in plant development including nodule formation. Although a few PC genes have been isolated from nodules, their functions are still unclear. Here, we performed a genome-wide identification of PC family in seven sequenced legume species (Medicago truncatula, Glycine max, Cicer arietinum, Cajanus cajan, Lotus japonicus, Vigna angularis and Phaseolus vulgaris) and found PCs experienced a remarkable expansion in M. truncatula and G. max. Further, we conducted an in-depth analysis of PC family in the model legume M. truncatula. Briefly, 82 MtPCs were divided into four subfamilies and clustered into seven clades, with a large proportion of tandem duplications and various cross-tissues expression patterns. Importantly, some PCs, such as MtPLC1, MtENODL27 and MtENODL28 were preferentially expressed in nodules. Further, RNA interference (RNAi) experiment revealed the knockdown of MtENDOL27 and MtENDOL28 impaired rhizobia infection, nodule numbers and nitrogenase activity. Moreover, in the MtENODL27-RNAi nodules, the infected cells were reduced and the symbiosomes did not reach the elongated stage, indicating MtENDOL27 is required for rhizobia infection and nodule development. In addition, co-expression analysis showed MtPLC1, MtENODL27 and MtENODL28 were grouped into two different functional modules and co-expressed with the known symbiotic nitrogen fixation-related genes, suggesting that they might participate in nodulation via different ways. In summary, this study provides a useful resource for future researches on the structure and function of PCs in nodulation.

Fine Tuning of Functional Features of the CuA Site by Loop-Directed Mutagenesis

Here we report the spectroscopic and electrochemical characterization of three novel chimeric Cu_A proteins in which either one or the three loops surrounding the metal ions in the Thermus thermophilus protein have been replaced by homologous human and plant sequences while preserving the set of coordinating amino acids. These conservative modifications mimic basic differences between Cu_A sites from different organisms and allow for fine tuning the energy gap between alternative electronic ground states of Cu_A.. This results in a systematic modulation of thermodynamic and kinetic electron transfer (ET) parameters and in the selection of one of two possible redox-active molecular orbitals, which differ in the ET reorganization energy by a factor of 2. Moreover, the ET mechanism is found to be frictionally controlled, and the modifications introduced into the different chimeras do not affect the frictional activation parameter.

Multiple Osmotic Stress Responses in Acidihalobacter prosperus Result in Tolerance to Chloride Ions

Extremely acidophilic microorganisms (pH optima for growth of ≤3) are utilized for the extraction of metals from sulfide minerals in the industrial biotechnology of “biomining.” A long term goal for biomining has been development of microbial consortia able to withstand increased chloride concentrations for use in regions where freshwater is scarce. However, when challenged by elevated salt, acidophiles experience both osmotic stress and an acidification of the cytoplasm due to a collapse of the inside positive membrane potential, leading to an influx of protons. In this study, we tested the ability of the halotolerant acidophile Acidihalobacter prosperus to grow and catalyze sulfide mineral dissolution in elevated concentrations of salt and identified chloride tolerance mechanisms in Ac. prosperus as well as the chloride susceptible species, Acidithiobacillus ferrooxidans. Ac. prosperus had optimum iron oxidation at 20 g L⁻¹ NaCl while At. ferrooxidans iron oxidation was inhibited in the presence of 6 g L⁻¹ NaCl. The tolerance to chloride in Ac. prosperus was consistent with electron microscopy, determination of cell viability, and bioleaching capability. The Ac. prosperus proteomic response to elevated chloride concentrations included the production of osmotic stress regulators that potentially induced production of the compatible solute, ectoine uptake protein, and increased iron oxidation resulting in heightened electron flow to drive proton export by the F₀F₁ ATPase. In contrast, At. ferrooxidans responded to low levels of Cl⁻ with a generalized stress response, decreased iron oxidation, and an increase in central carbon metabolism. One potential adaptation to high chloride in the Ac. prosperus Rus protein involved in ferrous iron oxidation was an increase in the negativity of the surface potential of Rus Form I (and Form II) that could help explain how it can be active under elevated chloride concentrations. These data have been used to create a model of chloride tolerance in the salt tolerant and susceptible species Ac. prosperus and At. ferrooxidans, respectively.

Spectroscopic characterization of a green copper site in a single-domain cupredoxin

Cupredoxins are widespread copper-binding proteins, mainly involved in electron transfer pathways. They display a typical rigid greek key motif consisting of an eight stranded β-sandwich. A fascinating feature of cupredoxins is the natural diversity of their copper center geometry. These geometry variations give rise to drastic changes in their color, such as blue, green, red or purple. Based on several spectroscopic and structural analyses, a connection between the geometry of their copper-binding site and their color has been proposed. However, little is known about the relationship between such diversity of copper center geometry in cupredoxins and possible implications for function. This has been difficult to assess, as only a few naturally occurring green and red copper sites have been described so far. We report herein the spectrocopic characterization of a novel kind of single domain cupredoxin of green color, involved in a respiratory pathway of the acidophilic organism Acidithiobacillus ferrooxidans. Biochemical and spectroscopic characterization coupled to bioinformatics analysis reveal the existence of some unusual features for this novel member of the green cupredoxin sub-family. This protein has the highest redox potential reported to date for a green-type cupredoxin. It has a constrained green copper site insensitive to pH or temperature variations. It is a green-type cupredoxin found for the first time in a respiratory pathway. These unique properties might be explained by a region of unknown function never found in other cupredoxins, and by an unusual length of the loop between the second and the fourth copper ligands. These discoveries will impact our knowledge on non-engineered green copper sites, whose involvement in respiratory chains seems more widespread than initially thought.

The sole tryptophan of amicyanin enhances its thermal stability but does not influence the electronic properties of the type 1 copper site

The cupredoxin amicyanin possesses a single tryptophan residue, Trp45. Its fluorescence is quenched when copper is bound even though it is separated by 10.1Å. Mutation of Trp45 to Ala, Phe, Leu and Lys resulted in undetectable protein expression. A W45Y amicyanin variant was isolated. The W45Y mutation did not alter the spectroscopic properties or intrinsic redox potential of amicyanin, but increased the pKa value for the pH-dependent redox potential by 0.5 units. This is due to a hydrogen-bond involving the His95 copper ligand which is present in reduced W45Y amicyanin but not in native amicyanin. The W45Y mutation significantly decreased the thermal stability of amicyanin, as determined by changes in the visible absorbance of oxidized amicyanin and in the circular dichroism spectra for oxidized, reduced and apo forms of amicyanin. Comparison of the crystal structures suggests that the decreased stability of W45Y amicyanin may be attributed to the loss of a strong interior hydrogen bond between Trp45 and Tyr90 in native amicyanin which links two of the β-sheets that comprise the overall structure of amicyanin. Thus, Trp45 is critical for stabilizing the structure of amicyanin but it does not influence the electronic properties of the copper which quenches its fluorescence.

Amino acid distribution rules predict protein fold

In the present article, we provide a brief overview of the main approaches to analysing the sequence-structure relationship of proteins and outline a novel method of structure prediction. The proposed method involves finding a set of rules that describes a correlation between the distribution of residues in a sequence and the essential structural characteristics of a protein structure. The residue distribution rules specify the 'favourable' residues that are required in certain positions of a polypeptide chain in order for it to assume a particular protein fold, and the 'unfavourable' residues incompatible with the given fold. Identification of amino acid distribution rules derives from examination of inter-residue contacts. We describe residue distribution rules for a large group of β-sandwich-like proteins characterized by a specific arrangement of strands in their two β-sheets. It was shown that this method has very high accuracy (approximately 85%). The advantage of the residue rule approach is that it makes possible prediction of protein folding even in polypeptide chains that have very low global sequence similarities, as low as 18%. Another potential benefit is that a better understanding of which residues play essential roles in a given protein fold may facilitate rational protein engineering design.

Comparative genomics in acid mine drainage biofilm communities reveals metabolic and structural differentiation of co-occurring archaea

Background

Metal sulfide mineral dissolution during bioleaching and acid mine drainage (AMD) formation creates an environment that is inhospitable to most life. Despite dominance by a small number of bacteria, AMD microbial biofilm communities contain a notable variety of coexisting and closely related Euryarchaea, most of which have defied cultivation efforts. For this reason, we used metagenomics to analyze variation in gene content that may contribute to niche differentiation among co-occurring AMD archaea. Our analyses targeted members of the Thermoplasmatales and related archaea. These results greatly expand genomic information available for this archaeal order.

Results

We reconstructed near-complete genomes for uncultivated, relatively low abundance organisms A-, E-, and Gplasma, members of Thermoplasmatales order, and for a novel organism, Iplasma. Genomic analyses of these organisms, as well as Ferroplasma type I and II, reveal that all are facultative aerobic heterotrophs with the ability to use many of the same carbon substrates, including methanol. Most of the genomes share genes for toxic metal resistance and surface-layer production. Only Aplasma and Eplasma have a full suite of flagellar genes whereas all but the Ferroplasma spp. have genes for pili production. Cryogenic-electron microscopy (cryo-EM) and tomography (cryo-ET) strengthen these metagenomics-based ultrastructural predictions. Notably, only Aplasma, Gplasma and the Ferroplasma spp. have predicted iron oxidation genes and Eplasma and Iplasma lack most genes for cobalamin, valine, (iso)leucine and histidine synthesis.

Conclusion

The Thermoplasmatales AMD archaea share a large number of metabolic capabilities. All of the uncultivated organisms studied here (A-, E-, G-, and Iplasma) are metabolically very similar to characterized Ferroplasma spp., differentiating themselves mainly in their genetic capabilities for biosynthesis, motility, and possibly iron oxidation. These results indicate that subtle, but important genomic differences, coupled with unknown differences in gene expression, distinguish these organisms enough to allow for co-existence. Overall this study reveals shared features of organisms from the Thermoplasmatales lineage and provides new insights into the functioning of AMD communities.

Evolutionarily consistent families in SCOP: sequence, structure and function

Background

SCOP is a hierarchical domain classification system for proteins of known structure. The superfamily level has a clear definition: Protein domains belong to the same superfamily if there is structural, functional and sequence evidence for a common evolutionary ancestor. Superfamilies are sub-classified into families, however, there is not such a clear basis for the family level groupings. Do SCOP families group together domains with sequence similarity, do they group domains with similar structure or by common function? It is these questions we answer, but most importantly, whether each family represents a distinct phylogenetic group within a superfamily.

Results

Several phylogenetic trees were generated for each superfamily: one derived from a multiple sequence alignment, one based on structural distances, and the final two from presence/absence of GO terms or EC numbers assigned to domains. The topologies of the resulting trees and confidence values were compared to the SCOP family classification.

Conclusions

We show that SCOP family groupings are evolutionarily consistent to a very high degree with respect to classical sequence phylogenetics. The trees built from (automatically generated) structural distances correlate well, but are not always consistent with SCOP (hand annotated) groupings. Trees derived from functional data are less consistent with the family level than those from structure or sequence, though the majority still agree. Much of GO and EC annotation applies directly to one family or subset of the family; relatively few terms apply at the superfamily level. Maximum sequence diversity within a family is on average 22% but close to zero for superfamilies.

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.

How the local geometry of the Cu-binding site determines the thermal stability of blue copper proteins

Identifying the factors that govern the thermal resistance of cupredoxins is essential for understanding their folding and stability, and for improving our ability to design highly stable enzymes with potential biotechnological applications. Here, we show that the thermal unfolding of plastocyanins from two cyanobacteria--the mesophilic Synechocystis and the thermophilic Phormidium--is closely related to the short-range structure around the copper center. Cu K-edge X-ray absorption spectroscopy shows that the bond length between Cu and the S atom from the cysteine ligand is a key structural factor that correlates with the thermal stability of the cupredoxins in both oxidized and reduced states. These findings were confirmed by an additional study of a site-directed mutant of Phormidium plastocyanin showing a reverse effect of the redox state on the thermal stability of the protein.

Modulation of copper site properties by remote residues determines the stability of plastocyanins

The metal cofactor determines the thermal stability in cupredoxins, but how the redox state of copper modulates their melting points remains unknown. The metal coordination environment is highly conserved in cyanobacterial plastocyanins. However, the oxidised form is more stable than the reduced one in thermophilic Phormidium, but the opposite occurs in mesophilic Synechocystis. We have performed neutral amino-acid substitutions at loops of Phormidium plastocyanin far from the copper site. Notably, mutation P49G/G50P confers a redox-dependent thermal stability similar to that of the mesophilic plastocyanin. Moreover, X-ray absorption spectroscopy reveals that P49G/G50P mutation makes the electron density distribution at the oxidised copper site shift towards that of Synechocystis plastocyanin.

Predicting structurally conserved contacts for homologous proteins using sequence conservation filters

The prediction of intramolecular contacts has a useful application in predicting the three-dimensional structures of proteins. The accuracy of the template-based contact prediction methods depends on the quality of the template structures. To reduce the false positive predictions associated with using the entire set of template-derived contacts, we develop selection filters that use sequence conservation information to predict subsets of contacts more likely to be structurally conserved between the template and the target. The method is developed specifically for protein families with few available templates such as the G protein-coupled receptor (GPCR) family. It is validated on a test set of 342 template-target pairs from three protein families, and applied to one template-target pair from the GPCR family. We find that the filter selection method increases the accuracy of contact prediction with sufficient coverage for structure prediction.

Metal-binding loop length and not sequence dictates structure

The C-terminal copper-binding loop in the beta-barrel fold of the cupredoxin azurin has been replaced with a range of sequences containing alanine, glycine, and valine residues to assess the importance of amino acid composition and the length of this region. The introduction of 2 and 4 alanines between the coordinating Cys, His, and Met results in loop structures matching those in naturally occurring proteins with the same loop lengths. A loop with 4 alanines between the Cys and His and 3 between the His and Met ligands has a structure identical to that of the WT protein, whose loop is the same length. Loop structure is dictated by length and not sequence allowing the properties of the main surface patch for interactions with partners, to which the loop is a major contributor, to be optimized. Loops with 2 amino acids between the ligands using glycine, alanine, and valine residues have been compared. An empirical relationship is found between copper site protection by the loop and reduction potential. A loop adorned with 4 methyl groups is sufficient to protect the copper ion, enabling most sequences to adequately perform this task. The mutant with 3 alanine residues between the ligands forms a strand-swapped dimer in the crystal structure, an arrangement that has not, to our knowledge, been seen previously for this family of proteins. Cupredoxins function as redox shuttles and are required to be monomeric; therefore, none have evolved with a metal-binding loop of this length.

Molecular dynamics of amicyanin reveals a conserved dynamical core for blue copper proteins

Molecular dynamics simulation has been carried out for the blue copper protein amicyanin from two different sources, Paracoccus denitrificans and Paraccocus versutus, to investigate the structural and dynamical properties common to the two molecules and to identify prominent features shared with proteins of the same family, the monomeric cupredoxins. The two amicyanins have almost identical secondary and tertiary structure. In the simulation, they differ for the number of hydrogen bonds in the main chain and the conformation of some beta-strands. However, they strictly maintain the arrangement of the portions of the beta-barrel that are conserved in the folding architecture of the blue copper proteins. Paracoccus versutus amicyanin equilibrates more rapidly, shows lower atomic deviation values, and is less rigid with respect to Paracoccus denitrificans amicyanin. Principal component analysis reveals that the conformational subspaces corresponding to eigenvectors with the same index for each of the two molecules are not necessarily equivalent. Nevertheless, a core scaffold with constrained dynamics exist for both amicyanins. In addition, two fairly flexible regions that are located on the opposite side with respect to the interaction sites with the partner molecules in the redox process have been evidenced in the protein structure. This description of amicyanin, with a few mobile regions remote from the active site and a rigid scaffold including most of the protein beta-barrel, has a close similarity with that of azurin and plastocyanin, two other cupredoxins previously investigated in simulation. Furthermore, similarities in the distribution of the atomic fluctuations indicate that amicyanin, azurin, and plastocyanin possess common dynamical features, in spite of differences in their structure. On the basis of these findings, we suggest that topological constraints imposed by the folding in correspondence of protein regions that are the most conserved determine the protein dynamics of the cupredoxin family. The dynamical properties of the cupredoxins might be controlled for functional advantages that include the binding mechanism with the biological partners and the collective inner motions of the protein matrix required for the electron transfer, whereas long-range conformational changes in the redox reaction should be excluded.

Contribution of A1 subunit residue Q316 in thrombin-activated factor VIII to A2 subunit dissociation

Blood coagulation factor VIII (fVIII) is activated by thrombin to form an A1/A2/A3-C1-C2 heterotrimer, which functions as a cofactor for factor IXa during intrinsic pathway factor X activation. Human thrombin-activated fVIII (fVIIIa) decays rapidly because of first-order dissociation of the A2 subunit, which may function to regulate the coagulation mechanism. The three fVIII A domains each consist of two cupredoxin-like subdomains. Substitution of the COOH-terminal A1 subdomain of porcine fVIIIa, which decays more slowly than human fVIIIa, reduces the dissociation rate constant for fVIIIa decay. Examination of a human fVIII A1-A2-A3 homology model [Pemberton, S., et al. (1997) Blood 89, 2413-2421) revealed a possible interaction between Q316 in the FG helix of the COOH-terminal A1 subdomain and M539 in the FG helix of the NH2-terminal A2 subdomain, which are sites where human and porcine fVIII differ. Decays of purified recombinant human and porcine fVIIIa and the human fVIIIa mutants Q316H, M539L and Q316H/M539L were compared at 23 and 37 degrees C. The decay rates of the Q316H and Q316H/M539L mutants, but not the M539L mutant, were significantly slower than human fVIIIa. These results indicate that the FG helix of the COOH-terminal A1 cupredoxin-like subdomain of fVIII may be under selective pressure by the requirements of hemostatic balance.

Beyond host–pathogen interactions: microbial defense strategy in the host environment

Many extracellular pathogenic bacteria colonize human or animal bodies through evasion of the host immune system, a process called host-pathogen interaction. What happens when other intruders try to invade the same host and try to establish themselves in the same niche is largely unknown. In one well-studied case, Pseudomonas aeruginosa is known to secrete the protein azurin as a weapon against such invaders as cancers, parasites and viruses. The production of such weapons by pathogenic bacteria could provide important insights into how a pathogen responds in the post-colonization state to impede other intruders for its own survival. Moreover, these molecules might find use in the pharmaceutical industry as next-generation therapeutics.

Structural, dynamical and functional aspects of the inner motions in the blue copper protein azurin

Molecular dynamics was applied to dissect out the internal motions of azurin, a copper protein performing electron transfer. Simulations of 16.5 ns were analyzed in search of coordinated displacements of amino acid residues that are important for the protein function. A region with high conformational instability was found in the 'southern' end of the molecule, far away from the copper site and the binding sites for the redox partners of azurin. By excluding the 'southern' region from the subsequent analysis, correlated motions were identified in the hydrophobic patch that surrounds the protein active site. The simulation results are in excellent agreement with recent NMR data on azurin in solution [A. V. Zhuravleva, D. M. Korzhnev, E. Kupce, A. S. Arseniev, M. Billeter, V. Y. Orekhov, Gated electron transfers and electron pathways in azurin: a NMR dynamic study at multiple fields and temperatures, J. Mol. Biol. 342 (2004) 1599-1611] and suggest a rationale for cooperative displacements of protein residues that are thought to be critical for the electron transfer process. A number of other structural and dynamic features of azurin are discussed in the context of the blue copper protein family and an explanation is proposed to account for the variability/conservation of some regions in the cupredoxins.

A new assay for nitric oxide reductase reveals two conserved glutamate residues form the entrance to a proton-conducting channel in the bacterial enzyme

A specific amperometric assay was developed for the membrane-bound NOR [NO (nitric oxide) reductase] from the model denitrifying bacterium Paracoccus denitrificans using its natural electron donor, pseudoazurin, as a co-substrate. The method allows the rapid and specific assay of NO reduction catalysed by recombinant NOR expressed in the cytoplasmic membranes of Escherichia coli. The effect on enzyme activity of substituting alanine, aspartate or glutamine for two highly conserved glutamate residues, which lie in a periplasmic facing loop between transmembrane helices III and IV in the catalytic subunit of NOR, was determined using this method. Three of the substitutions (E122A, E125A and E125D) lead to an almost complete loss of NOR activity. Some activity is retained when either Glu122 or Glu125 is substituted with a glutamine residue, but only replacement of Glu122 with an aspartate residue retains a high level of activity. These results are interpreted in terms of these residues forming the mouth of a channel that conducts substrate protons to the active site of NOR during turnover. This channel is also likely to be that responsible in the coupling of proton movement to electron transfer during the oxidation of fully reduced NOR with oxygen [U. Flock, N. J. Watmough and P. Adelroth (2005) Biochemistry 44, 10711-10719].

The atypical iron-coordination geometry of cytochrome f remains unchanged upon binding to plastocyanin, as inferred by XAS

The transient complex between cytochrome f and plastocyanin from the cyanobacterium Nostoc sp. PCC 7119 has been analysed by X-ray Absorption Spectroscopy in solution, using both proteins in their oxidized and reduced states. Fe K-edge data mainly shows that the atypical metal coordination geometry of cytochrome f, in which the N-terminal amino acid acts as an axial ligand of the heme group, remains unaltered upon binding to its redox partner, plastocyanin. This fact suggests that cytochrome f provides a stable binding site for plastocyanin and minimizes the reorganization energy required in the transient complex formation, which could facilitate the electron transfer between the two redox partners.

Divergence of interdomain geometry in two-domain proteins

For homologous protein chains composed of two domains, we have determined the extent to which they conserve (1) their interdomain geometry and (2) the molecular structure of the domain interface. This work was carried out on 128 unique two-domain architectures. Of the 128, we find 75 conserve their interdomain geometry and the structure of their domain interface; 5 conserve their interdomain geometry but not the structure of their interface; and 48 have variable geometries and divergent interface structure. We describe how different types of interface changes or the absence of an interface is responsible for these differences in geometry. Variable interdomain geometries can be found in homologous structures with high sequence identities (70%).

Internalization of bacterial redox protein azurin in mammalian cells: entry domain and specificity

Azurin is a member of a group of copper-containing redox proteins called cupredoxins. Different cupredoxins are produced by different aerobic bacteria as agents of electron transfer. Recently, we demonstrated that azurin enters into J774 and several types of cancer cells leading to the induction of apoptosis. We now demonstrate that azurin is internalized in J774 or cancer cells in a temperature-dependent manner. Azurin shows preferential entry into cancer compared with normal cells. An 28-amino-acid fragment of azurin fused to glutathione S-transferase (GST) or the green fluorescent protein (GFP), which are incapable of entering mammalian cells by themselves, can be internalized in J774 or human melanoma or breast cancer cells at 37 degrees C, but not at 4 degrees C. Competition experiments as well as studies with inhibitors such as cytochalasin D suggest that azurin may enter cells, at least in part, by a receptor-mediated endocytic process. The 28-amino-acid peptide therefore acts as a potential protein transduction domain (PTD), and can be used as a vehicle to transport cargo proteins such as GST and GST-GFP fusion proteins. Another member of the cupredoxin family, rusticyanin, that has also been shown to enter J774 and human cancer cells and exert cytotoxicity, does not demonstrate preferential entry for cancer cells and lacks the structural features characteristic of the azurin PTD.