Skip to main content

3D molecular structure generation for MD simulation

Project description

pygen-structures

Documentation StatusBuild Statuscodecovstatus

pygen-structures (pigeon structures) is a Python utility which allows for the generation of 3 dimensional molecular structures which can be used in molecular dynamics or Monte Carlo simulations. Molecules are generated from a list of residues and patches in the format of the CHARMM forcefield, and can be written out as valid PSF and PDB files. The package can be used as a command line utility, or as a Python library.

pygen-structures can be installed using pip (pip install pygen-structures), but relies upon RDKit, which is not pip-installable. RDKit can be installed in many ways, but the easiest way is to use the conda package manager. For full installation instructions, see the 'Installation' section of the readme. Python 3.6 and 3.7 are supported. To run the tests, pytest and OpenMM are required.

In essence, pygen-structures aims (eventually) to be a complete psfgen replacement with more autonomous functionality. At present, structures for small molecules can be generated. This should make it significantly easier to perform combinatorial searches on particular sequence lengths and linkages. This requires no manual intervention provided the molecules of interest are reasonably small (small enough that embedding coordinates is possible, and that the secondary structure is not vitally important) and the residue/patch definitions already exist in the forcefield.

Installation

There are other ways to install the required dependencies, but the easiest way by far is to use conda. Instructions, include the installation of test dependencies, are outlined below:

  1. Install the conda package manager. Make sure the conda executable is in your PATH.
  2. Set up a conda environment with the relevant dependencies (or install them in your base distribution). This can be done with the following command: conda create -n pygen-structures -c rdkit -c omnia 'python>=3.6' 'rdkit>=2018.3' numpy 'openmm>=7.4' pytest.
  3. Activate the conda environment: conda activate pygen-structures
  4. Use pip to install pygen-structures in this environment: pip install pygen-structures.
  5. Installation complete! You will have to activate this environment using conda activate pygen-structures each time you want to use it.
  6. Test the installation: pytest --pyargs pygen_structures

To install only the runtime dependencies, use the following command in step 2: conda create -n pygen-structures -c rdkit 'python>=3.6' 'rdkit>=2018.3' numpy

Command line usage

Command line usage for peptides is simple, and takes the following form:

pygen-structures SEQUENCE -o OUTPUT_PREFIX

Sequences are specified using the one letter protein code by default, and terminal patches can be supplied by using hyphens as delimiters (e.g. NNEU-AFK-CT2, note that both termini must be supplied). D-amino acids need only be preceded by a lowercase 'd'.

OUTPUT_PREFIX.psf and OUTPUT_PREFIX.pdb are created. If -o is not specified, the PDB file is written to stdout and no PSF file is generated.

The histidine form used can be set using --histidine, and defaults to HSE.

Patches can be supplied using --patches, the name of the patch, and the 0-based indices the patch is to be applied to (or 'FIRST'/'LAST').

To generate more complex structures, such as sugars, the residue names should be supplied (hyphen delimited) and the -u/--use-charmm-names option selected. Some example usage is given below.

--name and --segid control the names given in the COMPND record and the segment id respectively.

Examples

To produce a simple peptide sequence, the one letter code can be used. To produce the peptide HIS-GLU-TYR, creating HEY.psf and HEY.pdb:

pygen-structures HEY -o HEY

Supposing we think histidine should be protonated, we can change the protonation state of histidine by specifying a different histidine form:

pygen-structures HEY -o HEY --histidine HSP

Or we could use the three letter codes:

pygen-structures -u HSP-GLU-TYR -o HEY

For the trisaccharide raffinose, we must use the residue codes. The default segid is PROT, so we can specify a more specific segid using the --segid flag, and set the name given in the COMPND header using --name. The following command produces RAFF.psf and RAFF.pdb:

pygen-structures -u AGLC-BFRU-AGAL --patches RAFF 0 1 2 --segid RAFF --name Raffinose -o RAFF

We can also make glycopeptides. To link alpha-glucose to an arginine residue (in this case, from an ALA-ASN-ALA peptide), we can use the NGLA patch. Note that because the protein residue is not the last in the chain, we have to apply the C-terminus patch manually.

pygen-structures -u ALA-ASN-ALA-AGLC --patches CTER -2 NGLA 1 -1 -o ANA-NAGLC

By default, if parameters are missing then the files are not created and the missing parameters are written to stdout. Using the -v flag will disable verification.

$ pygen-structures AdP -o AdP
Missing parameters:
bonds {('CPD1', 'CC')}
$ pygen-structures -v AdP -o AdP
$

A different CHARMM distribution can be loaded using the -t option, with the path to the folder. pygen-structures ships with the latest CHARMM distribution (July 2019) at the time of writing, with some modifications to correct the D-amino acid parameters (these modifications are highlighted in the toppar README). The function which parses the folder will pick the latest versions of the parameter and topology files (36 over 27, 36m over 36), so if you plan on using an older version of the forcefield (this is not recommended) you will have to remove the newer versions and change the extensions to match the current conventions (.rtf, .prm).

Library usage

Information about classes and functions for usage as a library can be found on the project's ReadtheDocs page.

Contributing

Contributions are welcome! Please read our Code of Conduct.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

pygen-structures-0.2.5.tar.gz (7.8 MB view details)

Uploaded Source

Built Distribution

pygen_structures-0.2.5-py3-none-any.whl (8.1 MB view details)

Uploaded Python 3

File details

Details for the file pygen-structures-0.2.5.tar.gz.

File metadata

  • Download URL: pygen-structures-0.2.5.tar.gz
  • Upload date:
  • Size: 7.8 MB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/3.1.1 pkginfo/1.5.0.1 requests/2.22.0 setuptools/45.2.0 requests-toolbelt/0.9.1 tqdm/4.42.0 CPython/3.7.5

File hashes

Hashes for pygen-structures-0.2.5.tar.gz
Algorithm Hash digest
SHA256 6cd6c0e105ffd1862d88e28a3005997bd934b13cbe2941ead75641085e6c8bdb
MD5 84bc477bd63ac4f021343efe16038682
BLAKE2b-256 aaf85cf79d9255dcf2d7d14ac8444a752ef67cc54a50792908436eff1d49d5d5

See more details on using hashes here.

File details

Details for the file pygen_structures-0.2.5-py3-none-any.whl.

File metadata

  • Download URL: pygen_structures-0.2.5-py3-none-any.whl
  • Upload date:
  • Size: 8.1 MB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/3.1.1 pkginfo/1.5.0.1 requests/2.22.0 setuptools/45.2.0 requests-toolbelt/0.9.1 tqdm/4.42.0 CPython/3.7.5

File hashes

Hashes for pygen_structures-0.2.5-py3-none-any.whl
Algorithm Hash digest
SHA256 8421eb77e4a7106dc17b89649e78010b20870c2c77c2974c269ac8dbdd793d73
MD5 b3f06baebb68019f8fb23734ce51c341
BLAKE2b-256 5d345defda0bf3a928c600ecfdd63a83cb30784a2055f59475077225a3ca4504

See more details on using hashes here.

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page