FR 2871170 A1 20051209 FR 0406120 A 20040607 PROCEDE DE DETERMINATION DE LA CHARGE MUTATIONNELLE D'UNE BANQUE DE GENES OBTENUE PAR MUTAGENESE ALEATOIRE D'UN GENE D'INTERET ET LES MOYENS POUR SA MISE EN OEUVRE La présente invention se rapporte à un procédé permettant de déterminer la charge mutationnelle introduite dans une banque de gènes à l'issue d'une étape de mutagenèse aléatoire d'un gène d'intérêt. L'invention se rapporte également aux moyens pour la mise en oeuvre de ce procédé comme des kits et des supports destinés à être exploités par des ordinateurs. La mutagènese aléatoire des acides nucléiques permet d'introduire des modifications dans les séquences d'acides aminés des protéines. 1l s'agit d'une stratégie alternative performante à la mutagènese dirigée et au design rationnel des protéines. En effet, en l'absence de toutes données précises sur les caractéristiques d'une protéine comme sa structure, c'est à dire dans la majorité des cas, seule la mutagenèse aléatoire permet d'identifier rapidement des mutants ou des variants d'un gène ayant acquis de nouvelles propriétés. De manière générale, bien qu'il ait été possible d'isoler des variants améliorés dans des banques hyper-mutées (Zaccolo M. et al. (1999) Journal of Molecular Biology. 285:775-783), il est admis que l'introduction d'un trop grand nombre de mutations conduit à une majorité de protéines inactives et limite sévèrement l'identification de variants améliorés alors qu'un taux de mutation trop faible diminue la probabilité d'obtenir une mutation bénéfique. Ainsi plusieurs exemples d'évolution dirigée tendraient à montrer que le taux de mutation le plus efficace pour identifier des variants améliorés serait de 1-5 acides aminés par protéine, soit 2-7 nucléotides par gène (Cherry J.R. et al. (1999) Nature Biotechnology. 17:379-384. / Shafikhani S. et al. (1997) BioTechniquers 23:304-310. / Moore S.C. et al. (1997) Journal of Molecular Biology. 272:336-347.) Il est donc nécessaire de pouvoir contrôler le nombre de mutations introduites dans une banque de gènes par mutagenèse aléatoire. La charge mutationnelle (CM) est définie comme le nombre moyen de mutations par gène au sein d'un banque de gènes mutés. La méthode de référence de l'art antérieur pour déterminer précisément le taux de mutation dans une banque de gènes mutés reste aujourd'hui l'observation statistique par séquençage d'une fraction de gènes de la banque. Cette méthode reste encore lourde à mettre en oeuvre et surtout est très coûteuse. Comme alternative au séquençage, d'autres méthodes pour quantifier le taux de mutation ont été développées comme par exemple le Single Strand Conformation Polymorphism (SSCP) ou le Denaturing Gradient Gel Electrophoresis (DGGE) qui permettent de détecter des différences entre les séquences jusqu'à une base (Cariello et al. (1993) Mutational Research 288:103-112 / Brail et al, (1993) Mutational Research 303:171-175). Néanmoins ces techniques sont plutôt adaptées à une analyse qualitative comme la détection de mutation et la détermination de faible taux d'erreur qu'à une quantification de charge mutationnelle. Une autre méthode de l'art antérieur fait appel à l'analyse des changements dans le profil de restriction de séquences ayant intégrées ou perdues des sites de restrictions particuliers. Cette méthode s'applique donc uniquement à des séquences particulières avec une modification de sites de restriction représentatifs. 1l a aussi été proposé une méthode de suivi de la fréquence de mutation basée sur des tests phénotypiques où un gène rapporteur est soumis à une mutagenèse aléatoire avant clonage. Les fréquences de mutations sont alors simplement analysées en enregistrant la perte de fonctionnalité du produit du gène rapporteur due à l'apparition de mutation dans celui-ci. Cette dernière méthode est cependant limitée à des gènes pour lesquels l'introduction de mutations peut être observée facilement par un test phénotypique. Il est également décrit dans l'art antérieur des tables permettant d'associer à des conditions expérimentales une charge mutationnelle d'une banque de gènes mutés pour un gène d'intérêt donné en fonction de sa longueur. Ces tables sont construites sur la base d'un gène modèle soumis à différentes conditions expérimentales puis par le séquençage d'un nombre statistiquement représentatif de gènes de la banque de gènes mutés. Cependant, ces techniques ne tiennent compte que de la longueur dudit gène d'intérêt. De plus, aucune des méthodes de l'art antérieur n'utilisent de standard interne permettant de tenir compte des variabilités expérimentales susceptibles d'influencer la charge mutationnelle d'une banque. Ce qui signifie que ces méthodes ne sont pas suffisamment fiables et reproductibles pour qu'il soit possible de s'affranchir du séquençage afin de déterminer la charge mutationnelle obtenue dans une banque de gènes mutés. La présente invention a précisément pour but d'offrir dans le domaine de la mutagenèse aléatoire des méthodes fiables pour la prédiction de la charge mutationnelle de banques de gènes mutés. Ce but est atteint selon l'invention grâce à un procédé de détermination de la charge mutationnelle d'une banque de gènes obtenue par mutagenèse aléatoire d'un gène d'intérêt, caractérisé en ce qu'il comprend les étapes suivantes: a) la préparation d'un abaque liant la charge mutationnelle (CM) d'une banque de gènes obtenue par mutagenèse aléatoire d'un gène modèle à la fraction des gènes modèles mutés observées dans ladite banque; b) la mutagenèse aléatoire du gène modèle ayant servi pour la préparation de l'abaque de l'étape (a) et du gène d'intérêt pour obtenir les banques de gènes mutés correspondantes; c) la détermination de la charge mutationnelle (CM) de la banque de gènes obtenue à partir du gène modèle à l'étape (b) à partir de l'abaque établi à l'étape (a) ; d) l'application d'un facteur de correction (CF) à la charge mutationnelle (CM) de la banque de gènes modèles mutés déterminée à l'étape (c) pour déterminer la charge mutationnelle (CM) de la banque de gènes d'intérêt mutés de l'étape (b). Selon une forme toute préférée du procédé de l'invention, la technique de mutagenèse mise en oeuvre à l'étape (b) est identique pour le gène modèle et pour le gène d'intérêt et est également similaire à celle de l'étape (a) pour construire l'abaque à partir du gène modèle. Avantageusement, on entend par technique identique ou similaire, le fait que non seulement la technique est la même mais aussi les conditions de mise en oeuvre. La similarité ou l'identité de la technique de mutagenèse et de ses conditions n'excluent pas le fait que l'homme de métier puisse y apporter des aménagements qui ne soient pas de nature à modifier substantiellement les performances de la technique de mutagenèse. A titre d'exemple, les polymérases utilisées aux étapes (a) et (b) pour mettre en oeuvre une technique de PCR error-prone peuvent ne pas être identiques. De même, l'homme de métier est à même de déterminer un gène modèle pour l'étape b) différent de celui de l'étape a), mais présentant des propriétés substantiellement équivalente pour la mise en oeuvre desdites étapes. Le procédé de l'invention est remarquable en ce qu'il offre une précision suffisante pour constituer une alternative efficace au séquençage qui est très coûteux. Ainsi les étapes (c) et (d) peuvent s'affranchir du séquençage. Dans la suite de la description de l'invention qui suit, on utilisera indifféremment: banque de gènes obtenue par mutagenèse aléatoire d'un gène modèle ou banque de gènes modèles mutés ; et banque de gènes obtenue par mutagenèse aléatoire d'un gène d'intérêt ou banque de gènes d'intérêt mutés On entend par gène au sens de la présente invention, toute séquence d'acide nucléotidique, codante ou non. Une caractéristique essentielle de la présente invention réside dans le fait que les modifications de la séquence du gène modèle peuvent être observées après mutagenèse. On entend par gènes modèles mutés les gènes modèles dont la séquence nucléotidique est modifiée par rapport à celle du gène modèle avant mutagenèse. On entend par gène modèle muté observé un gène modèle dont la modification de séquence est observable. On entend par observable le fait que l'introduction aléatoire de mutation n'est pas analysée par séquençage mais par l'analyse directe ou indirecte de la modification d'une ou plusieurs propriétés des gènes modèles mutés par rapport au gène modèle avant mutagenèse. On entend par analyse indirecte le fait que la modification d'une ou plusieurs propriétés du gène muté est observé par une réaction intermédiaire. A titre d'exemple la modification d'un gène modèle peut être observée par la modification d'une ou plusieurs propriétés de la protéine codée par le gène modèle. La modification d'une ou plusieurs propriétés de la protéine peut consister en la perte ou le gain d'activité(s) caractérisant la protéine. Dans une forme préférée de réalisation de l'invention, l'observation des gènes modèles mutés est réalisée par un test phénotypique. On entend par test phénotypique, un test permettant de détecter la fonction de la protéine codée par le gène modèle, muté ou non. A titre d'exemple de gène dont la fonction peut être détectée par un test phénotypique, on peut citer un gène: - impliqué dans la résistances aux antibiotiques, - impliqué dans l'auxotrophie de la souche hôte, - codant pour une protéine présentant une activité aisément détectable (beta-galactosidase / xylanase / Green Fluorescent Protein (GFP...), codant pour une protéine impliquée dans une réaction immunologique de type antigène-anticorps. Le gène modèle est utilisé dans la méthode de l'invention pour: - établir un abaque tel que défini dans l'étape (a), et - être utilisé en parallèle du gène d'intérêt à l'étape (b) pour tenir compte des variabilités que peuvent avoir les conditions expérimentales sur les résultats de mutagenèse aléatoire, et ainsi servir de standard interne. Contrairement à un témoin classiquement utilisé en biologie, le standard interne défini dans le procédé de l'invention par le gène modèle, n'indique pas uniquement le bon déroulement ou le degré d'avancement d'une expérience, il donne en plus une mesure nécessaire à la réalisation du procédé. Dans une forme particulière de réalisation de l'invention, le gène modèle est un gène modèle dit universel, c'est à dire, qu'il possède des caractéristiques partagées par un grand nombre de gènes. A titre d'exemple, il présente une longueur comprise entre 500pb et 10kpb et un pourcentage en GC compris entre 20 et 80. Préférentiellement, le gène modèle universel présente un pourcentage en GC compris entre 40 et 60. A titre d'exemple d'un gène modèle universel, on peut citer le gène codant pour la protéine GFP, dont l'introduction de mutations peut être réalisée par un test phénotypique selon lequel la fluorescence des colonies sous ultraviolet reflète un phénotype GFP actif. Dans une autre forme préférée de réalisation de l'invention, le gène modèle présente des caractéristiques proches de celles du gène d'intérêt. Ces caractéristiques correspondent avantageusement à des propriétés intrinsèques. 1l est entendu par propriété intrinsèque d'un gène des propriétés inhérentes à la séquence de ce gène. Parmi les propriétés intrinsèques d'un gène envisagées ci-dessus, on peut citer à titre d'exemple et de manière non limitative: - sa longueur, - sa composition en base, la composition en base peut être avantageusement analysée sur un seul brin. - sa fréquence de bases répétées dans la séquence, - la présence de séquences symétriques dans le gène. Une fois le gène modèle sélectionné, une technique de mutagenèse aléatoire est choisie pour mettre en oeuvre la méthode de l'invention. Il est entendu par technique de mutagenèse aléatoire toute technique connue de l'homme de métier. Elle peut correspondre également à une combinaison de plusieurs techniques de mutagenèse. Elle peut être in vivo ou in vitro. Dans une forme de mise en oeuvre préférée, la méthode de l'invention se rapporte à des techniques de mutagenèse in vitro basée préférentiellement sur l'utilisation d'une polymérase ou sur la mise en oeuvre d'une technique d'amplification. A titre d'exemple, la technique de mutagenèse correspond à la PCR error-prone (Cadwell R.0 et al. (1992) PCR Methods Appl. 2:28-33. / Leung D.W. et al. (1989) Techniques 1:11-15.) Une fois le gène modèle sélectionné et la technique de mutagenèse choisie, un abaque permettant de lier la charge mutationnelle d'une banque de gènes obtenue par mutagenèse aléatoire d'un gène modèle à la fraction de gènes modèles mutés observés dans ladite banque est préparé (étape a). Il est entendu par abaque au sens de l'invention tout graphique, toute table ou toute équation permettant de lier de manière biunivoque la charge mutationnelle (CM) d'une banque de gènes obtenue par mutagenèse aléatoire d'un gène modèle à la fraction de gènes modèles mutés observés dans ladite banque. Un abaque est obtenu pour une technique de mutagenèse et un gène modèle donné. Ainsi, dans la forme préférée de l'invention où la modification de la séquence du gène modèle est observable par un test phénotypique, l'abaque obtenu à l'étape (a) permet de lier la charge mutationnelle d'une banque de gènes obtenue par mutagenèse aléatoire du gène modèle à la fraction de gènes présentant un phénotype mutant dans ladite banque. L'abaque peut avoir été réalisé au préalable par l'homme de l'art ou peut être obtenu de la manière suivante: i) Une technique de mutagenèse donnée est appliquée à un gène modèle préalablement sélectionné pour obtenir une banque de gènes modèles mutés. ii) La charge mutationnelle (CM) dans la banque de gènes modèles mutés est ensuite déterminée par séquençage. 1l est observé en parallèle la fraction de gènes modèles mutés dans ladite banque présentant une modification d'une ou plusieurs propriétés par rapport au gène sauvage. Dans une forme préférée de l'invention, un nombre de clones statistiquement représentatif est analysé et le gène modèle muté correspondant est séquencé. Dans une forme préférée de l'invention, entre 10 et 100 clones sont séquencés. Avantageusement, le séquençage est effectué sur 25 clones tandis que l'analyse de la fraction de clones présentant une modification de séquence est effectuée sur au moins 500 clones. Une liaison biunivoque peut ainsi être établie entre la charge mutationnelle d'une banque de gènes obtenue par mutagenèse aléatoire d'un gène modèle et la fraction de gènes modèles mutés observés dans ladite banque. Disposant de cette abaque, la technique de mutagenèse aléatoire utilisée pour construire l'abaque est appliquée en parallèle au gène modèle et à un gène d'intérêt à l'étape (b), le gène modèle de l'étape (b) étant de préférence le même que celui ayant servi à construire l'abaque à l'étape (a). A l'étape (b), la même technique de mutagènèse est appliquée simultanément aux gène modèle et au gène d'intérêt. Le fait d'utiliser le gène modèle avec le gène d'intérêt à l'étape (b) revient à utiliser le gène modèle comme standard interne. Le fait d'utiliser le gène modèle comme standard interne présente l'avantage: - de tenir compte des variations expérimentales qui peuvent avoir lieu lors d'une expérience effectuée à l'étape (b) sur le gène modèle par rapport à celles effectuées antérieurement pour obtenir la relation définie à l'étape (a) ; - de minimiser les variations expérimentales d'un échantillon à un autre, c'est à dire du traitement du gène modèle par rapport au gène d'intérêt lors de l'étape (b). Ces variations expérimentales sont non négligeables sur la charge mutationnelle observée. En effet, il a été démontré que des modifications même minimes des conditions expérimentales peuvent avoir des conséquences importantes sur la charge mutationnelle introduite dans une banque. Le tableau 1 ci-après présente les moyennes, les écart-types et les coefficients de variations C.V. (C.V.=écart- type/moyenne en %) des fractions de clones exprimant une GFP fonctionnelle à l'issue de six expériences répétées de PCR error-prone effectuées avec différentes concentrations de MgCl2 et MnClz (paramètres de contrôle du niveau de mutagènese). A partir d'une expérience répétée à l'identique, le résultat de la mutagènese error-proneest fortement variable avec des coefficients de variations sur les fractions de clones exprimant une GFP fonctionnelle, et donc du taux de mutation compris entre 10 et 40%. Tableau 1 Conditions réactionelles Fraction population clones gfp active [MgCl2] mM [MnCl2] mM Moyenne Ecart Type CV (%) 0 0,722 0,071 9,9 7 0,2 0,411 0,122 29,7 7 0,4 0,091 0,040 43,8 -------- ------- ---- 0,013 35 -5 x,036 ----Le tableau 1 ci-dessus indique la variabilité de la fraction de clones exprimant une GFP fonctionnelle dans des banques de gènes gfp mutés aléatoirement par PCR error-prone lors d'expériences répétées 6 fois à l'identiques. L'étape (b) est donc l'étape clé de la méthode de l'invention. Elle permet de tenir compte des variations expérimentales pouvant avoir lieu lors de la mise en oeuvre d'une technique de mutagenèse aléatoire sur un gène d'intérêt, lesdites variations pouvant influencer la définition de la charge mutationnelle de la banque de gènes ainsi obtenue. De ce fait, la méthode de l'invention permet de s'affranchir du séquençage au niveau des étapes (c) et (d) pour déterminer la charge mutationnelle d'une banque de gènes d'intérêt mutés. Dans une forme toute préférée de l'invention, le même mélange réactionnel est utilisé pour mettre en oeuvre en parallèle la technique de mutagenèse sur le gène modèle et sur le gène d'intérêt pour minimiser encore les variations expérimentales d'un échantillon à un autre. La banque de gènes modèles mutés est ensuite analysée pour déterminer la charge mutationnelle de la banque correspondante (étape c). A l'issue de l'étape (b), la fraction de gènes modèles mutés observés dans la banque obtenue après mutagenèse du gène modèle est déterminée. Grâce à l'abaque définie à l'étape (a), la fraction de gènes modèles mutés observés dans ladite banque permet de définir la charge mutationnelle de cette banque. A l'étape (d), la charge mutationnelle déterminée à l'étape (c) est corrigée en lui appliquant un facteur de correction (CF) pour déterminer la charge mutationnelle introduite dans la banque de gènes obtenue par mutagenèse aléatoire du gène d'intérêt. Ceci est réalisé selon l'équation suivante: CM (gène d'intérêt) = CM (gène modèle) x CF dans laquelle, CF représente le facteur de correction, et CM représente la charge mutationnelle. Le facteur de correction (CF) peut se décomposer en plusieurs sousfacteurs de correction suivant les paramètres sélectionnés comme influençant la charge mutationnelle d'une banque de gènes d'intérêt mutés par rapport à une banque de gènes modèle mutés. Un sous-facteur de correction permet de normaliser l'effet d'un paramètre sélectionné comme capable d'influencer la charge mutationnelle d'une banque par rapport à une autre. La construction du facteur de correction CF s'effectue par la multiplication des différents sous-facteurs que l'on souhaite prendre en compte pour la détermination de la charge mutationnelle de la banque de gènes d'intérêt mutés dans la méthode de l'invention. Les paramètres sélectionnés comme pouvant influencer la charge mutationnelle peuvent correspondre à titre d'exemple à un ou plusieurs des paramètres suivants: - paramètres intrinsèques - paramètres expérimentaux. Il est entendu par paramètres intrinsèques d'un gène ceux qui découlent de la séquence de ce gène comme indiqué précédemment. 1l est entendu par paramètre expérimental tout paramètre capable d'engendrer des différences dans l'introduction de mutations dans le gène modèle par rapport au gène d'intérêt lors de la mise en oeuvre de la technique de mutagenèse. La nature de ces paramètres expérimentaux sont dépendants de la technique de mutagenèse aléatoire utilisée. Ils ne peuvent être calculés qu'après avoir mis en oeuvre la technique de mutagenèse sur le gène modèle et sur le gène d'intérêt à l'étape (b). La considération de ce type de paramètre pour effectuer la correction à l'étape (d) s'ajoute à l'avantage d'utiliser un standard interne afin de tenir compte au mieux des variations expérimentales sur la détermination de la charge mutationnelle d'une banque de gènes. Lors de techniques de mutagenèse basées sur l'usage de la PCR (dont la PCR error-prone), les paramètres expérimentaux peuvent correspondre, de manière non limitative: - au nombre de multiplication du gène durant l'étape de mutagenèse aléatoire, - à la quantité de matrice initiale, au nombre de cycle d'amplification, - à la nature et à la quantité de polymérase utilisée - aux conditions expérimentales utilisées comme la concentration en cations divalents, la concentration et le ratio de concentration des dNTPs. La définition des différents sous-facteurs de correction doit s'effectuer suivant des lois mathématiques prenant en compte la nature des évènements biologiques qui seront normalisés. Par exemple, il est convenu dans l'art antérieur que le nombre de mutation introduites dans un gène est lié linéairement à la longueur du gène (paramètre intrinsèque). Ainsi CF[L], défini comme le sous-facteur 15 de correction lié à la longueur des gènes, est explicité comme suit: CF[L] = L(gène d'intérêt) / L(gène modèle) où L correspond à la longueur des gènes. De même, il est connu de l'art antérieur que les différentes méthodes de mutagenèse n'affectent pas de manière équivalente les différents nucléotides A, T, G ou C. Pour chaque méthode de mutagenèse on peut définir un sous-facteur de correction CF[base] qui tient compte des différences de composition en base entre le gène modèle et le gène d'intérêt vis à vis du biais d'introduction de mutation de la technique de mutagenèse. Pour établir un tel sous-facteur de correction il est nécessaire d'introduire la notion de 30 matrice cible des mutations. En effet, la matrice cible des mutations du gène modèle (Mtarget[X] où X représente le gène modèle) décomposent sur les quatre bases la probabilité de toucher un type de base particulier dans le gène modèle et permet de normaliser la composition en base par rapport au gène modèle utilisé. La matrice cible des mutations du gène modèle est déduite à partir de séquençages des gènes modèles mutés par une technique définie de mutagenèse aléatoire, et identification des mutations apparues. L'analyse des mutations apparues dans les séquences des gènes modèles mutés permet d'établir la matrice des mutations MutX où Mut, répertorie les fréquences des 12 types possibles de mutation, selon la matrice suivante: fA -, T fA --> c fA - Gl Mut A fT, c fA -* G fc,A fc-> T fc,G fG-4A fG-> T fG-*c dans laquelle, fXY représente la probabilité de muter la base X en base Y. La matrice cible de mutation du gène modèle X, Mtarget [X] est alors déduite de la matrice Mut, en sommant les trois fréquences affectant un type de base selon la matrice suivante: AMtarg[X] TMtarg[XJ Mtarget [X] = CMtarg[Xj GMtarg[X] ) dans laquelle, ZMtargpq Efz -aY ou Z est A, T, C ou G. Y=(A, T,C,G)-Z La matrice cible des mutations est dépendante de la nature du gène modèle. Plus un gène modèle sera riche dans un type de base particulier plus la probabilité d'affecter ce type de base avec une méthode définie de mutagenèse aléatoire sera élevée. La matrice cible des mutations est dépendante de la technique de mutagenèse utilisée car 2871170 16 cette dernière est fonction du biais de la méthode de mutagènese aléatoire. Ainsi CF[base], défini comme le sous- facteur de correction lié à la composition en base, est explicité selon l'équation suivante: C.F. [bases AM.gx)xA(g.interêt)+ TMt x]xT(g.interêt)+ CMt.dX]XC(g.interét)+ GMrarg X]XG(g.interêt) A(x) T(x) C(x) G(x) dans laquelle, A(q.interêt), Am, T(g.interêt), T(X), C(g, Cm, int erêt) G(g interét) et G(X) correspondent respectivement à la composition en base A, T, C et G du gène modèle X et du gène d'interêt. De plus, l'impact de certains mutagènes sur un nucléotide donné peut varier significativement selon le type de nucléotide qui lui est adjacent dans la séquence. C'est pourquoi le facteur correctif pourra intégrer lorsque ce sera pertinent un terme de comparaison entre les fréquences des dinucléotides dans le gène d'intérêt par rapport au gène modèle. La fréquence relative des trinucléotides entre gène modèle et gène d'intérêt pourra de même être prise en compte dans le facteur correctif. A titre d'exemple également, dans le cas d'une mutagènese in vitro basée sur la PCR, le taux de multiplication (paramètre expérimental) subit par la matrice va directement influencer le taux de mutation. En effet, l'efficacité de réplication est dépendante de la matrice ADN et du programme PCR utilisé. Dues à des différences de longueur, de composition en base, de la présence de structure secondaire, certains gènes sont plus difficiles à amplifier que d'autres. Ainsi, les éventuelles différences de rendement PCR entre la mutagènese du gène modèle et celle du gène d'intérêt pourront générer des différences dans le taux de 2871170 17 multiplication et ceci même en partant d'une quantité identique de matrice et en réalisant un même nombre de cycle PCR. Le taux de multiplication (d) est calculé d'après l'équation suivante: d = ln(Qf/Qi) /ln2 dans laquelle, Qi et Qf représentent respectivement les quantités molaires de matrice initiale et de produit à l'issue de la PCR. Les quantités initiale de matrice et finale de produit peuvent être déterminées par de nombreux moyens connus de l'homme de métier. Citons par exemple le dosage spectrophotométrique aux W ou le calibrage par rapport à une quantité connue d'ADN. Les paramètres expérimentaux peuvent également correspondre au nombre de générations d'une culture cellulaire lors d'une technique de mutagenèse in vivo . A titre d'exemple, dans le cas d'une mutagenèse in vivo , le nombre de mutations introduites dans les gènes va dépendre du nombre de générations de la culture. L'expérimentateur doit alors prendre soin de l'analyser, par exemple en mesurant la turbidité des cultures. Différentes formes de réalisation de l'invention peuvent être envisagées, où : - au moins un paramètre intrinsèque de chacun des gènes modèle et d'intérêt est considéré pour 30 calculer le facteur de correction, - au moins deux paramètres intrinsèques du gène modèle et du gène d'intérêt sont considérés pour calculer le facteur de correction selon l'équation suivante: CF = CF[ paramètre intrinsèque 1] x CF[ paramètre intrinsèque 2] 10 2871170 18 A titre d'exemple, ces paramètres intrinsèques correspondent à la longueur et à la composition en base des gènes, selon l'équation suivante: CF = CF[L] X CF [base] dans laquelle - CF[L] , le sous-facteurde correction est lié à la longueur des gènes, - CF[base], le sous-facteur de correction est lié à la composition en base des gènes, - au moins un paramètre intrinsèque et au moins un paramètre expérimental sont considérés pour établir le facteur de correction, avec C.F. déterminé selon l'équation suivante CF = CF[paramètre intrinsèque] X CF[paramètre expérimental] - au moins deux paramètres intrinsèques et au moins un paramètre expérimental sont considérés pour établir le facteur de correction, avec CF déterminé selon l'équation suivante CF = CF[paramètre intrinsèque 1] X CF[paramètre intrinsèque 2] x CF[paramètre expérimental] Dans le cas particulier où la technique de mutagenèse correspond à la PCR error-prone et le gène modèle correspond à la GFP (séquence de la Figure 1), différents sous-facteurs de correction correspondent pour: - paramètre intrinsèque 1 représente la longueur des gènes selon l'équation suivante: CF[L] = L(gène d'intérêt)/717 où L est en nombre de paire de base, paramètre intrinsèque 2 représente la composition en base sur les brins d'ADN selon l'équation suivante: 0,506xA(g.interêt) 0,346xT(g.sntet) 0,063x C(g.interêt) 0,085xG(g.interêt) C.F. [base]= + + + 0,34 0,26 0,197 0,203 ou A(g interêt) 1 (g.interêt) C(g.interêt) et G(g interêt) sont en en fraction, - paramètre expérimental représente le taux de multiplication selon l'équation suivante: CF[d] = d(gène d'interêt) / d(gfp) où d est le taux de multiplication calculé selon l'équation suivante: d = ln (Qf /Qi) /ln2 dans laquelle, Qi et Qf représentent respectivement les quantités molaires de la matrice initiale et des produits à l'issue de la PCR. L'invention concerne aussi les utilisations du procédé préalablement défini: 1) Le procédé de l'invention peut être utilisé pour obtenir et contrôler la charge mutationnelle désirée dans une banque de gènes d'intérêt mutés. L'expérimentateur, sur la base d'un abaque lui permettant de lier la charge mutationnelle d'une banque de gènes modèles mutés aux conditions expérimentales de mutagénèse aléatoire, va sélectionner les conditions expérimentales ou une gamme de conditions expérimentales lui permettant de se rapprocher d'une charge mutationnelle désirée pour sa banque obtenue après mutagénèse aléatoire d'un gène d'intérêt. L'expérimentateur va faire cette sélection en tenant compte au moins d'un paramètre intrinsèque des gènes modèles et d'intérêt. L'expérimentateur met alors en oeuvre les étapes (a) à (d) du procédé de l'invention préalablement décrite en utilisant à l'étape (b) les conditions expérimentales préalablement sélectionnées. Il utilise alors la méthode de l'invention pour vérifier que la charge mutationnelle introduite est bien celle désirée. Ainsi, l'utilisation de la méthode de l'invention pour obtenir et contrôler une charge mutationnelle désirée dans une banque de gènes d'intérêt mutés est caractérisée par la mise en oeuvre d'une étape supplémentaire de sélection des conditions expérimentales de mutagenèse aléatoire qui seront utilisées pour mettre en uvre l'étape (b). A titre d'exemple cette utilisation permet de contrôler la charge mutationnelle dans une banque de gènes d'intérêt mutés en préalable à une étude structure-fonction. L'expérimentateur réalise premièrement une banque de gène d'intérêt faiblement muté puis il vérifie, par la méthode de l'invention, que la charge mutationnelle introduite est bien de 1-2 mutations par gène, le nombre de mutation optimale pour ce type d'application (Vartanian J.P., Henry M., Wain-Hobson S. Hypermutagenic PCR involving all four transitions and a sizeable proportion of transversions. Nucleic Acids Res. (1996). 24:2627-2631). Le criblage des variants issus de la banque dans une fonction donnée et le séquençage des clones identifiés positifs permettent ensuite de lier sans ambiguïté la séquence à la fonction étudiée. 2) Le procédé de l'invention peut être utilisé pour contrôler la charge mutationnelle dans une banque de gènes d'intérêt mutés en préalable à une recherche de variants améliorés. L'expérimentateur réalise premièrement une banque de gène d'intérêt muté puis il vérifie, par la méthode de l'invention, que la charge mutationnelle introduite est bien de 3-7 mutations par gène, le nombre de mutation optimale pour ce type d'application (Wan L. et al. (1998) Proc. Natl. Acad. Sci. USA 95:1282512831. / You L et al. (1996)Protein Eng. 9:77-83). L'application d'un nombre optimal de mutation permet d'augmenter la probabilité d'identifier des gènes codant pour des protéines améliorées. 3) Le procédé de l'invention peut être utilisé pour étudier la plasticité protéique c'est à 2871170 21 dire la capacité à incorporer des mutations tout en conservant une activité présentée par la protéine sauvage (Shafikhani S. et al. (1997) BioTechniquers 23:304-310). L'expérimentateur réalise premièrement des banques de gènes d'intérêt plus ou moins mutés puis il détermine, par la méthode de l'invention, les charges mutationnelles associées aux différentes banques. L'analyse des fractions de clones produisant un variant protéique fonctionnel dans les différentes banques permet l'étude de la plasticité du produit du gène d'intérêt en liant ces dernières aux charges mutationnelles précédemment déterminées. 4) Le procédé de l'invention peut être utilisé pour étudier le taux d'évolution d'une protéine d'intérêt, c'est à dire la capacité à acquérir une nouvelle propriété avec une charge mutationnelle donnée (Miura T. et al. (2001) J. theor. Biol. 209:497-502). L'expérimentateur réalise premièrement des banques de gènes d'intérêt plus ou moins mutés puis il détermine, par la méthode de l'invention, les charges mutationnelles associées aux différentes banques. L'analyse des fractions de clones produisant un variant présentant une propriété originale par rapport à la protéine sauvage dans les différentes banques permet d'étudier le taux d'évolution en fonction de la charge mutationnelle. 5) L'invention peut être utilisée pour déterminer la charge mutationnelle optimale dans des banques de gène d'intérêt muté pour générer un lot de gènes parentaux avant la mise en oeuvre de méthodes de recombinaison moléculaire. L'expérimentateur réalise premièrement des banques de gène d'intérêt plus ou moins mutés puis il détermine, par la méthode de l'invention, les charges mutationnelles associées aux différentes banques. Deuxièmement la réalisation de criblages adéquats sur les banques de gènes d'intérêt mutés permettent d'identifier des clones présentant une/des propriété(s) améliorée(s) par rapport à celles présentées par la protéine sauvage. Chacune des banques de gènes d'intérêt mutés génèrent ainsi un lot de gènes optimisés, qui servent alors de gènes parentaux dans une étape de recombinaison. Les différentes banques de gène recombinés sont criblées pour identifier des protéines encore plus améliorées que celles produites par les gènes parentaux ou par le gène sauvage. L'analyse et la comparaison de ces banques, par exemple sur la fréquence d'apparition de clones améliorés ou sur la quantification des améliorations obtenues, permettent de déterminer la charge mutationnelle optimale à introduire initialement dans une banque de gène d'intérêt en préambule à une étape de recombinaison. Toutes les méthodes de recombinaison in vitro: DNA Shuffling (Stemmer W. (1994) Nature 370:389-391), L-Shuffling (WO0009679), StEP (Zhao H. et al. (1998) Nature Biotechnology 16:258-261), RACHiTT (Coco W. et al. (2001) Nature Biotechnology 19:354-359) et in vivo comme la recombinaison dans la levure (Pompon D. et al. (1989) Gene 83:15-24), ainsi que la chimerisation / réparation in vivo (Volkov A. et al. {1999) Nucleic Acids Research 27:e18) sont concernées par l'invention. 6/ Le procédé de l'invention peut être utilisé pour quantifier ou comparer la fidélité d'une (de) méthode(s) de polymérisation d'acides nucléiques. L'expérimentateur réalise une/des banque(s) de gènes modèles à partir des produits issus de polymérisation(s) réalisée(s) dans des conditions expérimentales données, composition en tampon, type de polymerase etc... Dans un deuxième temps il détermine, par la méthode de l'invention, les charges mutationnelles associées aux banques de gène modèle. Ces dernières, après la prise en compte éventuel du nombre de cycle de réplication, permettent d'estimer la fidélité de la méthode de polymérisation, dans les conditions expérimentales utilisées. La présente invention concerne également les kits pour la mise en oeuvre du procédé et de ses différentes applications ou utilisations. Un tel kit comprend: - au moins un gène modèle par exemple sous la forme d'un plasmide; - au moins un abaque permettant de lier la charge mutationnelle d'une banque de gènes obtenue par mutagenèse aléatoire du gène modèle à la fraction de gènes modèles mutés observés dans ladite banque; - les réactifs nécessaires à la mise en oeuvre de la technique de mutagenèse aléatoire comme par exemple une ADN polymérase thermostable, dNTP, magnésium, manganèse, tampon de PCR dans le cas où la technique de mutagenèse correspond à la PCR errorprone; - les réactifs, si nécessaires, pour observer des gènes modèle mutant; - une notice indiquant les sous-facteurs de correction définis pour un paramètre sélectionné comme influençant la charge mutationnelle. Il y a autant d'abaques que de gènes modèles dans ledit kit. Dans une forme de mise en oeuvre où la technique de mutagenèse utilise une étape d'amplification, le kit peut posséder les amorces pour l'amplification du gène modèle. Lorsque le kit est mis en oeuvre pour obtenir et contrôler une charge mutationnelle désirée dans une banque de gènes mutés, il contient en outre au moins un abaque lui permettant de lier la charge mutationnelle (CM) d'une banque de gènes obtenue par mutagénèse aléatoire d'un gène modèle aux conditions expérimentales de mutagénèse aléatoire. Les sous-facteurs de correction fournis 5 seront complétés par une ou plusieurs des données suivantes telles que et de façon non limitative - la longueur du gène d'intérêt, - la composition en base A du gène d'intérêt, d'intérêt, d'intérêt, 15 d'intérêt,la composition en base C du gène - la composition en base G du gène - la composition en base T du gène - le nombre de multiplication de la matrice du gène modèle - le nombre de multiplication de la matrice du gène d'intérêt, dont les valeurs correspondantes sont intégrées par l'expérimentateur, en fonction du gène d'intérêt. Le choix de ces données se fera selon que l'on souhaite tenir compte d'un ou plusieurs paramètres intrinsèque et ou expérimental. Dans le cas où le gène modèle correspond à la gfp et la technique de mutagenèse correspond à la PCR error-prone, les différents sous-facteurs de correction sont ceux indiqués précédemment. Dans une forme préférée, le kit possède plusieurs gènes modèles et plusieurs abaques pour chaque gène modèle de manière à choisir le gène modèle qui ait les propriétés intrinsèques les plus proches du 35 gène d'intérêt. Selon une forme particulière de mise en oeuvre du kit de l'invention, les sous-facteurs de correction sont enregistrés dans un fichier numérique sauvegardé sur un support approprié (CD-Rom, disquette, etc.). De même dans un kit selon l'invention, l'abaque peut être enregistré sous forme d'une base de données numériques sauvegardées sur un support approprié. L'invention se rapporte enfin a un procédé dans lequel l'une au moins des étapes (a), (c) et (d) est mise en oeuvre par des moyens informatiques. Lesdits moyens informatiques peuvent: - lier la charge mutationnelle d'une banque de gènes obtenue par mutagenèse aléatoire d'un gène modèle à la fraction de gènes modèles mutés observés dans ladite banque. être incrémentés par chaque nouvel abaque créé. Ce procédé intègre la mise en oeuvre d'une base de données de différents abaques, chacun défini pour une technique de mutagenèse et un gène modèle donné. - posséder les données pour établir au moins un facteur de correction pour la détermination de la charge mutationnelle d'une banque de gènes d'intérêt mutés. L'expérimentateur devra alors préciser, en fonction de ses travaux, des données numériques relatives à l'une au moins des caractéristiques suivantes, que lesdits moyens informatiques exploitent le gène modèle et ses propriétés intrinsèques - la technique de mutagenèse 10 15 20 25 30 le ou les paramètres sélectionnés comme capables d'influencer la charge mutationnelle d'une banque de gène par rapport à une autre pour définir les sous-facteurs de correction. A titre d'exemple, l'expérimentateur entrera une ou plusieurs des données déjà indiquées précédemment. D'autres avantages et caractéristiques de l'invention apparaîtront dans les exemples qui suivent, où il sera fait référence aux dessins en annexe dans lesquels: La figure 1 représente la séquence du gène 15 de la gfp. La figure 2 représente la relation entre la fraction de clones exprimant une GFP fonctionelle et la charge mutationnelle dans la banque de gènes gfp mutés. La figure 3 représente la matrice des mutations MUtGFP pour le gène modèle gfp et pour la technique de mutagènese aléatoire par PCR error- prone. La figure 4 représente la matrice cible de mutation Marger [gfp] sur le gène gfp pour la technique de mutagenèse aléatoire par PCR error-prone. La figure 5 représente la platicité de TEM- 1 P-lactamase, soit le taux de survie des protéines mutées en fonction de la charge mutationnelle des banques de gènes TEM-1 mutés. La figure 6 représente l'évolution de TEM-1 (3-lactamase, soit la fréquence d'apparition de clones résistant à la cefotaxime en fonction de la charge mutationnelle des banques de gènes TEM-1 mutés. Exemple 1: Détermination de la charge mutationnelle introduite dans une banque de crènes de xvlanase de Thermotoga maritima (xvnA) à l'issue d'une étape de mutagenèse aléatoire grâce à la mise en oeuvre des étapes suivantes: 1) Sélection du gène modèle Dans l'exemple développé, le gène de la Green Fluorescent Protein (gfp) a été choisi comme gène modèle grâce à la facile mise en oeuvre d'un test phénotypique basé sur l'expression du produit du gène gfp: les colonies exprimant une protéine correspondant à la Green Fluorescent Protein ou GFP fonctionnelle présentent une couleur verte aux ultraviolets. 2) Sélection de la méthode de mutagenèse aléatoire Du fait de sa facilité de mise en oeuvre, les mutagenèses aléatoires ont été menées par PCR error-prone comme décrit par Cadwell (Cadwell R.C., Joyce G.F. (1992) PCR Methods Appl. 2:28-33). 3) Obtention d'un abaque permettant de lier la charge mutationnelle d'une banque de gènes modèle qfp mutés aléatoirement à la fraction de gènes modèle mutés observés dans ladite banque Pour déterminer la correspondance entre la charge mutationnelle dans la banque de gènes modèle gfp mutés et la fraction de clones exprimant une GFP fonctionelle, trois expériences de mutagènese aléatoire (A,B et C) par PCR error-prone ont été effectuées. Les réactions de PCR error-prone ont été réalisées dans des volumes de 100uL contenant 10mM de Tris-HC1 pH9, 50rnM de KC1, 0.1% de Triton X-100, 0.2mg/mL de BSA, 20pmol d'amorces universelles pET5' et pET3', 2,5U de Taq polymerase (QBiogen). Les concentrations en cations divalents sont de 4mM MgCl2 pour l'expérience A, 7mM MgCl2 plus 0,3mM MnC12 pour l'expérience B et 7mM de MgC12 plus 0,5mM MnC12 pour l'expérience C. Les concentrations en dNTP utilisées sont 0, 2mM de dATP, 0, 2mM de dGTP, 1mM de dTTP et 1mM de dCTP tandis que lfmol de vecteur a servi comme matrice ADN. Un programme PCR de 94 C pendant 5min; (91 C pendant 30s, 60 C pendant 30s, 72 C pendant 1min.) 25 fois puis 72 C pendant 10 min en guise d'étape d'élongation finale, a été appliqué dans un thermocycler MJ Research PTC-200. Les produits PCR ont alors été purifiés en utilisant le kit de purification sur colonne QIAquick (Qiagen), vérifiés sur gel et les quantités d'ADN obtenues ont été quantifiées par le kit PicoGreen (Molecular Probes) comme décrit par le fournisseur. Les produits de PCR error-prone (lpg) ont été digérés par Eco RI et Nde I (New England Biolabs) pendant 1h à 37 C dans le tampon de réaction adéquat comme décrit par le fournisseur puis purifiés sur colonne QlAquick et insérés dans le vecteur pET26b+ (Stratagene) préalablement digéré avec Eco RI et Nde I. Les cellules E. coui MC1061(DE3) ont été transformées par choc thermique comme décrit par Maniatis (Sambrook J. et al. (1989) Molecular cloning: a laboratory manual Cold Spring Harbor Laboratory Press -New York) avec 1pL des réactions de ligation. Les transformants ont été sélectionnés sur boites LB agar complémentées avec 0,6mg/rnL de kanamycine. Après clonage et transformation, les trois banques sont analysées pour la fraction de clones exprimant une GFP fonctionelle mais également par séquençage pour déterminer les charges mutationnelles dans les trois banques de gènes gfp mutés à l'issu des expériences de mutagenèse A, B et C. Le criblage des clones exprimant une protéine GFP fonctionnelle a été réalisé en identifiant les colonies présentant une fluorescence verte sous banc UV (longueur d'onde d'excitation à 355nm). Parallèlement, le rendement d'insertion des gènes gfp mutés dans les différentes banques a été analysé par PCR sur des colonies choisies aléatoirement. Par banque, le criblage pour la recherche de clone présentant une GFP fonctionnelle a été réalisé sur environ 1000 clones et le rendement d'insertion (nombre de clones présentant un insert par nombre total de clones) a été estimé à partir de 96 PCR sur colonies, ceci de façon à être statistiquement pertinent. Les banques de gènes gfp mutés sont ainsi analysées sur le nombre total de colonie (cfu_total), le nombre de colonies exprimant une GFP fonctionnelle (cfu_fluo) et le rendement d'insertion des gènes gfp mutés dans les banques (Ri). La fraction de clone exprimant une GFP fonctionnelle dans une banque donnée (Fa[GFP]) est alors déduites de l'équation suivante: cfu total x Ri 15 Pour déterminer les charges mutationnelles dans les banques A, B ou C de gènes gfp mutés, 24 clones par banque contenant un insert gfp ont été sélectionnés pour une amplification PCR avec les amorces pET5'/pET3'. Les produits de digestion ont été purifiés sur colonne QIAquick puis séquencés en totalité. Toutes les mutations (modification de base, insertion et délétion) détectées dans les séquences par rapport à celle présentée par le gène gfp sauvage ont été vérifiées sur les chromatogrammes. Pour chacune des trois banques, le séquençage a porté sur 13623 à 15774 pb. Le séquençage d'un nombre suffisamment important de base permet d'obtenir une mesure statistiquement fiable des charges mutationnelles, même dans le cas de la banque contenant les gènes les moins mutés (banque A). 1l est à noter que le ratio entre les clones exprimant une GFP fonctionnelle et ceux exprimant une GFP inactive a été conservé pour la sélection des 24 clones dédiés au Fa[GFP] = cfu fluo séquençage. En effet, il a été décrit qu'à l'issue d'une étape de mutagènese aléatoire les gènes codant pour une protéine active pourraient avoir une charge mutationnelle plus faible que des gènes non sélectionnés (Daugherty P et al. (2000) Proc. Natl. Acad. Sci. USA. 97:2029-2034). Ainsi, la sélection aléatoire pour le séquençage d'un nombre statistiquement faible de clones aurait pu biaiser l'observation de la charge mutationnelle. Le tableau 2 ci-après représente l'analyse de la fraction de clone exprimant une GFP fonctionnelle en fonction de la charge mutationnelle des banques de gfp mutés. Il résume les résultats obtenus sur l'analyse du test phénotypique GFP, tandis que la figure 2 présente le logarithme de la fraction de GFP fonctionnelle en fonction de la charge mutationnelle dans les banques de gfp mutés. Tableau 2 Expérience Conditions réactionelles Nombre de Nombre de Charge mutationeile Fraction GFP PCR-ep [MgCl2) mM [MnCl2j pb séquencées mutations (mutations/gfp) fonctionelle (%) A 4 0 15057 33 1, 57 75,1 ---------------- -------- ------------------- 212 --------------------- -- - 0 15774 9 16 5 ---.3,64 --- 6 - 7 0,5 13623 285 15 _ ^ C 7 A partir des données générées, il est possible d'établir une courbe de régression polynomiale d'ordre 2, représentant l'abaque selon l'équation suivante: ln(Fa[GFP]) = -0, 0083XCM[gfp]2 - 0,1111xCM[gfp] liant la charge mutationnelle d'une banque de gènes obtenue par PCR error-prone du gène modèle gfp (C.M.[gfp]) à la fraction de gènes gfp mutés observés dans ladite banque lors du test phénotypique (Fa[GFP]). La valeur du coefficient de régression (R2) étant de 0,999 (cf figure 2), l'abaque précédemment déterminé permet d'estimer de façon précise la charge mutationnelle dans une banque de gènes gfp mutés à l'issu d'une étape de mutagènese aléatoire par PCR error-prone à partir de la fraction de clone exprimant une GFP fonctionelle. 2871170 31 L'analyse des séquençages sur les gène gfp mutés a permis également d'analyser la matrice de mutation MUtcF,P (figure 3) et d'en déduire la matrice cible de mutation sur le gène gfp, Mtarge,[gfp] (figure 4). 4) Application en parallèle de la mutacenèse aléatoire au crène modèle gfp et au gène d'intérêt xvnA de Thermotoga maritima Le gène xynA de xylanase issue de Thermotoga maritima (numéro d'accession EMBL Z46264) a été soumis à mutagenèse aléatoire par PCR error-prone avec 7rnM MgC12 et 0,3mM MnC12 (conditions B). Le protocole expérimental est identique à celui précédemment décrit, excepté que l'étape d'élongation PCR a été fixée à 3 minutes au lieu de 1 min. Il est à noter que le même mélange réactionnel a servi pour la PCR error- prone sur le gène gfp et pour la PCR errorprone sur le gène xynA. L'utilisation du gène gfp comme standard interne dans le cadre de l'invention doit permettre de s'affranchir des variations expérimentales sur les charges mutationnelles dans les banques de gènes modèle gfp mutés. Les rendements des expériences de PCR error-prone ont été déterminés par comparaison avec une quantité connue de standard ADN après migration de lpL de réaction sur gel agarose 1%. La quantité initiale de vecteur étant de lfmole dans chacune des deux expériences menées, le nombre de multiplication (d) a pu être déterminé à 15 et 12 pour les PCR errorprone gfp et xynA respectivement. 5) Détermination de la charge mutationnelle de la banQue obtenue à partir du crène modèle gfp A l'issue de la PCR error-prone réalisée dans les conditions B, la fraction de gènes gfp codant 2871170 32 pour une protéine fonctionnelle dans la banque a été déterminée à 0.165. En utilisant l'abaque donné par l'équation suivante (équation 4) . ln(Fa[GFP]) = -0,0083xCM[gfp]2 - 0,1111xCM[gfp] liant la charge mutationnelle d'une banque de gènes modèle gfp mutésà la fraction de gènes modèle mutés observés, la charge mutationnelle dans la banque de gfp mutés est simulée à 9,5 mutations par gène. Il est à noter que la charge mutationnelle dans cette même banque a été déterminée par séquençage à 9,64 mutations par gène gfp soit une erreur de 1,5% entre la charge mutationnelle donnée par l'abaque et celle déterminée par séquençage. 6) Détermination de la charge mutationnelle introduite dans la banque de gènes xvnA mutés obtenue à partir de la charge mutationnelle dans la banque de gènes modèle qfp mutés La charge mutationnelle introduite dans la banque de gènes xynA mutés à l'issue de l'expérience de PCR errorprone est déterminée en corrigeant la charge mutationnelle calculée dans la banque de gènes modèle gfp mutés avec le facteur de correction C.F. définie selon l'équation suivante: C.M. (gène d'intérêt) = C.M. (gène modèle) x C.F. dans laquelle, C.F. représente le facteur de correction, et C.M. représente la charge mutationnelle. a) En tenant compte d'au moins un paramètre intrinsèque Dans cette analyse, le paramètre intrinsèque retenu est la longueur des gènes, ce dernier étant reconnu par l'homme de métier comme ayant une forte influence sur la charge mutationnelle des banques de gènes mutés. Le gène modèle correspondant à la gfp et la technique de mutagenèses correspondant à la PCR errorprone, L'unique sous-facteur de correction correspond pour: paramètre intrinsèque 1 = longueur des gènes Avec C.F.[L] = L(gène d'intérêt)/717 où L est en nombre de paire de base. Le tableau 3 c1-après reprend les principaux paramètres intrinsèques, notamment la longueur du gène modèle gfp et du gène d'intérêt xynA. Tableau 3 Composition en base (%) Gène L(pb) A T C G gfp 717 34 26 19,7 20,3 xynA 3120 32 19 23 26 Ainsi la résolution donne: C.F. = C.F.[L] = 3120/717 = 4,35 A l'issu de la simulation avec la longueur comme unique paramètre intrinsèque, la charge mutationnelle calculée dans la banque de gène xynA est de. C.M.[xynA]siml = 4,35 x 9,5 = 41,3 mutations/gène xynA b) En tenant compte d'au moins deux paramètres intrinsèques et d'un paramètre expérimental Dans cette analyse, plus poussée que la précédente, les paramètres intrinsèques retenus sont la longueur des gènes et la composition en base. Le paramètre expérimental retenu est le nombre de multiplication de la matrice ADN lors de la PCR errorprone. Ces paramètres sont connus de l'homme de métier comme pouvant influencer la charge mutationnelle dans les banques de gènes à l'issu d'une mutagenèse par PCR error-prone. Le gène modèle correspondant à la gfp et la technique de mutagenèse correspondant à la PCR errorprone, différents sous-facteurs de correction correspondent pour: - paramètre intrinsèque 1 représente la longueur des gènes selon l'équation suivante: C.F.[L] = L(gène d'intérêt)/717 où L est en nombre de paire de base. - paramètre intrinsèque 2 représente la composition en base sur les brins codant d'ADN selon l'équation suivante: (0,506x Ag.mter,, + 0,346x + 0,063 x Cg.i.è, + 0,085 x Gg.,tre l 0,34 0,26 0,197 0,203, ou Ag.interêt, Tg.interêt, Cg-interêt et Gg. interêt sont en fractions. - paramètre expérimental représente le taux de multiplication selon l'équation suivante: C.F.[d] = d(gène d'interêt) / d(gfp) Le tableau 3 reprend les principaux paramètres intrinsèques, longueur et composition en base, du gène modèle gfp et du gène d'intérêt xynA. Par ailleurs, à l'issu des PCR error-prone réalisées sur le gène modèle gfp et le gène d'intérêt xynA, les taux de multiplication sont de 15 et 12 respectivement pour gfp et xynA. La résolution donne: C.F.[L] = 3120/717 = 4,35 0,506 x 0,32 + 0,346x 0,19 + 0,063 x 0,23 + 0,085 x 0,26 0,34 0,26 0, 197 0,203 C.F.[d] = 12/15 = 0,80 La construction du facteur de correction 30 C.F. s'effectuant par la multiplication des différents sous-facteurs pris en compte pour la détermination de Avec C.F. [base] = C.F. [base] _ = 0,91 la charge mutationnelle de la banque de gènes d'intérêt mutés: C.F. = C.F.[L]X C.F.[base]X C.F.[d] Soit. C.F. = 4,35x0,91x0,8=3,17 A l'issu de la simulation avec la longueur et la composition en base comme paramètres intrinsèques et un le taux de multiplication comme paramètre expérimental, la charge mutationnelle calculée dans la banque de gènes xynA mutés est de: C.M.[xynA]sim2 = 3, 17x 9,5 = 30,1 mutations/gène xynA La charge mutationnelle dans la banque de gènes xynA mutés a également été déterminée par séquençage pour valider les méthodes de détermination de charge mutationnelle découlant de l'invention. Douze clones recombinant xynA Thermotoga maritima mutés ont été sélectionnés pour séquençage. Les gènes xylanases mutés ont été entièrement séquencés comme décrit dans l'exemple 1, (3). Parmi les 38160pb séquencées, 381 mutations ont été détectées dans les gènes mutés xynA. Ainsi la charge mutationnelle observée dans la banque de gènes xynA mutés est de 31,8 mutations / gène xynA. L'erreur constatée entre la charge mutationnelle observée et la charge mutationnelle simulée dans la banque de gène xynA est de 29,9% pour la simulation 1 (réalisée avec la seule longueur comme paramètre intrinsèque) , et de 5,3% pour la simulation 2 (réalisée avec la longueur et la composition en base comme paramètres intrinsèques et le taux de multiplication comme paramètre expérimental). La prise en compte, dans les simulations, de la composition en base et du taux de multiplication du gène modèle gfp et du gène d'intérêt xynA permet de diminuer l'erreur sur la charge mutationnelle d'environ un facteur six. Un tel résultat souligne l'intérêt, présenté par l'invention, d'intégrer dans le facteur de correction C.F. non pas la seule longueur des gènes mais le plus grand nombre de paramètres intrinsèques et expérimentaux identifiés comme susceptibles d'influencer la charge mutationnelle. Par ailleurs la justesse du résultat de la charge mutationnelle obtenu dans le cadre de la simulation 2 par rapport à celle observée lors d'un séquencage (moins de 6% d'erreur) démontre la pertinence de l'invention à pouvoir remplacer de manière fiable l'étude de la charge mutationnelle dans une banque de gènes mutés par séquençage. L'ensemble des résultats sur les charges mutationnelles simulées ou observée dans la banque de gènes xynA mutés est présenté dans le tableau 4 ci-après. Ce tableau 4 récapitule les résultats de charge mutationnelle sur une banque de gènes xynA mutés, observée par séquençage (C.M.[xynA]obs), simulée en tenant compte uniquement de la longueur des gènes (simulation A) ou simulée en tenant compte de la longueur des gènes, de la composition en base et du nombre de multiplication de la matrice ADN (simulation B). Tableau 4 Simulation A B C.M.[xynA]sim Erreur sim./obs. 41,3 29,9% 30,1 5,3% C.M.[xynA]obs 31,8 Exemple 2. Analyse de la plasticité protéique de TEM-1 (3-lactamase Cette analyse est effectuée selon les étapes suivantes: 1) sélection du gène modèle; 2) sélection de la méthode de mutagènese aléatoire; 3) obtention d'un abaque permettant de lier la charge mutationnelle d'une banque de gènes modèle gfp mutés aléatoirement à la fraction de gènes modèle mutés observés dans ladite banque; Ces étapes sont identiques à celles développées dans l'exemple 1. Par contre cet exemple comprend les étapes 4 à 7 décrites ci-après. 4) Application en parallèle de la mutagenèse aléatoire par PCR errorprone au gène modèle qfp et au gène d'intérêt TEM-1 Le gène modèle gfp et le gène d'intérêt TEM-1 (numéro d'accession EMBL AAR25033) ont été soumis à quatres expériences de mutagenèse aléatoire par PCR error-prone avec 4mM MgCl2 (condition I), 7mM MgCl, plus 0,15mM MnCl2 (condition II), 7mM MgC12 plus 0,3mM MnCl2 (condition III) et 7mM MgC12 plus 0,5mM MnCl2(condition IV). Le protocole expérimental est identique à celui décrit dans l'exemple 1. Les rendements des expériences de PCR error-prone ont été déterminés comme décrit dans l'exemple 1. Le nombre de multiplication a ainsi pu être déterminé à 15 et 14 pour les PCR error-prone gfp et TEM-1 respectivement. 5) Détermination des charges mutationnelles dans les banques obtenues à partir du gène modèle qfp A l'issue des quatre expériences de PCR errorprone, les fractions des banques gfp mutés codant cfu total x Ri Fa[GFP] = cfu fluo pour une protéine GFP fonctionnelle ont été déterminées avec l'équation suivante (équation 3) . La détermination des charges mutationnelles dans les banques de gènes gfp mutés est effectuée en utilisant l'abaque décrit dans l'équation suivante (équation 4) . ln(Fa[GFP]) = -0,0083xCM[gfp]2 - 0,1111xCM[gfp] à partir de la fraction de gène gfp codant pour une protéine GFP active (Fa[GFP]). Le tableau 5 ci-après reprend les résultats obtenus à l'issus de cette analyse. 1l donne les charges mutationnelles dans les banques de gènes gfp mutés à partir de l'analyse du test phénotypique mis en place sur la protéine GFP. Ri est le rendement d'insertion des gènes gfp mutés dans les différentes banques déterminé par PCR/colonie, cfu_totales et cfu_fluo sont respectivement le nombre de clones analysés et le nombre de clones fluorescents vert aux U.V., Fa[GFP] est la fraction fonctionnelle de clones GFP tandis que C.M.[gfp] est la charge mutationnelle dans les différentes banques de gènes gfp mutés. Tableau 5 Conditions PCR-ep cfu_totales cfufluo Ri Fa[GFP] C.M.[gfp] I 270 226 0, 958 0,874 1,1 II 340 179 0,958 0,550 4,1 III 775 135 0,958 0,182 9,1 IV 1160 52 0,958 0,047 13,6 6) Détermination des charges mutationnelles introduites dans les banques de gènes TEM-1 mutés obtenues à partir des charges mutationnelles dans les banques de gènes modèle qfp mutés Comme précédemment, les charges mutationnelles introduites dans les banques de gènes TEM 1 mutés à l'issue des expériences de PCR errorprone sont déterminées en corrigeant les charges mutationnelles calculées dans les banques de gènes modèle gfp mutés avec le facteur de correction C.F. selon l'équation suivante (équation 1) C.M. (gène d'intérêt) = C.M. (gène modèle) x C.F. où C.F. représente le facteur de correction et C.M. représente la charge mutationnelle. Pour obtenir un maximum de précision sur la détermination des charges mutationnelles dans les banques de gènes TEM-1 mutés, le facteur de correction est calculé en tenant compte de deux paramètres intrinsèques (la longueur et la composition en base des gènes) et d'un paramètre expérimental (le taux de multiplication des matrices ADN lors de la PCR errorprone). Le gène modèle correspondant à la gfp et la technique de mutagenèse correspondant à la PCR error- prone, différents sous-facteurs de correction correspondent pour: - paramètre intrinsèque 1 représente la longueur des gènes selon l'équation suivante: C.F.[L] = L(gène d'intérêt)/717 où L est en nombre de paire de base. - paramètre intrinsèque 2 représente la composition en base sur les brins codant d'ADN selon l'équation suivante: 0,506 x Ag.mter& 0,346 x 0,063 x Cg...,er t 0,085 x Gg.ni ti Avec C.F. [base] _ + + + 0,34 0,26 0,197 0,203 où Ag interêt, Tg interêt, Cg.interêt et Gg interêt sont en fraction. - paramètre expérimental représente le taux de multiplication selon l'équation suivante: C.F.[d] = d(gène d'interêt) / d(gfp) Le tableau 6 reprend les principaux paramètres intrinsèques, longueur et composition en base du gène modèle gfp et du gène d'intérêt TEM-1. Par ailleurs, à l'issu des PCR error-prone réalisées sur le gène modèle gfp et le gène d'intérêt xynA, les taux de multiplication sont de 15 et 14 pour gfp et TEM-1 respectivement. Composition en base (%) Gène L (pb) A T C G gfp 717 34 26 19,7 20,3 TEM- 1 861 25,9 24,7 23,6 25,8 Tableau 6. Paramètres intrinsèques du gène modèle gfp et du gène d'intérêt TEM-1. La résolution donne: C.F.[L] = 861/717 = 1,2 (0,506x 0,259 + 0,346 x0,247 + 0,063 x 0,236 + 0,085 x 0,258 0,34 0,26 0, 197 0,203 C.F.[d] = 14/15 = 0,93 Soit. C.F. = 1,2x0,9x0,93=1,01 Dans la détermination du facteur de correction C. F., la plus grande longueur du gène TEM-1 est compensée par un plus fort pourcentage en base C et G et une diminution du rendement de PCR errorprone dans le cas du gène TEM-1 par rapport au gène modèle gfp. A l'issu de la simulation avec deux paramètres intrinsèques et un paramètre expérimental, les charges mutationnelles calculées dans les banques de gène TEM-1 sont de: - Pour la condition I: C.M.[TEM-1]sim = 1,01 x 1,1 = 1,1 mutations/gène TEM-1 - Pour la condition II: C.M.[TEM-1]sim = 1,01 x 4,1 = 4,1 mutations/gène TEM-1 - Pour la condition III: C.M.[TEM-1]sim = 1,01 x 9,1 = 9,2 mutations/gène TEM-1 C.F.[base] _ = 0,9 - Pour la condition IV: C.M.[TEM-1]sim = 1,01 X 13,6 = 13,7 mutations/gène TEM-1 6) Analyse de la plasticité protéicue de TEM-1 (3-lactamase L'étude de la plasticité protéique nécessite de pouvoir lier le taux de survie des produits du gène d'intérêt muté, sous un crible donné, dans une banque de gènes d'intérêt muté en fonction de la charge mutationnelle dans ladite banque (Shafikhani S. et al. {1997) BioTechniques 23:304-310 / Daugherty P. et al. (2000)Proc. Natl. Acad. Sci. USA 97:2029-2034). Pour analyser le taux de survie de TEM-1 13-lactamase dans les quatres banques de PCR error-prone précédemment générées, les produits de PCR (lpg) ont été digérés par Eco RI et Nde I (New England Biolabs) pendant 1h à 37 C dans le tampon de réaction adéquat comme décrit par le fournisseur puis purifiés sur colonne QlAquick et insérés dans le vecteur pET26b+ (Stratagene) préalablement digéré avec Eco RI et Nde I. Les cellules E. coli MC1061(DE3) ont été transformées par choc thermique comme décrit par Maniatis (Sambrook J. et al. (1989) Molecular cloning: a laboratory manual Cold Spring Harbor Laboratory Press -New York) avec 5uL des réactions de ligation. Le nombre de transformants dans les différentes banques a été estimé en étalant une fraction de la banque sur boites LB agar complémentées avec 0,6mg/mL de kanamycine tandis que le nombre de clones exprimant une TEM-1 (3-lactamase fonctionnelle a été estimé en étalant le reste de la banque sur boites LB agar complémentées avec 0,6mg/mL de kanamycine et 10ug/mL d'ampicilline. Le rendement d'insertion des gènes mutés TEK-1 a été estimé par PCR sur colonie comme précédemment décrit. Le tableau 7 reprend les résultats obtenus sur les taux de survie de TEM1 0-lactamase mutées dans les 4 banques d'error-prone. La fraction fonctionnelle de clones TEM-1 (Fa[TEM]), est déterminée à partir du rendement d'insertion des gènes TEM-1 mutés (Ri), du nombre de clones analysés (cfu_totales) et du nombre de clones résistants sur ampicilline (cfu_R). Conditions cfu_totales cfu_R Ri Fa ITEM] I 761 539 0,958 0,739 II 746 265 0,958 0,371 III 1215 55 0,938 0,048 IV 1143 7 0,938 0,0065 Tableau 7: Taux de survie de la protéine TEM-1 13-lactamase dans les 4 banques de PCR errorprone. La figure 5 présente l'analyse de la plasticité protéique de TEM-1 Betalactamase, à savoir le taux de survie des protéines TEM-1 Beta-lactamase mutés en fonction de la charge mutationnelle dans les banques de gènes TEM-1 mutés. En conclusion, la présente invention permet d'étudier, dans une forme particulière de mise en oeuvre, la plasticité protéique sans recours au séquençage. Exemple 3: Analyse de l'évolution de TEM-1 Beta- lactamase vers la résistance à l'antibiotique cefotaxime en fonction de la charge mutationnelle des banques de PCR error-prone. Cette analyse comprend les étapes suivantes: 1) sélection du gène modèle. 2) sélection de la méthode de mutagènese aléatoire. 3) obtention d'un abaque permettant de lier la charge mutationnelle d'une banque de gènes modèle gfp mutés aléatoirement à la fraction de gènes modèle mutés observés dans ladite banque. Ces étapes sont identiques à celles développées dans l'exemple 1 et 2. Par contre, cette analyse comprend aussi les étapes 4-7 décrites ciaprès. 4) Application en parallèle de la mutagenèse aléatoire par PCR errorprone au crène modèle gfp et au gène d'intérêt TEM-1 Le gène modèle gfp et le gène d'intérêt TEM-1 ont été soumis à cinq expériences de mutagenèse aléatoire par PCR error-prone avec 5mM MgCl2 (condition 1), 7mM MgCl, plus 0,1mM MnCl2 (condition 2), 7mM MgCl, plus 0,2mM MnCl2 (condition 3), 7mM MgCl2 plus 0,3mM MnCl2 (condition 4) et 7mM MgCl2 plus 0,5mM MnCl2 (condition 5). Le protocole expérimental est identique à celui décrit dans l'exemple 2. Les rendements des expériences de PCR error-prone ont été déterminés comme décrit dans l'exemple 1 et sont identiques à ceux trouvés dans l'exemple 2. Le taux de multiplication est de 15 et de 14 pour les PCR error-prone gfp et TEM-1 respectivement. 5) Détermination des charqes mutationnelles dans les banques obtenues à partir du crène modèle gfp A l'issue des cinq expériences de PCR errorprone, les fractions des banques gfp mutés codant pour une protéine fonctionnelle ont été déterminées avec l'équation suivante (équation 3): Fa[GFP] = du flux cfu total x Ri La détermination des charges mutationnelles dans les banques de gènes gfp mutés est effectuée en utilisant l'abaque décrit dans l'équation suivante: ln(Fa[GFP]) = -0, 0083xCM[gfp]2 - 0,1111xCM[gfp] à partir de la fraction de gènes gfp mutés codant pour une protéine GFP active (Fa[GFP]). Le tableau 8 ci-après reprend les résultats obtenus à l'issue de cette analyse. Ce tableau rapporte la détermination des charges mutationnelles dans les banques de gènes gfp mutés à partir de l'analyse du test phénotypique GFP. Ri est le rendement d'insertion des gènes gfp mutés dans les différentes banques, cfu_totales et cfu_fluo sont respectivement le nombre de clones analysés et le nombre de clones fluorescents aux U.V., Fa[GFP] est la fraction fonctionnelle de clones GFP tandis que C.M.[gfp] est la charge mutationnelle dans les différentes banques de gènes gfp mutés. Tableau 8 Conditions PCR-ep cfu_totales cfu_fluo Ri Fa[GFP] C.M.[gfp] 1 425 305 0,96 0,749 2,2 2 450 268 0,96 0,621 3,4 3 568 212 0,96 0,389 5,9 4 842 100 0,96 0,124 10,5 1093 26 0,75 0,032 14,8 20 6) Détermination des charges mutationnelles introduites dans les banques de gènes TEK-1 mutés obtenues à partir des charges mutationnelles dans les banques de gènes modèle qfp mutés Le facteur de correction est calculé en tenant compte de deux paramètres intrinsèques, la longueur et la composition en base des gènes, et d'un paramètre expérimental, le taux de multiplication des matrices ADN lors de la PCR error-prone. Le gène modèle gfp et le gène d'intérêt TEM-1 étant identique à celui de l'exemple 2 et les taux de multiplication des matrices gfp et TEM-1 étant également inchangés entre les exemple 1 et 2, le facteur de correction C. F. est identique à celui déterminé dans l'exemple 2. Soit C.F. = 1,01 A l'issu de la simulation avec deux paramètres intrinsèques et un paramètres expérimental, les charges mutationnelles calculées dans les cinq banques de gènes TEM-1 mutés sont de: - Pour la condition 1 C.M.[TEM-1]sim = 1,01 x 2,2 = 2,2 mutations/gène TEM-1 - Pour la condition 2 C.M.[TEM-1]sim = 1,01 x 3,4 = 3,4 mutations/gène TEM-1 Pour la condition 3 C.M.[TEM-1]sim = 1,01 x 5,9 = 6,0 mutations/gène TEM-1 - Pour la condition 4 C.M.[TEM-1]sim = 1,01 X 10,5 = 10,6 mutations/gène TEM-1 - Pour la condition 5 C.M.[TEM-1]sim = 1,01 x 14,8 = 14,9 mutations/gène TEM-1 7) Analyse de l'acquisition de la résistance à la cefotaxime des clones recombinant TEM-1 Ji-lactamase mutés Le système d'évolution choisi décrit l'influence de la charge mutationnelle de la banque de gène TEM-1 mutés sur la fréquence d'apparition de clones résistant à la cefotaxime grâce à l'évolution du gène TEM-1. Pour analyser la fréquence d'apparition de clones résistants, les produits PCR sont clonés dans le vecteur pET26b+ comme décrit précédemment pour former 5 banques de gènes TEM-1 mutés. Les banques de clones sont réalisées par électrotransformation de cellules E. coli MC1061(DE3) avec 10x2pL de ligation. Le nombre de transformants dans chacune des cinq banques a été estimé en étalant 1/25eme de la banque sur boites LB agar complémentées avec 0,6mg/mL de kanamycine tandis que les clones résistant à la cefotaxime ont été isolés en étalant le reste de la banque (24/25eme) sur boites LB agar complémentées avec 0,6mg/mL de kanamycine et 15ng/mL de cefotaxime. Le rendement d'insertion des gènes TEM-1 mutés a été estimé par PCR sur colonie. Le tableau 9 ci-après reprend les résultats obtenus sur la fréquence d'apparition de clones, résistant à la cefotaxime, dans les 5 banques de TEM-mutés. La fréquence de clones améliorés dans les différentes banques est calculée à partir du nombre de clones résistant sur cefotaxime (cfu_cefoR) et du nombre de clones recombinant (Nb_recornb.), lui même déduit du nombre de clones dans la banque (cfu_totales) et du rendement d'insertion des gènes TEK-1 (Ri) mutés. Tableau 9 Conditions cfu_totales Ri Nb recomb. cfu cefoR Fréq. améliorés 0/00 1 103200 0,96 98900 51 0,516 2 53750 0,96 51510 41 0,796 3 189000 0,96 181125 50 0,276 4 158000 0,94 148125 18 0,122 223500 0,92 204875 11 0,054 La figure 6 montre l'influence de la charge mutationnelle des banques de gènes TEM-1 mutés sur la fréquence d'apparition d'un variant TEM-1 améliorés. La fréquence d'apparition de clones résistant à la cefotaxime en fonction de la charge mutationnelle suit une courbe en cloche passant par une charge mutationnelle optimale. La méthode de l'invention peut ainsi être utilisée pour analyser les taux d'évolution des banques de gènes mutés mais également pour se placer à la charge mutationnelle optimale et ceci sans recourir au séquençage. En conclusion, cet exemple montre également l'importance de contrôler précisément la charge mutationnelle dans les banques de gènes mutés puisque la fréquence d'apparition de variants TEM-1 résistant à la céfotaxime est diminuée de trois fois pour une charge mutationnelle passant de 3,4 à 6 mutations par gène TEM-1. La présente invention permet, sans séquençage, un contrôle de la charge mutationnelle beaucoup plus fin que les méthodes proposés par l'art antérieur et répond parfaitement à ce type de besoin. Exemple 4: comparaison et estimation de la fidélité d'un jeux de ADN polymerases dans une réaction de PCR donnée. Cette mise en oeuvre particulière l'invention comprend les étapes suivantes: 1) sélection du gène modèle. 2) sélection de la méthode de mutagènese aléatoire. 3) obtention d'un abaque permettant de lier la charge mutationnelle d'une banque de gènes modèle gfp mutés aléatoirement à la fraction de gènes modèle mutés observés dans ladite banque. Ces étapes sont identiques à celles développées dans les exemple 1 à 3. L'analyse comprend en outre les étapes 4-6 décrites ci-après. 4) Application en parallèle de la PCR au gène modèle qfp avec un -jeux de ADN polymerases Le gène modèle gfp a été soumis à trois expériences de PCR contenant 7mM de MgCl, et 0,2mM de MnC12. Le protocole expérimental est identique à celui décrit dans les exemples précédents avec un temps d'élongation de 1 minute, excepté que 2,5U de Taq polymerase (Q-Biogen), de VentR (New England Biolabs) ou de Pfu Turbo (Stratagène) ont été utilisée comme ADN polymérase. Les rendements des expériences de PCR ont été déterminés comme décrit précédemment. Le nombre de multiplication est de dA, dB et d, pour les PCR effectuées avec la Taq DNA polymerase, la Vent}, DNA polymerase et la Pfu DNA polymerase respectivement. 5) Détermination des charges mutationnelles dans les banques obtenues à partir du crène modèle qfp A l'issue des expériences de PCR, les fractions des banques gfp mutés codant pour une GFP fonctionnelle ont été déterminées avec l'équation suivante (équation 3) . Fa[GFP] = cfu_fluo cfu_total x Ri Les charges mutationnelles dans les banques de gènes gfp mutés sont déterminées en utilisant l'abaque décrit dans l'équation suivante (équation 4) : ln(Fa[GFP]) = -0,0083xCM[gfp]2 - 0,1111xCM[gfp] à partir de la fraction de gènes gfp mutés codant pour une protéine GFP active (Fa[GFP]). Les charges mutationnelles sont de A, B et C pour les PCR effectuées avec la Taq DNA polymerase, la VentR DNA polymerase et la Pfu DNA polymerase respectivement. 6) Comparaison de la fidélité des trois polymerases utilisées dans les expériences de PCR Les taux d'erreur des différentes polymerases dans les conditions réactionnelles 2871170 49 utilisées sont déduits des charges mutationnelles dans les banques de gfp mutés en normalisant par la longueur du gène modèle gfp et les taux de multiplication des PCR. La longueur du gène modèle gfp est de 717 paires de base et les taux de multiplication à l'issu des PCR sont de dA, dB et d, avec la Taq DNA polymerase, la VentR DNA polymerase et la Pfu DNA polymerase respectivement. Ainsi Taux erreur [Taq DNA polymerase] _ A/(717xdA) erreur /pb/ duplication Taux erreur [VentR DNA polymerase] _ B/(717xdB) erreur /pb/ duplication Taux erreur [Pfu DNA polymerase] _ C/(717xd,) erreur /pb/ duplication En conclusion, la méthode de l'invention permet de déterminer précisément les fidélités d'un jeux de ADN polymerase dans un mélange réactionnel donné sans recourir au séquençage, et ceci en s'affranchissant des variations de conditions expérimentales. 15 REVENDICATIONS 1) Procédé de détermination de la charge mutationnelle d'une banque de gènes obtenue par mutagenèse aléatoire d'un gène d'intérêt, caractérisé en ce qu'il comprend les étapes suivantes: a) la préparation d'un abaque liant la charge mutationnelle (CM) d'une banque de gènes obtenue par mutagenèse aléatoire d'un gène modèle à la fraction des gènes modèles mutés observée dans ladite banque; b) la mutagenèse aléatoire du gène modèle ayant servi pour la préparation de l'abaque de l'étape (a) et du gène d'intérêt pour obtenir les banques de gènes mutés correspondantes; c) la détermination de la charge mutationnelle (CM) de la banque de gènes obtenue à partir du gène modèle à l'étape (b) à partir de l'abaque établi à l'étape (a) ; d) l'application d'un facteur de correction (CF) à la charge mutationnelle (CM) de la banque de gènes modèle mutés déterminée à l'étape (c) pour déterminer la charge mutationnelle (CM) de la banque de gènes d'intérêt mutés de l'étape (b). 2) Procédé selon la revendication 1, caractérisé en ce que la technique de mutagenèse mise en oeuvre à l'étape (b) est identique pour le gène modèle et pour le gène d'intérêt et est également similaire à celle de l'étape (a) pour construire l'abaque à partir du gène modèle. 3) Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que les étapes (c) et (d) ne comprennent pas de séquençage. 2871170 51 4) Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que pour la préparation de l'abaque de l'étape (a) l'observation des gènes modèle mutés dans ladite banque est réalisée en analysant directement ou indirectement une ou plusieurs modifications d'une ou plusieurs propriétés de la séquence nucléotidique des gènes modèles mutés par rapport au gène modèle avant mutagenèse. 5) Procédé selon la revendication 4, caractérisé en ce que ladite observation consiste à observer la modification d'une ou plusieurs propriétés de la protéine codée par le gène modèle. 6) Procédé selon la revendication 5, caractérisé en ce que l'observation de la modification d'une ou plusieurs propriétés de la protéine codée par le gène modèle est réalisée par un test phénotypique. 7) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le gène modèle possède des caractéristiques partagées par un grand nombre de gènes. 8) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le gène modèle code la protéine GFP. 9) Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que le gène modèle possède des caractéristiques proches de celles du gène d'intérêt. 10) Procédé selon l'une quelconque des revendications 7 ou 9, caractérisé en ce que les caractéristiques du gène modèle comprennent - sa longueur, sa composition en base, - sa fréquence en bases répétées dans la séquence, - la présence de séquences symétriques dans le gène. 11) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la mutagenèse aléatoire des étapes (a) et (b) met en oeuvre la polymérisation. 12) Procédé selon la revendication 11, caractérisé en ce que la polymérisation est la PCR error-prone. 13) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'à l'étape (b) le même mélange réactionnel est utilisé pour la mutagenèse du gène modèle et du gène d'intérêt. 14) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'à l'étape (d), la charge mutationnelle déterminée à l'étape (c) est corrigée en appliquant un facteur de correction (C.F.) qui se décompose en plusieurs sous- facteurs définis suivant les paramètres sélectionnés comme influençant la charge mutationnelle d'une banque de gènes d'intérêt mutés par rapport à une banque de gènes modèle mutés. 15) Procédé selon la revendication 14, caractérisé en' ce que les paramètres sélectionnés pouvant influencer la charge mutationnelle sont choisis parmi - au moins un et de préférence au moins deux paramètres intrinsèques, - au moins un paramètre expérimental. 16) Procédé selon la revendication 15, caractérisé en ce que les paramètres intrinsèques sont choisis parmi les caractéristiques suivantes d'un gène - sa longueur, - sa composition en base - sa fréquence en bases répétées dans la séquence, - la présence de séquences symétriques dans le gène. 17) Procédé selon la revendication 15, caractérisé en ce que, lorsque la mutagenèse est basée sur une amplification en chaîne utilisant une polymérase, les paramètres expérimentaux correspondent à un ou plusieurs des facteurs suivants: - le taux d'amplification du gène durant l'étape de mutagenèse, - la quantité de matrice initiale, - le nombre de cycle d'amplification, - à la nature et à la quantité de polymérase utilisée aux conditions expérimentales utilisées comme la concentration en cations divalents, la concentration et le ratio de concentration des dNTPs. 2871170 54 18) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que au moins un paramètre intrinsèque et au moins un paramètre expérimental sont considérés pour établir le facteur de correction. 19) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que lorsque la mutagenèse est réalisée par PCR errorprone, et le gène modèle est la gfp, les différents sous-facteurs de correction sont choisis parmi: - le paramètre intrinsèque 1 représentant la longueur des gènes selon l'équation suivante: C.F.[L] = L(gène d'intérêt)/717 où L est en nombre de paire de base, - le paramètre intrinsèque 2 représentant C.F.[base] = ( 0,506 x A(gintetët) + 0,346 x T(g.intera + 0,063 x C(gsntee) + 0,085 x G(g.intert) 0,34 0,26 0,197 0,203 > la composition en base selon l'équation suivante: Ou A(g.interêt)/ T (g.interêt), C(g.interêt) et G(g.interêt) sont en fraction, - le paramètre expérimental représente le taux de multiplication selon l'équation suivante: C.F.[d] = d(gène d'interêt) / d(gfp) où d est le taux de multiplication calculé selon l'équation suivante: d = ln(Qf/Qi) /1n2 dans laquelle, Qi et Qf représentent respectivement les quantités molaires de la matrice et des produits à l'issue de la PCR. 2871170 55 20) Utilisation d'un procédé selon l'une quelconque des revendications 1 à 19 pour contrôler la charge mutationnelle dans une banque de gènes mutés caractérisé en ce qu'il comprend une étape supplémentaire de sélection des conditions expérimentales de mutagenèse aléatoire qui seront utilisées pour mettre en oeuvre l'étape (b). 21) Utilisation d'un procédé selon l'une quelconque des revendications 1 à 19 pour étudier la capacité d'une protéine à incorporer des mutations tout en conservant une activité présentée par la protéine sauvage. 22) Utilisation d'un procédé selon l'une quelconque des revendications 1 à 19 pour étudier la capacité d'une protéine à acquérir une nouvelle propriété avec une charge mutationnelle donnée. 23) Utilisation d'un procédé selon l'une quelconque des revendications 1 à 19 pour déterminer la charge mutationnelle optimale dans des banques de gènes d'intérêt mutés pour générer un lot de gènes parentaux avant la mise en oeuvre de méthodes de recombinaison moléculaire. 24) Utilisation d'un procédé selon l'une quelconque des revendications 1 à 19 pour quantifier ou comparer la fidélité d'une (de) méthode(s) de polymérisation d'acides nucléiques. 25) Kit susceptible d'être mis en oeuvre dans un procédé pour déterminer la charge mutationnelle 2871170 56 d'une banque obtenue à l'issu de la mutagenèse aléatoire d'un gène d'intérêt selon l'une quelconque des revendications 1 à 19 ou dans une utilisation selon les revendications 20 à 24, caractérisé en ce qu'il comprend: - au moins un gène modèle par exemple sous la forme d'un plasmide; - au moins un abaque permettant de lier la charge mutationnelle d'une banque de gènes obtenue par mutagenèse aléatoire du gène modèle à la fraction de gènes modèles mutés observés dans ladite banque; - les réactifs nécessaires à la mise en oeuvre de la technique de mutagenèse aléatoire comme par exemple une ADN polymérase thermostable, dNTP, magnésium, manganèse, tampon de PCR dans le cas où la technique de mutagenèse correspond à la PCR errorprone; - les réactifs, si nécessaires, pour observer des gènes modèle mutant; - une notice d'utilisation indiquant les sous-facteurs de correction définis pour un paramètre sélectionné comme influençant la charge mutationnelle. 26) Kit susceptible d'être mis en oeuvre dans une utilisation selon la revendication 20 caractérisé en ce qu'il comprend en plus des composant du kit selon la revendication 25 au moins un abaque permettant de lier la charge mutationnelle (CM) d'une banque de gènes obtenue parmutagénèse aléatoire d'un gène modèle aux conditions expérimentales de mutagénèse aléatoire. 27) Kit selon la revendication 25 ou 26, caractérisé en ce qu'il possède les amorces pour 2871170 57 l'amplification du gène modèle dans le cas où la technique de mutagenèse est basée sur une technique d'amplification. 28) Kit selon l'une quelconque des revendications 25 à 27, caractérisé en ce que les sous-facteurs de correction fournis doivent être complétés par une ou plusieurs des données comme celles choisies notamment dans le groupe suivant: - la longueur du gène d'intérêt, - la composition en base A du gène d'intérêt, d'intérêt, la composition en base C du gène - la composition en base G du gène d'intérêt, - la composition en base T du gène d'intérêt, - le nombre de multiplication de la matrice 20 du gène modèle - le nombre de multiplication de la matrice du gène d'intérêt, dont les valeurs correspondantes sont intégrées par l'expérimentateur, en fonction du gène 25 d'intérêt et du résultat de mutagenèse aléatoire. 29) Kit selon l'une quelconque des revendications 25 à 28, caractérisé en ce que le gène modèle est celui de la gfp. 30) Kit selon la revendication précédente caractérisé, en ce que la mutagenèse est réalisée par PCR error-prone, le gène modèle est celui de la gfp, 15 2871170 58 les différents sous-facteurs de correction sont choisis parmi: - le paramètre intrinsèque 1 représentant la longueur des gènes selon l'équation suivante: C.F.[L] = L(gène d'intérêt)/717 où L est en nombre de paire de base, - le paramètre intrinsèque 2 représentant la composition en base sur les brins codant d'ADN selon C.F.[base] 0,506 x A(g.interet) + 0,346 x T(g. interêt) + 0,063 x C(g.interêt) + 0,085 x G(g.interêt) 0,34 0,26 0,197 0,203 l'équation suivante: où A(g.interêt)i T(g.interêt)/ C (g.interêt) et G(g.interêt) sont en fraction, - le paramètre expérimental représente le taux de multiplication selon l'équation suivante: C.F.[d] = d(gène d'interêt) / d(gfp) où d est le taux de multiplication calculé selon l'équation suivante: d = ln(Qf/Qi) /ln2 dans laquelle, Qi et Qf représentent respectivement les quantités molaires de la matrice 20 initiale et des produits à l'issue de la PCR. 31) Kit selon l'une quelconque des revendications 25 à 30, caractérisé en ce qu'il possède plusieurs gènes modèles et plusieurs abaques pour chaque gène modèle de manière 'à choisir le gène modèle ayant des caractéristiques "intrinsèques les plus proches du gène d'intérêt. 32) Procédé selon l'une quelconque des 30 revendications 1 à 19, caractérisé en ce que l'une au 2871170 59 moins des étapes (a), (c) et (d) est mise en oeuvre par des moyens informatiques. 33) Procédé selon la revendication 32, caractérisé en ce que lesdits moyens informatiques exploitent des données numériques relatives à l'une au moins des caractéristiques suivantes: - le gène modèle sélectionné - la technique de mutagenèse choisie - le ou les paramètres sélectionnés comme capables d'influencer la charge mutationnelle d'une banque de gènes mutés par rapport à une autre pour définir les sous-facteurs de correction. 34) Procédé selon l'une des revendications 32 ou 33, caractérisé en ce que lesdits moyens informatiques exploitent dés données numériques relatives à l'une au moins des caractéristiques suivantes: - la longueur du gène d'intérêt, - la composition en base A du gène d'intérêt, - la composition en base C du gène d'intérêt, - la composition en base G du gène d'intérêt, la composition en base T du gène d'intérêt, - le nombre de multiplication de la matrice du gène modèle - le nombre de multiplication de la matrice du gène d'intérêt, dont les valeurs correspondantes sont intégrées par l'expérimentateur, en fonction du gène d'intérêt et du résultat de mutagenèsemutagenèse aléatoire.