Efecto isotópico en el proceso de transferencia de doble protones de porficeno investigado por el método QM/MM mejorado

Resumen

Aquí se presenta un protocolo que utiliza el método QM/MM mejorado para investigar el efecto isotópico en el proceso de transferencia de doble protones en porficeno.

Resumen

La sustitución de deuterio único en porficeno conduce a una geometría molecular asimétrica, que puede afectar al proceso de transferencia de doble protones en la molécula de porficeno. En este estudio, aplicamos un método mejorado QM/MM llamado SITS-QM/MM para investigar los efectos del isótopo de hidrógeno/deuterio (H/D) en la transferencia de doble protones en porficéo. Los cambios de distancia en las simulaciones de dinámica molecular SITS-QM/MM sugirieron que el deuterio sustituyó al porficeno adoptó el mecanismo de transferencia de doble protones escalonado. El análisis estructural y los cambios de energía libre del proceso de transferencia de doble protones indicaron que la sustitución isotópica asimétrica comprimió sutilmente los enlaces de hidrógeno covalentes y puede alterar la ubicación original del estado de transición.

Introducción

El proceso de transferencia de protones en porficéos tiene aplicaciones potenciales en el desarrollo de interruptores moleculares, transistores y dispositivos de almacenamiento de información^1,2. En particular, la tautomerización en porficéos a través del proceso de transferenciade doble protones ha atraído un amplio interés en los campos de la espectroscopia y la fotofísica 2. Los átomos de hidrógeno interno de porficeno pueden migrar de un isómero trans al otro isómero trans equivalente a través del proceso de transferencia de doble protones como se muestra en la Figura 1. Se han propuesto dos mecanismos para el proceso de transferencia de doble protones: el mecanismo concertado y el mecanismo escalonado^3,4. En el proceso concertado de transferencia de doble protones, ambos átomos de protones se mueven al estado de transición sincrónicamente de una manera simétrica, mientras que un protón completa la transferencia antes que el otro protón en un proceso escalonado. Dos átomos de hidrógeno pueden transferirse simultáneamente o escalonadamente dependiendo de la fuerza de correlación entre dos átomos de hidrógeno⁵.

La sustitución isotópica se ha utilizado para detectar las propiedades estructurales de las moléculas y velocidad constantes de la cinética de reacción⁶. La sustitución de deuterio único en el hidrógeno interno del porficeno conduce a una forma asimétrica de la molécula. El enlace de hidrógeno puede expandirse o contraerse debido a la diferencia de masa entre los átomos de hidrógeno y deuterio. La sustitución isotópica introduce una perturbación en el andamio del porficeno. Se plantea la cuestión de si la estructura asimétrica afectaría el proceso de transferencia de protones. Limbach y compañeros de trabajo informaron que la sustitución del hidrógeno por deuterio comprimirá tanto los enlaces de hidrógeno, como el acoplamiento cooperativo de dos enlaces de hidrógeno en porficeno puede favorecer el mecanismo concertado^7,mientras que Yoshikawa declaró que el deuteración haría que el mecanismo escalonado contribuyera más que el mecanismo concertado⁸. Se han desarrollado técnicas experimentales, como la espectroscopia de fuerza, paracapturar detalles de tautomerización en un solo porficeno 9. Sin embargo, sigue siendo difícil determinar los detalles atómicos de la transferencia de protones experimentalmente debido a su naturaleza transitoria.

Los cálculos teóricos y las simulaciones pueden actuar como herramientas complementarias para esclarecer los mecanismos de reacción de la transferencia de protones. Entre los diferentes métodos teóricos, las simulaciones de dinámica molecular (MD) pueden monitorear los movimientos dinámicos de cada átomo, y se ha utilizado ampliamente para revelar mecanismos complejos en reacciones químicas y enzimáticas. Sin embargo, las simulaciones MD regulares tienden a sufrir problemas de muestreo insuficientes, especialmente cuando existe una barrera de energía alta en el proceso de interés. Por lo tanto, se han desarrollado métodos de muestreo mejorados, que incluyen el muestreo de ruta de transición^10,11, muestreo paraguas (EE.UU.)^12,13, y muestreo de templado integrado (ITS)^{14, 15}. La combinación de diferentes métodos de muestreo mejorados puede aumentar aún más la eficiencia de muestreo^16,17,18. Para aprovechar los algoritmos de muestreo mejorados en la simulación de reacciones químicas, hemos implementado el método de muestreo de templado integrado selectivo (SITS) con potencial mecánico cuántico y mecánico molecular (QM/MM) recientemente^19. El método SITS-QM/MM propuesto combina las ventajas de ambos métodos: el método SITS acelera el muestreo y puede explorar todos los canales de reacción posibles sin conocimiento previo del mecanismo de reacción, y QM/MM proporciona una descripción más precisa de el proceso de formación de unión y ruptura de unión, que no puede ser simulado por los métodos MM únicamente. El enfoque SITS-QM/MM implementado ha descubierto con éxito la transferencia concertada de doble protones, mecanismo de transferencia de doble protones no correlacionado y correlacionado en diferentes sistemas, sin predefinir las coordenadas de reacción¹⁹. En el caso del porficeno, se ha notificado el carácter de transferencia de protones escalonapero pero correlacionado¹⁹. El método híbrido SITS-QM/MM se utilizó para investigar el efecto isotópico en el porficeno en nuestro estudio, y a continuación se muestran las descripciones detalladas del algoritmo y el protocolo de nuestro método.

Hemos implementado el método SITS con potencial sin QM/MM híbrido. Se definió el potencial efectivo de SITS para incluir la energía potencial a diferentes temperaturas con los factores de ponderación n_k para cubrir rangos de temperatura más amplios,

donde, N es el número de términos canónicos, la temperatura inversa es la temperatura inversa, y n_k es el factor de ponderación correspondiente para cada componente canónico. U_E (R) y U_N(R) representan los términos mejorados y no mejorados en SITS y se definen como,

U _s, U_se y U_e son la energía potencial del subsistema, la interacción entre el subsistema y el medio ambiente, y la energía potencial del medio ambiente. El potencial QM/MM se expresa como una suma híbrida de tres componentes,

donde U_qm, U_qm/mmy U_mm son el término de energía interna del subsistema QM, la energía de interacción entre las regiones QM y MM y la energía de interacción dentro del subsistema MM, respectivamente. El término U_qm/mm se puede dividir en tres componentes, que incluyen los términos electrostáticos, van der Waals y energía de interacción covalente entre los átomos QM y MM,

Asignamos, y en un término de U_s en SITS,

El potencial total del sistema se descompuso entonces en la energía del subsistema U s, la energía de interacción entre el subsistema y el medio ambiente U_se, y la energía del medio ambiente U_e. Por ejemplo, en el sistema de la obra actual, el subsistema es la porficencia, y el medio ambiente el agua.

El perfil PMF a lo largo de una variable colectiva (R) se deriva como,

Las coordenadas de reacción generalmente utilizadas para cada transferencia de hidrógeno de N1-H1- N2 son q₁ o (r₁ar₂)/2 y q₂ a r₁ + r₂, donde se encuentran₁ es la distancia de N1-H1 y r₂ es la distancia de H1-N2.

El método se ha implementado en el paquete de simulación QM/MM MD QM⁴D²⁰. El código fuente completo y la documentación se pueden encontrar aquí: http://www.qm4d.info/.

Generalmente, las simulaciones SITS-QM/MM MD implican cuatro pasos: pre-equilibrio (pre-sits); optimización n_k (opt-sits); simulación de producción y análisis de datos.

Protocolo

1. Modelo de construcción

1. Construir la estructura del porficeno: Abra el software GaussView haciendo doble clic en el ratón. A continuación, haga clic en el botón Fragmento de elemento en el menú de GaussView para elegir los elementos necesarios. Construir pórfido. A continuación, haga clic en el botón Archivo para guardar como el archivo pdb.
2. Solvate el modelo: Solvate porficene en una caja de agua cúbica TIP3P²¹ con una longitud de borde de 38o mediante la emisión del comando en el sistema operativo linux: genbox_d -cp prp-vac.pdb -cs spc216.gro -o solv.pdb -maxsol 1484 -box 3.8.
3. Construir deuterio porficeno: Ejecute el siguiente comando para generar el archivo de topología: cns < ppi_solv.inp. A continuación, abra prp-wat.psf con el comando vi y cambie la masa de H1 de 1.00800 a 2.01600 para reemplazar un átomo de hidrógeno intramolecular de porficeno con deuterio para construir un solo deuterio sustituido por el porficeno.
4. Establezca los parámetros normales de simulación de MD: Método de entrada scctb, integral 0.5 fsy cutoff 12 en el archivo de entrada MD abriéndolo con el comando vi.
   NOTA: Adoptar una distancia de corte de 12o para calcular las interacciones tanto vdW como electrostáticas. Simular la molécula de porficeno con el método DFTB/MIO²². Establezca el paso de tiempo de integración como 0,5 fs para simulaciones MD. Mantenga la temperatura del sistema de simulación a 300 K con el termostato Langevin. A continuación, realice las simulaciones con el software QM⁴D siguiendo los pasos a continuación.

2. Pre-sits

1. Configurar parámetros de temperatura: Entrada templow 260, temphigh 1100 y ntemp 160 en el archivo de entrada.
   NOTA: El rango de temperatura de 260K a 1100K fue distribuido por el software QM4D a 160 puntos de temperatura durante las simulaciones MD. Los archivos de entrada de plantilla se incluyen en Archivos suplementarios.
2. Inicie las pre-sits: Establezca runtype 100 y step 120,000 en el archivo de entrada. A continuación, ejecute el siguiente comando: $PATH/qm4d $INPUTFILE > $OUTPUTFILE.
   NOTA: El paso total es de 120.000, pero se puede ajustar en función de la necesidad específica. Los parámetros sugeridos de las simulaciones MD se guardan en el $INPUTFILE. El mismo comando también se utiliza en los siguientes pasos de opt-sits y simulación de producción, con el archivo de entrada modificado en consecuencia.
3. Cálculo de las energías descompuestas
   1. Extraer cambios de energía: Durante la etapa de pre-sits, monitoree la energía de cada término para calcular los valores medios, como se muestra en la Figura1. Utilice el comando grep Linux para extraer la energía de la siguiente manera:
      grep 'SITS-ener0' $INPUTFILE á awk ''a+'$3;b+'$4;c+'$5'END'print a/NR,b/NR,c/NR'.
   2. Modificar energías medias en el archivo de entrada MD: Calcule las energías medias basadas en la salida de la línea de comandos anterior, y modifique la línea vshift0 -30801.95; vshift1 -26.88; vshift2 -13888.28 en el archivo de entrada con los promedios recién generados.
      NOTA: Los números -30801.85, -26.88 y -13888.28 son las energías promedio en el sistema modelo actual. Modifique los valores en función de los sistemas específicos.

3. Opt-sits

1. Inicie opt-sits: Establezca runtype 0 en el archivo de entrada. A continuación, inicie el programa QM⁴D escribiendo el comando como se muestra en el paso 2.2 para iniciar el paso de optimización.
2. Monitoreo de los cambios energéticos y n_k valores.
   1. Trazar la propagación de energía con el programa de "gracia" y asegurarse de que la fluctuación de energía puede cubrir los extremos más bajos y más altos del rango de temperatura.
   2. Después de la optimización, guarde los valores n_k finales del paso opt-sits en un nuevo archivo, que se denomina nk.dat en este protocolo.

4. Ejecución de simulaciones de producción

1. Preparar el archivo de entrada MD: Establezca runtype 1 en el nuevo archivo de entrada para iniciar el paso de simulación de producción. Especifique el nombre de archivo del archivo n_k almacenado como nkfile nk.dat en el archivo de entrada. El número de pasos de tiempo se estableció como 6.400.000 en los sistemas actuales.
2. Iniciar simulación MD de producción: Ejecute el siguiente comando para iniciar simulaciones MD: $INPUTFILE $PATH/qm4d > $OUTPUTFILE.
   NOTA: Asegúrese de que el software QM⁴D puede leer los valores n_k. El tiempo de simulación depende del sistema, así que cambie el paso de simulación en función de demandas específicas. Seleccione un número adecuado de pasos de tiempo para garantizar suficiente tiempo de simulación para su propio sistema. Este paso probablemente consume mucho tiempo, así que guarde los archivos de reinicio para evitar reiniciar la producción desde el principio una vez que se interrumpa.

5 Análisis de datos

1. Seguimiento de los cambios de distancia
   1. Supervise el proceso de formación y rotura de enlaces durante la fase de producción, utilice el comando grep para comprobar los cambios de distancia de H1-N1 y H1-N2 a lo largo del tiempo de simulación. La misma operación se puede llevar a cabo para H2-N3 y H2-N4. A continuación, trace la propagación de distancia utilizando el valor de distancia acumulado durante las simulaciones de producción.
2. Extraer coordenadas de reacción
   1. Extraiga las coordenadas de reacción y los términos de energía del archivo de salida de producción generado por el comando QM⁴D por grep:
      grep 'dist 1' $OUTPUTFILE awk ''imprimir $5'' > distancia1;
      grep 'ener0' $OUTPUTFILE > ener0.
   2. Organice los datos en cuatro columnas: q₁, q₂, U₀ y U' (U₀ y U' son la energía normal de salida y la energía ponderada), y escríbalos en el archivo de datos en cada período de tiempo.
3. Cálculo de la energía libre
   1. Calcule la energía libre emitiendo el siguiente comando:
      sits-pmf 300 $INPUTFILE PMF2 [hist_minx hist_maxx num_binsx] [hist_miny hist_maxy num_binsy] > $OUTPUTFILE.
      NOTA: sits-pmf es el método de análisis basado en histograma. [hist_minx hist_maxx num_binsx] define el intervalo y el número de ubicaciones para la primera coordenada de reacción. Las coordenadas de reacción de la segunda se pueden establecer mediante [hist_miny hist_maxy num_binsy].
   2. Para proyectar la energía libre en el paisaje bidimensional, escriba el siguiente comando:
      sits-pmf 300 h1-2d.dat PMF2 -0.6 0.6 24 2.45 4.25 36 > sits-pmf.out.
      NOTA: Utilice un total de 24 bins y 36 bins para cubrir los cambios de distancia en dos coordenadas de reacción seleccionadas, q₁ y q₂, respectivamente. Guarde los datos PMF 2D para cada hidrógeno/deuterio en el archivo sits-pmf.out.

Resultados

El efecto de sustitución de deuterio único en el proceso de transferencia de doble protones en porficeno se examinó en el protocolo actual (Figura1). Se comprobó la energía potencial del subsistema QM y el agua durante el paso de preequilibrio y optimización para asegurarse de que la energía se ha ampliado a un rango de energía más amplio (Figura2). Los cambios representativos de distancia y ángulo (Figura3 y

Discusión

La estructura del porficeno se mostró en la Figura1. El potencial híbrido QM/MM de incrustación electrostática con el método SITS se utilizó para describir las reacciones químicas en el agua^23,24. La transferencia de protones se produce dentro del porficeno³ y por lo tanto el porficeno se establece como región QM y el agua de recordatorio se establece como región MM. Aquí adoptamos DFTB/MIO como nues...

Materiales

Name	Company	Catalog Number	Comments
operating system	CentOS Linux release 6.0
QM4D software	http://www.qm4d.info/		in-house program
Computer desktop	HP