Closed j-wags closed 4 years ago
Summary: Only three unique bonds, ~10 unique angles and impropers, but lots of proper torsions
(summary: All errors under 3e-3 kJ/mol)
Ah, this may have been prematurely merged. But it's actually looking OK. I'll open an Issue for the remaining discrepancies.
Closes #19
This PR-in-progress adds
ConvertTripeptideParameters.py
: Takes all residue permutations generated byGenerateTripeptides.sh
, and finds all valence parameters which have at least one tagged atom in both residues 1 and 2. Then, makes a SMARTS for each parameter, including only the tagged atoms and immediate neighbors. The vast majority of these parameters end up being redundant, and are thus automatically deduplicated. Since it's short, I'll post the resulting OFFXML in comments as we work on this.test_parameterize_tripeptides.py
: Constructs a ForceField by first loading the "backbone parameters" (ConvertTripeptideParameters.py
/test_backbone.offxml
) and then the residue parameters from our original workflow (ConvertParameters.py
/test.offxml
). Appending them in this way makes a "safe" hierarchy, where the per-residue parameters are "more specific" and tehrefore have precedence over the potentially-promiscuous backbone parameters.