Al examinar un sistema computacionalmente, puede interesarle saber cómo cambia la energía libre en función de alguna coordenada intermolecular o intramolecular (como la distancia entre dos átomos o un ángulo de torsión). La superficie de energía libre a lo largo de la coordenada elegida se denomina potencial de fuerza media (PMF). Si el sistema de interés está en un solvente, entonces el PMF también incorpora los efectos del solvente. [1]
El PMF se puede obtener en simulaciones de Monte Carlo o de dinámica molecular para examinar cómo cambia la energía de un sistema en función de algún parámetro de coordenadas de reacción específico. Por ejemplo, puede examinar cómo cambia la energía del sistema en función de la distancia entre dos residuos, o a medida que una proteína es arrastrada a través de una bicapa lipídica. Puede ser una coordenada geométrica o una coordenada energética (disolvente) más general. A menudo, las simulaciones PMF se utilizan junto con el muestreo paraguas , porque normalmente la simulación PMF no logrará muestrear adecuadamente el espacio del sistema a medida que avanza. [2]
El Potencial de Fuerza Media [3] de un sistema con N partículas es por construcción el potencial que da la fuerza media sobre todas las configuraciones de todas las n+1...N partículas que actúan sobre una partícula j en cualquier configuración fija manteniendo fijo un conjunto de partículas 1...n
Arriba, es la fuerza promedio, es decir, la "fuerza media" sobre la partícula j . Y es el llamado potencial de fuerza media. Para , es el trabajo promedio necesario para llevar las dos partículas desde una separación infinita a una distancia . También está relacionada con la función de distribución radial del sistema, , por: [4]
El potencial de fuerza media se aplica generalmente en el método de inversión de Boltzmann como una primera estimación del potencial de interacción de pares efectivos que debería reproducir la función de distribución radial correcta en una simulación mesoscópica. [5] Lemkul et al. han utilizado simulaciones de dinámica molecular dirigida para calcular el potencial de fuerza media para evaluar la estabilidad de las protofibrillas amiloides de Alzheimer. [6] Gosai et al. también han utilizado simulaciones de muestreo de paraguas para demostrar que el potencial de fuerza media disminuye entre la trombina y su aptámero (un complejo proteína-ligando) bajo el efecto de campos eléctricos. [7]