Navigation

Latest Release

v1.4.2 is a bugfix release. v1.4 series come with Python 3 support and new sequence analysis features. See changes or go to downloads.

How to Cite ProDy

Bakan A, Meireles LM, Bahar I ProDy: Protein Dynamics Inferred from Theory and Experiments
Bioinformatics 2011 27(11):1575-1577. See all.

Table Of Contents

Previous topic

Bibliography

Next topic

License

Changes

Release 1.4.3 (in development)

Changes:

Release 1.4.2 (April 19, 2013)

Improvements:

Changes:

  • MSE is included in the definition of non-standard amino acids, i.e. nonstdaa.

Bugfixes:

  • Atom selection problems related to using all and none in composite selections, e.g. 'calpha and all', is fixed by defining these keywords as Atom flags.
  • Fasta files with sequence labels using multiple pipe characters would cause C parser (and so parseMSA()) to fail. This issue is fixed by completely disregarding pipe characters.
  • Empty chain identifiers for PDB hits would cause a problem in parsing XML results file and blastPDB() would throw an exception. This case is handled by slicing the chain identifier string.
  • A problem in viewNMDinVMD() related to module imports is fixed.
  • A problem with handling weights in loadEnsemble() is fixed.

Release 1.4.1 (Dec 16, 2012)

New Features:

Improvements:

  • refineMSA() function and evol refine application can perform MSA refinements by removing similar sequences.
  • writePDB() function takes beta and occupancy arguments to be outputted in corresponding columns.
  • MSA indexing and slicing are revised and improved.
  • parseMSA() is improved to handle indexing of sequences that have the same label in an MSA file, e.g. domains repeated in a protein.
  • prody anm, prody gnm, and prody pca applications can write heatmap files for visualization using NMWiz and Heatmapper plugins.
  • Several improvements made to handling sequence labels in Pfam MSA files. Files that contain sequence parts with same protein UniProt ID are handled delicately.

Changes:

  • ProDy will not emit a warning message when a wwPDB server is not set using wwPDBServer(), and use the default US server.
  • Indexing MSA returns Sequence instances.
  • Iterating over MSA and MSAFile yields Sequence instances.

Bugfixes:

  • Fixed a syntax problem that prevented running ProDy using Python 2.6.
  • Fixed NMA indexing problem that was introduced in v1.4.

NMWiz

  • NMWiz can visualize heatmaps linked to structural view via Heatmapper. Clicking on the heatmap will highlight atom or residue pairs.
  • ProDy interface has the option to write and load cross-correlations.
  • NMWiz can determined whether a model is an extended model. For extended models plotting mobility has been improved. Only a single value per residue will be plotted, and clicking on the plot will highlight all of the residue atoms.

Release 1.4 (Dec 2, 2012)

New Features:

Python 3 Support

  • ProDy has been refactored to support Python 3. Windows installers for Python 2.6, 2.7, 3.1, and 3.2 are available in Getting ProDy.
  • Unit tests are compatible with Python 2.7 and 3.2, and running them with other versions gives errors due to unavailability of some unittest features.

Sequence Analysis

Bugfix:

Release 1.3.1 (Nov 6, 2012)

New Features:

Improvements:

  • Improved gunzip() function to handle .gz extensions and string buffers.

Changes:

Bugfix:

  • Fixed a bug in HierView that would cause wrong assignment of residue/chain indices to atoms when residue or chain atoms are separated by atoms of other entities. This would also caused problems when making keyword selections, such as protein.
  • Added dummy atom check in Ensemble.setAtoms() and Trajectory.setAtoms() methods to avoid indexing problems.

Release 1.3 (Sep 30, 2012)

Improvements:

  • select module and its documentation are completely rewritten. Select class uses simplest possible parser to evaluate selection strings and achieves more than 25% speed-up on average.
  • Atom Selections become more forgiving of small typos, but will issue warning messages when they are detected via SelectionWarning. These messages can be turned of using confProDy()
  • Functions used in ProDy Applications have been refactored to allow for using them directly. See apps for their documentation.

Bugfix:

Release 1.2.1 (Sep 6, 2012)

If you are upgrading from ProDy v1.1, see also the below changes introduced in v1.2.

Bugfix:

  • A problem in select module regarding Numpy numeric types is fixed. Problem would emerge on platforms which do not offer some numeric types, e.f. np.float16.
  • Fixed problems in prody anm, prody gnm, and prody fetch related to writing output files.

Changes:

  • The way that prody fetch command handles files containing PDB identifiers has changed.

Release 1.2 (Aug 30, 2012)

Important Changes:

Package folder prody is moved into lib folder to prevent exceptions related to importing compiled packages from the installation folder.

Some changes in Trajectory and Ensemble methods related to linking, setting, and selecting atoms were made to make the interface more intuitive. These changes, which may break your code, are as follows:

For usage examples see Trajectory analysis, Trajectory analysis II, Atoms and Frames, and Trajectory output.

New Features:

  • Atom flags, that are used in Atom Selections, is implemented. See its documentation for handy usage examples.
  • sortAtoms() function is implemented.
  • pickCentralConf() function is implemented to pick the conformation or the active coordinate set that is closest to the average of coordinate sets.
  • writePSF(), a simple PSF file writer, is implemented.
  • glob() utility function is implemented.
  • iterPDBFilenames() function is implemented, which can be used to iterate over all PDB files stored in a local mirror of Protein Data Bank.
  • findPDBFiles() function is implemented, which can be used to access PDB files in a path.

Improvements:

Changes:

  • showOverlapTable() displays first set of modes along x axis of the plot.
  • AtomGroup.setData() does not accept arrays with boolean data type, use AtomGroup.setFlags() instead.
  • writePDB() function argument model is changed to csets that indicates the coordinate set index of atoms argument.
  • PackageLogger.timing() does not return elapsed time, only logs this information.
  • PackageLogger.startLogfile() is deprecated for removal in v1.3, use PackageLogger.start() instead.
  • PackageLogger.closeLogfile() is deprecated for removal in v1.3, use PackageLogger.close() instead.
  • from prody.utilities import * will not work anymore due to potential name conflicts with Python standard library functions. Import required functions explicitly.
  • writePDB() appends .pdb extension to filename when it is not present
  • prody select command positional argument order is changed to allow for handling multiple PDBs at a time. Old older will be supported until v1.4, but a warning message will be issued.
  • select argument in alignCoordsets() is removed, make selection outside of the function instead.

Deprecations:

  • AtomGroup.getHeteros() method has been deprecated for removal in v1.3, use getFlags('hetatm') instead.
  • AtomMap.getMappedFlags() and AtomMap.getDummyFlags() methods have been deprecated for removal in v1.3, use getFlags('mapped') and getFlags('dummy') instead.
  • getVerbosity() and setVerbosity() are deprecated for removal in v1.3, use confProDy() instead which save changes permanently.
  • NMA.getModes() and ModeSet.getModes() methods are deprecated for removal in v1.3, use list(), e.g. list(model), instead.

Bugfixes:

  • Fixed a bug in prody contacts command that arose problems when when selecting a subset of the target atoms.

NMWiz

Improvements:

  • ProDy Interface shows the size of the trajectory output file for PCA calculations.
  • Mode Graphics Options allows for copying arrows settings from one mode to another.
  • Color scale method and midpoint for protein coloring based on mobility and bfactors can be adjusted from Protein Graphics Options panel.

Release 1.1 (June 1, 2012)

New Features:

Improvements:

  • moveAtoms() function is improved to move atoms to a specified location.
  • DCDFile and parseDCD() take astype keyword argument for automatic type recasting for coordinate arrays. This option can be used to convert 32-bit coordinate arrays to 64-bit automatically for higher precision calculations.
  • Commands prody anm, prody gnm, and prody pca can extend a coarse grained model to backbone or all atoms of the residues. See their documentation pages.

Changes:

  • Color scale used by showOverlapTable() is normalized by default.
  • tools module is depracated for removal, use utilities instead.
  • array argument in moveAtoms() is replaced with by keyword argument.
  • which argument in AtomGroup.copy() method is deprecated for removal in version 1.2.
  • DCDFile does not log information for most common type of DCD file, i.e. 32-bit CHARMM format.
  • Trajectory.getNextIndex() method is deprecated for removal in v1.2, use nextIndex() instead.

Bugfixes:

  • Fixed several problems in iterNeighbors() function and Contacts class that were introduced after transition to new KDTree interface.
  • Fixed a problem in setting selection strings of fragments identified using findFragments().
  • Fixed a problem in calcCenter() related to weighted center calculation.
  • Fixed a problem of in copying AtomMap instances, which would emerge when bond information was present in unusual mappings, such as when atom orders are changed or an atom is present multiple times in the mapping.

NMWiz

Improvements:

  • Mode scaling options are improved.
  • Options added for extending coarse grained NMA models to residue backbone or all atoms.

Release 1.0.4 (May 2, 2012)

Bugfixes:

  • Fixed a problem in calcPhi() function that raised a name error.
  • Fixed a problem in KDTree.getDistances() method that raised a name error when unitcell is provided.
  • Fixed a problem in buildDistMatrix() and calcDistance() functions causing miscalculations when unitcell is given.
  • Revised KDTree methods dealing with to handle special cases where unitcell might have some dimensions zero.

Changes:

  • buildKDTree() method is removed, earlier than planned due to unexpected bugfix releases.

Release 1.0.3 (May 1, 2012)

Bugfixes:

  • Fixed kdtree import problem.

New Features:

  • buildDistMatrix() function that can take periodic boundary conditions is implemented.

Improvements:

  • calcDistance() function is improved to take periodic boundary conditions into account when provided by the users.

Release 1.0.2 (May 1, 2012)

New Features:

  • Methods to deal with connected subsets of atoms are implemented, see AtomGroup.iterFragments() and AtomGroup.numFragments().
  • pickCentral() method is implemented for picking the atom that is closest to the centroid of a group or subset of atoms.
  • ProDy configuration option auto_secondary is implemented to allow for parsing and assigning secondary structure information from PDB file header data automatically. See assignSecstr() and confProDy() for usage details.
  • prody align makes use of --select when aligning multiple structures. See usage examples: prody align
  • printRMSD() function that prints minimum, maximum, and mean RMSD values when comparing multiple coordinate sets is implemented.
  • findFragments() function that identifies fragments in atom subsets, e.g. Selection, is implemented.
  • A new KDTree interface with coherent method names and capability to handle periodic boundary conditions is implemented.

Improvements:

Changes:

  • There is a change in the behavior of addition operation on instances of AtomGroup. When operands do not have same number of coordinate sets, the result will have one coordinate set that is concatenation of the active coordinate sets of operands.
  • buildKDTree() function is deprecated for removal, use the new KDTree class instead.

Bugfixes:

  • A problem in building hierarchical views when making selections using resindex, chindex, and segindex keywords is fixed.
  • A problem in Chain and Residue selection strings that would emerge when a HierView is build using a selection is fixed.
  • A problem with copying AtomGroup instances whose coordinates are not set is fixed.
  • AtomGroup fragment detection algorithm is rewritten to avoid the problem of reaching maximum recursion depth for large molecules with the old recursive algorithm.
  • A problem with picking central atom of AtomGroup instances in pickCentral() function is fixed.
  • A problem in Select class that caused exceptions when evaluating complex macro definitions is fixed.
  • Fixed a problem in handling multiple trajectory files. The problem would emerge when a file was added (addFile()) to a Trajectory after atoms were set (setAtoms()). Newly added file would not be associated with the atoms and coordinates parsed from this file would not be set for the AtomGroup instance.

Release 1.0.1 (Apr 6, 2012)

New Features:

  • ProDy can be configured to automatically check for updates on a regular basis, see checkUpdates() and confProDy() functions for details.
  • alignPDBEnsemble() function is implemented to align PDB files using transformations calculated in ensemble analysis. See usage example in PCA from a Blast search example.
  • PDBConformation.getTransformation() is implemented to return the transformation that was used to superpose conformation onto reference coordinates. This transformation can be used to superpose the original PDB file onto the reference PDB file.
  • Amino acid sequences with regular expressions can be used to make atom selections, e.g. 'sequence "C..C"'. See Atom Selections for usage details.
  • calcCrossProjection() function is implemented.

Improvements:

  • Select class raises a SelectionError when potential typos are detected in a selection string, e.g. 'chain AB' is a grammatically correct selection string that will return None since no atoms have chain identifier 'AB'. In such cases, an exception noting that values exceed maximum number of characters is raised.
  • prody align command accepts percent sequence identity and overlap parameters used when matching chains from given multiple structures.
  • When using prody align command to align multiple structure, all models in NMR structures are aligned onto the reference structure.
  • prody catdcd command accepts --align SELSTR argument that can be used to align frames when concatenating files.
  • showProjection() and showCrossProjection() functions are improved to evaluate list of markers, color, labels, and texts. See usage example in PCA of X-ray structures: Plotting.
  • Trajectory instances can be used for calculating and plotting projections using calcProjection(), showProjection(), calcCrossProjection(), and showCrossProjection() functions.

Changes:

  • Phosphorylated amino acids, phosphothreonine (TPO), O-phosphotyrosine (PTR), and phosphoserine (SEP), are recognized as acidic protein residues. This prevents having breaks in protein chains which contains phosphorylated residues. See Atom Selections for definitions of protein and acidic keywords.
  • Hit dictionaries from PDBBlastRecord will use percent_overlap instead of percent_coverage. Older key will be removed in v1.1.
  • Transformation.get4x4Matrix() method is deprecated for removal in v1.1, use Transformation.getMatrix() method instead.

Bugfixes:

  • A bug in some ProDy Applications is fixed. The bug would emerge when invalid arguments were passed to effected commands and throw an unrelated exception hiding the error message related to the arguments.
  • A bug in 'bonded to ...' is fixed that emerged when '...' selected nothing.
  • A bug in 'not' selections using . operator is fixed.

Release 1.0 (Mar 7, 2012)

Improvements:

  • ANM.buildHessian() method is not using a KDTree by default, since with some code optimization the version not using KDTree is running faster. Same optimization has gone into GNM.buildKirchhoff() too, but for Kirchoff matrix, version using KDTree is faster and is the default. Both methods have kdtree argument to choose whether to use it or not.
  • prody script is updated. Importing Prody and Numpy libraries are avoided. Script responses to help queries faster. See ProDy Applications for script usage details.
  • Added bonded to ... selection method that expands a selection to immediately bound atoms. See Atom Selections for its description.
  • fetchPDBLigand() parses bond data from the XML file.
  • fetchPDBLigand() can optionally save compressed XML files into ProDy package folder so that frequent access to same files will be more rapid. See confProDy() function for setting this option.
  • Select class is revised. All exceptions are handled delicately to increase the stability of the class.
  • Distance based atom selection is 10 to 15% faster for atom groups with more than 5K atoms.
  • Added uncompressed file saving option to prody blast command.

Changes:

Bugfixes:

  • A problem in distance based atom selection which would could cause problems when a distance based selection is made from a selection is fixed.
  • Changed prody blast so that when a path for downloading files are given files are not save to local PDB folder.

Release 0.9.4 (Feb 4, 2012)

Changes:

  • setAtomGroup() and getAtomGroup() methods are renamed as Ensemble.setAtoms() and Ensemble.getAtoms().
  • AtomGroup class trajectory methods, i.e. AtomGroup.setTrajectory(), AtomGroup.getTrajectory(), AtomGroup.nextFrame(), AtomGroup.nextFrame(), and AtomGroup.gotoFrame() methods are deprecated. Version 1.0 will feature a better integration of AtomGroup and Trajectory classes.

Bugfixes:

Release 0.9.3 (Feb 1, 2012)

New Features:

  • DBRef class is implemented for storing references to sequence databases parsed from PDB header records.

  • Methods for storing coordinate set labels in AtomGroup instances are implemented: getACSLabel(), and getACSLabel().

  • calcCenter() and moveAtoms() functions are implemented for dealing with coordinate translation.

  • Hierarchical view, HierView, is completely redesigned. PDB files that contain non-empty segment name column (or when such information is parsed from a PSF file), new design delicately handles this information to identify distinct chains and residues. This prevents merging distinct chains in different segments but with same identifiers and residues in those with same numbers. New design is also using ordered dictionaries collections.OrderedDict and lists so that chain and residue iterations yield them in the order they are parsed from file. These improvements also bring modest improvements in speed.

  • Segment class is implemented for handling segments of atoms defined in molecular dynamics simulations setup, using psfgen for example.

  • Context manager methods are added to trajectory classes. A trajectory file can be opened as follows:

    with Trajectory('mdm2.dcd') as traj:
    for frame in traj:
        calcGyradius(frame)

  • Chain slicing is implemented:

    p38 = parsePDB('1p38')
chA = p38['A']
res_4to10 = chA[4:11]
res_100toLAST = chA[100:]

  • Some support for bonds is implemented to AtomGroup class. Bonds can be set using setBonds() method. All bonds must be set at once. iterBonds() or iterBonds() methods can be used to iterate over bonds in an AtomGroup or an Atom.

  • parsePSF() parses bond information and sets to the atom group.

  • Selection.update() method is implemented, which may be useful to update a distance based selection after coordinate changes.

  • buildKDTree() and iterNeighbors() methods are implemented for facilitating identification of pairs of atoms that are proximal.

  • iterAtoms() method is implemented to all atomic classes to provide uniformity for atom iterations.

  • calcAngle(), calcDihedral(), calcPhi(), calcPsi(), and calcOmega() methods are implemented.

Improvements:

Changes:

  • Residue methods getNumber(), setNumber(), getName(), setName() methods are deprecated and will be removed in v1.0.
  • Chain methods getIdentifier() and setIdentifier() methods are deprecated and will be removed in v1.0.
  • Polymer attribute identifier is renamed as chid.
  • Chemical attribute identifier is renamed as resname.
  • getACSI() and setACSI() are renamed as getACSIndex() and setACSIndex(), respectively.
  • calcRadiusOfGyration() is deprecated and will be removed in v1.0. Use calcGyradius() instead.

Bugfixes:

  • Fixed a problem in parsePDB() that caused loosing existing coordinate sets in an AtomGroup when passed as ag argument.
  • Fixed a problem with "same ... as ..." argument of Select that selected atoms when followed by an incorrect atom selection.
  • Fixed another problem with "same ... as ..." which result in selecting multiple chains when same chain identifier is found in multiple segments or multiple residues when same residue number is found in multiple segments.
  • Improved handling of negative integers in indexing AtomGroup instances.

Release 0.9.2 (Jan 11, 2012)

New Features:

  • prody catdcd command is implemented for concatenating and/or slicing .dcd files. See prody catdcd for usage examples.
  • DCDFile can be opened in write or append mode, and coordinate sets can be added using write() method.
  • getReservedWords() can be used to get a list of words that cannot be used to label user data.
  • confProDy() function is added for configuring ProDy.
  • ProDy can optionally backup existing files with .BAK (or another) extension instead of overwriting them. This behavior can be activated using confProDy() function.

Improvements:

  • writeDCD() file accepts AtomGroup or other Atomic instances as trajectory argument.
  • prody align command can be used to align multiple PDB structures.
  • prody pca command allows atom selections for DCD files that are accompanied with a PDB or PSF file.

Changes:

  • DCDFile instances, when closed, raise exception, similar to behavior of file objects in Python.
  • Title of AtomGroup instances resulting from copying an Atomic instances does not start with ‘Copy of’.
  • changeVerbosity() and getVerbosityLevel() are renamed as setVerbosity() and getVerbosity(), respectively. Old names will be removed in v1.0.
  • ProDy applications (commands) module is rewritten to use new argparse module. See ProDy Applications for details of changes.
  • argparse module is added to the package for Python versions 2.6 and older.

Bugfixes:

Release 0.9.1 (Nov 9, 2011)

Bug Fixes:

  • Fixed problems with reading and writing configuration files.
  • Fixed problem with importing nose for testing.

Release 0.9 (Nov 8, 2011)

New Features:

  • PDBML and mmCIF files can be retrieved using fetchPDB() function.
  • getPDBLocalFolder() and setPDBLocalFolder() functions are implemented for local PDB folder management.
  • parsePDBHeader() is implemented for convenient parsing of header data from .pdb files.
  • showProtein() is implemented to allow taking a quick look at protein structure.
  • Chemical and Polymer classes are implemented for storing chemical and polymer component data parsed from PDB header records.

Changes:

Warning

This release introduces numerous changes in method and function names all aiming to improve the interactive usage experience. All changes are listed below. Currently these functions and methods are present in both old and new names, so code using ProDy must not be affected. Old function names will be removed from version 1.0, which is expected to happen late in the first quarter of 2012.

Old function names are marked as deprecated, but ProDy will not issue any warnings until the end of 2011. In 2012, ProDy will automatically start issuing DeprecationWarning upon calls using old names to remind the user of the name change.

For deprecated methods that are present in multiple classes, only the affected modules are listed for brevity.

Note

When modifying code using ProDy to adjust the name changes, turning on deprecation warnings may help locating all use cases of the deprecated names. See turnonDepracationWarnings() for this purpose.

Functions:

The following function name changes are mainly to reduce the length of the name in order to make them more suitable for interactive sessions:

Old name New name
applyBiomolecularTransformations() buildBiomolecules()
assignSecondaryStructure() assignSecstr()
scanPerturbationResponse() calcPerturbResponse()
calcCrossCorrelations() calcCrossCorr()
calcCumulativeOverlap() calcCumulOverlap()
calcCovarianceOverlap() calcCovOverlap()
showFractOfVariances() showFractVars()
showCumFractOfVariances() showCumulFractVars()
showCrossCorrelations() showCrossCorr()
showCumulativeOverlap() showCumulOverlap()
deform() deformAtoms()
calcSumOfWeights() calcOccupancies()
showSumOfWeights() showOccupancies()
trimEnsemble() trimPDBEnsemble()
getKeywordResidueNames() getKeywordResnames()
setKeywordResidueNames() setKeywordResnames()
getPairwiseAlignmentMethod() getAlignmentMethod()
setPairwiseAlignmentMethod() setAlignmentMethod()
getPairwiseMatchScore() getMatchScore()
setPairwiseMatchScore() setMatchScore()
getPairwiseMismatchScore() getMismatchScore()
setPairwiseMismatchScore() setMismatchScore()
getPairwiseGapOpeningPenalty() getGapPenalty()
setPairwiseGapOpeningPenalty() setGapPenalty()
getPairwiseGapExtensionPenalty() getGapExtPenalty()
setPairwiseGapExtensionPenalty() setGapExtPenalty()

Coordinate methods:

All getCoordinates() and setCoordinates() methods in atomic and ensemble classes are renamed as getCoords() and setCoords(), respectively.

getNumOf methods:

All method names starting with getNumOf now start with num. This change brings two advantages: method names (i) are considerably shorter, and (ii) do not suggest that there might also be corresponding set methods.

Old name New name Affected modules
getNumOfAtoms() numAtoms() atomic, ensemble, dynamics
getNumOfChains() numChains() atomic
getNumOfConfs() numConfs() ensemble
getNumOfCoordsets() numCoordsets() atomic, ensemble
getNumOfDegOfFreedom() numDOF() dynamics
getNumOfFixed() numFixed() ensemble
getNumOfFrames() numFrames() ensemble
getNumOfResidues() numResidues() atomic
getNumOfMapped() numMapped() atomic
getNumOfModes() numModes() dynamics
getNumOfSelected() numSelected() ensemble
getNumOfUnmapped() numUnmapped() atomic

getName method:

getName() methods are renamed as getTitle() to avoid confusions that might arise from changes in atomic method names listed below. All classes in atomic, ensemble, and dynamics are affected from this change.

In line with this change, parsePDB() and parsePQR() name arguments are changed to title, but name argument will also work until release 1.0.

This name change conflicted with DCDFile.getTitle() method. The conflict is resolved in favor of the general getTitle() method. An alternative method will be implemented to handle title strings in DCD files.

get/set methods of atomic classes:

Names of get and set methods allowing access to atomic data are all shortened as follows:

Old name New name
getAtomNames() getNames()
getAtomTypes() getTypes()
getAltLocIndicators() getAltlocs()
getAnisoTempFactors() getAnisos()
getAnisoStdDevs() getAnistds()
getChainIdentifiers() getChains()
getElementSymbols() getElements()
getHeteroFlags() getHeteros()
getInsertionCodes() getIcodes()
getResidueNames() getResnames()
getResidueNumbers() getResnums()
getSecondaryStrs() getSecstrs()
getSegmentNames() getSegnames()
getSerialNumbers() getSerials()
getTempFactors() getBetas()

This change affects all atomic classes, AtomGroup, Atom, Chain, Residue, Selection and AtomMap.

Other changes in atomic methods:

  • getSelectionString() renamed as getSelstr()

Methods handling user data (which was previously called attribute) are renamed as follows:

Old name New name
getAttribute() getData()
getAttrNames() getDataLabels()
getAttrType() getDataType()
delAttribute() delData()
isAttribute() isData()
setAttribute() setData()

To be removed:

Finally, the following methods will be removed, but other suitable methods are overloaded to perform their action:

  • removed AtomGroup.getBySerialRange(), overloaded AtomGroup.getBySerial()
  • removed getProteinResidueNames(), overloaded getKeywordResnames()
  • removed setProteinResidueNames(), overloaded setKeywordResnames()

Scripts:

The way ProDy scripts work has changed. See ProDy Applications for details. Using older scripts will start issuing deprecation warnings in 2012.

Bug Fixes:

  • Bugs in execDSSP() and execSTRIDE() functions that caused exceptions when compressed files were passed is fixed.
  • A problem in scripts for PCA of DCD files is fixed.

NMWiz

Development of NMWiz is finalized and it will not be distributed in the ProDy installation package anymore. See NMWiz pages for instructions on installing it.

Release 0.8.3 (Oct 16, 2011)

New Features:

  • Functions to read and write PQR files: parsePQR() and writePQR().
  • Added PDBEnsemble.getIdentifiers() method that returns identifiers of all conformations in the ensemble.
  • ProDy tests are incorporated to the package installer. If you are using Python version 2.7, you can run the tests by calling prody.test().

Improvements:

  • blastPDB() function and PDBBlastRecord class are rewritten to use faster and more compact code.
  • New PackageLogger function is implemented to unify logging and reporting task progression.
  • Improvements in PDB ensemble support functions, e.g. trimPDBEnsemble(), are made.
  • Improvements in ensemble concatenations are made.

Bug Fixes:

  • Bugfixes in PDBEnsemble() slicing operation. This may have affected users when slicing a PDB ensemble for plotting projections in color for different forms of the protein.

Release 0.8.2 (Oct 14, 2011)

New Features:

  • fetchPDBClusters(), loadPDBClusters(), and getPDBCluster() functions are implemented for handling PDB sequence cluster data. These functions can be used instead of blastPDB() function for fast access to structures of the same protein (at 95% sequence identity level) or similar proteins.
  • Perturbation response scanning method described in [CA09] is implemented as scanPerturbationResponse() based on the code provided by Ying Liu.

Changes:

  • fetchPDBLigand() returns the URL of the XML file in the ligand data dictionary.
  • Name of the ProDy configuration file in user home directory is renamed as .prodyrc (used to be .prody).
  • applyBiomolecularTransformations() and assignSecondaryStructure() functions raise ValueError when the function fails to perform its action due to missing data in header dictionary.
  • fetchPDB() decompresses PDB files found in the working directory when user asks for decompressed files.
  • parsePDB() appends chain and subset arguments to AtomGroup() name.
  • chain argument is added to PDBBlastRecord.getHits().

Improvements:

  • Atom selection class Select is completely redesigned to prevent breaking of the parser when evaluating invalid selection strings.
  • Improved type checking in parsePDB() function.

Bug Fixes:

  • Bugfixes in parseDSSP(): one emerged problems in lines indicating chain breaks, another did not parse bridge-partners correctly. Both fixes are contributed by Kian Ho.
  • Bugfix in parsePDB() function. When only header is desired (header=True, model=0), would return a tuple containing an empty atom group and the header.

Developmental:

  • Unit tests for proteins and select modules are developed.

Release 0.8.1 (Sep 16, 2011)

New Features:

  • fetchLigandData() is implemented for fetching ligand data from Ligand Expo.
  • parsePSF() function is implemented for parsing X-PLOR format PSF files.

Changes:

  • __slots__ is used in AtomGroup and Atomic classes. This change prevents user from assigning new variables to instances of all classes derived from the base Atomic.
  • pyparsing is updated to version 1.5.6.

Bug Fixes:

  • A bug in AtomGroup.copy() method is fixed. When AtomGroup instance itself is copied, deep copies of data arrays were not made.
  • A bug in Select class raising exceptions when negative residue number values are present is fixed.
  • Another bug in Select class misinterpreting same residue as ... statement when specific chains are involved is fixed.
  • A bug in AtomGroup.addCoordset() method duplicating coordinates when no coordinate sets are present in the instance is fixed.

NMWiz

Changes:

  • Version number in main window is iterated.
  • Mode graphics material is stored for individual modes.
  • Mode scaling factor is printed when active mode or RMSD is changed.
  • All selections are deleted to avoid memory leaks.

Release 0.8 (Aug 24, 2011)

Note

After installing v0.8, you may need to make a small change in your existing scripts. If you are using Ensemble class for analyzing PDB structures, rename it as PDBEnsemble. See the other changes that may affect your work below and the class documentation for more information.

New Features:

  • DCDFile is implemented for handling DCD files.
  • Trajectory is implemented for handling multiple trajectory files.
  • writeDCD() is implemented for writing DCD files.
  • Trajectory analysis example to illustrate usage of new classes for handling DCD files. EDA of MD Trajectories example is updated to use new ProDy classes.
  • PCA supports Trajectory and DCDFile instances.
  • Ensemble and PDBEnsemble classes can be associated with AtomGroup instances. This allows selecting and evaluating coordinates of subset of atoms. See setAtomGroup(), select(), getAtomGroup(), and getSelection() methods.
  • execDSSP(), parseDSSP(), and performDSSP() functions are implemented for executing and parsing DSSP calculations.
  • execSTRIDE(), parseSTRIDE(), and performSTRIDE() functions are implemented for executing and parsing DSSP calculations.
  • parsePDB() function parses atom serial numbers. Atoms can be retrieved from an AtomGroup instance by their serial numbers using getBySerial() and getBySerialRange() methods.
  • calcADPs() function can be used to calculate anisotropic displacement parameters for atoms with anisotropic temperature factor data.
  • getRMSFs() is implemented for calculating root mean square fluctuations.
  • AtomGroup and Mode or Vector additions are supported. This adds a new coordinate set to the AtomGroup instance.
  • getAttrNames() is implemented for listing user set attribute names.

Improvements:

Changes:

  • Ensemble does not store conformation names. Instead, newly implemented PDBEnsemble class stores identifiers for individual conformations (PDB IDs). This class should be used in cases where source of individual conformations is important.
  • calcProjection(), showProjection(), and showCrossProjection() function calculate/display root mean square deviations, by default.
  • Oxidized cysteine residue abbreviation CSO is added to the definition of protein keyword.
  • getMSF() method is renamed as getMSFs().
  • parseDCD() function returns Ensemble instances.

Bug Fixes:

  • A bug in select module causing exceptions when regular expressions are used is fixed.
  • Another bug in select module raising exception when “(not ..,” is passed is fixed.
  • Various bugfixes in ensemble module.
  • Problem in prody fetch that occurred when a file is found in a local mirror is fixed.
  • Bugfix in AtomPointer.copy() method.

NMWiz

New Features:

  • NMWiz can be used to compare two structures by calculating and depicting structural changes.
  • Arrow graphics is scaled based on a user specified RMSD value.

Improvements:

  • NMWiz writes DCD format trajectories for PCA using ProDy. This provides significant speed up in cases where IO rate is the bottleneck.

Changes:

  • Help is provided in a text window to provide a cleaner GUI.

Release 0.7.2 (Jun 21, 2011)

New Features:

  • parseDCD() is implemented for parsing coordinate sets from DCD files.

Improvements:

  • parsePDB() parses SEQRES records in header sections.

Changes:

  • Major classes can be instantiated without passing a name argument.
  • Default selection in NMWiz ProDy interface is changed to ensure selection only protein Cα atoms.

Bug Fixes:

  • A bug in writeNMD() function causing problems when writing a single mode is fixeed.
  • Other bugfixes in dynamics module functions.

Release 0.7.1 (Apr 28, 2011)

Highlights:

  • Atomic __getattribute__() is overloaded to interpret atomic selections following the dot operator. For example, atoms.calpha is interpreted as atoms.select('calpha'). See :ref:`` for more details.
  • AtomGroup class is integrated with HierView class. Atom group instances now can be indexed to get chains or residues and number of chains/residues can be retrieved. A hierarchical view is generated and updated when needed. See :ref:`` for more details.

New Features:

  • matchAlign() is implemented for quick alignment of protein structures. See Extract ligands usage example.
  • setAttribute(), getAttribute(), delAttribute(), and isAttribute() functions are implemented for AtomGroup class to facilitate storing user provided atomic data. See Storing data in AtomGroup example.
  • saveAtoms() and loadAtoms() functions are implemented to allow for saving atomic data and loading it This saves custom atomic attributes and much faster than parsing data from PDB files.
  • calcCollectivity() function is implemented to allow for calculating collectivity of deformation vectors.

Improvements:

  • parsePDB() can optionally return biomolecule when biomol=True keyword argument is passed.
  • parsePDB() can optionally make secondary structure assignments when secondary=True keyword argument is passed.
  • calcSqFlucts() function is changed to accept Vector instances, e.g. deformation vectors.

Changes:

  • Changes were made in calcADPAxes() function to follow the conventions in analysis ADPs. See its documentation.

Bug Fixes:

Release 0.7 (Apr 4, 2011)

New Features:

  • Regular expressions can be used in atom selections. See select module for details.
  • User can define selection macros using defSelectionMacro() function. Macros are saved in ProDy configuration and loaded in later sessions. See select module for other related functions.
  • parseSparseMatrix() function is implemented for parsing matrices in sparse format. See the usage example in Using an External Matrices.
  • deform() function is implemented for deforming coordinate sets along a normal mode or linear combination of multiple modes.
  • sliceModel() function is implemented for slicing normal mode data to be used with functions calculating atomic properties using normal modes.

Improvements:

  • Atom selections using bare keyword arguments is optimized. New keyword definitions are added. See select module for the complete list.
  • A new keyword argument for calcADPAxes() allows for comparing largest axis to the second largest one.

Changes:

  • There are changes in function used to alter definitions of selection keywords. See select for details.
  • assignSecondaryStructure() function assigns SS identifiers to all atoms in a residue. Residues with no SS information specified is assigned coil conformation.
  • When Ensemble and NMA classes are instantiated with an empty string, instances are called “Unnamed”.
  • sliceMode(), sliceVector() and reduceModel() functions return the atom selection in addition to the sliced vector/mode/model instance.

Bug Fixes:

  • Default selection for calcGNM() function is set to “calpha”.

NMWiz

New Features:

  • NMWiz supports GNM data and can use ProDy for GNM calculations.
  • NMWiz can gather normal mode data from molecules loaded into VMD. This allows NMWiz to support all formats supported by VMD.
  • User can write data loaded into NMWiz in NMD format.
  • An Arrow Graphics option allows the user to draw arrows in both directions.
  • User can select Licorice representation for the protein if model is an all atom mode.
  • User can select Custom as the representation of the protein to prevent NMWiz from chancing a user set representation.
  • Trace is added as a protein backbone representation option.

Improvements:

  • NMWiz remembers all adjustments on arrow graphics for all modes.
  • Plotting Clear button clears only atom labels that are associated with the dataset.
  • Removing a dataset removes all associated molecule objects.
  • Selected atom representations are turned on based on atom index.
  • Padding around interface button has been standardized to provide a uniform experience between different platforms.

Release 0.6.2 (Mar 16, 2011)

New Features:

  • performSVD() function is implemented for faster and more memory efficient principal compoment analysis.
  • extrapolateModel() function is implemented for extrapolating a coarse-grained model to an all atom model. See the usage example Extend a coarse-grained model.
  • plog() is implemented for enabling users to make log entries.

Improvements:

  • compare functions are improved to handle insertion codes.
  • HierView allows for indexing using chain identifier and residue numbers. See usage example Hierarchical view.
  • Chain allows for indexing using residue number and insertion code. See usage example Hierarchical view.
  • addCoordset() function accepts Atomic and Ensemble instances as coords argument.
  • New method HierView.getAtoms() is implemented.
  • AtomGroup set functions check the correctness of dimension of data arrays to prevent runtime problems.
  • prody pca script is updated to use the faster PCA method that uses SVD.

Changes:

  • “backbone” definition now includes the backbone hydrogen atom (Thanks to Nahren Mascarenhas for pointing to this discrepancy in the keyword definition).

Bug Fixes:

  • A bug in PCA allowed calculating covariance matrix for less than 3 coordinate sets is fixed.
  • A bug in mapOntoChain() function that caused problems when mapping all atoms is fixed.

Release 0.6.1 (Mar 2, 2011)

New Features:

  • setWWPDBFTPServer() and getWWPDBFTPServer() functions allow user to change or learn the WWPDB FTP server that ProDy uses to download PDB files. Default server is RCSB PDB in USA. User can change the default server to one in Europe or Japan.
  • setPDBMirrorPath() and getPDBMirrorPath() functions allow user to specify or learn the path to a local PDB mirror. When specified, a local PDB mirror is preferred for accessing PDB files, over downloading them from FTP servers.
  • mapOntoChain() function is improved to map backbone or all atoms.

Improvements:

  • WWPDB_PDBFetcher can download PDB files from different WWPDB FTP servers.
  • WWPDB_PDBFetcher can also use local PDB mirrors for accessing PDB files.

Changes:

  • RCSB_PDBFetcher is renamed as WWPDB_PDBFetcher.
  • mapOntoChain() and matchChains() functions accept "ca" and "bb" as subset arguments.
  • Definition of selection keyword “protein” is updated to include some non-standard amino acid abbreviations.

Bug Fixes:

  • A bug in WWPDB_PDBFetcher causing exceptions when non-string items passed in a list is fixed.
  • An important bug in parsePDB() is fixed. When parsing backbone or Cα atoms, residue names were not checked and this caused parsing water atoms with name "O" or calcium ions with name "CA".

Release 0.6 (Feb 22, 2011)

New Features:

  • Biopython module pairwise2 and packages KDTree and Blast are incorporated in ProDy package to make installation easier. Only NumPy needs to be installed before ProDy can be used. For plotting, Matplotlib is still required.
  • NMWiz is distributed with ProDy source. On Linux, if VMD is installed, ProDy installer locates VMD plugins folder and installs NMWiz. On Windows, user needs to follow a separate set of instructions (see NMWiz).
  • Gamma class is implemented for facilitating use of force constants based on atom type, residue type, or property. An example derived classes are GammaStructureBased and GammaVariableCutoff.
  • calcTempFactors() function is implemented to calculate theoretical temperature factors.
  • 5 new ProDy Applications are implemented, and existing scripts are improved to output figures.
  • getModel() method is implemented to make function development easier.
  • resetTicks() function is implemented to change X and/or Y axis ticks in plots when there are discontinuities in the plotted data.

Improvements:

  • ANM.buildHessian() and GNM.buildKirchhoff() classes are improved to accept Gamma instances or other custom function as gamma argument. See also Developing Gamma Functions.
  • Select class is changed to treat single word keywords differently, e.g. “backbone” or “protein”. They are interpreted 10 times faster and in use achieve much higher speed-ups when compared to composite selections. For example, using the keyword “calpha” instead of the name CA and protein, which returns the same selection, works >20 times faster.
  • Optimizations in Select class to increase performance (Thanks to Paul McGuire for providing several Pythonic tips and Pyparsing specific advice).
  • applyBiomolecularTransformations() function is improved to handle large biomolecular assemblies.
  • Performance optimizations in parsePDB() and other functions.
  • Ensemble class accepts Atomic instances and automatically adds coordinate sets to the ensemble.

Changes:

  • PDBlastRecord is renamed as PDBBlastRecord.
  • NMA instances can be index using a list or tuple of integers, e.g. anm[1,3,5].
  • “ca”, “bb”, and “sc” keywords are defined as short-hands for “calpha”, “backbone”, and “sidechain”, respectively.
  • Behavior of calcANM() and calcGNM() functions have changed. They return the atoms used for calculation as well.

Bug Fixes:

NMWiz

  • NMWiz can be used as a graphical interface to ProDy. ANM or PCA calculations can be performed for molecules that are loaded in VMD.
  • User can set default color for arrow graphics and paths to ANM and PCA scripts.
  • Optionally, NMWiz can preserve the current view in VMD display window when loading a new dataset. Check the box in the NMWiz GUI main window.
  • A bug that prevented selecting residues from plot window is fixed.

Release 0.5.3 (Feb 11, 2011)

New Features:

  • Membership, equality, and non-equality test operation are defined for all atomic classes. See Operations on Selections.
  • Two functions are implemented for dealing with anisotropic temperature factors: calcADPAxes() and buildADPMatrix().
  • NMA.setEigens() and NMA.addEigenpair() methods are implemented to assist analysis of normal modes calculated using external software.
  • parseNMD() is implemented for parsing NMD files.
  • parseModes() is implemented for parsing normal mode data.
  • parseArray() is implementing for reading numeric data, particularly normal mode data calculated using other software for analysis using ProDy.
  • The method in [BH02] to calculate overlap between covariance matrices is implemented as calcCovOverlap() function.
  • trimEnsemble() to trim Ensemble instances is implemented.
  • checkUpdates() to check for ProDy updates is implemented.

Changes:

  • Change in default behavior of parsePDB() function. When alternate locations exist, those indicated by A are parsed. For parsing all alternate locations user needs to pass altloc=True argument.
  • getSumOfWeights() is renamed as calcSumOfWeights().
  • mapAtomsToChain() is renamed as mapOntoChain().
  • ProDyStartLogFile() is renamed as startLogfile().
  • ProDyCloseLogFile() is renamed as closeLogfile().
  • ProDySetVerbosity() is renamed as changeVerbosity().

Improvements:

  • A few bugs in ensemble and dynamics classes are fixed.
  • Improvements in RCSB_PDBFetcher allow it not to miss a PDB file if it exists in the target folder.
  • writeNMD() is fixed to output B-factors (Thanks to Dan Holloway for pointing it out).

Release 0.5.2 (Jan 12, 2011)

Bug Fixes:

  • An important fix in sampleModes() function was made (Thanks to Alberto Perez for finding the bug and suggesting a solution).

Improvements:

  • Improvements in ANM.calcModes(), GNM.calcModes(), and PCA.calcModes() methods prevent Numpy/Scipy throwing an exception when more than available modes are requested by the user.
  • Improvements in blastPDB() enable ProDy throw an exception when no internet connection is found, and warn user when downloads fail due to restriction in network regulations (Thanks to Serkan Apaydin for helping identify these improvements).
  • New example Write PDB file.

Release 0.5.1 (Dec 31, 2010)

Changes in dependencies:

  • Scipy (linear algebra module) is not required package anymore. When available it replaces Numpy (linear algebra module) for greater flexibility and efficiency. A warning message is printed when Scipy is not found.
  • Biopython KDTree module is not required for ENM calculations (specifically for building Hessian (ANM) or Kirchoff (GNM) matrices). When available it is used to increase the performance. A warning message is printed when KDTree is not found.

Release 0.5 (Dec 21, 2010)

New Features:

Bug Fixes:

  • A bug in parsePDBStream() that caused unidentified errors when a model in a multiple model file did not have the same number of atoms is fixed.

Changes:

  • Iterating over a Chain instance yields Residue instances.
  • Vector instantiation requires an array only. name is an optional argument.
  • Functions starting with get and performing a calculations are renamed to start with calc, e.g. getRMSD() is now calcRMSD().

Release 0.2 (Nov 16, 2010)

Important Changes:

  • Single word keywords not followed by “and” logical operator are not accepted, e.g. “protein within 5 of water” will raise a SelectionError, use “protein and within 5 of water” instead.
  • findMatchingChains() is renamed to matchChains().
  • showOverlapMatrix() is renamed to showOverlapTable().
  • Modules are reorganized.

New Features:

Improvements:

  • proteins.compare functions are improved. Now they perform sequence alignment if simple residue number/identity based matchin does not work, or if user passes pwalign=True argument. This impacts the speed of X-ray ensemble analysis.
  • Select can cache data optionally. This results in speeds up from 2 to 50 folds depending on number of atoms and selection operations.
  • Implementation of showProjection() is completed.

NMWiz

Release 0.2.3

  • For each mode a molecule for drawing arrows and a molecule for showing animation is formed in VMD on demand. NMWiz remembers a color associated with a mode.
  • Deselecting a residue by clicking on a plot is possible.
  • A bug causing incorrect parsing of NMD files from ANM server is fixed.

Release 0.2.2

  • Selection string option allows user to show a subset of arrows matching a VMD selection string. Optionally, this selection string may affect protein and animation representations.
  • A bug that caused problems when over plotting modes is removed.
  • A bug affecting line width changes in plots is removed.
  • Selected residue representations are colored according to the color of the plot.

Release 0.2.1

  • Usability improvements.
  • Loading the same data file more than once is prevented.
  • If a GUI window for a dataset is closed, it can be reloaded from the main window.
  • A dataset and GUI can be deleted from the VMD session via the main window.

Release 0.2

  • Instant documentation is improved.
  • Problem with clearing selections is fixed.
  • Plotting options frame is populated.
  • Multiple modes can be plotted on the same canvas.

Release 0.1.2 (Nov 9, 2010)

  • Important bug fixes and improvements in NMA helper and plotting functions.
  • Documentation updates and improvements.

Release 0.1.1 (Nov 8, 2010)

  • Important bug fixes and improvements in chain comparison functions.
  • Bug fixes.
  • Source clean up.
  • Documentation improvements.

Release 0.1 (Nov 7, 2010)

  • First release.

Navigation

© Copyright 2010-2013, Ahmet Bakan. Last updated on May 13, 2013. Created using Sphinx 1.1.3.