A General Computational Framework for Sequence-Level Thermodynamic Analysis of Modeled RNA–DNA and DNA–DNA Duplex States: A Proof-of-Concept Study in HIV-1

Inside the microscopic world of a cell, genetic material is not always a static double helix. Sometimes, the strands of DNA that usually pair together with each other are interrupted when a strand of RNA comes along and binds to one of them, pushing the other DNA strand aside. This three-part structure, where RNA and DNA hold hands while a third strand waits in the wings, is known as an R-loop. These formations are natural features of biology, playing roles in how genes are read and regulated, but they can also cause instability if they appear in the wrong place or at the wrong time. Understanding where and why these structures might form is a key question for scientists studying viruses, because the genetic makeup of a virus can influence how it interacts with the human cells it infects. The challenge lies in distinguishing between a structure that is merely possible according to the laws of chemistry and one that actually happens inside a living cell.

A new study takes a different approach to this problem by focusing entirely on the chemical energy required to build these structures. The researchers did not look inside a cell or measure biological activity directly. Instead, they built a computational framework to calculate the thermodynamic stability of specific viral sequences. In this context, thermodynamics refers to the energy changes that occur when molecules come together; a more negative energy value generally means a bond is stronger and more likely to form. The team focused on the Human Immunodeficiency Virus type 1, or HIV-1, using its most common reference genetic sequence as a map. They selected sixteen distinct regions along the virus's genetic code, ranging from short segments to longer stretches, and subjected each one to a rigorous two-step calculation.

First, the researchers modeled what would happen if the RNA version of a specific viral region tried to bind with a complementary DNA strand. They used a specialized software tool designed to predict how these RNA-DNA hybrids behave. Next, they modeled what would happen if that same DNA sequence tried to bind with its own matching DNA partner, essentially asking how stable the natural DNA double helix is in that same spot. By comparing the energy of the RNA-DNA pairing against the energy of the DNA-DNA pairing, they created a specific contrast value for each region. This comparison allowed them to see if the virus's genetic code has a natural chemical tendency to form the three-stranded R-loop structure more easily than it holds onto its own double helix.

The results revealed a clear and consistent pattern across the virus's genome. Out of the sixteen regions tested, fifteen showed a negative contrast value, meaning the RNA-DNA hybrid was calculated to be more stable than the corresponding DNA-DNA structure. The average difference was eight kilocalories per mole, a significant margin in the world of molecular chemistry. Only one region, a segment known as U5, showed a positive value, indicating that in that specific spot, the DNA preferred to stay paired with itself. The statistical analysis confirmed that this trend was not a random occurrence; the likelihood of seeing this many negative results by chance was extremely low. The researchers also checked to ensure the result wasn't skewed by overlapping sections of the genome, and the pattern held steady even when those overlapping parts were removed from the calculation.

However, the study is careful to define the limits of what these numbers mean. The findings are strictly a prediction based on chemical models and do not prove that these R-loops actually form inside a human body or that they cause the virus to integrate into human DNA. The calculations did not account for the complex environment of a cell, such as the presence of proteins that might block these interactions, the physical twisting of the DNA, or the concentration of molecules. The author describes their work as a way to prioritize which parts of the viral genome are most chemically likely to form these structures, serving as a guide for future experiments rather than a final answer. In essence, the study provides a detailed map of where the virus's genetic code is most eager to make a specific type of bond, offering a testable hypothesis for scientists to verify in the laboratory.

The research also explored whether the length of the genetic segment or the specific mix of chemical building blocks within it influenced the results. While there were some hints that regions with certain chemical compositions showed stronger tendencies, the small number of regions tested meant these observations were preliminary. The primary conclusion remains focused on the robustness of the energy difference itself. By establishing a reproducible method to compare RNA-DNA and DNA-DNA stability across a viral genome, the study offers a new tool for understanding the chemical potential of viral sequences. It suggests that for most of the HIV-1 genome, the conditions are chemically favorable for the formation of these three-stranded structures, but it leaves the question of whether nature actually uses this potential to the next generation of biological experiments.