Automated simulation-based membrane protein refinement into cryo-EM data.

The resolution revolution has increasingly enabled single-particle cryogenic electron microscopy (cryo-EM) reconstructions of previously inaccessible systems, including membrane proteins-a category that constitutes a disproportionate share of drug targets. We present a protocol for using density-guided molecular dynamics simulations to automatically refine atomistic models into membrane protein cryo-EM maps. Using adaptive force density-guided simulations as implemented in the GROMACS molecular dynamics package, we show how automated model refinement of a membrane protein is achieved without the need to manually tune the fitting force ad hoc. We also present selection criteria to choose the best-fit model that balances stereochemistry and goodness of fit. The proposed protocol was used to refine models into a new cryo-EM density of the membrane protein maltoporin, either in a lipid bilayer or detergent micelle, and we found that results do not substantially differ from fitting in solution. Fitted structures satisfied classical model-quality metrics and improved the quality and the model-to-map correlation of the x-ray starting structure. Additionally, the density-guided fitting in combination with generalized orientation-dependent all-atom potential was used to correct the pixel-size estimation of the experimental cryo-EM density map. This work demonstrates the applicability of a straightforward automated approach to fitting membrane protein cryo-EM densities. Such computational approaches promise to facilitate rapid refinement of proteins under different conditions or with various ligands present, including targets in the highly relevant superfamily of membrane proteins.

Biophysical characterization of calcium-binding and modulatory-domain dynamics in a pentameric ligand-gated ion channel.

Pentameric ligand-gated ion channels (pLGICs) perform electrochemical signal transduction in organisms ranging from bacteria to humans. Among the prokaryotic pLGICs, there is architectural diversity involving N-terminal domains (NTDs) not found in eukaryotic relatives, exemplified by the calcium-sensitive channel (DeCLIC) from a Desulfofustis deltaproteobacterium, which has an NTD in addition to the canonical pLGIC structure. Here, we have characterized the structure and dynamics of DeCLIC through cryoelectron microscopy (cryo-EM), small-angle neutron scattering (SANS), and molecular dynamics (MD) simulations. In the presence and absence of calcium, cryo-EM yielded structures with alternative conformations of the calcium-binding site. SANS profiles further revealed conformational diversity at room temperature beyond that observed in static structures, shown through MD to be largely attributable to rigid-body motions of the NTD relative to the protein core, with expanded and asymmetric conformations improving the fit of the SANS data. This work reveals the range of motion available to the DeCLIC NTD and calcium-binding site, expanding the conformational landscape of the pLGIC family. Further, these findings demonstrate the power of combining low-resolution scattering, high-resolution structural, and MD simulation data to elucidate interfacial interactions that are highly conserved in the pLGIC family.

Probing solution structure of the pentameric ligand-gated ion channel GLIC by small-angle neutron scattering.

Pentameric ligand-gated ion channels undergo subtle conformational cycling to control electrochemical signal transduction in many kingdoms of life. Several crystal structures have now been reported in this family, but the functional relevance of such models remains unclear. Here, we used small-angle neutron scattering (SANS) to probe ambient solution-phase properties of the pH-gated bacterial ion channel GLIC under resting and activating conditions. Data collection was optimized by inline paused-flow size-exclusion chromatography, and exchanging into deuterated detergent to hide the micelle contribution. Resting-state GLIC was the best-fit crystal structure to SANS curves, with no evidence for divergent mechanisms. Moreover, enhanced-sampling molecular-dynamics simulations enabled differential modeling in resting versus activating conditions, with the latter corresponding to an intermediate ensemble of both the extracellular and transmembrane domains. This work demonstrates state-dependent changes in a pentameric ion channel by SANS, an increasingly accessible method for macromolecular characterization with the coming generation of neutron sources.

Allosteric potentiation of a ligand-gated ion channel is mediated by access to a deep membrane-facing cavity.

Theories of general anesthesia have shifted in focus from bulk lipid effects to specific interactions with membrane proteins. Target receptors include several subtypes of pentameric ligand-gated ion channels; however, structures of physiologically relevant proteins in this family have yet to define anesthetic binding at high resolution. Recent cocrystal structures of the bacterial protein GLIC provide snapshots of state-dependent binding sites for the common surgical agent propofol (PFL), offering a detailed model system for anesthetic modulation. Here, we combine molecular dynamics and oocyte electrophysiology to reveal differential motion and modulation upon modification of a transmembrane binding site within each GLIC subunit. WT channels exhibited net inhibition by PFL, and a contraction of the cavity away from the pore-lining M2 helix in the absence of drug. Conversely, in GLIC variants exhibiting net PFL potentiation, the cavity was persistently expanded and proximal to M2. Mutations designed to favor this deepened site enabled sensitivity even to subclinical concentrations of PFL, and a uniquely prolonged mode of potentiation evident up to ∼30 min after washout. Dependence of these prolonged effects on exposure time implicated the membrane as a reservoir for a lipid-accessible binding site. However, at the highest measured concentrations, potentiation appeared to be masked by an acute inhibitory effect, consistent with the presence of a discrete, water-accessible site of inhibition. These results support a multisite model of transmembrane allosteric modulation, including a possible link between lipid- and receptor-based theories that could inform the development of new anesthetics.