Ho provato a costruire un modello 3D sia mediante allineamento pairwise che mediante allineamento multiple templates. Nel secondo caso però le performance sono inferiori. C'è un modo per ottimizzare la struttura? Come o dove trovo i valori di RMSD del modello? grazieeeeee!!!
Quando ho avuto a che fare con degli allineamenti multipli io ho imposto (sulla base del profilo di predizione di struttura secondaria della sequenza target) costraints per le porzioni di struttura secondaria predette.
Per l'RMSD del modello non ho capito bene cosa intendi... l'RMSD si calcola tra coppie di strutture, tu usi + di un template...
provo a interpretare...
Nel caso degli allineamenti multipli Modeller costruisce un set di coordinate atomiche del backbone basata sulla media delle coordinate di ogni templato (dopo averli sovrapposti tutti tra loro), e poi le trasferisce alla sequenza target (con tutte le valutazioni energetiche del caso). Ora, se vuoi ottenere la struttura del modello iniziale da cui è partito Modeller dovresti avviare un job in cui riduci al minimo lo step di ottimizzazione del modello (poca minimizzazione, niente dinamica molecolare); per far questo dovresti aggiungere nel tuo file .py le seguenti righe:
A questo punto forse avrebbe senso fare RMSD incrociati tra: - modello pairwise ottimizzato - modello pairwise "grezzo" - modello multiple templates ottimizzato - modello multiple templates "grezzo"
Per ottimizzare la struttura io proverei prima di tutto a riconsiderare gli allineamenti, prima di giocare sulla costruzione del modello. Hai provato a fare threading per ogni coppia target/templato? Come rankano gli allineamenti prodotti in termini di energia di threading?