BGPT: Paper Review: The Linked Conservation of Structure and Function in a Family of High Diversity

Explore by Goal

Quick Explanation Copied

Bottom line: Gough & Chothia (2004) present a careful structure-first analysis showing a conserved ~35–57 residue core around the buried Cu site in monomeric cupredoxins despite large sequence diversity; methods, high‑resolution structures, and sequence counts support the claim, but the sample (seven domains + homologs) and exclusion of multidomain forms limit generality and functional tests

Long Explanation

Visual Paper Analysis — Gough & Chothia (2004): Monomeric Cupredoxins

One-sentence summary

High-resolution structural superposition of seven monomeric cupredoxins reveals a small, highly conserved structural core (35–57 residues) centred on the buried copper binding site; sequence identity is low but structural/functional constraints preserve geometry and residue-class patterns around the Cu ligand.

Primary source:

Quantitative view of the dataset (from Table 1)

Interpretation (visual-first, concise)

High-resolution X-ray structures (resolutions 1.33–1.9 Å) form a robust basis for structural comparisons; PDB set spans 96–155 residues and shows aligned segments substantially smaller than full length in many cases, indicating variable loops/insertions outside the core ().
Conserved core: two related structural sets (I = five cupredoxins; II = two phytocyanins). Set I: 57 residues with conserved conformation; intersection of both sets reduces to ~35 residues focused on Cu site and immediate supporting strands ().
Functional constraint: Cu coordination geometry is described as a distorted tetrahedron (Cys + Met + 2 His) with a frequently weak fifth backbone-oxygen interaction (noted in azurin) — a structural necessity for electron transfer chemistry ().

Critical appraisal — strengths and limitations

Strengths

Structure-first approach using high-resolution X-ray data and explicit hydrogen-bond/contact mapping (PINQ) — robust for identifying geometric conservation around functional sites ().
Clear mapping of sequence-class constraints (surface/neutral/buried classes) and the observation that residue-class patterns (s/n/b) are conserved even when exact residues vary — a useful evolutionary insight.

Limitations and blindspots

Sampling: seven monomeric structures (plus homologous sequences) are informative but small by modern standards; multidomain cupredoxins and membrane contexts were explicitly excluded — limits generality to monomeric soluble domains ().
Functional inference is structural/sequence-based; there are few direct biochemical assays in the paper linking small structural variations to measured electron-transfer rates or redox potentials (authors cite prior functional studies). Structural conservation implies functional constraint but does not prove identical kinetics across divergent members ().
Evolutionary breadth: homolog counts and sequence alignments are used but modern large-scale MSA/coevolution tools (DCA, BPPS/SIPRIS) were not available at the time; reanalysis with current sequence databases and coevolution/contact predictors could validate or refine the 35/57-residue core boundaries.

Where conclusions are robust vs. tentative

Robust: conservation of the core geometry that places the Cu ligands and supporting β-strands is strongly supported by high-resolution structural superpositions and PINQ analyses ().
Tentative: functional equivalence across diverse sequences — the structural evidence strongly suggests conserved electron-transfer function, but small geometric or side-chain changes (e.g., Met -> Gln in stellacyanin) can alter redox potentials; explicit biochemical confirmation across broader family members is needed to quantify effect sizes ().

Possible modern follow-ups (concise)

Recompute multiple sequence alignments across large nr/UniRef databases; run DCA/EVcouplings and AlphaFold/AF2-based structural predictions for many homologs to test whether the 35-residue core and coevolving pairs replicate with far larger sampling.
Measure redox potentials and electron-transfer kinetics for representative divergent cupredoxins (including auracyanin, phytocyanins, multidomain constructs) and map functional shifts to measured structural differences.
Use mutational scans (site-saturation at conserved core positions) to quantify which positions are strictly functionally constrained vs. structurally tolerated — combine with deep mutational scanning for energetic maps.

Reproducibility and data availability

Authors used published PDB coordinates (Table 1) and Arthur Lesk's PINQ program for hydrogen-bond/contact analyses; methods are described sufficiently to reproduce the structural superpositions using the same PDBs and PINQ or modern equivalents (DSSP, STRIDE, in-house scripts). However, raw alignment files and PINQ outputs are not fully archived in the paper supplement, so re-running requires rebuilding the pipeline from the methods and PDB files ().

Summary judgement (concise)

Takeaway: This is a solid, methodical structural analysis (2004 era) that convincingly shows conservation of a small structural/functional core in monomeric cupredoxins despite wide sequence divergence; claims are well-supported for the monomeric family members analysed but should be framed as domain-specific (monomeric) and validated further across broader datasets and functional assays.

Author reviews

Citations used in this review are the primary paper and its internal tables/figures; additional contextual references were limited to those cited by the authors in-text (e.g., Baker 1988, Guss et al. 1996) but the review above only directly cites the primary paper per your instruction.

Feedback:

Updated: March 17, 2026