Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

DNA-RNA hybrids #72

Merged
merged 69 commits into from
Oct 9, 2023
Merged

DNA-RNA hybrids #72

merged 69 commits into from
Oct 9, 2023

Conversation

lorenzo-rovigatti
Copy link
Owner

Add the possibility of simulating DNA-RNA hybrids (by @eryykr )

eryykr and others added 22 commits May 23, 2023 16:15
Fixed bugs caused by double inheritance (i think), updated function for checking input sanity to accommodate for both DNA and RNA in the system.
Interaction type with modified backbone potential for relaxing systems containing both DNA and RNA.
Adding the new topology format to the hybrid model
Trying to fix (suspected) undefined behaviour
Version of hybrid model which uses inheritance
Add the DNA-RNA hybrid model
This commit adds 6 more tests to the `run` level. These tests check that the DNA, RNA and DRH interactions compute the right energies for a short nicked duplex
The member was used to check in the RNA and DNA classes whether the instance of the interaction was actually a DRH interaction. We now check this with a local variable since we need this information only at the very beginning
src/Interactions/DNAInteraction.cpp Outdated Show resolved Hide resolved
src/Interactions/DRHInteraction_relax.cpp Show resolved Hide resolved
src/defs.h Outdated Show resolved Hide resolved
Comment on lines 150 to 152

/// Needed for DNA-RNA interactions
char acid_type;
Copy link
Owner Author

@lorenzo-rovigatti lorenzo-rovigatti Oct 8, 2023

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@eryykr Can we do this bookkeeping in the DRH interaction rather than here? Do you foresee any issues?

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think we could so long as there is a way of distinguishing between DNA and RNA nucleotides without acid_type

Copy link
Owner Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

OK! Have a look at my changes in b695623 and let me know if you agree with them.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good!

* Remove some unusued define's
* Turn a char[] into a std::string
* Move the information about the type of each nucleotide (R or D) from BaseParticle to the DRHInteraction
@lorenzo-rovigatti
Copy link
Owner Author

lorenzo-rovigatti commented Oct 9, 2023

While working on this branch, yesterday I have also added a few more tests that check the values of the energy terms present in a nicked double strand as computed by the DNA2, RNA2 and DNA_RNA force fields. These are run automatically every time a push is performed, but you can also run them manually with make test_run.
Pinging @ErikPoppleton @sulcgroup: if you collaborate with someone who works on the code you can tell them to run the tests locally and check that their changes do not impact the core code.

These tests will end up in master as soon as we merge the branch (which will happen soon, hopefully)

@lorenzo-rovigatti lorenzo-rovigatti merged commit 7c97406 into master Oct 9, 2023
3 checks passed
@lorenzo-rovigatti lorenzo-rovigatti deleted the dna_rna branch October 9, 2023 15:02
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.

2 participants