Integrated NMR/Molecular Dynamics Determination of the Ensemble Conformation of a Thermodynamically Stable CUUG RNA Tetraloop.

Both experimental and theoretical structure determinations of RNAs have remained challenging due to the intrinsic dynamics of RNAs. We report here an integrated nuclear magnetic resonance/molecular dynamics (NMR/MD) structure determination approach to describe the dynamic structure of the CUUG tetraloop. We show that the tetraloop undergoes substantial dynamics, leading to averaging of the experimental data. These dynamics are particularly linked to the temperature-dependent presence of a hydrogen bond within the tetraloop. Interpreting the NMR data by a single structure represents the low-temperature structure well but fails to capture all conformational states occurring at a higher temperature. We integrate MD simulations, starting from structures of CUUG tetraloops within the Protein Data Bank, with an extensive set of NMR data, and provide a structural ensemble that describes the dynamic nature of the tetraloop and the experimental NMR data well. We thus show that one of the most stable and frequently found RNA tetraloops displays substantial dynamics, warranting such an integrated structural approach.

Double Mutant of Chymotrypsin Inhibitor 2 Stabilized through Increased Conformational Entropy.

The conformational heterogeneity of a folded protein can affect not only its function but also stability and folding. We recently discovered and characterized a stabilized double mutant (L49I/I57V) of the protein CI2 and showed that state-of-the-art prediction methods could not predict the increased stability relative to the wild-type protein. Here, we have examined whether changed native-state dynamics, and resulting entropy changes, can explain the stability changes in the double mutant protein, as well as the two single mutant forms. We have combined NMR relaxation measurements of the ps-ns dynamics of amide groups in the backbone and the methyl groups in the side chains with molecular dynamics simulations to quantify the native-state dynamics. The NMR experiments reveal that the mutations have different effects on the conformational flexibility of CI2: a reduction in conformational dynamics (and entropy estimated from this) of the native state of the L49I variant correlates with its decreased stability, while increased dynamics of the I57V and L49I/I57V variants correlates with their increased stability. These findings suggest that explicitly accounting for changes in native-state entropy might be needed to improve the predictions of the effect of mutations on protein stability.

Fitting Force Field Parameters to NMR Relaxation Data.

We present an approach to optimize force field parameters using time-dependent data from NMR relaxation experiments. To do so, we scan parameters in the dihedral angle potential energy terms describing the rotation of the methyl groups in proteins and compare NMR relaxation rates calculated from molecular dynamics simulations with the modified force fields to deuterium relaxation measurements of T4 lysozyme. We find that a small modification of Cγ methyl groups improves the agreement with experiments both for the protein used to optimize the force field and when validating using simulations of CI2 and ubiquitin. We also show that these improvements enable a more effective a posteriori reweighting of the MD trajectories. The resulting force field thus enables more direct comparison between simulations and side-chain NMR relaxation data and makes it possible to construct ensembles that better represent the dynamics of proteins in solution.

Activation and substrate specificity of the human P4-ATPase ATP8B1.

Asymmetric distribution of phospholipids in eukaryotic membranes is essential for cell integrity, signaling pathways, and vesicular trafficking. P4-ATPases, also known as flippases, participate in creating and maintaining this asymmetry through active transport of phospholipids from the exoplasmic to the cytosolic leaflet. Here, we present a total of nine cryo-electron microscopy structures of the human flippase ATP8B1-CDC50A complex at 2.4 to 3.1 Å overall resolution, along with functional and computational studies, addressing the autophosphorylation steps from ATP, substrate recognition and occlusion, as well as a phosphoinositide binding site. We find that the P4-ATPase transport site is occupied by water upon phosphorylation from ATP. Additionally, we identify two different autoinhibited states, a closed and an outward-open conformation. Furthermore, we identify and characterize the PI(3,4,5)P3 binding site of ATP8B1 in an electropositive pocket between transmembrane segments 5, 7, 8, and 10. Our study also highlights the structural basis of a broad lipid specificity of ATP8B1 and adds phosphatidylinositol as a transport substrate for ATP8B1. We report a critical role of the sn-2 ester bond of glycerophospholipids in substrate recognition by ATP8B1 through conserved S403. These findings provide fundamental insights into ATP8B1 catalytic cycle and regulation, and substrate recognition in P4-ATPases.

A distinct mechanism of C-type inactivation in the Kv-like KcsA mutant E71V.

C-type inactivation is of great physiological importance in voltage-activated K+ channels (Kv), but its structural basis remains unresolved. Knowledge about C-type inactivation has been largely deduced from the bacterial K+ channel KcsA, whose selectivity filter constricts under inactivating conditions. However, the filter is highly sensitive to its molecular environment, which is different in Kv channels than in KcsA. In particular, a glutamic acid residue at position 71 along the pore helix in KcsA is substituted by a valine conserved in most Kv channels, suggesting that this side chain is a molecular determinant of function. Here, a combination of X-ray crystallography, solid-state NMR and MD simulations of the E71V KcsA mutant is undertaken to explore inactivation in this Kv-like construct. X-ray and ssNMR data show that the filter of the Kv-like mutant does not constrict under inactivating conditions. Rather, the filter adopts a conformation that is slightly narrowed and rigidified. On the other hand, MD simulations indicate that the constricted conformation can nonetheless be stably established in the mutant channel. Together, these findings suggest that the Kv-like KcsA mutant may be associated with different modes of C-type inactivation, showing that distinct filter environments entail distinct C-type inactivation mechanisms.

Fitting Side-Chain NMR Relaxation Data Using Molecular Simulations.

Proteins display a wealth of dynamical motions that can be probed using both experiments and simulations. We present an approach to integrate side-chain NMR relaxation measurements with molecular dynamics simulations to study the structure and dynamics of these motions. The approach, which we term ABSURDer (average block selection using relaxation data with entropy restraints), can be used to find a set of trajectories that are in agreement with relaxation measurements. We apply the method to deuterium relaxation measurements in T4 lysozyme and show how it can be used to integrate the accuracy of the NMR measurements with the molecular models of protein dynamics afforded by the simulations. We show how fitting of dynamic quantities leads to improved agreement with static properties and highlight areas needed for further improvements of the approach.

Shifts in the selectivity filter dynamics cause modal gating in K+ channels.

Spontaneous activity shifts at constant experimental conditions represent a widespread regulatory mechanism in ion channels. The molecular origins of these modal gating shifts are poorly understood. In the K+ channel KcsA, a multitude of fast activity shifts that emulate the native modal gating behaviour can be triggered by point-mutations in the hydrogen bonding network that controls the selectivity filter. Using solid-state NMR and molecular dynamics simulations in a variety of KcsA mutants, here we show that modal gating shifts in K+ channels are associated with important changes in the channel dynamics that strongly perturb the selectivity filter equilibrium conformation. Furthermore, our study reveals a drastically different motional and conformational selectivity filter landscape in a mutant that mimics voltage-gated K+ channels, which provides a foundation for an improved understanding of eukaryotic K+ channels. Altogether, our results provide a high-resolution perspective on some of the complex functional behaviour of K+ channels.