Simulation et échantillonnage basés sur la structure des mouvements des protéines du facteur de transcription le long de l’ADN, du passage à l’échelle atomique à la diffusion à grain grossier

Résumé

L’objectif de ce protocole est de révéler la dynamique structurelle de la diffusion unidimensionnelle des protéines le long de l’ADN, en utilisant une protéine du domaine WRKY du facteur de transcription végétale comme système exemplaire. Pour ce faire, des simulations de dynamique moléculaire atomistique et à gros grains ainsi que des échantillonnages informatiques approfondis ont été mis en œuvre.

Résumé

Le glissement unidimensionnel (1-D) de la protéine du facteur de transcription (TF) le long de l’ADN est essentiel pour faciliter la diffusion du TF afin de localiser le site d’ADN cible pour la régulation génétique. La détection de la résolution de la paire de bases (pb) du glissement tf ou du pas sur l’ADN est encore expérimentalement difficile. Nous avons récemment effectué des simulations de dynamique moléculaire (MD) entièrement atomiques capturant un pas spontané de 1 pb d’une petite protéine TF de domaine WRKY le long de l’ADN. Basé sur le chemin de pas WRKY de 10 μs obtenu à partir de telles simulations, le protocole montre ici comment effectuer des échantillonnages conformationnels plus étendus des systèmes TF-ADN, en construisant le modèle d’état de Markov (MSM) pour le pas de protéine 1-bp, avec différents nombres de micro- et macro-états testés pour la construction MSM. Afin d’examiner la recherche par diffusion 1D processive de la protéine TF le long de l’ADN avec la base structurelle, le protocole montre en outre comment mener des simulations MD à grains grossiers (CG) pour échantillonner la dynamique à longue échelle du système. De telles modélisations et simulations CG sont particulièrement utiles pour révéler les impacts électrostatiques protéine-ADN sur les mouvements de diffusion processifs de la protéine TF au-dessus de dizaines de microsecondes, par rapport aux mouvements de pas de protéines submicrosecondes à microsecondes révélés par les simulations de tous les atomes.

Introduction

Les facteurs de transcription (TF) recherchent l’ADN cible pour lier et réguler la transcription des gènes et les activités connexes¹. Mis à part la diffusion tridimensionnelle (3D), la diffusion facilitée de TF a été suggérée comme étant essentielle pour la recherche d’ADN cible, dans laquelle les protéines peuvent également glisser ou sauter le long de l’ADN unidimensionnel (1D), ou sauter avec transfert intersegmental sur l’ADN 2,3,4,5,6,7.

Dans une étude récente, nous avons mené des dizaines de microsecondes (μs) de simulations de dynamique moléculaire d’équilibre (MD) d’équilibre tout atome sur une plante TF - la protéine du domaine WRKY sur l’ADN⁸. Un pas complet de 1 bp de WRKY sur l’ADN poly-A en quelques microsecondes a été capturé. Les mouvements de la protéine le long du sillon de l’ADN et la dynamique de rupture-reformage des liaisons hydrogène (HB) ont été observés. Bien qu’une telle trajectoire représente un chemin échantillonné, un paysage global de pas de protéines fait encore défaut. Ici, nous montrons comment étendre les échantillonnages informatiques autour de la voie d’étape protéique initialement capturée avec le modèle d’état de Markov (MSM) construit, qui ont été largement mis en œuvre pour simuler une variété de systèmes biomoléculaires impliquant des changements conformationnels substantiels et une séparation à l’échelle de temps 9,10,11,12,13,14,15,16, 17,18,19. Le but est de révéler l’ensemble conformationnel et les états méta-stables de la diffusion de la protéine TF le long de l’ADN pour une étape cyclique.

Alors que la simulation MD ci-dessus révèle la résolution atomique des mouvements des protéines pour 1 pb sur l’ADN, la dynamique structurelle de la diffusion processive à long terme du TF le long de l’ADN à la même haute résolution est difficilement accessible. La réalisation de simulations de DM à grains grossiers (CG) au niveau des résidus est toutefois techniquement accessible. L’échelle de temps de simulation CG peut être effectivement étendue à des dizaines ou des centaines de fois plus que les simulations atomiques ^{20,21,22,23,24,25,26,27,28,29}. Ici, nous montrons les simulations CG réalisées en implémentant le logiciel CafeMol développé par Takada lab³⁰.

Dans le protocole actuel, nous présentons d’abord les simulations atomiques de la protéine du domaine WRKY le long de l’ADN poly-A et la construction MSM, qui se concentrent sur l’échantillonnage des mouvements de pas de la protéine pour seulement 1 pb le long de l’ADN. Ensuite, nous présentons la modélisation CG et les simulations du même système protéine-ADN, qui étendent l’échantillonnage informatique à la diffusion processive protéique sur des dizaines de bps le long de l’ADN.

Ici, nous utilisons le logiciel GROMACS^31,32,33 pour effectuer des simulations MD et MSMbuilder³⁴ pour construire le MSM pour des instantanés conformationnels échantillonnés, ainsi que pour utiliser VMD³⁵ pour visualiser les biomolécules. Le protocole exige que l’utilisateur puisse installer et implémenter le logiciel ci-dessus. L’installation et la mise en œuvre du logiciel CafeMol³⁰ sont alors nécessaires pour réaliser les simulations CG MD. D’autres analyses des trajectoires et de la visualisation sont également effectuées dans VMD.

Protocole

1. Construction du modèle d’état de Markov (MSM) à partir de simulations md atomiques

1. Voie spontanée de pas protéique et collecte des structures initiales
   1. Utilisez une trajectoire MD⁸ de 10 μs entièrement atome obtenue précédemment pour extraire 10000 images uniformément d’un chemin de pas « vers l’avant » de 1 bp (c’est-à-dire une trame pour chaque nanoseconde). Le nombre total de trames doit être suffisamment grand pour inclure toutes les conformations représentatives.
   2. Préparez le chemin de transition avec 10000 images dans VMD en cliquant sur Fichier > Enregistrer les coordonnées, tapez protéine ou nucléique dans la zone Atomes sélectionnés et choisissez des images dans la zone Cadres, cliquez sur Enregistrer pour obtenir les images nécessaires.
      REMARQUE: Une trajectoire de simulation MD de 10 μs tout atome précédemment obtenue (appelée « trajectoire de pas en avant » ici) pour wrky marchant à 1 bp sur un ADN poly-A homogène de 34 pb⁸ a été utilisée comme voie initiale pour lancer d’autres échantillonnages conformationnels. Notez que dans la plupart des pratiques, cependant, un chemin initial est construit, en effectuant des simulations MD pilotées ou ciblées, ou en mettant en œuvre des méthodes générales de génération de chemins, etc. ^36,37,38,39.
   3. Alignez l’axe long de l’ADN de référence (de la structure cristalline) sur l’axe des x et définissez le centre de masse initial (COM) de l’ADN complet de 34 pb à l’origine de l’espace de coordonnées pour faciliter une analyse plus approfondie des données. Pour ce faire, cliquez sur Extensions > console Tk dans VMD, puis tapez dans la fenêtre de commande de la console Tk :
      source rotate.tcl
      Le script tcl se trouve dans le fichier supplémentaire 3.
   4. Calculez ensuite la distance quadratique moyenne (RMSD) de l’épine dorsale protéique en alignant l’ADN central de 10 bp (A 14 à 23 et T 14' à 23') à celui de la structure cristalline⁴⁰, et le RMSD représente des mesures géométriques des systèmes (voir figure 1A). Pour ce faire, cliquez sur VMD > Extensions > Analysis >'outil de trajectoire RMSD et tapez les noyaux et les résidus 14 à 23 et 46 à 55 dans la zone de sélection d’atomes, cliquez sur Aligner , puis sur RMSD pour calculer les valeurs RMSD.
   5. Calculer le degré de rotation de la protéine autour de l’ADN Θ(t) sur le plan y-z dans MATLAB en tapant la commande
      rad2deg(atan(z/y))
      avec le positionnement angulaire initial défini comme Θ(0)=0, comme conduit précédemment⁸.
   6. Tapez la commande suivante dans MATLAB⁴¹ pour utiliser les méthodes K-means 42,43,44 et classez les 10000 structures en 25 clusters en tapant :
      [idx, C]=kmeans( X, 25)
      ici X est une matrice 2D de RMSD et d’angle de rotation de WRKY sur l’ADN. Rassemblez les structures de ces 25 centres de cluster pour d’autres simulations MD.
      REMARQUE: Étant donné que la protéine RMSD échantillonnée par rapport à l’ADN couvre une plage d’environ 25 Å, nous choisissons 25 grappes pour avoir une grappe par angstrom.
2. Réalisation de^{la 1 ère} série de simulations MD et des paramètres de simulation
   1. Construisez des systèmes atomistiques pour les 25 structures en utilisant le logiciel GROMACS 5.1.2³² sous le champ de force parmbsc1⁴⁵ et en utilisant le fichier buildsystem.sh du fichier supplémentaire 2 dans le shell.
   2. Effectuez des simulations MD 60-ns pour ces 25 systèmes sous ensemble NPT avec un pas de temps de 2 fs en tapant la commande suivante dans shell :
      gmx_mpi grompp -f md.mdp -c npt.gro -p topol.top -o md.tpr
      gmx_mpi mdrun -deffnm md
3. Regroupement du 1^St trajectoires autour de MD
   1. Supprimez les 10 premières ns de chaque trajectoire de simulation en tapant shell :
      gmx_mpi trjcat -f md.xtc -b 10000 -e 600000 -o newtraj.xtc
      et recueillir les conformations des trajectoires de 25 × 50 ns pour le regroupement afin de préparer les structures d’entrée pour les échantillonnages ultérieurs plus étendus (simulations MD 2^ème tour).
      REMARQUE: Pour réduire l’impact de la trajectoire initiale et permettre l’équilibrage local, 10-ns de la période initiale de simulations ont été supprimés.
   2. Choisissez des paires de distance entre la protéine et l’ADN comme paramètres d’entrée pour la projection de l’analyse en composantes indépendante du temps (tICA)^46,47,48. Utilisez la commande make_ndx dans GROMACS pour ce faire :
      gmx_mpi make_ndx entrée -f.pdb -o index.ndx
      REMARQUE: Ici, les atomes de protéine CA et les atomes lourds (NH1, NH2, OH, NZ, NE2, ND2) du résidu Y119, K122, K125, R131, Y133, Q146, K144, R135, W116, R117, Y134, K118, Q121 qui peuvent former des liaisons hydrogène (HB) avec le nucléotide d’ADN ont été sélectionnés, qui s’associent aux atomes O1P O2P et N6 du nucléotide d’ADN (A14-20, T19-23). Les acides aminés sélectionnés peuvent former des HB stables ou des ponts salins avec l’ADN.
   3. Copiez l’index atom sélectionné ci-dessus à partir du fichier index.ndx dans un nouveau fichier texte (index.dat). Obtenez les informations de paire entre ces atomes par le script python à partir du fichier supplémentaire 1 generate_atom_indices.py et tapez:
      Index generate_atom_indices.py python2.6.dat > AtomIndices.txt
      Cela génère les 415 paires de distance entre la protéine et l’ADN.
   4. Calculez les 415 paires de distance de chaque trajectoire en tapant la commande suivante dans la fenêtre de commande MSMbuilder :
      msmb AtomPairsFeaturizer -out pair_features --pair_indices AtomIndices.txt --top references.pdb --trjs « trajectories/*.xtc » --transformed pair_features --stride 5
   5. Effectuer tICA pour réduire la dimension des données sur les 2 premiers composants ou vecteurs indépendants du temps (tIC) en tapant :
      msmb tICA -i .. /tica_rc_a/tmp/ -o tica_results --n_components 2 --lag_time 10 --gamma 0,05 -t tica_results,h5
      REMARQUE: tICA est une méthode de réduction de dimension qui calcule la valeur propre de la matrice de corrélation décalée dans le temps pour déterminer les degrés de liberté de relaxation les plus lents du système de simulation par l’équation:
      
      où X_i(t) est la valeur de la i-ième coordonnée de réaction au temps t, et X_j(t+Δt) est la valeur de la j-ième coordonnée de réaction au temps t+Δt. est la valeur attendue du produit des trajectoires de simulation globales X_i(t) et X_j(t + Δt). Les directions le long des degrés de liberté de relaxation les plus lents correspondent aux valeurs propres les plus grandes de la matrice de corrélation décalée dans le temps ci-dessus. Ici, 2 tIC semblent être un ensemble minimal pour différencier trois macro-états sur notre construction MSM (abordé plus tard). On peut également calculer le score de quotient de Rayleigh à matrice généralisée (GMRQ)⁴⁹, par exemple, pour explorer un ensemble optimal de composants à utiliser.
   6. Utilisez la commande dans MSMbuilder pour regrouper les jeux de données projetés en 100 clusters par la méthode K-center^43,44 (voir Figure 1B) :
      msmb KCenters -i ./tica_results.h5 -o kcenters_output -t kcenters_output --n_clusters 100.
      Sélectionnez la structure centrale de chaque cluster comme structure initiale pour la 2^ème série de simulations MD. Conservez les informations de simulation des 100 structures simulées, y compris les positions, les températures, les pressions, etc., à l’exception des vitesses.
      REMARQUE: Après la première série de 25 simulations, la mémoire du chemin initial a été réduite, de sorte que nous générons plus de grappes, par exemple 100 grappes, dans la deuxième ronde, pour élargir considérablement les échantillons conformationnels.
4. Réalisation des simulations approfondies du^{2 e} tour de MD
   1. Effectuez des simulations MD 60-ns à partir de ces 100 structures initiales après avoir imposé des vitesses initiales aléatoires à tous les atomes. Ajoutez les vitesses initiales aléatoires en activant la génération de vitesse dans le fichier mdp, c’est-à-dire en changeant le fichier md.mdp gen_vel = no en gen_vel = yes.
   2. Supprimez les 10 premiers ns de chaque simulation comme décrit à l’étape 1.3.1, collectez 2 500 000 instantanés des trajectoires de 100 × 50 ns uniformément pour construire le MSM.
      REMARQUE: Notez que dans la construction ultérieure des macro-états, un petit nombre d’états hors trajectoire avec une population particulièrement faible (~ 0,2%, au bas du plan X-Θ) ont été trouvés. Ces états hors chemin sont classés comme un macro-état lorsque le nombre total de macro-états est défini sur 3 à 6 (Figure 2B). Étant donné qu’un macro-état de population aussi faible ne comprend que 3 trajectoires, qui ont été supprimées à la fin, les résultats montrés dans ce protocole ont en effet été obtenus à partir de 97 × 50 ns trajectoires, avec un total de 2 425 000 images ou instantanés.
5. Regroupement des 2^èmes trajectoires autour de MD
   1. Effectuer tICA pour les trajectoires du 2^e tour comme fait précédemment. Tapez MSMbuilder:
      msmb tICA -i .. /tica_rc_a/tmp/ -o tica_results --n_components 2 --lag_time 10 --gamma 0,05 -t tica_results,h5
   2. Calculer l’échelle de temps implicite pour valider les paramètres du temps de retard de corrélation Δt et les nombres de micro-états (voir Figure 1C),
      
      où τ représente le temps de latence utilisé pour construire la matrice de probabilité de transition (TPM); μ_k(τ) représente la kème valeur propre du TPM sous un temps de décalage de τ. Utilisez le script python du fichier supplémentaire 1 pour ce python BuildMSMsAsVaryLagTime.py -d .. / -f .. /trajlist_num -i 50 -m 1000 -t 10 -n 20 -s 500.
   3. Variez le temps de latence τ et le nombre de micro-états en modifiant les paramètres utilisés ci-dessus :
      python BuildMSMsAsVaryLagTime.py -d .. / -f .. /trajlist_num -i 50 -m 1000 -t 5 10 20 30 40 -n 20 -s 20 200 400 500 800 2000
      REMARQUE: Le système est considéré comme markovien lorsque les courbes d’échelle de temps implicites commencent à se stabiliser avec la séparation d’échelle de temps. Ensuite, choisissez le Dt comme temps de retard de corrélation, et le τ le temps de décalage où l’échelle de temps implicite commence à se stabiliser pour construire MSM.
   4. En conséquence, choisissez un nombre relativement grand (mais pas trop grand) d’états, N = 500, et un temps de retard de corrélation relativement court Δt = 10 ns. Le temps de latence s’est avéré être τ = 10 ns pour construire MSM.
   5. Classez les conformations en 500 clusters (voir Figure 1D) à l’aide de la commande :
      msmb KCenters -i ./tica_results.h5 -o kcenters_output -t kcenters_output --n_clusters 500
6. Construction de VSS
   1. Regroupez les 500 micro-états en 3 à 6 macro-états pour déterminer le nombre de macro-états qui conviennent le mieux selon l’algorithme PCCA+⁵⁰ dans MSMbuilder, en utilisant le script python dans le fichier supplémentaire 1 python msm_lumping_usingPCCAplus.py. Identifier un réseau cinétique réduit de modèles pour les changements conformationnels les plus essentiels des biomolécules, en construisant un petit nombre de macro-états, c’est-à-dire en regroupant cinétiquement des centaines de micro-états comme décrit ci-dessous^17,51.
   2. Cartographier les conformations de haute dimension au X (mouvement de la protéine le long de l’axe long de l’ADN) et l’angle de rotation de la protéine le long de l’ADN pour chaque macro-état, comme décrit aux étapes 1.1.3 et 1.1.4 (par exemple, aucun état avec une population trop faible < 1%; voir la figure 2C). Trouvez ensuite les 3 macro-états qui représentent le mieux le système (Figure 1E). Voir la figure 2D pour des instantanés du mouvement des protéines le long de l’ADN et de l’angle de rotation des protéines autour de l’ADN.
      REMARQUE: Dans des travaux antérieurs générant la voie de pas en avant de la protéine spontanée de 10 μs, nous avons également effectué des simulations MD d’équilibre de 5 x 4 μs pour élargir modérément les échantillons. Nous avons montré la cartographie de la trajectoire avant d’origine (voir la figure 2A à gauche) et d’autres trajectoires d’échantillonnage de 4 μs sur la trajectoire avant effectuée précédemment (voir la figure 2A à droite)⁸. La cartographie des trajectoires originales de 100 × 50 ns (voir figure 2B à gauche)⁸ et des trajectoires de 97 × 50 ns utilisées dans ce travail est montrée (voir la figure 2B à droite).
7. Calcul des temps moyens de premier passage (MFPT)
   1. Effectuez cinq trajectoires Monte Carlo (MC) de 10 ms basées sur le TPM du MSM à 500 micro-états avec le temps de latence de 10 ns défini comme pas de temps de MC. Calculez MFPT⁵² entre chaque paire de macro-états (Figure 3) par le script python dans le fichier supplémentaire 1 python mfpt_msm3.py.
   2. Calculez l’erreur moyenne et l’erreur type du MFPT à l’aide du fichier bash dans le fichier supplémentaire 2, tapez :
      sh mfpt_analysis.bash

2. Effectuer une simulation à gros grains (CG) pour échantillonner la dynamique à long terme

1. Effectuez des simulations CG à l’aide du logiciel CafeMol 3.0³⁰. Consultez les paramètres de simulation CG spécifiés dans le fichier de configuration d’entrée avec une extension .inp, y compris les structures d’entrée, les paramètres de simulation, les fichiers de sortie, etc. Tapez la commande suivante sur le terminal pour exécuter la simulation CG :
   cafemol XXX.inp
2. Spécifiez les blocs suivants dans le fichier d’entrée, chaque bloc commençant par l’étiquette < and ending with >>>>.
   1. Définissez le bloc des noms de fichiers (obligatoire) pour spécifier les répertoires de travail et le chemin d’accès au magasin de fichiers d’entrée/sortie. Tapez ce qui suit pour le bloc de noms de fichiers pour ces simulations :
      <<<< noms de fichiers
      chemin = XXXXX (chemin d’accès)
      filename = wrky (les noms des fichiers de sortie)
      SORTIE psf pdb film dcd rst
      path_pdb = XXXXX (chemin de structure natif d’entrée)
      path_ini = XXXXX (chemin de structure initial d’entrée)
      path_natinfo = XXXXX (chemin d’accès au fichier d’informations natif)
      path_para = XXXXX (chemin des fichiers de paramètres)
      >>>>
      REMARQUE: Comme le modèle Go⁵³ est utilisé dans la modélisation CG, c’est-à-dire que la protéine sera biaisée vers la conformation native, il faut donc définir la structure modélisée comme conformation native. Ici, la structure cristalline d’entrée a été définie comme la conformation native.
   2. Définissez le bloc de contrôle de tâche (obligatoire) pour définir le mode d’exécution des simulations. Tapez la commande suivante :
      <<<< job_cntl
      i_run_mode = 2 (= 2 la simulation de température constante)
      i_simulate_type = 1 (=1 dynamique de Langevin)
      i_initial_state = 2 (=2 signifie que la configuration initiale est une configuration native)
      >>>>
      Sélectionnez les simulations de dynamique Langevin à température constante.
   3. Définissez l’unité et le bloc d’état (obligatoire) pour définir les informations pour les structures d’entrée. Tapez la commande suivante :
      <<<< unit_and_state
      i_seq_read_style = 1 (=1 signifie des séquences de lecture à partir du fichier PDB)
      i_go_native_read_style = 1 (=1 signifie que la structure native provient du fichier PDB)
      1 protéine .pdb (unité et état molecular_type native_structure)
      2-3 ADN ADN.pdb (unité et état molecular_type native_structure)
      >>>>
      REMARQUE: Les fichiers de structure d’entrée initiaux (protéine.pdb et ADN.pdb ici) sont nécessaires. Les structures sont écrites au format pdb. Deux fichiers pdb sont nécessaires ici: l’un est le fichier de structure protéique contenant les coordonnées des atomes lourds de WRKY (unité 1), et l’autre est les coordonnées de l’ADN double brin (ds) de 200 bp (unité 2-3). La protéine est initialement placée à 15 Å de l’ADN.
   4. Définissez le bloc de fonction d’énergie (obligatoire) défini dans le bloc energy_function. Tapez la commande suivante :
      <<<< energy_function
      LOCAL(1) L_GO
      LOCAL(2-3) L_DNA2
      NLOCAL(1/1) ALLER EXV ELE
      ADN NLOCAL(2-3/2-3) ELE
      NLOCAL(1/2-3) EXV ELE
      i_use_atom_protein = 0
      i_use_atom_dna = 0
      i_para_from_ninfo = 1
      i_triple_angle_term = 2
      >>>>
      REMARQUE: Dans les simulations CG, la protéine est grossièrement grainée par le modèle Go⁵³ avec chaque acide aminé représenté par une particule CG placée à sa position Cα. La conformation protéique sera alors biaisée vers la structure native, ou structure cristalline ici, sous le potentiel Go (Figure 4A à gauche). L’ADN est décrit par le modèle 3SPN.2⁵⁴, dans lequel chaque nucléotide est représenté par 3 particules CG S, P, N, qui correspondent respectivement au sucre, au phosphate et à la base azotée (Figure 4A à droite). Les interactions électrostatiques et vdW sont considérées entre différentes chaînes. Les interactions électrostatiques entre la protéine et l’ADN dans la simulation CG sont approchées par le potentiel de Debye-Hückel⁵⁵. L’énergie répulsive vdW prend la même forme que dans le modèle Go.
   5. Définissez le bloc md_information (obligatoire) pour définir les informations de simulation. Tapez la commande suivante :
      <<<< md_information
      n_step_sim = 1
      n_tstep(1) = 500000000
      tstep_size = 0,1
      n_step_save = 1000
      n_step_neighbor = 100
      i_com_zeroing = 0
      i_no_trans_rot = 0
      tempk = 300,0
      n_seed = -1
      >>>>
      La n_tstep est l’étape de simulation. Définissez la tstep_size comme la durée de chaque étape MD, chaque pas de temps CG Cafemol est d’environ 200 fs³⁰, donc chaque étape MD ici est de 200 × 0,1 fs en principe. Mettez à jour la liste des voisins toutes les 100 étapes MD (n_step_neighbor = 100). Réglez la température de simulation sur 300 K. Contrôlez la température en utilisant l’algorithme Verlet de type vitesse pour mettre à jour la structure des protéines avec le thermostat Berendsen⁵⁶.
      REMARQUE : Le n_step_sim est le numéro de bassin du potentiel basé sur le modèle Go, ou le nombre minimal local de la courbe d’énergie. Un potentiel à plusieurs bassins permet à la conformation des protéines d’être biaisée en différentes conformations, de sorte que la conformation des protéines peut passer d’un minimum local à un autre. Ici, seul le modèle Go à bassin unique est utilisé, ce qui signifie une seule conformation biaisée (structure cristalline) pour les protéines dans les simulations. Pendant ce temps, comme il n’y a pas d’interaction de liaison hydrogène protéine-ADN, etc. modélisée dans le contexte CG, les mouvements moléculaires peuvent être échantillonnés encore plus rapidement, c’est-à-dire > 10 fois que dans les simulations atomiques.
   6. Définissez le bloc électrostatique (requis uniquement lorsque l’interaction électrostatique est utilisée) car l’interaction électrostatique est considérée entre différentes chaînes, utilisez donc ce bloc pour définir les paramètres de l’interaction électrostatique en tapant:
      <<<< électrostatique
      cutoff_ele = 10,0
      ionic_strength = 0,15
      >>>>
      Réglez la longueur de Debye dans l’interaction électrostatique à 10 Å, correspondant à l’état de la solution. Réglez la force ionique à 0,15 M, comme à l’état physiologique.

Résultats

Glissement couplé à la rotation ou pas de 1 bp de WRKY de la construction MSM
Toutes les conformations protéiques sur l’ADN sont mappées au mouvement longitudinal X et à l’angle de rotation de la protéine COM le long de l’ADN (voir la figure 3A). Le couplage linéaire de ces deux degrés indique un pas couplé à la rotation de la protéine du domaine WRKY sur l’ADN. Les conformations peuvent être regroupées en 3 macro-états (S1, S2 et S3) dans le MSM. Le...

Discussion

Ce travail porte sur la façon de mener des simulations et des échantillonnages informatiques basés sur la structure pour révéler un facteur de transcription ou une protéine TF se déplaçant le long de l’ADN, non seulement au détail atomique du pas, mais aussi dans la diffusion processive, ce qui est essentiel pour la diffusion facilitée de TF dans la recherche de cible d’ADN. Pour ce faire, le modèle d’état de Markov ou MSM d’une petite protéine du domaine TF WRKY marchant pour 1-bp le long de l’ADN...

matériels

Name	Company	Catalog Number	Comments
CafeMol	Kyoto University		coarse-grained (CG) simulations
GROMACS	University of Groningen Royal Institute of Technology Uppsala University		molecular dynamics simulations software
Matlab	MathWorks		Numerical calculation software
MSMbuilder	Stanford University		build MSM
VMD	UNIVERSITY OF ILLINOIS AT URBANA-CHAMPAIGN		molecular visualization program