Team:Aalto-Helsinki/Model Setup
Computational methods and simulation set-up
Molecular modeling methods ultimately aim for direct comparison with experimental measurements. As such, a good model of molecular interaction is essential. Quantum chemistry based Ab initio molecular dynamics methods aim to reduce the amount of �tting and guesswork required for accurate modeling of molecular interactions. However, such approaches are generally limited to small systems and short timescales due to the added computational demand[15]. In classical molecular dynamics, molecules are described using stick-and-ball models: spherical atoms are connected by springs that represent bonds. As such, internal molecular forces can be described by simple mathematical models. For example, Hooke's law can be used to describe bonded interactions, while non-bonded interactions can be described by Lennard-Jones potential.[13, 18]
Molecular dynamics is based on numerical, step-by-step, evaluation of Newton's equations of motion.
Due to the many-body nature of the problem, Newton's equations of motion are discretized and solved
numerically. MD trajectories, that describe the time evolution of the system, consist of both position
and velocity vectors of the particles in the system. The position and velocity vectors are reevaluated
according to �nite time interval by using numerical integrators. The position vectors de�ne the geometric
con�guration of the system while the velocity vectors de�ne the kinetic energy and temperature of the
system.[13, 18, 15]
We investigated the behavior of DCD-1L in aqueous solvents at varying salt concentrations and tem- peratures using molecular dynamics. the simulations make use of the GROMOS force�eld parameter set 53a6[17]. GROMOS53A6 is considered a united-tom force�eld, which maps a carbon and its associated apolar hydrogens as one interaction center. The parametrization of the force�eld is based on accurate reproduction of free enthalpies of hydration and apolar solvation of a wide range of compounds. The relative free enthalpy of solvation between apolar and polar environments is an important property in many biological phenomena, including protein folding and membrane formation.
System set-up
We probed the behavior of a single DCD-1L solvated in water with added 2 Na+ counter-ions to neutralize the net negative charge of the peptide, as well as added salt (NaCl), with the salt concentrations varying from 20 mM to 500 mM. Additionally we probed the temperature dependency of DCD-1L structure by simulating a single peptide in water with 2 Na+ counter-ions at temperatures of 290 K, 300 K, 310 K and 320 K. All the simulation trajectories were calculated using the Gromacs v.4.6.7 simulation package[19, 20].
As starting structure for our simulations, chose the helical crystal structure of dermcidin derived by Song et al.[6] available in the RCSB protein Data Bank (entry 2YMK). Due to the presence of missing atoms in the original model by Song et al., we constructed a a homology model based on the DCD-1L sequence used in our laboratory constructs using the SWISS-MODEL automated protein structure homology-modeling server[21, 22, 23, 24]. The resulting structure was then mapped to the GROMOS53A6 force�eld using the Gromacs 4.6.7 tool pdb2gmx.
In all simulations, the dimensions of the simulation box were cubic with a side length of 11.2 nm. The peptide was inserted and centered into the simulation box and solvated in water. For all simulations, the simple point-charge (SPC) water model was used. After solvation, two of the water molecules were replaced by 2 Na+ ions to act as counter-ions and e�ectively neutralize the charge on the simulated system. Additionally, for simulations involving added salt, an adequate number of water molecules were replaced by an equal number of Na+ and Cl ions to mimic the desired salt concentration of either 20 mM, 50 mM, 70 mM, 100 mM, 120 mM, 150 mM, 200 mM, 300 mM or 500 mM. Visualizations of the simulated system set-up con�gurations are presented in �gure X.
