Skip to content

fixed 2'-hydroxyl dihedral calculation#2231

Merged
orbeckst merged 3 commits intodevelopfrom
issue-2218-2-OH
Apr 3, 2019
Merged

fixed 2'-hydroxyl dihedral calculation#2231
orbeckst merged 3 commits intodevelopfrom
issue-2218-2-OH

Conversation

@orbeckst
Copy link
Member

@orbeckst orbeckst commented Apr 2, 2019

Fixes #2218

Changes made in this Pull Request:

  • fixed hydrogen selection (H2'' --> H2')
  • updated nucleic acid test file: now contains properly equilibrated fully solvated MD system (CHARMM36m)
  • updated all regression tests to correspond to the new RNA test system
  • removed variable datafiles.NUCL; instead use datafiles.RNA_PSF and datafiles.RNA_PDB

PR Checklist

  • Tests?
  • n/a Docs?
  • CHANGELOG updated?
  • Issue raised/referenced?

- fixes #2218
- RNA 2-hydroxyl selection now agrees with our docs:
  H2'' (wrong) --> H2' (correct)
@orbeckst orbeckst requested a review from richardjgowers April 2, 2019 06:29
- replaced 1k5i.pdb (which had the 2'-OH H in strange positions) with a simulation
  system of model 1 of 1k5i in CHARMM36m (generated by CHARMM-GUI and equilibrated
  for 1 ns with Gromacs 2018.6). Octahedral box, K+ counter ions, TIP3P water.
  - topology datafiles.RNA_PSF data/analysis/1k5i_c36.psf
  - coordinates datafiles.RNA_PDB data/analysis/1k5i_c36.pdb
- added data files as gzipped datafiles (instead of bz2 because of #2232)
- removed datafiles.NUCL = data/1k5i
- updated all regression tests in test_nucleic.py with numbers appropriate for
  the new test file; manually confirmed the values for the 2'-hydroxyl dihedrals
  in the test by manually selecting and calculating in VMD
- replaced NUCL in atomselections test
  (also removed superfluous/outdated utils/test_nuclinfo)
- updated CHANGELOG
@codecov
Copy link

codecov bot commented Apr 2, 2019

Codecov Report

Merging #2231 into develop will not change coverage.
The diff coverage is n/a.

Impacted file tree graph

@@           Coverage Diff            @@
##           develop    #2231   +/-   ##
========================================
  Coverage    90.71%   90.71%           
========================================
  Files           15       15           
  Lines         1982     1982           
========================================
  Hits          1798     1798           
  Misses         184      184

Continue to review full report at Codecov.

Legend - Click here to learn more
Δ = absolute <relative> (impact), ø = not affected, ? = missing data
Powered by Codecov. Last update 9428866...4f8dc17. Read the comment docs.

I got a weird failure in PR #2231 in Travis job
https://travis-ci.com/MDAnalysis/mdanalysis/jobs/189553940 :

   Exception occurred:
     File "/home/travis/miniconda/envs/test/lib/python3.5/site-packages/sphinx/ext/napoleon/docstring.py", line 703, in _parse_raises_section
       m = self._name_rgx.match(_type).groupdict()
   AttributeError: 'NoneType' object has no attribute 'groupdict'

and the docs build find locally (under Python 3.6)
@orbeckst
Copy link
Member Author

orbeckst commented Apr 2, 2019

CI was failing the doc build in https://travis-ci.com/MDAnalysis/mdanalysis/jobs/189553940 :

 Exception occurred:
   File "/home/travis/miniconda/envs/test/lib/python3.5/site-packages/sphinx/ext/napoleon/docstring.py", line 703, in _parse_raises_section
     m = self._name_rgx.match(_type).groupdict()
 AttributeError: 'NoneType' object has no attribute 'groupdict'

Locally under 3.6 the docs build just fine. I am trying bumping our doc build environment to the currently "best", i.e., 3.7.

@orbeckst orbeckst merged commit 6a60c8a into develop Apr 3, 2019
@orbeckst orbeckst deleted the issue-2218-2-OH branch April 3, 2019 01:30
@orbeckst orbeckst mentioned this pull request Apr 3, 2019
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

nucleic acid analysis 2-hydroxyl dihedral bug

2 participants