Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations.

To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χOL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 µs, for a total of 552 µs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.

Lipids Alter Rhodopsin Function via Ligand-like and Solvent-like Interactions.

Rhodopsin, a prototypical G protein-coupled receptor, is a membrane protein that can sense dim light. This highly effective photoreceptor is known to be sensitive to the composition of its lipidic environment, but the molecular mechanisms underlying this fine-tuned modulation of the receptor's function and structural stability are not fully understood. There are two competing hypotheses to explain how this occurs: 1) lipid modulation occurs via solvent-like interactions, where lipid composition controls membrane properties like hydrophobic thickness, which in turn modulate the protein's conformational equilibrium; or 2) protein-lipid interactions are ligand-like, with specific hot spots and long-lived binding events. By analyzing an ensemble of all-atom molecular dynamics simulations of five different states of rhodopsin, we show that a local ordering effect takes place in the membrane upon receptor activation. Likewise, docosahexaenoic acid acyl tails and phosphatidylethanolamine headgroups behave like weak ligands, preferentially binding to the receptor in inactive-like conformations and inducing subtle but significant structural changes.

A nuclear export sequence in GPN-loop GTPase 1, an essential protein for nuclear targeting of RNA polymerase II, is necessary and sufficient for nuclear export.

XAB1/Gpn1 is a GTPase that associates with RNA polymerase II (RNAPII) in a GTP-dependent manner. Although XAB1/Gpn1 is essential for nuclear accumulation of RNAPII, the underlying mechanism is not known. A XAB1/Gpn1-EYFP fluorescent protein, like endogenous XAB1/Gpn1, localized to the cytoplasm but it rapidly accumulated in the cell nucleus in the presence of leptomycin B, a chemical inhibitor of the nuclear transport receptor Crm1. Crm1 recognizes short peptides in substrate proteins called nuclear export sequences (NES). Here, we employed site-directed mutagenesis and fluorescence microscopy to assess the functionality of all six putative NESs in XAB1/Gpn1. Mutating five of the six putative NESs did not alter the cytoplasmic localization of XAB1/Gpn1-EYFP. However, a V302A/L304A double mutant XAB1/Gpn1-EYFP protein was clearly accumulated in the cell nucleus, indicating the disruption of a functional NES. This functional XAB1/Gpn1 NES displays all features present in most common and potent NESs, including, in addition to Φ1-Φ4, a critical fifth hydrophobic amino acid Φ0. Therefore, in human Gpn1 this NES spans amino acids 292-LERLRKDMGSVAL-304. XAB1/Gpn1 NES is remarkably conserved during evolution. XAB1/Gpn1 NES was sufficient for nuclear export activity, as it caused a complete exclusion of EYFP from the cell nucleus. Molecular modeling of XAB1/Gpn1 provided a mechanistic reason for NES selection, as functionality correlated with accessibility, and it also suggested a mechanism for NES inhibition by intramolecular masking. In conclusion, we have identified a highly active, evolutionarily conserved NES in XAB1/Gpn1 that is critical for nucleo-cytoplasmic shuttling and steady-state cytoplasmic localization of XAB1/Gpn1.