Schrodinger software was used to prepare the Panx2 model, including building the missing loop (amino acids 185–187, 267–272, and 374–376) via the Prime
60 (link), adding protons to all amino acids via PROPKA
61 (link), and energy-minimizing the system with the OPLS4 force field
62 (link) to ensure that no positional conflicts were in the model. Subsequently, a web-based platform for generating the inputs for molecular dynamics named CHARMM-GUI
63 (link) was applied to build the system. The N- and C-termini of the model were treated as amino termini with a positive charge and carboxyl termini with a negative charge, respectively, as they were freely exposed to the solvation. Two intramolecular disulfide bonds (C81–C279, C99–C259) identified in the experimental structures were created for each chain. The state of the protein after the above processes was regarded as its initial geometry in our simulations. The orientation of the protein and the position of lipid bilayers were determined by the Positioning of Proteins in Membranes (PPM) 2.0 server
64 (link) and checked manually. The lipid bilayers with the 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) model and water molecules with the TIP3P model were generated, resulting in a 120 × 120 × 165 Å box with 200 POPC molecules. The pore water molecules (water in the ion channel) were also added based on the channel geometry. Finally, the entire system was neutralized with 150 mM NaCl.
GPU-accelerated Gromacs 2021
65 (link) was used to perform the molecular dynamics simulation and CHARMM36m force field
66 (link) for the molecules. The MD simulation consisted of energy minimization, pre-equilibration and production simulations. The system was first energy minimized with the steepest descent algorithm while keeping 4000 kJ/(mol nm
2) force constant on backbone atoms and ligand atoms, 2000 kJ/(mol nm
2) force constant on side-chain atoms and 1000 kJ/(mol nm
2) force restraint on lipid atoms. Then, six-step pre-equilibration simulations (0.6 ns, 0.6 ns ps, 1 ns, 1 ns, 1 ns, and 1 ns) were carried out, where restraint was reduced slowly (4000, 2000, 1000, 500, 300, 0 kJ/(mol nm
2) on the backbone and ligand atoms, 2000, 1000, 500, 200, 50, 0 kJ/(mol nm
2) on side-chain atoms, and 1000, 400, 400, 200, 40, 0 kJ/(mol nm
2) on lipid atoms) to relax the system. Finally, a production simulation was performed for 100 ns using Langevin thermostat, with a constant temperature of 310 K and a constant pressure of 1 atm. Periodic boundary conditions (PBCs) were introduced during the all molecular dynamics simulations. In all steps, the time step was 2 fs, and atomic coordinates were written every 5 ps. After the MD simulation, root mean squared deviations, root mean squared fluctuations and distances between atoms were analyzed by Gromacs. Details of MD simulations have also been deposited in github at [
https://github.com/shiyu-wangbyte/panx2-simulation].
Zhang H., Wang S., Zhang Z., Hou M., Du C., Zhao Z., Vogel H., Li Z., Yan K., Zhang X., Lu J., Liang Y., Yuan S., Wang D, & Zhang H. (2023). Cryo-EM structure of human heptameric pannexin 2 channel. Nature Communications, 14, 1118.