Molecular simulations of DEAH-box helicases reveal control of domain flexibility by ligands: RNA, ATP, ADP, and G-patch proteins.

DEAH-box helicases use the energy from ATP hydrolysis to translocate along RNA strands. They are composed of tandem RecA-like domains and a C-terminal domain connected by flexible linkers, and the activity of several DEAH-box helicases is regulated by cofactors called G-patch proteins. We used all-atom molecular dynamics simulations of the helicases Prp43, Prp22, and DHX15 in various liganded states to investigate how RNA, ADP, ATP, or G-patch proteins influence their conformational dynamics. The simulations suggest that apo helicases are highly flexible, whereas binding of RNA renders the helicases more rigid. ATP and ADP control the stability of the RecA1-RecA2 interface, but they have only a smaller effect on domain flexibility in absence of a RecA1-RecA2 interface. Binding of a G-patch protein to DHX15 imposes a more structured conformational ensemble, characterized by more defined relative domain arrangements and by an increased conformational stability of the RNA tunnel. However, the effect of the G-patch protein on domain dynamics is far more subtle as compared to the effects of RNA or ATP/ADP. The simulations characterize DEAH-box helicase as dynamic machines whose conformational ensembles are strongly defined by the presence of RNA, ATP, or ADP and only fine-tuned by the presence of G-patch proteins.

Continuous millisecond conformational cycle of a DEAH box helicase reveals control of domain motions by atomic-scale transitions.

Helicases are motor enzymes found in every living organism and viruses, where they maintain the stability of the genome and control against false recombination. The DEAH-box helicase Prp43 plays a crucial role in pre-mRNA splicing in unicellular organisms by translocating single-stranded RNA. The molecular mechanisms and conformational transitions of helicases are not understood at the atomic level. We present a complete conformational cycle of RNA translocation by Prp43 in atomic detail based on molecular dynamics simulations. To enable the sampling of such complex transition on the millisecond timescale, we combined two enhanced sampling techniques, namely simulated tempering and adaptive sampling guided by crystallographic data. During RNA translocation, the center-of-mass motions of the RecA-like domains followed the established inchworm model, whereas the domains crawled along the RNA in a caterpillar-like movement, suggesting an inchworm/caterpillar model. However, this crawling required a complex sequence of atomic-scale transitions involving the release of an arginine finger from the ATP pocket, stepping of the hook-loop and hook-turn motifs along the RNA backbone, and several others. These findings highlight that large-scale domain dynamics may be controlled by complex sequences of atomic-scale transitions.

Simulated Tempering-Enhanced Umbrella Sampling Improves Convergence of Free Energy Calculations of Drug Membrane Permeation.

Molecular dynamics simulations have been widely used to study solute permeation across biological membranes. The potential of mean force (PMF) for solute permeation is typically computed using enhanced sampling techniques such as umbrella sampling (US). For bulky drug-like permeants, however, obtaining converged PMFs remains challenging and often requires long simulation times, resulting in an unacceptable computational cost. Here, we augmented US with simulated tempering (ST), an extended-ensemble technique that consists in varying the temperature of the system along a pre-defined temperature ladder. Simulated tempering-enhanced US (STeUS) was employed to improve the convergence of PMF calculations for the permeation of methanol and three common drug molecules. To obtain sufficient sampling of the umbrella histograms, which were computed only from the ground temperature, we modified the simulation time fraction spent at the ground temperature between 1/K and 50%, where K is the number of ST temperature states. We found that STeUS accelerates convergence, when compared to standard US, and that the benefit of STeUS is system-dependent. For bulky molecules, for which standard US poorly converged, the application of ST was highly successful, leading to a more than fivefold accelerated convergence of the PMFs. For the small methanol solute, for which conventional US converges moderately, the application of ST is only beneficial if 50% of the STeUS simulation time is spent at the ground temperature. This study establishes STeUS as an efficient and simple method for PMF calculations, thereby strongly reducing the computational cost of routine high-throughput studies of drug permeability.

GABAA Alpha 2,3 Modulation Improves Select Phenotypes in a Mouse Model of Fragile X Syndrome.

Fragile X syndrome (FXS) is the most common cause of inherited intellectual disability. FXS is caused by functional loss of the Fragile X Protein (FXP), also known as Fragile X Mental Retardation Protein (FMRP). In humans and animal models, loss of FXP leads to sensory hypersensitivity, increased susceptibility to seizures and cortical hyperactivity. Several components of the GABAergic system, the major inhibitory system in the brain, are dysregulated in FXS, and thus modulation of GABAergic transmission was suggested and tested as a treatment strategy. However, so far, clinical trials using broad spectrum GABAA or GABAB receptor-specific agonists have not yielded broad improvement of FXS phenotypes in humans. Here, we tested a more selective strategy in Fmr1 knockout (KO) mice using the experimental drug BAER-101, which is a selective GABAA α2/α3 agonist. Our results suggest that BAER-101 reduces hyperexcitability of cortical circuits, partially corrects increased frequency-specific baseline cortical EEG power, reduces susceptibility to audiogenic seizures and improves novel object memory. Other Fmr1 KO-specific phenotypes were not improved by the drug, such as increased hippocampal dendritic spine density, open field activity and marble burying. Overall, this work shows that BAER-101 improves select phenotypes in Fmr1 KO mice and encourages further studies into the efficacy of GABAA-receptor subunit-selective agonists for the treatment of FXS.

The structure of Prp2 bound to RNA and ADP-BeF3- reveals structural features important for RNA unwinding by DEAH-box ATPases.

Noncoding intron sequences present in precursor mRNAs need to be removed prior to translation, and they are excised via the spliceosome, a multimegadalton molecular machine composed of numerous protein and RNA components. The DEAH-box ATPase Prp2 plays a crucial role during pre-mRNA splicing as it ensures the catalytic activation of the spliceosome. Despite high structural similarity to other spliceosomal DEAH-box helicases, Prp2 does not seem to function as an RNA helicase, but rather as an RNA-dependent ribonucleoprotein particle-modifying ATPase. Recent crystal structures of the spliceosomal DEAH-box ATPases Prp43 and Prp22, as well as of the related RNA helicase MLE, in complex with RNA have contributed to a better understanding of how RNA binding and processivity might be achieved in this helicase family. In order to shed light onto the divergent manner of function of Prp2, an N-terminally truncated construct of Chaetomium thermophilum Prp2 was crystallized in the presence of ADP-BeF3- and a poly-U12 RNA. The refined structure revealed a virtually identical conformation of the helicase core compared with the ADP-BeF3-- and RNA-bound structure of Prp43, and only a minor shift of the C-terminal domains. However, Prp2 and Prp43 differ in the hook-loop and a loop of the helix-bundle domain, which interacts with the hook-loop and evokes a different RNA conformation immediately after the 3' stack. On replacing these loop residues in Prp43 by the Prp2 sequence, the unwinding activity of Prp43 was abolished. Furthermore, a putative exit tunnel for the γ-phosphate after ATP hydrolysis could be identified in one of the Prp2 structures.