# Champ de force (chimie)

Dans le cadre de la mécanique moléculaire, un champ de force est un ensemble de potentiels et de paramètres permettant de décrire la structure de l'énergie potentielle d'un système de particules (typiquement, des atomes, mais non exclusivement). L'usage de l'expression champ de force en chimie et biologie numériques diffère ainsi de celui de la physique, où il indique en général un gradient négatif d'un potentiel scalaire. Les potentiels et paramètres d'un champ de force donné sont déterminés à partir de résultats expérimentaux et de calculs précis en mécanique quantique.
Il existe, en fonction des systèmes étudiés, de nombreux champs de force différents. Les champs de force tout-atome (all atom) donnent des paramètres pour tous les atomes du système, y compris l'hydrogène, alors que les champs de force atomes unifiés (unified atom) ne prennent pas en compte les atomes d'hydrogène sauf si ces derniers sont polaires (H est traité comme un centre d'interaction). Les champs de force basés sur les atomes lourds ne prennent pas du tout en compte les atomes d'hydrogène. Certains champs de force sont basés sur des modèles simplifiés de tout ou partie d'un ou plusieurs objets comme les champs à gros grains (coarse grained) fréquemment utilisés dans les simulations longues de protéines, ou les champs basés sur des grilles et réseaux (in lattice).
Le choix d'un champ, spécifique d'un système et plus ou moins efficace, dépend des études que l'on doit produire et des propriétés recherchées.

## Forme générale

La forme fonctionnelle de base d'un champ de force comprend les termes de liaisons reliés aux atomes liés par des liaisons covalentes, et les termes d'interactions (termes dit non-liés ou non-covalents) qui décrivent les interactions à longue portée (électrostatique et de Van der Waals). La décomposition spécifique en chacun de ces termes dépend du champ de force, mais la forme générale de l'énergie totale dans un champ de force additif peut être écrite :
    $\ E_{totale}=E_{liee}+E_{non-liee}$
où les composants des contributions covalentes et non-covalentes sont données par les sommes suivantes :
    $\ E_{liee}=E_{liaison}+E_{angle}+E_{angle\ diedre}$
    $\ E_{non-liee}=E_{electrostatique}+E_{van\ der\ Waals}$
Les termes de liaison et d'angle sont habituellement modélisés comme des oscillateurs harmoniques dans les champs de force ne permettant pas la formation/destruction de liaison. Une description plus réaliste d'une liaison covalente soumise à un étirement important peut être donnée par un potentiel de Morse, plus couteux. La forme fonctionnelle du reste des termes liés est très variable. Des potentiels dièdres sont également habituellement inclus. De plus, des termes de « torsion impropre » peuvent être ajoutés pour forcer la planéité des cycles aromatiques et autres systèmes conjugués, ainsi que des termes « croisés » décrivant le couplage entre les différentes variables internes, comme les angles et les longueurs de liaison. Certains champs de force incluent des termes explicites pour les liaisons hydrogène.
Les termes non-liés sont plus couteux en temps de calcul car ils doivent prendre en compte bien plus d'interactions par atome. Un des choix les plus courants est de limiter les interactions aux énergies de paire. Le terme de van der Waals est en général calculé à l'aide d'un potentiel de Lennard-Jones, et le terme électrostatique par la loi de Coulomb, bien que ces deux termes puissent être atténués ou augmentés par un facteur constant afin de tenir compte de la polarisabilité électronique, et produire ainsi un meilleur accord avec l'expérience.

## Paramétrage

En plus de la forme fonctionnelle des potentiels, un champ de force définit un ensemble de paramètres pour chacun des types d'atomes. Ainsi par exemple, un champ de force comprendra des paramètres distincts pour un atome d'oxygène dans un groupe carbonyle et dans un groupe hydroxyle. Un ensemble typique de paramètres comprend des valeurs pour la masse atomique, le rayon de van der Waals et la charge partielle d'atomes individuels, et des valeurs d'équilibre de longueurs de liaison, de mesures d'angles plans et dièdres pour des paires, triplets et quadruplets d'atomes liées, et les valeurs de constante de force pour chaque potentiel. Les champs de force les plus courants utilisent un modèle de charge fixe dans lequel chaque atome se voit affecter une charge électrostatique qui n'est pas affectée par l'environnement local. Les développements proposés pour les champs de force de prochaine génération incluent des modèles pour la polarisabilité dans lesquels les charges des particules sont influencées par celles de leurs voisines. La polarisabilité peut, par exemple, être approchée par l'introduction de dipôles induits ; elle peut aussi être représentée par des particules de Drude, ou des sites porteurs de charges virtuels et sans masses rattachés par un potentiel harmonique de type ressort à chaque atome polarisable. L'introduction de la polarisabilité dans les champs de force pour un usage courant est ralentie par le cout important associé au calcul du champ électrique local.
Bien que de nombreuses simulations moléculaires portent sur des macromolécules biologiques comme les protéines, l'ADN ou l'ARN, les paramètres utilisés pour un type d'atomes donné proviennent généralement des observations de molécules organiques, plus facilement appréhendables que cela soit pour des études expérimentales ou pour des calculs quantiques. Des champs de force différents peuvent être issus de types données expérimentales distinctes, comme l'enthalpie de vaporisation (comme dans le champ OPLS), l'enthalpie de sublimation (comme dans le champ CFF), les moments dipolaires, ou des paramètres spectroscopiques variés (CFF).
Les ensembles de paramètres et les formes fonctionnelles sont définies par les développeurs de champs de force afin d'être auto-cohérents. Parce que les formes fonctionnelles des potentiels varient de manière extensive y compris pour des champs de force très proches (ou entre versions successives d'un même champ de force), les paramètres d'un champ de force ne doivent pas être utilisés avec les potentiels définis pour un autre champ.

## Défauts et manques

Tous les champs de force sont fondés sur de nombreuses approximations et construits à partir de différents types de données expérimentales. Ils sont en conséquence appelés empiriques. Certains champs ne tiennent habituellement pas compte de la polarisation électronique de l'environnement, ce qui a pour effet de réduire de manière significative les interactions électrostatiques des charges atomiques partielles. Ce problème fut considéré au moyen du développement des « champs de force polarisables »$^,$ ou par utilisation d'une constante diélectrique macroscopique. Cependant, l'utilisation d'une seule valeur pour la constante diélectrique est sujette à caution dans les environnements hautement hétérogènes des protéines ou des membranes biologiques, et la nature du diélectrique dépend du modèle utilisé.
Les forces de van der Waals sont aussi très fortement dépendantes de leur environnement, car résultant de l'interaction entre dipôles induits et « instantanés » (voir force intermoléculaire). La théorie initiale de Fritz London sur ces forces n'est valable que dans le vide. Une théorie plus générale sur les forces van der Waals dans un milieu dense fut développée par A.D. McLachlan en 1963, pour laquelle la théorie de Fritz-London est considérée comme un cas spécifique. Cette théorie indique que les interactions de van der Waals dans un milieu sont plus faibles que dans le vide et suivent la loi « comme s'il y avait dissolution », ce qui signifie que les atomes de types différents interagissent plus faiblement que les atomes de même type. Cette théorie contraste avec les lois de combinaison ou l'équation de Slater-Kirkwood utilisées pour le développement des champs de force classiques. Les lois de combinaison indiquent que l'énergie d'interaction de deux atomes distincts (comme C et N) est une moyenne des énergies d'interaction des paires d'atomes correspondantes (soit C-C et N-N). Selon la théorie de McLachlan, les interactions entre particules dans un milieu peuvent aller jusqu'à être complètement répulsives, comme observé dans l'hélium liquide. Les conclusions de la théorie de McLachlan sont étayées par des mesures directes des forces d'attraction entre différents matériaux (constante d'Hamaker), comme indiqué par Jacob Israelachvili dans son livre Intermolecular and surface forces. Il y était indiqué que « l'interaction entre hydrocarbures dans l'eau est environ 10 % de celle dans le vide ». De tels effets ne sont pas pris en compte en mécanique moléculaire classique.
Un autre type de critique vient des applications pratiques, comme la détermination de structure protéique. Ainsi, aucun des participants au CASP n'essaya de raffiner ses modèles afin d'éviter « un embarras majeur pour la mécanique moléculaire, c'est-à-dire le fait que la minimisation de l'énergie ou la dynamique moléculaire conduisent généralement à un modèle peu en accord avec la structure expérimentale ». En fait, les champs de force ont été appliqués avec succès dans les déterminations de structures protéiques dans différentes applications en diffractométrie de rayons X et en résonance magnétique nucléaire, par exemple avec le code XPLOR. Cependant, ces déterminations sont fondées sur des contraintes expérimentales, alors que les champs de force sont principalement destinés à supprimer les interférences atomiques. Les résultats des calculs sont pratiquement les mêmes avec des potentiels de sphères durs (par exemple dans le code DYANA, pour des calculs fondés sur les données RMN), ou avec des programmes de raffinement cristallographique qui n'utilisent aucune fonction de l'énergie. Les lacunes des champs de force restent un goulet d'étranglement pour la modélisation par homologie des protéines. Cette situation a donné naissance au développement de fonctions empiriques de notation spécifiques pour l'assemblage de ligands, le repliement de protéine$^,$$^,$, la construction de protéines numérique$^,$$^,$, ou la modélisation des protéines dans les membranes.
Il existe également une opinion selon laquelle la mécanique moléculaire peut opérer avec une énergie non pertinente pour le repliement des protéines ou la liaison de ligand. Les paramètres des champs de force typiques permettent de reproduire l'enthalpie de sublimation. Cependant, le repliement des protéines ou la liaison de ligand sont thermodynamiquement très semblables à une cristallisation ou à des transitions liquide-solide, car ces processus constituent un « gel » de molécules mobiles dans un milieu condensé$^,$$^,$. Ainsi, les variations de l'énergie libre durant ces processus sont supposées représenter une combinaison d'une énergie de type enthalpie de fusion, d'une contribution d'entropie conformationnelle et d'énergie libre de solvatation. L'enthalpie de fusion est significativement plus petite que l'enthalpie de sublimation. Par conséquent, les potentiels décrivant le repliement des protéines ou la liaison de ligand doivent être plus faibles que les potentiels en mécanique moléculaire. En fait, les énergies des liaisons hydrogène dans les protéines sont de ~ -1,5 kcal/mol lorsqu'elles sont estimées à partir de l'ingénierie protéique ou des données sur la transition entre hélice alpha et pelote$^,$, mais si l'on estime les mêmes quantités à partir des enthalpies de sublimation des cristaux moléculaires, elles sont de l'ordre de -4 à -6 kcal/mol. La profondeur des potentiels de Lennard-Jones modifiés dérivés des données de l'ingénierie protéique sont aussi plus petites que dans un champ de force classique et respectent une loi « comme s'il y avait dissolution », comme indiqué par la théorie de McLachlan.

## Champs de force « courants »

Différents champs de force ont été construits afin d'être appliqués à différents problèmes.
Le champ MM2 a ainsi été développé par Norman Allinger dans un premier temps pour l'analyse conformationnelle des hydrocarbures et autres petites molécules organiques. Il a été conçu afin de reproduire les géométries « covalentes » d'équilibre aussi précisément que possible. Il comprend un grand ensemble de paramètres continument redéfinis et mis à jour pour différentes classes de composés organiques (MM3 et MM4)$^,$$^,$$^,$$^,$.
Le champ CFF a été développé par Wrshel, Lifson et collaborateurs comme une méthode générale d'unification des études sur les énergies, les structures et vibrations des molécules et cristaux moléculaires. Le programme CFF, développé par Levitt et Warshell, est basé sur une représentation cartésienne de tous les atomes, et sert de base à de nombreux programmes de simulations plus importants.
Le champ ECEPP a été développé de manière spécifique pour la modélisation des peptides et protéines. Il utilise des géométries fixés de résidus d'acides aminés afin de simplifier la surface d'énergie potentielle. La minimisation d'énergie est ainsi menée dans l'espace des angles de torsion des protéines. Les champs MM2 et ECEPP incluent tous deux des potentiels pour les liaisons hydrogène et des potentiels de torsion afin de décrire les rotations autour des liaisons simples. ECEPP/3 a été intégré dans Internal Coordinate Mechanics et FANTOM, par exemple.
Les champs AMBER, CHARMM et GROMOS ont été développés en premier lieu pour la dynamique moléculaire des macromolécules, bien qu'ils soient communément employés pour la minimisation d'énergie. Les coordonnées de tous les atomes y sont considérées comme des variables libres.

### Champs de force classiques

Voici quelques champs de force « classiques » :
* AMBER (Assisted Model Building and Energy Refinement). Très utilisé pour la modélisation des protéines, ADN et ARN.
* CHARMM (Chemistry at HARvard Molecular Mechanics). Développé initialement à Harvard, très utilisé pour les petites molécules ou les macromolécules. CHARMm est sa version commerciale, proposée par Dassault Systemes BIOVIA.
* CVFF également utilisé massivement pour les petites molécules ou les macromolécules.
* COSMOS-NMR. Champ de force hybride QM/MM (quantique/classique) adapté à de nombreux composés inorganiques, organiques ou macromolécules biologiques. Il inclut le calcul semi-empirique des charges atomiques et propriétés de RMN. Il est d'ailleurs optimisé pour la résolution des structures basée sur la RMN et est introduit dans l'ensemble COSMOS de modélisation moléculaire.
* GROMACS. Champ de force optimisé pour l'ensemble du même nom.
* GROMOS. Champ de force intégré dans l'ensemble généraliste GROMOS (GROningen MOlecular Simulation package), dédié à la simulation dynamique des systèmes biomoléculaires. Le champ de force GROMOS (version A) a été développé pour être appliqué aux solutions aqueuses ou apolaires de protéines, nucléotides ou sucres. Cependant, une version pour simulation en phase gaz (version B) de molécules isolées est disponible.
* OPLS (Optimized Potential for Liquid Simulations). Champ de force développé par William L. Jorgensen au département de chimie de Yale. Il en existe plusieurs évolutions comme OPLS-AA, OPLS-UA, OPLS-2001, ou OPLS-2005.
* ECEPP/2. Premier champ de force dédié aux molécules polypeptidiques, développé par F.A.Momany, H.A.Scheraga et collaborateurs.
* QCFF/PI. Champ de force dédié aux molécules conjugées$^,$.

### Champs de force de deuxième génération

* CFF est une famille de champs de force adaptés à une large gamme de composés organiques, comme pour les polymères, les métaux, etc.
* MMFF (Merck Molecular Force Field) développé pour une large gamme de molécules.
* MM2, MM3, MM4, développés par Norman Allinger, paramétrée pour une large gamme de molécules.
* QVBMM développé par Vernon G.S. Box, paramétré pour traiter les biomolécules et une large gamme de molécules organiques, et inclus dans STR3DI32.

### Champs de force polarisables basés sur des dipôles induits

* CFF/ind et ENZYMIX. Le premier champ de force polarisable utilisé de manière importante dans de nombreuses applications aux systèmes biologiques.
* DRF90 développé par P. Th. van Duijnen et collaborateurs.

### Champs de force polarisables basés sur des charges ponctuelles

* PFF (Polarizable Force Field) développé par Richard A. Friesner et collaborateurs.
* Approche SP-basis Chemical Potential Equalization (CPE) développée par R. Chelli et P. Procacci.
* Champ de force polarisable CHARMM développé par S. Patel (Université du Delaware) et C.L. Brooks III (Université du Michigan).
* Champ de force polarisable AMBER développé par Jim Caldwell et collaborateurs.

### Champs de force polarisables basés sur des multipôles distribués

* SIBFA (Sum of Interactions Between Fragments Ab initio computed). Champ de force pour petites molécules et protéines flexibles, développé par Nohad Gresh (Université Paris V René Descartes) et Jean-Philip Piquemal (Université Paris VI Pierre et Marie Curie). SIBFA est un procédé de mécanique moléculaire formulé et calibré sur la base de calculs ab initio sur des supermolécules. Son but est de permettre des calculs fiables et simultanés à la fois sur les énergies conformationnelles et intermoléculaires qui induisent les propriétés de liaison spécifiques des molécules d'intérêt biologique et pharmaceutique. Cette procédure permet un traitement précis des métaux de transition. L'inclusion d'une contribution de champ de ligand permet des calculs sur les protéines métalliques « couches ouvertes ».
* AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications). Champ de force développé par Pengyu Ren (Université du Texas à Austin) et Jay W. Ponder (Université de Washington).
* ORIENT procedure. Champ de force développé par Anthony J. Stone (Université de Cambridge) et collaborateurs.
* Non-Empirical Molecular Orbital (NEMO) procédé développé par Gunnar Karlström et collaborateurs à l'Université de Lund (Suède).

### Champs de force polarisables basés sur la densité

* Gaussian Electrostatic Model (GEM)$^,$$^,$. Champ de force polarisable développé par Thomas A. Darden et G. Andrés Cisneros au NIEHS et Jean-Philip Piquemal (Université Paris VI Pierre et Marie Curie).
* Procédure de polarisation basée sur l'approche Kim-Gordon développée par Jürg Hutter et collaborateurs (Université de Zurich).

### Champs de force polarisables basés sur la théorie de la polarisation de liaison

* COSMOS-NMR. Ce champ de force hybride QM/MM permet des calculs en mécanique quantique explicites des propriétés électrostatiques en utilisant les orbitales de liaison localisées avec le formalisme de la théorie de la polarisation de liaison rapide. La fluctuation de charge atomique est possible pour chaque pas de dynamique moléculaire.

### Champs de force réactifs

* ReaxFF. Champ de force réactif développé par Adri van Duin, William Goddard et collaborateurs. Il est rapide, transférable et constitue une méthode numérique intéressante pour les simulations à l'échelle atomique de réactions chimiques.
* EVB (empirical valence bond). Ce champ de force réactif, introduit par Warshel et collaborateurs, est probablement l'une des méthodes les plus fiables et physiquement cohérente d'utilisation des champs de force de modéliser les réactions chimiques des différents environnements. L'EVB facilite les calculs des énergies libres d'activation dans les phases condensées et les enzymes.
* RWFF. Champ de force réactif pour l'eau développé par Detlef W. M. Hofmann, Liudmila N. Kuleshova et Bruno D'Aguanno. Le champ de force rapide, reproduisant les données expérimentales de diffusion de neutron de manière précise, et permet de simuler la formation ou la rupture de liaison de l'eau et des acides.

### Champs à gros grains

* Virtual atom molecular mechanics (VAMM). Champ de force développé par Korkut et Hendrickson. VAMM est un champ de force à gros grains spécifiquement construit pour des simulations de mécanique moléculaire comme des transitions conformationnelles à large échelle basées sur des interactions d'atomes Cα. C'est un champ de force basé sur des préconnaissances et formulé afin d'identifier les caractéristiques dépendant de la structure secondaire et sur l'information de contact spécifique aux résidus dans les protéines.

### Autres

* VALBOND. Fonction construite pour la modification d'angle basée sur la théorie de la liaison de valence et marche pour de grandes distorsions d'angles, les molécules hypervalentes et les complexes de métaux de transition. Elle peut être introduite dans d'autres champs de force (comme CHARMM et UFF).

### Modélisation de l'eau

L'ensemble des paramètres utilisé pour modéliser l'eau ou les solutions aqueuses (de manière basique un champ de force pour l'eau) est appelé modèle d'eau. L'eau a beaucoup attiré l'attention en raison de ses propriétés inhabituelles et de son rôle majeur comme solvant. De nombreux modèles d'eau ont été proposés, comme : TIP3P, TIP4P, SPC, Flexible SPC ou ST2.