openmolecules.org

 
Home » DataWarrior » Bug Reports » Generate conformers - write to file does nothing
Generate conformers - write to file does nothing [message #984] Wed, 08 July 2020 21:01 Go to next message
Paul is currently offline  Paul
Messages: 16
Registered: May 2019
Location: United States
Junior Member
Hello,

When I try to write conformers into file, a 1KB dwar file is generated, which contains only what appear to be the file markup tags, i.e <datawarrior-fileinfo> . . . </datawarrior properties>.

Conformers are generated successfully and get a new 3D molecule view pane when "write to File" is not selected, and Max conformer count per stereoisomer is set to 1.

When Max conformer count is set to 2 or larger, a single conformer is displayed. The appropriate number of energies is given in the new Energy column, but only 1 conformer is provided.

Best Regards,

Paul
Re: Generate conformers - write to file does nothing [message #987 is a reply to message #984] Thu, 09 July 2020 00:58 Go to previous messageGo to next message
nbehrnd is currently offline  nbehrnd
Messages: 224
Registered: June 2019
Senior Member
Hello Paul,

the generation of conformers may, but need not to, include the application of a force field. DW
includes MMFF94s, MMFF94s+, and Indorsia; if used (which might be your default) some of the
conformers may coalesce. The fourth option next to minimize energy is «don't minimize»; then
retaining better the results of the parameters «Algorithm» and «Initial torsions» in display
and file written.

Norwid
Re: Generate conformers - write to file does nothing [message #988 is a reply to message #987] Thu, 09 July 2020 02:46 Go to previous messageGo to next message
Paul is currently offline  Paul
Messages: 16
Registered: May 2019
Location: United States
Junior Member
Thanks - I have been able to use the various force fields with success, but only when generating a single conformer as new columns in the DWAR document.

Regardless of the force field, or choosing not to minimize, the output from generate conformers → write to file is a blank DWAR file.

Paul
Re: Generate conformers - write to file does nothing [message #991 is a reply to message #988] Thu, 09 July 2020 10:10 Go to previous messageGo to next message
nbehrnd is currently offline  nbehrnd
Messages: 224
Registered: June 2019
Senior Member
Hello Paul,

starting with the liquid_crystal.dwar file attached below, and using the parameters
shown in the screen photo below

index.php?t=getfile&id=241&private=0

the newly written conformers_test.dwar contains multiple models. Not 16 which was
set as the maximum number, but at least four. May you share a minimal example,
including the parameters you use, leading to the empty .dwar file?

Norwid
Re: Generate conformers - write to file does nothing [message #992 is a reply to message #991] Thu, 09 July 2020 18:23 Go to previous messageGo to next message
Paul is currently offline  Paul
Messages: 16
Registered: May 2019
Location: United States
Junior Member
Thank you - your files helped me to determine the reason my attempts to write conformers to files have failed.

After I downloaded your files, I was able to write conformers to file from those dwar files.

Immediately noticeable was that your file contained only an option "Structure" in the Structure Column dropdown menu in the generate conformers dialog.

My files contained only an option "Structure of SMILES" in the same dialog, even though they also contain 3D structures (data column 3-D structure of SMILES) that are clearly visible in Datawarrior's 3D-Structure pane. With this option, writing conformers to file generates an empty DWAR container.

When I Delete the SMILES and/or structure of SMILES column(s), the image in the 3D-Structure pane disappears. This suggested that the conformer displayed in the 3D viewer pane was somehow dynamic, but when I edit the SMILES, neither the structure of SMILES nor the previously generated conformer update to the new SMILES.

Since you had shown me that writing conformers to file could work, I then continued to try different approaches to obtaining a file without a Structure of SMILES column, but instead a Structure column.

I saved my file as SDF, choosing 3D structure if available. When I opened the new SDF, I then had the desired "Structure" column (not "Structure of SMILES). The SDF still contained a SMILES column, but DWAR used the 3D structure instead of generating structures from SMILES. Writing conformers to file now works!

However, if a DWAR file contains both SMILES and 3D structures, Datawarrior uses the SMILES > Structure of SMILES and will not write conformers to file.

At last, I successfully wrote conformers to file (either from the SDF or after saving it back to DWAR, so long as the SDF did not contain a SMILES column)! The SMILES column must be deleted after opening the SDF, not from the DWAR file, then saved to a DWAR file.

So then, I am thinking that there must be a better way to accomplish these ends without (1) generating a conformer from structure of smiles in Datawarrior, (2) Saving it to SDF, and (3) Using the new SDF to generate conformers (or deleting the smiles column then saving it back to DWAR).

On further exploration, I determined that I can write conformers to file without the SDF intermediate if I generate 2D atom coordinates (from Structure of SMILES, of course).

The generate conformers dialog does not change - the only structure available in the dropdown menu is still Structure of SMILES (even though I now have columns for 3D-structure of smiles, Structure of SMILES [2d coordinates], and Structure of SMILES).


Conclusions:

Datawarrior's Generate Conformers > Structure Column dropdown menu lists only one structure column, even when more are available.

Datawarrior will generate conformers from Structure of SMILES within a current document (SDF or DWAR), but not to a new document. Attempts to write to file generate empty DWAR containers.

Datawarrior will write conformers to file when the DWAR-generated Structure of SMILES [2D coordinates] data column exists, even though there is no apparent difference between the generate conformers dialog (it indicates conformers are generated from Structure of SMILES).


Thank you very much for helping me work this out.


Now I'm wondering - does the Structure dropdown menu ever populate with more than one option?


Best Regards,

Paul

Re: Generate conformers - write to file does nothing [message #996 is a reply to message #984] Fri, 10 July 2020 13:22 Go to previous message
thomas is currently offline  thomas
Messages: 706
Registered: June 2014
Senior Member
I cannot reproduce this. If I generate 10 random molecules with otherwise default settings, and if I then generate conformers (Max conformer count: 2, otherwise default), I get two conformers for most molecules:

index.php?t=getfile&id=244&private=0

I tried with V5.2.1 as well as with the most recent dev version. Writing 16 conformers per structure into a file also worked. Maybe it somehow is related to your input structures? Could you share a problematic structure?
  • Attachment: t.png
    (Size: 40.25KB, Downloaded 637 times)
Previous Topic: Uncaught Exception:/ by zero error when using structure view.
Next Topic: Dropbox dead links
Goto Forum:
  


Current Time: Wed Oct 30 21:02:16 CET 2024

Total time taken to generate the page: 0.00687 seconds