Integrazione della flessibilità e della dinamica della struttura della proteina bersaglio nella scoperta di farmaci computazionali utilizzando l'analisi di docking basata su ensemble

La nostra ricerca si concentra sull'applicazione di tecniche computazionali per progettare farmaci più efficaci, con l'obiettivo di accelerare la scoperta di farmaci e, in ultima analisi, migliorare i risultati del trattamento e la qualità della vita del paziente.

L'attuale progettazione di farmaci assistita da computer spesso trascura la flessibilità delle proteine bersaglio. Il protocollo discusso affronta questa lacuna incorporando più conformazioni proteiche derivate dalla simulazione dinamica molecolare. La progettazione di farmaci basata su ensemble migliora l'accuratezza considerando la flessibilità delle proteine. Il nostro obiettivo è integrare l'intelligenza artificiale per una scoperta di farmaci più rapida, personalizzata ed efficace in futuro.

[Presentatore] Per iniziare, avvia il software Avogadro su un sistema informatico. Per l'analisi dei cluster, digitare il comando fornito sullo schermo. Quando richiesto, digitare 1 per il gruppo proteico per calcolare l'adattamento dei minimi quadrati e la deviazione quadratica media, o RMSD. Quindi digitare nuovamente 1 per l'output del sistema. Aprire il file cluster-size.xvg. Se il numero di cluster è basso, aumentare il valore limite RMSD. In alternativa, se il numero è alto, diminuire il limite. Eseguite l'analisi del taglio della tolleranza con i comandi visualizzati sullo schermo utilizzando valori di cutoff RMSD diversi in base alle esigenze. Ora apri il software Chimera e cerca cluster.pdb. Clicca su Regali e Pubblicazione 1, silhouette, nastro arrotondato. per la rappresentazione visiva. Fare clic in sequenza su Seleziona, Catena, No ID, quindi su cluster.pdb, #10, quindi selezionare e invertire i modelli selezionati. Vai su Azioni, quindi premi Atomi/Legami e fai clic su Elimina per isolare la catena. Ora seleziona File. Fare clic su Salva PDB. Assegnare al file il nome cluster1.pdb e premere Salva per salvare il file. Per l'ancoraggio basato su ensemble, avviare il software dello strumento Autodock per aprirlo. Inserire i file cluster1.pdb e ligand.pdb in una nuova cartella. Ora fai clic su File, Preferenze e Imposta. Nel popup, incolla l'indirizzo della cartella nel campo della directory di avvio e fai clic su Imposta. Fare clic sull'icona della cartella blu, scegliere cluster1.pdb e fare clic su Apri. Vai su Modifica, quindi premi Cariche, Aggiungi cariche di Kollman e fai clic su OK. Fai clic su Griglia, Macromolecole, Scegli. Selezionare il cluster 1 nella casella Scegli macromolecole, quindi premere Seleziona molecole. Fare clic su OK per generare un file di macromolecola AutoDock4 modificato. Salvarlo come cluster1.pdbqt. Quindi, svuotare l'area di lavoro facendo clic su Modifica, quindi premere Elimina ed eliminare tutte le molecole prima di fare clic su Continua, quindi premere Ligando, Input, Apri. Quando viene visualizzato un file Ligand per la cartella Autodock4, selezionare ligand.pdb e fare clic su Apri, quindi su OK. Ora scegli Ligando, Albero di torsione e Rileva radice per definire la flessibilità torsionale del ligando. Vai su Ligand, Output, Salva come PDBQT e salva la cartella Formatted Autotors Molecules come ligand.pdbqt. Dopo aver svuotato l'area di lavoro, aprire il file cluster1.pdbqt facendo clic su Griglia, Macromolecole e Apri, quindi premere Sì e OK. Passare nuovamente a Griglia, premere Imposta tipi di mappa e scegliere Apri ligando. Selezionare e aprire ligand.pdbqt. Ora vai all'opzione Grid Box sotto Grid. Nella casella Opzioni griglia, impostare il numero di punti nelle dimensioni X, Y e Z su 120 e la spaziatura su 0,375 angstrom. Lasciare le impostazioni centrali come predefinite. Quindi fare clic su File e chiudere Salvataggio corrente. Vai su Griglia, Output e premi Salva GPF. Quando viene visualizzato il file di output dei parametri della griglia, immettere grid.gpf come nome del file e fare clic su Salva. Quindi, fai clic su Esegui ed Esegui AutoGrid. Nella scheda Nome file parametro, fare clic su Sfoglia. Aprire il file grid.gpf. Ora sfoglia il percorso del programma. Cerca autogrid4.exe e fai clic su Apri e avvia. In sequenza, fare clic su Docking, seguito da Macromolecole e Imposta nomi file rigidi. Quando viene visualizzato un file di macromolecole PDBQT, selezionare cluster1.pdbqt e fare clic su Apri. Scegliere il ligando dal menu Docking. Quando viene visualizzata la casella Scegli ligandi, selezionare Ligando e fare clic su Seleziona ligando, quindi premere Accetta. Ora vai su Algoritmo genetico da Docking. Quando viene visualizzata la casella Parametri algoritmo genetico, impostare GA Runs su 100 e fare clic su Accetta. Fare clic su Docking, Output, Lamarckian GA 4.2. Quando viene visualizzato un file di output dei parametri di ancoraggio GALS di Autodock4.2, denominarlo docking.dpf e fare clic su Salva. Ora premi Esegui ed Esegui Autodoc. Verrà visualizzata la casella Esegui Autodoc. In Nome file di parametri, fare clic su Sfoglia. Quando viene visualizzato un file di parametri autodock4, selezionare docking.dpf e fare clic su Apri. In Nome percorso programma, fare clic su Sfoglia. Verrà visualizzato un file autodock4. Cerca autodock4.exe e fai clic su Apri, quindi su Avvia. Elimina tutte le molecole come dimostrato in precedenza e ripeti il processo per tutti i file cluster. Sono state ottenute la struttura chimica e la rappresentazione strutturale 3D della flavokawaina B e del lisozima allo stato iniziale prima della simulazione dinamica molecolare. L'energia totale della struttura proteica è risultata stabile durante la simulazione e la deviazione quadratica media si è stabilizzata dopo 20 nanosecondi. La fluttuazione quadratica media ha rivelato un'elevata flessibilità nelle regioni comprese tra i residui da 40 a 50, da 60 a 80 e da 100 alla fine. Un totale di 15 cluster strutturali sono stati ottenuti dal clustering basato sulla deviazione quadratica media di 10.001 sistemi di traiettoria, con il cluster più grande contenente 5.818 membri. Le conformazioni sovrapposte di tutti gli ammassi hanno mostrato variazioni strutturali visibili tra le traiettorie. L'aggancio molecolare della flavokawaina B con le strutture rappresentative dei primi quattro cluster ha mostrato un legame coerente nello stesso sito in tutte le conformazioni, con il cluster 2 che ha mostrato l'energia di legame più bassa di -29,37 kJ/mol. La mappatura elettrostatica della superficie ha confermato siti di legame identici in tutte le conformazioni di cluster con flavokawaina B annidata nella stessa regione tascabile. L'analisi dettagliata dell'interazione ha mostrato che il legame con la flavokawaina B era stabilizzato da diversi residui circostanti, tra cui alanina 31, glutammina 35, leucina 56, acido gamma-carbossiglutammico 57, isoleucina 58, alanina 95, isoleucina 98, triptofano 108, valina 109, alanina 110, triptofano 111 e arginina 114.