Extending the MARTINI 3 Coarse-Grained Forcefield to Polypeptoids

In dit werk wordt de eerste MARTINI 3-compatible coarse-grained krachtveld voor peptoiden ontwikkeld, dat 19 veelvoorkomende residu-types omvat, geoptimaliseerd is met behulp van all-atom simulaties en geïntegreerd is in de martinize2-tool om efficiënte mesoscale simulaties van peptoid-structuur en zelfassemblage mogelijk te maken.

De kern

De "Vereenvoudigde Bouwplaat" voor Kunstmatige Eiwit-Nabootsers

Stel je voor dat je een gigantische, ingewikkelde machine wilt bouwen, zoals een robot die medicijnen door je lichaam kan vervoeren. In de natuur zijn dit soort machines vaak eiwitten. Maar wetenschappers hebben een speciale, kunstmatige versie van deze machines bedacht die ze peptoiden noemen. Ze lijken op eiwitten, maar ze zijn sterker, blijven langer stabiel en zijn makkelijker te ontwerpen.

Het probleem is echter: deze machines zijn zo complex dat ze als een enorme, verwarrende lade vol losse schroefjes, boutjes en veertjes voelen. Als je ze in een computer wilt simuleren om te zien hoe ze werken, moet je elke individuele schroef (elk atoom) berekenen. Dat is voor computers net zo zwaar als het proberen te voorspellen van het weer op elke individuele steen in een stad. Het duurt te lang en kost te veel energie.

De Oplossing: De "MARTINI 3" Bouwplaat

In dit artikel vertellen de onderzoekers hoe ze een nieuwe, slimme manier hebben bedacht om deze peptoiden te simuleren. Ze gebruiken een methode die MARTINI 3 heet.

• De Analogie: Stel je voor dat je in plaats van elke losse schroef te tekenen, je de machine in grote blokken verdeelt. Een blokje "motor" staat voor 10 schroeven, een blokje "wiel" staat voor 20 boutjes. Je ziet nog steeds hoe de machine werkt en hoe hij beweegt, maar je hebt 57 keer minder details om te rekenen.
• Het Nieuwe: Vroeger waren er geen goede "bouwplaten" voor deze specifieke kunstmatige peptoiden die werkten met de moderne MARTINI 3-systeem. De onderzoekers hebben nu de eerste complete set bouwplaten gemaakt voor 19 verschillende soorten peptoid-blokjes.

Hoe hebben ze dit gedaan?

1. De "Tijdsreismachine" (PBMetaD): Omdat de peptoiden heel flexibel zijn en soms vastlopen in een bepaalde vorm (net als een deur die vastzit in de scharnieren), hebben de onderzoekers een speciale techniek gebruikt. Dit is alsof je de deur even met een hefboom open duwt om te zien wat erachter zit, en dan de deur weer laat vallen. Zo hebben ze alle mogelijke vormen van de peptoiden in kaart gebracht.
2. Het "Vertalen" (Boltzmann Inversie): Ze hebben gekeken naar de gedetailleerde "schroefjes-versie" (de atomaire simulatie) en gezegd: "Als de schroeven zich zo gedragen, hoe moet ons grote blokje dan bewegen?" Ze hebben de regels vertaald van de fijne details naar de grote blokken.
3. De "Standaardset": Ze hebben ervoor gezorgd dat hun nieuwe blokjes perfect passen in het bestaande MARTINI 3-systeem. Je kunt nu peptoiden mengen met vetten, water en andere moleculen zonder dat het systeem in de war raakt.

Wat levert dit op?

• Snelheid: De simulaties zijn nu 57 keer sneller. Wat vroeger maanden duurde, gaat nu in dagen of uren.
• Betrouwbaarheid: Ze hebben getest of de nieuwe bouwplaten kloppen. Ze hebben gekeken of de peptoiden zich in het water gedragen zoals ze dat in de echte wereld zouden doen (bijvoorbeeld: klitten ze samen of drijven ze los?). Het antwoord is ja: de nieuwe modellen zijn heel nauwkeurig.
• Toekomst: Met deze nieuwe tool kunnen wetenschappers nu makkelijker nieuwe medicijnen ontwerpen of slimme materialen bouwen die zichzelf kunnen assembleren, zoals nanobots die ziektes bestrijden.

Kortom:
De onderzoekers hebben een "verkleinde versie" gemaakt van complexe kunstmatige eiwitten. Hierdoor kunnen computers deze moleculen veel sneller bestuderen, waardoor we sneller nieuwe medicijnen en materialen kunnen uitvinden. Het is alsof ze van een ingewikkeld legpuzzel met 10.000 stukjes zijn gegaan naar een set van 200 grote, logische blokken die precies hetzelfde resultaat geven, maar veel sneller te bouwen zijn.

Technische samenvatting

Probleemstelling

Polypeptiden (poly-N-gesubstitueerde glycines) zijn synthetische peptidomimetic polymers waarbij de zijketens aan het stikstofatoom van het amide-ruggegraat zijn gebonden in plaats van aan de $\alpha$-koolstof, zoals bij natuurlijke peptiden. Deze structuurvoet geeft polypeptiden unieke eigenschappen, zoals een gebrek aan ruggegraat-waterstofbruggen en een trage cis/trans-isomerisatie van de $\omega$-dihedrale hoek. Hoewel er groeiende interesse is in polypeptiden voor biomaterialen, ontbreekt er tot nu toe een robuust coarse-grained (CG) model dat compatibel is met het moderne MARTINI 3-kader.

Bestaande modellen zijn vaak beperkt tot MARTINI 2 of zijn op maat gemaakt (bespoke), wat de overdraagbaarheid beperkt. Bovendien is de cis/trans-isomerisatie een zeldzame gebeurtenis met een hoge activeringsbarrière (~63 kJ/mol), wat het moeilijk maakt om deze te convergeren in conventionele all-atom (AA) simulaties binnen haalbare tijdschalen. Dit beperkt de mesoscale simulatie van polypeptidenstructuur, zelfassemblage en interacties met membranen.

Methodologie

De auteurs hebben een nieuw, bottom-up CG-krachtenveld ontwikkeld dat volledig compatibel is met MARTINI 3 en 19 veelgebruikte residuen dekt. De aanpak omvatte de volgende stappen:

1. All-Atom (AA) Referentie Simulaties:

   • Er werden uitgebreide AA-simulaties uitgevoerd met GROMACS en het MoSiC-CGenFF-NTOID-krachtenveld.
   • Om de trage cis/trans-isomerisatie van de $\omega$-dihedrale hoek grondig te sample, werd Parallel Bias Metadynamics (PBMetaD) toegepast. Dit zorgde voor geconvergeerde steekproeven van de conformationele ruimte.
   • Systemen varieerden in ketenlengte (n = 4 tot 20) en omvatten homopolymeer en heterotrimers.

2. CG Mapping Strategie:

   • Ruggegraat: Er werd een "center-of-geometry" (COG) mapping gebruikt in plaats van center-of-mass, om de moleculaire volume te behouden. De ruggegraat werd gemapt naar een uniforme P2-bead (onafhankelijk van secundaire structuur), wat compatibel is met MARTINI 3. Verschillende mapping-schemata werden getest; Mapping 1 (alleen zware atomen) bleek het beste te werken voor unimodale verdelingen.
   • Zijketens: De zijketens werden gemapt volgens de standaard MARTINI 3 bead-definities (polaire, apolaire, geladen, aromatische types). Voor een subset van residuen werden aanpassingen gedaan aan de waterstofatoom-representatie in de SC1-bead om asymmetrische verschuivingen te voorkomen.

3. Parameterisatie:

   • Gebonden interacties: Bonden, hoeken en dihedrale hoeken werden afgeleid via Direct Boltzmann Inversion (DBI) van de verdelingsfuncties uit de PBMetaD-simulaties.
   • Unificatie: De meeste ruggegraat-gebonden parameters werden gemiddeld over de 19 residuen om een uniforme ruggegraat-representatie te creëren. Een uitzondering was de BB–BB+–SC1+ hoek, die sterk residu-afhankelijk was. Voor deze hoek werd k-means clustering gebruikt om de 18 residuen in 5 clusters te verdelen, elk met specifieke parameters.
   • Niet-gebonden interacties: Deze werden grotendeels overgenomen uit de standaard MARTINI 3-parameterbibliotheek, zonder aanpassing, omdat de mapping alleen standaard bead-types gebruikte.

4. Validatie:

   • Het model werd gevalideerd tegen AA-referenties voor:
     • Polymer-schaling (straal van gyratie, $R_g$).
     • Vrije-energieoppervlakken van ruggegraat-dihedrale hoeken.
     • Potentiaal van gemiddelde kracht (PMF) voor dimerisatie.
     • Vroege-stadium aggregatie van heterotrimers (bijv. FXF, FKF).

Belangrijkste Bijdragen

• Eerste MARTINI 3 Polypeptiden Model: Dit is het eerste CG-krachtenveld dat polypeptiden volledig integreert in het MARTINI 3-ecosysteem, wat interoperabiliteit met lipiden, eiwitten en andere biomoleculen mogelijk maakt.
• Geautomatiseerde Workflow: Alle parameters en workflows zijn geïntegreerd in de open-source martinize2-tool (via GitHub), wat geautomatiseerde generatie van topologieën en structuren voor de gemeenschap mogelijk maakt.
• Efficiëntie: Het model biedt een 57-voudige toename in computationele efficiëntie ten opzichte van all-atom simulaties, waardoor simulaties op nanometer- en microseconden-schalen haalbaar worden.
• Residu-specifieke Clustering: De innovatieve toepassing van k-means clustering voor de BB–BB+–SC1+ hoek balanceert transferabiliteit met de noodzaak om residu-specifiek gedrag (vooral gerelateerd aan cis/trans-isomerisatie) vast te houden.

Resultaten

• Structuur en Schaling: Het CG-model reproduceerde de straal van gyratie ($R_g$) van korte ketens (N ≤ 10) goed (afwijking < 30%) voor de meeste residuen. Voor langere ketens en geladen residuen (zoals NAsp) werden kleine afwijkingen waargenomen, wat toegeschreven wordt aan de vereenvoudigde ladingsrepresentatie in CG-modellen.
• Conformationele Vrijheid: De vrije-energieoppervlakken van de ruggegraat-dihedrale hoeken ($\psi_1, \psi_2$) in het CG-model kwamen sterk overeen met de geconvergeerde PBMetaD AA-referenties, wat aantoont dat het model de complexe cis/trans-dynamiek correct vastlegt.
• Interacties en Aggregatie:
  • De dimerisatie-PMF voor hydrofobe, polaire, kationische en aromatische tetramers toonde goede overeenstemming met AA-resultaten.
  • Bij simulaties van 50 heterotrimers (FXF, FKF, KFF, XFF) reproduceerde het model de experimenteel waargenomen trend waarbij FXF lineaire nanodraden vormt, terwijl de anderen globulaire structuren vormen. Hoewel het CG-model systematisch iets grotere clusters voorspelde (een bekend effect van MARTINI), werd de sequentie-afhankelijke aggregatietrend correct vastgelegd.

Significantie

Deze studie legt een fundamentele basis voor het simuleren van grote polypeptiden-systemen op mesoscale. Door de compatibiliteit met MARTINI 3 kunnen onderzoekers nu:

• De zelfassemblage van polypeptiden in complexe omgevingen (zoals membranen) bestuderen.
• De rational design van next-generation functionele materialen versnellen.
• De dynamiek van sequentie-gedefinieerde polymeren onderzoeken op tijdschalen die voorheen ontoegankelijk waren voor all-atom simulaties.

Het werk verlegt de grenzen van computationele materialwetenschap voor synthetische biomimetica en biedt een robuust, openbaar beschikbaar platform voor de verdere ontwikkeling van polypeptiden-toepassingen in geneesmiddelenlevering en nanotechnologie.