Atomic Density Distributions in Proteins: Structural and Functional Implications

Technische Samenvatting: Atoomdichtheidsdistributies in Eiwitten: Structurele en Functionele Implicaties

Probleemstelling
Atoompakkingsgraad is een fundamentele metriek voor het karakteriseren van eiwitstructuren, wat invloed heeft op stabiliteit, evolutionaire snelheden en functionele rollen. Hoewel het vaststaat dat de binnenkant van eiwitten over het algemeen compacter is dan de oppervlakte en dat pakkingsdefecten (holtes) de stabiliteit kunnen compromitteren, is de distributie van de atoomdichtheid niet uniform over alle eiwitten. Eerdere studies hebben grotendeels vertrouwd op gemiddelde statistische eigenschappen of specifieke structurele metrieken (bijv. contactgetallen) om verpakking te analyseren. Echter, een uitgebreide analyse waarbij individuele atoomdichtheidsprofielen over een enorme dataset van eiwitten worden vergeleken om systematische structurele en functionele patronen te identificeren, ontbrak. Deze studie adresseert de behoefte om eiwitdichtheid te karakteriseren, niet door middel van gemiddelden, maar door direct individuele dichtheidsdistributies afgeleid van een grote, niet-redundante set eiwitstructuren te vergelijken.

Methodologie
De auteurs analyseerden een dataset van 21.255 niet-redundante eiwitstructuren (afkomstig van een initiële set van 22.478, gefilterd om structuren met >80.000 atomen uit te sluiten) verkregen uit de Protein Data Bank (PDB). De selectiecriteria omvatten een resolutieplafond van 2,2 Å, een R-factor plafond van 0,25 en een sequentie-identiteit plafond van ≤50%.

De kernmethodologie bestond uit de volgende stappen:

1. Systeemvoorbereiding: Biologische assemblage-eenheden werden gebruikt voor alle berekeningen. Ontbrekende waterstofatomen werden toegevoegd met behulp van OpenBabel, en structuren werden gesolvateerd met vooraf geëquilibreerde watermoleculen met behulp van het Solvate-programma om in vitro-condities na te bootsen, terwijl het toevoegen van water in begraven holtes werd vermeden.
2. Dichtheidsberekening: Voor elk atoom in een structuur werd de lokale atoomdichtheid (in Da/ų) berekend binnen een sfeer gecentreerd op dat atoom. De massa van alle atomen die binnen de sfeer vielen, werd opgeteld en gedeeld door het volume van de sfeer. Dit werd recursief uitgevoerd voor elk atoom in het eiwit.
3. Parametertesten: Drie sfereradii (5 Å, 6 Å en 7 Å) werden getest om bias door een arbitraire radiusselectie te mitigeren. Radii >7 Å bleken onvoldoende resolutie te bieden door excessieve middeling.
4. Datalverwerking: Dichtheidswaarden voor elk eiwit werden omgezet in histogrammen met 100 bins. Deze distributies werden genormaliseerd door het totaal aantal atomen om frequentiedistributies te genereren.
5. Statistische Analyse:
   • Principal Component Analysis (PCA): Gebruikt om de distributie van de data te visualiseren en uitschieters te identificeren.
   • Hiërarchische Clustering: Euclidische afstanden tussen alle paren dichtheidsdistributies werden berekend om afstandsmatrices te genereren. Deze matrices werden onderworpen aan hiërarchische clustering (met behulp van de "complete" linkage methode) om eiwitten met vergelijkbare dichtheidsprofielen te groeperen.
   • Correlatieanalyse: Relaties tussen dichtheidclusters en structurele kenmerken (grootte, secundaire structuur, B-factoren, waterovervloed) en functionele classificaties (EC-families) werden onderzocht.

Belangrijkste Resultaten

• Significante Variabiliteit: Cumulatieve dichtheidsdistributies toonden aan dat hoewel de meeste eiwitten rond een gemiddelde clusteren, een statistisch significant aantal structuren afwijkt als uitschieters (sommigen door >3σ).
• Clustering van Distributies: Hiërarchische clustering van de 6 Å radius-data (gekozen als de optimale balans van informatiegehalte) identificeerde 12 duidelijke clusters. Verdere analyse van de vijf meest onderscheidende clusters (exclusief kleine groepen) onthulde systematische verschuivingen in dichtheid:
  • Lage-dichtheidclusters: Bevatten eiwitten zoals guanidinobutyrase.
  • Hoge-dichtheidclusters: Bevatten cytokromen en coiled-coil structuren.
• Structurele Correlaties:
  • Eiwitgrootte: Een negatieve correlatie werd waargenomen tussen de eiwitgrootte (aantal residuen) en de mediane atoomdichtheid. Grotere eiwitten neigen compacter gepakt te zijn en vertonen een nauwer bereik van dichtheidswaarden vergeleken met kleinere eiwitten, die een grotere variabiliteit vertonen.
  • Stabiliteitsindicatoren: Hoge pakkingsdichtheid correleerde positief met de abundantie van geordende kristallografische watermoleculen en negatief met atomaire B-factoren (temperatuurfactoren), wat suggereert dat compact gepakte structuren rigider en stabieler zijn.
  • Secundaire Structuur: Dichte clusters vertoonden een hoger percentage α-helices en een lager percentage β-strands. Specifiek waren "coiled-coil" motieven significant oververtegenwoordigd in de dichtste clusters, consistent met hun "knobs-into-holes" verpakkingsvereiste voor mechanische stabiliteit.
• Functionele Correlaties:
  • Cytokromen: Deze elektronentransporteiwitten waren significant oververtegenwoordigd in de dichtste clusters, wat de hypothese ondersteunt dat een compacte verpakking efficiënt elektronentransport faciliteert.
  • Enzymfamilies: Hydrolasen toonden een voorkeur voor dichtere clusters, terwijl transferasen en ligasen de voorkeur gaven aan lossere clusters. De auteurs suggereren dat dit verband kan houden met de "stabiliteit-activiteit"-trade-off en de mechanistische verschillen in het omgaan met water (abundant in de cellulaire omgeving) versus andere substraten.

Betekenis en Claims
De auteurs beweren dat hun benadering van het vergelijken van individuele dichtheidsprofielen een robuuste methode biedt voor het kwantitatief karakteriseren van structurele en functionele patronen in eiwitten. De studie valideert en breidt eerdere bevindingen over de relatie tussen verpakking, stabiliteit en eiwitgrootte uit.

Kernclaims zijn onder meer:

1. Methodologische Nuttigheid: De berekende atoomdichtheidsdistributies dienen als kwaliteitsmetrieken die het verpakkingsniveau van eiwitstructuren kunnen classificeren en uitschieters met unieke structurele eigenschappen kunnen identificeren.
2. Structuur-Functie Koppeling: Specifieke structurele motieven (coiled-coils) en functionele families (cytokromen, hydrolasen) vertonen duidelijke verpakkingsvoorkeuren. Compacte verpakking lijkt een functionele vereiste voor mechanische stabiliteit en elektronentransport, terwijl lossere verpakking geassocieerd kan zijn met specifieke enzymatische mechanismen (bijv. transferasen).
3. Evolutionair Perspectief: De waargenomen variabiliteit in atoomdichtheid wordt geïnterpreteerd, niet noodzakelijkerwijs als het resultaat van actieve, huidige selectie voor specifieke dichtheden, maar potentieel als een restant van moleculaire evolutieprocessen. De auteurs suggereren dat de trade-off tussen stabiliteit en nieuwe enzymatische functies de waargenomen distributie van verpakkingsdichtheden over verschillende eiwitfamilies kan drijven.

Het artikel concludeert dat hoewel systematische patronen in atoomdichtheid bestaan, deze niet in een vaste, één-op-één manier overeenkomen met structurele of functionele kenmerken, wat de complexe interactie tussen evolutionaire geschiedenis en functionele beperkingen weerspiegelt.