Merge branch 'main' of github.com:prody/ProDy into hpb

Author: jamesmkrieger

CONTRIBUTING.md

@@ -1,4 +1,4 @@
-This is a first draft, which is based on the ProDy website page http://prody.csb.pitt.edu/manual/devel/develop.html. Please see that page for updates.
+This is a first draft, which is based on the ProDy website page http://www.bahargroup.org/prody/manual/devel/develop.html. Please see that page for updates.
 
 
 Install Git and a GUI

PKG-INFO

@@ -32,23 +32,23 @@ Description: .. image:: https://secure.travis-ci.org/prody/ProDy.png?branch=mast
         You can run ProDy on all major platforms.  For download and installation
         instructions see:
 
-        * http://prody.csb.pitt.edu/downloads
+        * http://www.bahargroup.org/prody/downloads
 
 
         DOCUMENTATION
         -------------
 
-        * Homepage: http://prody.csb.pitt.edu/
+        * Homepage: http://www.bahargroup.org/prody/
 
-        * Tutorials: http://prody.csb.pitt.edu/tutorials
+        * Tutorials: http://www.bahargroup.org/prody/tutorials
 
-        * Reference: http://prody.csb.pitt.edu/manual
+        * Reference: http://www.bahargroup.org/prody/manual
 
-        * Applications: http://prody.csb.pitt.edu/manual/apps
+        * Applications: http://www.bahargroup.org/prody/manual/apps
 
-        * NMWiz GUI: http://prody.csb.pitt.edu/nmwiz
+        * NMWiz GUI: http://www.bahargroup.org/prody/nmwiz
 
-        * Changes: http://prody.csb.pitt.edu/manual/release
+        * Changes: http://www.bahargroup.org/prody/manual/release
 
 
         SOURCE CODE

docs/apps/prody/index.rst

@@ -40,4 +40,4 @@ the output files are prefixed with :file:`p38_anm`::
 
 The output file :file:`p38_anm.nmd` can be visualized using `NMWiz`_.
 
-.. _NMWiz: http://prody.csb.pitt.edu/nmwiz/
+.. _NMWiz: http://www.bahargroup.org/prody/nmwiz/

docs/conf.py

@@ -158,13 +158,13 @@
     'numpy': ('http://docs.scipy.org/doc/numpy/', None),
     'scipy': ('http://docs.scipy.org/doc/scipy/reference/', None),
     'matplotlib': ('http://matplotlib.sourceforge.net/', None),
-    'prodywebsite': ('http://prody.csb.pitt.edu/', None),
+    'prodywebsite': ('http://www.bahargroup.org/prody/', None),
 }
 
 rst_epilog = u"""
 
-.. _ProDy: http://prody.csb.pitt.edu
-.. _Tutorials: http://prody.csb.pitt.edu/tutorials
+.. _ProDy: http://www.bahargroup.org/prody
+.. _Tutorials: http://www.bahargroup.org/prody/tutorials
 .. _NMWiz: http://csb.pitt.edu/NMWiz
 .. _VMD: http://www.ks.uiuc.edu/Research/vmd
 .. _PDB: http://www.pdb.org

docs/index.rst

@@ -11,7 +11,7 @@ ProDy Manual
 .. only:: rtd and html
 
     This is a partial copy of ProDy documentation.  Please visit
-    `ProDy Homepage <http://prody.csb.pitt.edu>`_ for complete
+    `ProDy Homepage <http://www.bahargroup.org/prody>`_ for complete
     documentation with tutorials.
 
 

prody/compounds/pdbligands.py

@@ -8,6 +8,7 @@
 from prody import LOGGER, SETTINGS, getPackagePath, PY3K
 from prody.atomic import AtomGroup, ATOMIC_FIELDS
 from prody.utilities import openFile, makePath, openURL
+from .ccd import parseCCD
 
 __all__ = ['PDBLigandRecord', 'fetchPDBLigand', 'parsePDBLigand']
 
@@ -18,13 +19,20 @@ def __init__(self, data):
         self._rawdata = data
 
     def getCanonicalSMILES(self):
-        return self._rawdata['CACTVS_SMILES_CANONICAL']
+        canonical = None
+        for row in self._rawdata[2].data:
+            if row['_pdbx_chem_comp_descriptor.type'] == 'SMILES_CANONICAL':
+                canonical = row['_pdbx_chem_comp_descriptor.descriptor']
+        return canonical
 
 
 def fetchPDBLigand(cci, filename=None):
-    """Fetch PDB ligand data from PDB_ for chemical component *cci*.
+    """Handle PDB ligand data from PDB_ for chemical component *cci*.
     *cci* may be 3-letter chemical component identifier or a valid XML
-    filename.  If *filename* is given, XML file will be saved with that name.
+    filename. If *filename* is given, XML file will be saved with that name.
+
+    This function may not work as _PDB is not hosting ligand XML files anymore.
+    It is kept for use with existing ligand XML files.
 
     If you query ligand data frequently, you may configure ProDy to save XML
     files in your computer.  Set ``ligand_xml_save`` option **True**, i.e.
@@ -50,6 +58,7 @@ def fetchPDBLigand(cci, filename=None):
     ideal (energy minimized) coordinate sets:
 
     .. ipython:: python
+       :okexcept:
 
        from prody import *
        ligand_data = fetchPDBLigand('STI')
@@ -233,8 +242,8 @@ def fetchPDBLigand(cci, filename=None):
     return dict_
 
 
-def parsePDBLigand(cci, filename=None):
+def parsePDBLigand(cci):
     """See :func:`.fetchPDBLigand`"""
-    lig_dict = fetchPDBLigand(cci, filename)
+    lig_dict = parseCCD(cci)
     return PDBLigandRecord(lig_dict)
 

prody/ensemble/functions.py

@@ -411,7 +411,7 @@ def buildPDBEnsemble(atomics, ref=None, title='Unknown', labels=None, atommaps=N
 
     if 'mapping_func' in kwargs:
         raise DeprecationWarning('mapping_func is deprecated. Please see release notes for '
-                                 'more details: http://prody.csb.pitt.edu/manual/release/v1.11_series.html')
+                                 'more details: http://www.bahargroup.org/prody/manual/release/v1.11_series.html')
     start = time.time()
 
     if not isListLike(atomics):

prody/proteins/pdbfile.py

@@ -490,11 +490,10 @@ def _parsePDBLines(atomgroup, lines, split, model, chain, subset,
         elements = np.zeros(asize, dtype=ATOMIC_FIELDS['element'].dtype)
         bfactors = np.zeros(asize, dtype=ATOMIC_FIELDS['beta'].dtype)
         occupancies = np.zeros(asize, dtype=ATOMIC_FIELDS['occupancy'].dtype)
-        anisou = None
         siguij = None
     else:
         radii = np.zeros(asize, dtype=ATOMIC_FIELDS['radius'].dtype)
-
+    anisou = None
     asize = 2000 # increase array length by this much when needed
 
     start = split
@@ -537,19 +536,7 @@ def _parsePDBLines(atomgroup, lines, split, model, chain, subset,
     dec = True
     while i < stop:
         line = lines[i]
-        if not isPDB:
-            fields = line.split()
-            if len(fields) == 10:
-                fields.insert(4, '')
-            elif len(fields) != 11:
-                LOGGER.warn('wrong number of fields for PQR format at line %d'%i)
-                i += 1
-                continue
-
-        if isPDB:
-            startswith = line[0:6].strip()
-        else:
-            startswith = fields[0]
+        startswith = line[0:6].strip()
 
         if startswith == 'ATOM' or startswith == 'HETATM':
             if isPDB:
@@ -560,6 +547,17 @@ def _parsePDBLines(atomgroup, lines, split, model, chain, subset,
                 else:
                     resname = line[17:20].strip()
             else:
+                fields = line.split()
+                if fields[5].find('.') != -1:
+                    # coords too early as no chid
+                    fields.insert(4, '')
+                if len(fields) != 11:
+                    try:
+                        fields = fields[:6] + [line[30:38].strip(), line[38:46].strip(), line[46:54].strip()] + line[54:].split()
+                    except:
+                        LOGGER.warn('wrong number of fields for PQR format at line %d'%i)
+                        i += 1
+                        continue
                 atomname= fields[2]
                 resname = fields[3]
 
@@ -1689,7 +1687,10 @@ def writePQRStream(stream, atoms, **kwargs):
     write = stream.write
 
     calphas = atoms.ca
-    ssa = calphas.getSecstrs()
+    if calphas is not None:
+        ssa = calphas.getSecstrs()
+    else:
+        ssa = None
     helix = []
     sheet = []
     if ssa is not None:

prody/proteins/waterbridges.py

@@ -638,9 +638,9 @@ def analyseFrame(i, interactions_all):
 
             if max_proc == 1:
                 interactions_all = []
-                for j0, frame0 in enumerate(traj, start=start_frame):
+                for i in range(len(atoms.getCoordsets()[start_frame:stop_frame])):
                     interactions_all.append([])
-                    analyseFrame(j0, start_frame, frame0, interactions_all)
+                    analyseFrame(i, interactions_all)
             else:
                 with mp.Manager() as manager:
                     interactions_all = manager.list()

prody/tests/datafiles/__init__.py

@@ -381,6 +381,26 @@
         'file': 'pdb7pbl.pdb',
         'atoms': 31396,
         'nucleoside': 252
+    },
+    'pqrUnknown': {
+        'file': 'pqr_snippet1.pqr',
+        'atoms': 5,
+        'models': 1
+    },
+    'pqrTranscomp': {
+        'file': 'pqr_snippet2_transcomp.pqr',
+        'atoms': 5,
+        'models': 1
+    },
+    'pqrFpocket': {
+        'file': 'pqr_snippet3_fpocket.pqr',
+        'atoms': 5,
+        'models': 1
+    },
+    'pqrPymol': {
+        'file': 'pqr_snippet4_pymol.pqr',
+        'atoms': 5,
+        'models': 1
     }
 }