exercises:2015_cecam_tutorial:urea
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
exercises:2015_cecam_tutorial:urea [2015/08/19 14:22] – fix file links and some markup tmueller | exercises:2015_cecam_tutorial:urea [2020/08/21 10:15] (current) – external edit 127.0.0.1 | ||
---|---|---|---|
Line 2: | Line 2: | ||
Problem: QM/MM study of the Urea Zwitterion in water by means of a QM/MM Hamiltonian. | Problem: QM/MM study of the Urea Zwitterion in water by means of a QM/MM Hamiltonian. | ||
+ | |||
+ | * Original author: Marcella Iannuzzi | ||
+ | * Complete source and output files: [[http:// | ||
===== Introduction ===== | ===== Introduction ===== | ||
Line 25: | Line 28: | ||
* '' | * '' | ||
* '' | * '' | ||
- | * '' | + | * '' |
The tasks we will complete in this tutorial exercise are: | The tasks we will complete in this tutorial exercise are: | ||
Line 32: | Line 35: | ||
* Based on the averages of the NPT we will equilibrate at the NVT level with an average simulation box | * Based on the averages of the NPT we will equilibrate at the NVT level with an average simulation box | ||
* One the system is equilibrated at the classical Hamiltonian level we will switch to a QM/MM Hamiltonian, | * One the system is equilibrated at the classical Hamiltonian level we will switch to a QM/MM Hamiltonian, | ||
- | * Study the chemical reactivity of the Zwitterion in solution: by inspecting the possibility of having a reversal reaction with formation of neutral urea or alternatively the elimination reaction, with formation of NH$_3$ | + | * Study the chemical reactivity of the Zwitterion in solution: by inspecting the possibility of having a reversal reaction with formation of neutral urea or alternatively the elimination reaction, with formation of < |
===== Theoretical Background ===== | ===== Theoretical Background ===== | ||
Line 38: | Line 41: | ||
Urea is formed in large quantities as a product of catabolism of nitrogen-containing compounds. Owing to its resonance stabilization, | Urea is formed in large quantities as a product of catabolism of nitrogen-containing compounds. Owing to its resonance stabilization, | ||
- | Cyanate ion further readily undergoes conversion to CO2 and ammonia. In contrast, when catalyzed by ureases, urea is generally believed to undergo hydrolysis rather then ammonia elimination producing either HCO3- and NH4+ or ammonium carbamate, depending on the buffer system. Activation energies for urea decomposition in water at different pH have been obtained experimentally. For neutral pH, the reported activation energy ranges from 28.4 Kcal/mol to 32.4 Kcal/mol. There have been also numerous theoretical investigations of the decomposition of urea and related systems. In all of them the explicit representation of the solvent was found to be essential for detailed resolution of the mechanism, identification of the rate determining step and evaluation of the barrier. In particular in a very recent paper by Jorgensen et al. the hydrogen bonded water molecules were found to act as hydrogen shuttle for the first step of the elimination reaction. The forming zwitterionic intermediate, | + | Cyanate ion further readily undergoes conversion to < |
The goal of this exercise will be to inspect the chemical reactivity of the zwitterion surrounded by water molecules. In this exercise, in a very simplicistic way, only urea will be considered QM while the rest of the system will be described with a classical Hamiltonian. By some high level calculations, | The goal of this exercise will be to inspect the chemical reactivity of the zwitterion surrounded by water molecules. In this exercise, in a very simplicistic way, only urea will be considered QM while the rest of the system will be described with a classical Hamiltonian. By some high level calculations, | ||
Line 60: | Line 63: | ||
&END CELL | &END CELL | ||
& | & | ||
- | CONN_FILE_NAME | + | CONN_FILE_NAME ${ROOT}/ |
CONNECTIVITY AMBER | CONNECTIVITY AMBER | ||
- | COORD_FILE_NAME | + | COORD_FILE_NAME ${ROOT}/ |
COORDINATE CRD | COORDINATE CRD | ||
&END TOPOLOGY | &END TOPOLOGY | ||
Line 74: | Line 77: | ||
&MM | &MM | ||
& | & | ||
- | parm_file_name | + | parm_file_name ${ROOT}/ |
parmtype AMBER | parmtype AMBER | ||
&spline | &spline | ||
Line 90: | Line 93: | ||
&END | &END | ||
</ | </ | ||
- | In the FORCE_FIELD | + | In the '' |
- | The core of the evaluation in a classical run, is the evaluation of the electrostatic. We can adjust these parameters in the &POISSON section (similarly to the DFT calculations). For classical runs we can employ either standard EWALD summations, Particle-Mesh Ewald (PME) sums or Smooth-Particle-Mesh Ewald ones (SPME). | + | The core of the evaluation in a classical run, is the evaluation of the electrostatic. We can adjust these parameters in the '' |
For this exercise we emply the SPME with a grid mesh of 54 for all 3 dimensions and the $\alpha$ parameter for the reciprocal space contributions is equal to 0.4. | For this exercise we emply the SPME with a grid mesh of 54 for all 3 dimensions and the $\alpha$ parameter for the reciprocal space contributions is equal to 0.4. | ||
- | The control of the NPT equilibration is specified instead by the MD section: | + | The control of the NPT equilibration is specified instead by the '' |
< | < | ||
Line 137: | Line 140: | ||
===== Second task: MM isothermal ensemble ===== | ===== Second task: MM isothermal ensemble ===== | ||
- | Using the average parameters of the cell lattice, as determined in the previous run, we setup an input file to run an NVT equilibration, | + | Using the average parameters of the cell lattice, as determined in the previous run, we setup an input file to run an NVT equilibration, |
< | < | ||
Line 159: | Line 162: | ||
Starting from the MM system, equilibrated at the right pressure and temperature, | Starting from the MM system, equilibrated at the right pressure and temperature, | ||
- | All the informations about a QM/MM run are specified in the QMMM section being part of the FORCE_EVAL%QMMM. In particular, we need to specify first the QM CELL. This is mandatory and important for DFT calculations (performance, | + | All the informations about a QM/MM run are specified in the QMMM section being part of the '' |
< | < | ||
Line 190: | Line 193: | ||
For this specific example, we need also to introduce an additional modification in the force-field. In fact, in the classical force-field, | For this specific example, we need also to introduce an additional modification in the force-field. In fact, in the classical force-field, | ||
- | We can do that with an additional FORCEFIELD section inside the QMMM one: | + | We can do that with an additional |
< | < | ||
Line 224: | Line 227: | ||
==== Homeworks ==== | ==== Homeworks ==== | ||
- | Try to convert this input to use GPW. Hints: when setting-up a correct | + | Try to convert this input to use GPW. Hints: when setting-up a correct |
===== Fourth task: QM/MM Metadyanamics simulations ===== | ===== Fourth task: QM/MM Metadyanamics simulations ===== | ||
Line 288: | Line 291: | ||
==== Elimination ==== | ==== Elimination ==== | ||
- | The elimination reaction is sampled along the CV representing the bond between the NH$_3$ | + | The elimination reaction is sampled along the CV representing the bond between the < |
< | < | ||
Line 301: | Line 304: | ||
===== Questions ===== | ===== Questions ===== | ||
- | Evaluate the free-energy for both processes: from Zwitterionic to Neutral form and for the elimination pathway. How do these numbers compare with the 5 kcal/mol predicted in several published works? Inspect carefully the metadynamics trajectory in order to find a solution (what is the first attempt of the hydrogen of NH3(+) before moving towards the NH(-) group? what would happen if the nearby water molecules would be treated QM?). | + | Evaluate the free-energy for both processes: from Zwitterionic to Neutral form and for the elimination pathway. How do these numbers compare with the 5 kcal/mol predicted in several published works? Inspect carefully the metadynamics trajectory in order to find a solution (what is the first attempt of the hydrogen of < |
===== Homeworks ===== | ===== Homeworks ===== | ||
Take into account a primary solvation shell of water molecules as a part of the QM subsystem, using a FLEXIBLE_PARTITIONING scheme to prevent the diffusion of the QM water molecules. Re-run the equilibration steps and perform the Zwitterionic-Neutral metadynamics. Do you see any change in the barrier energy? Why? | Take into account a primary solvation shell of water molecules as a part of the QM subsystem, using a FLEXIBLE_PARTITIONING scheme to prevent the diffusion of the QM water molecules. Re-run the equilibration steps and perform the Zwitterionic-Neutral metadynamics. Do you see any change in the barrier energy? Why? |
exercises/2015_cecam_tutorial/urea.1439994133.txt.gz · Last modified: 2020/08/21 10:14 (external edit)