2021-06-16 Meeting notes

Date

Jun 16, 2021

Participants

• @Jeffrey Wagner

• @Connor Davel

• @Owen Madin

Discussion topics

Item	Notes

Item	Notes
Detecting proton migration	CD – Worked on other things CD + OM – We should prioritize making a database of this
Non-geomeTRIC fixes:	CD found some sqm keywords that may help -- maxcyc and SHAKE, and one other CD – I’ve tested this out.
GeomeTRIC background and starting points	How to start a restriction scheme?
Database of failing cases/molecules	Proton vs. heavy atom migrations How to make a database of these? Initial test case Zwitterions from protein structures Non-protein cases Non-zwitterion case CD – Diversity in functional groups? JW – Possibly add N1NC(=N)C(=N)1) CD – Which molecule datasets do we have available? Minidrugbank.sdf https://github.com/openforcefield/openff-toolkit/tree/master/openff/toolkit/data/molecules Are perception differences here a major issue? JW – No, in most cases both ways to interpret an ambiguous molecule are valid. So we’ll just want to make sure to regenerate the conformers to place hydrogens consistently/accoding to the graph CD – I already do generate new conformers. So I’d load an SDF, generate new conformers (disregarding the roiginal), output those conformers to NEW sdfs, and then run those with antechamber. JW – Perfect. LOW VALUE, POSSIBLY HARD: DrugBank_X.mol2, AlhEthOH files, and FreeSolv Antechamber can read mol2 directly. But there are three “sub specifications” of mol2. So you can try to chunk out the mol2 files from the link above and see if antechamber can read any of them. burn-in-sdf Open
What’s the “gold standard”?	Is OpenEye AM1 always the goal to match? Or should we look at something like RESP, or physical property accuracy?

Reproducing Case

Two conformers of the following molecule were tested as a proof-of-concept for an AM1-restriction scheme:

Conformer 0 (no proton rearrangement)	Conformer 1 (proton rearrangement from the highlighted atom to the oxygen)

Conformer 0 (no proton rearrangement)	Conformer 1 (proton rearrangement from the highlighted atom to the oxygen)
Open Before AM1 minimization	Open Before AM1 minimization

Conformer 0 does not undergo any proton rearrangement whereas conformer 1 sees proton migration from a nitrogen to an oxygen. Since we are concerned with the affect of this proton migration on partial charge assignment, the resulting partial charges were plotted using both amber and openeye. The black dots represent openeye’s charge assignments, which where very consistent between conformer 0 and 1 and show up as a single point. In the amber case, notice differences between the charge assignment of conformer 0 and 1.

For reference, atom 15 is the highlighted nitrogen, atoms 28 and 29 are the carboxylate oxygens, and atoms 16 and 17 are the nitrogen protons (one of which migrates). Rearrangement of one of the conformers results in a large charge discrepancy between the two conformers. This is most apparent on the nitrogen and oxygens, the atoms to/from which the proton migrates.

Solutions

The following section quickly explores 3 solutions that could be easily implemented without much change in how sqm is run/called. Each solution will be compared for internal consistency (do both conformers give the ~same partial charges?) and for consistency with other methods (in this first case, openeye). While interesting and good to know, these initial solution were ultimately ineffective at producing “good” partial charges. The solutions are:

1. Restricting any and all geometry changes with amber’s “-ek maxcyc=0” flag, which restricts the number of iterations in which the geometry changes to 0.

2. Restricting atom-proton bonds with amber’s “qmshake=1” flag, which holds all atom-proton bonds fixed and so does not allow any proton migrations to occur.

3. First minimizing the geometry with MD (forcefield openff-1.3.0.offxml) and then running AM1 with “-ek maxcyc=0”. This intends to first find a global energy minimum before running AM1, hopefully eliminating the problems with method 1.

To help visualize, the resulting molecules of each simulation are tabulated below. Only Conformer 1 (the second conformer) is shown.

Each of these solutions “solved” the problem in their own way. That is to say, protons 16 and 17 are still attached to the nitrogen. Two things to note are: the folding of the molecule in Method 2 (it seems that the protons are still somehow very strong attracted to the carboxylate oxygens during simulation!) and the very slight (almost unnoticeable) change in molecule geometry in method 3. It is clear the MD simulations cannot very accurately mimic the geometry changes that occur in qm simulations.

In a similar way as before, the partial charges were compared to openeye charges.

Method 1 (maxcyc=0)

Method 2 (qmshake=1)

Method 3 (MD geometry minimization followed by AM1 with maxcyc=0)

Two things are clear at this point:

1. A new metric of success is needed. Openeye is a great first metric to test our results against. Still, in the future, it may be useful to use other qm simulations to compare the success of our methods. Another option could be physical property tests.

2. A restriction scheme is the next step. It would be useful if we could apply restraints to the starting internal geometry of the molecule, or at least of certain parts of the molecule. This is likely where this project is heading.

Action items

Decisions