A Practical Quantum Mechanics Molecular Mechanics Method for the Dynamical Study of Reactions in Biomolecules

Two-step ATP-driven opening of cohesin head

The cohesin ring is a protein complex composed of four core subunits: Smc1A, Smc3, Rad21 and Stag1/2. It is involved in chromosome segregation, DNA repair, chromatin organization and transcription regulation. Opening of the ring occurs at the "head" structure, formed of the ATPase domains of Smc1A and Smc3 and Rad21. We investigate the mechanisms of the cohesin ring opening using techniques of free molecular dynamics (MD), steered MD and quantum mechanics/molecular mechanics MD (QM/MM MD). The study allows the thorough analysis of the opening events at the atomic scale: i) ATP hydrolysis at the Smc1A site, evaluating the role of the carboxy-terminal domain of Rad21 in the process; ii) the activation of the Smc3 site potentially mediated by the movement of specific amino acids; and iii) opening of the head domains after the two ATP hydrolysis events. Our study suggests that the cohesin ring opening is triggered by a sequential activation of the ATP sites in which ATP hydrolysis at the Smc1A site induces ATPase activity at the Smc3 site. Our analysis also provides an explanation for the effect of pathogenic variants related to cohesinopathies and cancer.

Quantum Mechanics/Molecular Mechanics Free Energy Maps and Nonadiabatic Simulations for a Photochemical Reaction in DNA: Cyclobutane Thymine Dimer

The absorption of ultraviolet radiation by DNA may result in harmful genetic lesions that affect DNA replication and transcription, ultimately causing mutations, cancer, and/or cell death. We analyze the most abundant photochemical reaction in DNA, the cyclobutane thymine dimer, using hybrid quantum mechanics/molecular mechanics (QM/MM) techniques and QM/MM nonadiabatic molecular dynamics. We find that, due to its double helix structure, DNA presents a free energy barrier between nonreactive and reactive conformations leading to the photolesion. Moreover, our nonadiabatic simulations show that most of the photoexcited reactive conformations return to standard B-DNA conformations after an ultrafast nonradiative decay to the ground state. This work highlights the importance of dynamical effects (free energy, excited-state dynamics) for the study of photochemical reactions in biological systems.

MEPSA: minimum energy pathway analysis for energy landscapes

From conformational studies to atomistic descriptions of enzymatic reactions, potential and free energy landscapes can be used to describe biomolecular systems in detail. However, extracting the relevant data of complex 3D energy surfaces can sometimes be laborious. In this article, we present MEPSA (Minimum Energy Path Surface Analysis), a cross-platform user friendly tool for the analysis of energy landscapes from a transition state theory perspective. Some of its most relevant features are: identification of all the barriers and minima of the landscape at once, description of maxima edge profiles, detection of the lowest energy path connecting two minima and generation of transition state theory diagrams along these paths. In addition to a built-in plotting system, MEPSA can save most of the generated data into easily parseable text files, allowing more versatile uses of MEPSA's output such as the generation of molecular dynamics restraints from a calculated path.

AVAILABILITY AND IMPLEMENTATION

MEPSA is freely available (under GPLv3 license) at: http://bioweb.cbm.uam.es/software/MEPSA/ CONTACT: pagomez@cbm.csic.es.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.