Analysis scripts for molecular dynamics simulations of stratum corneum lipid multilayers.
The package can be installed via pip
- Clone the repository
git clone https://github.com/uppittu11/analysis.git- Install the package
cd analysis
pip install -e .- Run tests for the package
pytest -v analysisI suggest you install this package in a separate conda environment so the analyze command does not clobber anything in your default namespace
To calculate basic structural properties (tilt angle, repeat distance, nematic order, and area per lipid) using a parallel multiprocessing pool, use the bash command analyze.
-f trajectory file (needs to be a format loadable by MDTraj)
-c topology file (needs to be a format loadable by MDTraj)
-o output directory
-n number of leaflets
--cg use this flag to denote a coarse-grain system
--reload use this flag to ignore cached trajectory (typically located in the output directory)
--min minimum z position of lipid headgroups to include in analysis
--max maximum z position of lipid headgroups to include in analysis. Use --min and --max to specify one or more layers in a multilayer system
analysis will output a pickled dictionary file containing structural properties calculated for each frame of the trajectory (default: outputDirectory/results.p. This output file can be loaded and analyzed using functions in the analysis.data module.
- Tilt Angles
- Bilayer/Multilayer Height
- Nematic Order Parameter
- Area per Lipid
- Area per Tail