Powered by OpenAIRE graph
Found an issue? Give us feedback
image/svg+xml art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos Open Access logo, converted into svg, designed by PLoS. This version with transparent background. http://commons.wikimedia.org/wiki/File:Open_Access_logo_PLoS_white.svg art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos http://www.plos.org/ ZENODOarrow_drop_down
image/svg+xml art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos Open Access logo, converted into svg, designed by PLoS. This version with transparent background. http://commons.wikimedia.org/wiki/File:Open_Access_logo_PLoS_white.svg art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos http://www.plos.org/
ZENODO
Dataset . 2025
License: CC BY
Data sources: ZENODO
ZENODO
Dataset . 2025
License: CC BY
Data sources: Datacite
ZENODO
Dataset . 2025
License: CC BY
Data sources: Datacite
versions View all 2 versions
addClaim

MD simulation run 4 of TARP gamma-7 from set 1

Authors: Krieger, James;

MD simulation run 4 of TARP gamma-7 from set 1

Abstract

Molecular dynamics simulations of pig TARP γ-7 were set up an AlphaFold 3 model for UniProt ID F2Z4Y2. This run belongs to a first set of 4 runs from a first γ-7 simulation system that was prepared directly from the AlphaFold model without any modification. It is run 4 of 4 and has a duration of 100 ns. The CHARMM-GUI Membrane Builder was used to prepare the simulation systems, including trimming the C-terminus at residue 206, neutralising the truncated C-terminus with N-methylamide cap (C-terminal patch CT3), orienting and positioning the protein in a bilayer with PPM, building a 75 Å by 75 Å POPC bilayer around the protein using the pseudoatom replacement method, adding water and 0.15 M Na+ and Cl- ions, performing initial energy minimisation in CHARMM, and preparing GROMACS input files. 4 MD simulations were run from the system prepared with CHARMM-GUI using GROMACS 2021.5 with the CHARMM36 all-atom force field, including the CHARMM36m update for proteins and the CHARMM TIP3P water model and standard parameters, such as particle mesh Ewald electrostatics and force switching for van der Waals with a cutoff of 1.2 nm and switching from 1.0 to 1.2 nm. The simulation protocol in GROMACS was the standard one created by CHARMM-GUI, run on a computer cluster using the scripts gromsubmit_cpu.sh and gromsubmit_cpu_prod_run1.sh based on the provided README file, which included slurm submission parameters, and a modified version of the provided step7_production.mdp file with a larger number of time steps, corresponding to 100 ns. Energy minimisation was run for 5000 steepest descent steps with position and dihedral restraints on the protein and lipid to relax the solvent around them. Harmonic position restraints were used for protein backbone and side-chain heavy atoms with force constants of 4000 and 2000 kJ·mol⁻¹·nm⁻² were in all directions, respectively. Lipids were treated with planar position restraints were used in the z-direction with force constants of 1000 kJ·mol⁻¹·nm⁻² to maintain lipid head groups in the planes of the membrane leaflets and harmonic dihedral restraints with force constants of 1000 kJ·mol⁻¹·nm⁻² to maintain correct orientation of phosphate groups and lipid tails relative to the core choline head group. These are close to the previously published values of 10.0, 5.0 and 2.5 kcal·mol⁻¹·Å⁻², followed by several rounds of equilibration of different duration, gradually reducing the restraints as previously described [ref 2]. All steps after energy minimisation used the Bussi-Donadio-Parrinello stochastic velocity rescaling thermostat with three coupling groups corresponding to protein, lipid and solvent atoms and a time constant of 0.1 ps and a target temperature of 303.15 K. The first two equilibration steps were in the NVT (constant number of atoms, volume and temperature) ensemble, only equilibrating the temperature. This was followed by four rounds of NPT (constant number of atoms, pressure and temperature) equilibration, adding semi-isotropic pressure equilibration with the Bernetti-Bussi stochastic cell rescaling barostat with a time constant of 5.0 ps and compressibility of 4.5 x 10-5 bar-1, equilibrating the system to a pressure of 1.0 bar, which was also maintained for all steps after NVT equilibration. Finally, the restraints were removed and production MD simulations were run in the NPT ensemble with the same thermostat and barostat with the same parameters. The first three equilibration steps used 1-fs time steps, after which 2-fs time steps where used. Bonds containing hydrogen were constrained with LINCS for all steps, including minimisation. 

Related Organizations
  • BIP!
    Impact byBIP!
    selected citations
    These citations are derived from selected sources.
    This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically).
    0
    popularity
    This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network.
    Average
    influence
    This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically).
    Average
    impulse
    This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network.
    Average
Powered by OpenAIRE graph
Found an issue? Give us feedback
selected citations
These citations are derived from selected sources.
This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically).
BIP!Citations provided by BIP!
popularity
This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network.
BIP!Popularity provided by BIP!
influence
This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically).
BIP!Influence provided by BIP!
impulse
This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network.
BIP!Impulse provided by BIP!
0
Average
Average
Average