26 settembre 2025
Questo articolo presenta un protocollo per la costruzione di un modello farmacoforico di consenso integrando caratteristiche molecolari da più ligandi. Questo metodo è applicabile alle attività di scoperta di farmaci mirati a qualsiasi bersaglio biologico con conformazioni note legate al ligando, consentendo l'identificazione di caratteristiche di interazione chiave per lo screening virtuale e la progettazione razionale di farmaci.
La nostra ricerca esplora la modellazione farmacoforica con ampie librerie di ligandi. In questo studio, ci concentriamo su SARS-CoV-2 Mpro per identificare i modelli strutturali che guidano la scoperta di farmaci. I recenti progressi nella modellazione farmacoforica utilizzano set di dati più ampi e diversificati, che migliorano l'accuratezza nell'identificazione dei principali alimentatori molecolari per la scoperta di farmaci.
Le simulazioni dinamiche molecolari e l'apprendimento automatico sono attualmente utilizzati per far progredire la modellazione dei farmacofori e migliorare i risultati della scoperta di farmaci. Una sfida chiave è la gestione di set di dati di ligandi ampi e diversificati, migliorando al contempo l'accuratezza dell'estrazione futura e riducendo i costi computazionali. Abbiamo dimostrato che le costose librerie di ligandi migliorano l'accuratezza futura dei farmacofori, consentendo previsioni più affidabili e una migliore identificazione di potenziali farmaci candidati.
Per iniziare, avvia il software PyMOL. Allineare tutti i complessi di ligandi proteici sul software. Per ogni complesso allineato, estrarre il conformatore del ligando e salvarlo come file separato in formato SDF.
Carica ogni file ligando individualmente per coltivarlo, utilizzando l'opzione Carica funzionalità. Dopo il caricamento, utilizzare l'opzione Salva sessione per scaricare il file del farmacoforo in formato JSON. Archivia tutti i file JSON del farmacoforo scaricati in un'unica cartella per un utilizzo successivo nell'ambiente Google Colab.
Avvia Google Colab in un browser web e crea un nuovo taccuino. Installa Conda e PyMOL, quindi esegui la cella facendo clic sull'icona di riproduzione o premendo Maiusc e Invio per verificare la corretta esecuzione. Se completata correttamente, verrà visualizzata una barra verde sopra la cella.
Ora, installa il pacchetto Python CANFAR e importa i moduli richiesti. Esegui la cella facendo clic sull'icona di riproduzione per confermare l'avvenuta installazione. Al termine dell'esecuzione, verrà visualizzato un messaggio di conferma.
Creare quindi una nuova cartella per l'archiviazione dei file JSON del farmacoforo. Apri la cartella, fai clic con il pulsante destro del mouse al suo interno e seleziona Carica per caricare i file JSON. Per analizzare e consolidare le funzionalità del farmacoforo, estrarre le funzionalità da tutti i file caricati.
Fare clic sull'icona di riproduzione per eseguire la cella e generare un frame di dati consolidato contenente le caratteristiche di ogni singolo ligando. Visualizzare tutti i descrittori del farmacoforo, quindi eseguire la cella per visualizzare i modelli di clustering. Ora, salva il modello del farmacoforo nel formato PyMOL.
Esegui la cella per generare un file pse, quindi salva il modello come file JSON eseguendo la cella di codice corrispondente. Genera output e dendogrammi raggruppati di funzionalità utilizzando il codice fornito. Fare clic su Riproduci ed esportare i dati del farmacoforo di consenso in un file CSV eseguendo il comando di esportazione.
Utilizzare il modello di farmacoforo di consenso generato per lo screening virtuale. 100 complessi principali di proteasi o Mpro co-cristallizzati con diversi inibitori non covalenti sono stati allineati per visualizzare la diversità del ligando all'interno della tasca di legame. Ogni ligando è stato estratto individualmente dai complessi allineati per generare file separati per lo screening virtuale.
I ligandi estratti sono stati caricati sul server pharmit e utilizzati per generare file JSON strutturati per la modellazione di farmacofori. I file JSON hanno prodotto 1.450 caratteristiche farmacoforiche, raggruppate in 110 cluster, tra cui 23 aromatici, 30 accettori di legami idrogeno, 16 donatori di legami idrogeno, 36 cluster idrofobici e cinque anionici. Per costruire il modello di consenso, sono stati mantenuti solo i cluster più densamente popolati, mentre i cluster di anioni sono stati esclusi, a causa di membri insufficienti.
Il modello farmacoforico di consenso finale consisteva in 11 caratteristiche, comprendenti tre aromatici, per accettori di legami idrogeno, due donatori di legami idrogeno e due caratteristiche idrofobiche. Dopo aver rimosso la caratteristica aromatica meno rappresentativa, circa uno, sono stati identificati due composti candidati. La struttura chimica bidimensionale del composto 101 267 741 è stata visualizzata per confermarne la compatibilità con il modello farmacoforico.
La mescola 101 267 741 si inserisce profondamente nella tasca di attacco Mpro, occupando le tasche secondarie S1 e S2 in modo più efficace rispetto al ligando di riferimento 38a, che si diffonde verso la tasca S1. Il composto 101 267 741 ha formato 11 interazioni intermolecolari, sette legami idrogeno e quattro contatti idrofobici.
Questo articolo presenta un protocollo per la costruzione di un modello farmacoforo consensuale, integrando caratteristiche molecolari provenienti da più ligandi, con un focus specifico sulla Mpro di SARS-CoV-2. Questo metodo potenzia gli sforzi di scoperta di farmaci identificando caratteristiche chiave delle interazioni per lo screening virtuale e la progettazione razionale di farmaci.
La modellizzazione del farmacoforo consensuale mediante l'uso di librerie estese di ligandi consente una maggiore affidabilità predittiva nella fase iniziale della scoperta di farmaci, in particolare per bersagli come la Mpro di SARS-CoV-2. L'integrazione di conformazioni legate a ligandi diversi riduce il bias del modello e supporta uno screening virtuale robusto, influenzando direttamente l'identificazione dei composti promettenti e la selezione del portafoglio. Questo approccio semplifica la selezione razionale di candidati per bersagli dotati di dati strutturali ricchi, migliorando l'efficienza della ricerca e sviluppo aziendale.
Questo protocollo di modellizzazione del farmacoforo consensuale è applicabile dalla fase iniziale della scoperta fino all'identificazione del composto leader, sfruttando la biologia strutturale e l'informatica per guidare lo screening virtuale e la priorizzazione dei candidati.