Structure, mechanism, and inhibition of Hedgehog acyltransferase.

The Sonic Hedgehog (SHH) morphogen pathway is fundamental for embryonic development and stem cell maintenance and is implicated in various cancers. A key step in signaling is transfer of a palmitate group to the SHH N terminus, catalyzed by the multi-pass transmembrane enzyme Hedgehog acyltransferase (HHAT). We present the high-resolution cryo-EM structure of HHAT bound to substrate analog palmityl-coenzyme A and a SHH-mimetic megabody, revealing a heme group bound to HHAT that is essential for HHAT function. A structure of HHAT bound to potent small-molecule inhibitor IMP-1575 revealed conformational changes in the active site that occlude substrate binding. Our multidisciplinary analysis provides a detailed view of the mechanism by which HHAT adapts the membrane environment to transfer an acyl chain across the endoplasmic reticulum membrane. This structure of a membrane-bound O-acyltransferase (MBOAT) superfamily member provides a blueprint for other protein-substrate MBOATs and a template for future drug discovery.

Molecular mechanism of anion permeation through aquaporin 6.

Aquaporins (AQPs) are recognized as transmembrane water channels that facilitate selective water permeation through their monomeric pores. Among the AQP family, AQP6 has an intriguing characteristic as an anion channel, which is allosterically controlled by pH conditions and is eliminated by a single amino acid mutation. However, the molecular mechanism of anion permeation through AQP6 remains unclear. Using molecular dynamics simulations in the presence of a transmembrane voltage utilizing an ion concentration gradient, we show that chloride ions permeate through the pore corresponding to the central axis of the AQP6 homotetramer. Under low pH conditions, a subtle opening of the hydrophobic selectivity filter (SF), located near the extracellular part of the central pore, becomes wetted and enables anion permeation. Our simulations also indicate that a single mutation (N63G) in human AQP6, located at the central pore, significantly reduces anion conduction, consistent with experimental data. Moreover, we demonstrate that the pH-sensing mechanism in which the protonation of H184 and H189 under low pH conditions allosterically triggers the gating of the SF region. These results suggest a unique pH-dependent allosteric anion permeation mechanism in AQP6 and could clarify the role of the central pore in some of the AQP tetramers.

Cryo-EM reveals two distinct serotonin-bound conformations of full-length 5-HT3A receptor.

The 5-HT3A serotonin receptor1, a cationic pentameric ligand-gated ion channel (pLGIC), is the clinical target for management of nausea and vomiting associated with radiation and chemotherapies2. Upon binding, serotonin induces a global conformational change that encompasses the ligand-binding extracellular domain (ECD), the transmembrane domain (TMD) and the intracellular domain (ICD), the molecular details of which are unclear. Here we present two serotonin-bound structures of the full-length 5-HT3A receptor in distinct conformations at 3.32 Å and 3.89 Å resolution that reveal the mechanism underlying channel activation. In comparison to the apo 5-HT3A receptor, serotonin-bound states underwent a large twisting motion in the ECD and TMD, leading to the opening of a 165 Å permeation pathway. Notably, this motion results in the creation of lateral portals for ion permeation at the interface of the TMD and ICD. Combined with molecular dynamics simulations, these structures provide novel insights into conformational coupling across domains and functional modulation.

PtdIns(4,5)P2 stabilizes active states of GPCRs and enhances selectivity of G-protein coupling.

G-protein-coupled receptors (GPCRs) are involved in many physiological processes and are therefore key drug targets1. Although detailed structural information is available for GPCRs, the effects of lipids on the receptors, and on downstream coupling of GPCRs to G proteins are largely unknown. Here we use native mass spectrometry to identify endogenous lipids bound to three class A GPCRs. We observed preferential binding of phosphatidylinositol-4,5-bisphosphate (PtdIns(4,5)P2) over related lipids and confirm that the intracellular surface of the receptors contain hotspots for PtdIns(4,5)P2 binding. Endogenous lipids were also observed bound directly to the trimeric Gαsßγ protein complex of the adenosine A2A receptor (A2AR) in the gas phase. Using engineered Gα subunits (mini-Gαs, mini-Gαi and mini-Gα12)2, we demonstrate that the complex of mini-Gαs with the ß1 adrenergic receptor (ß1AR) is stabilized by the binding of two PtdIns(4,5)P2 molecules. By contrast, PtdIns(4,5)P2 does not stabilize coupling between ß1AR and other Gα subunits (mini-Gαi or mini-Gα12) or a high-affinity nanobody. Other endogenous lipids that bind to these receptors have no effect on coupling, highlighting the specificity of PtdIns(4,5)P2. Calculations of potential of mean force and increased GTP turnover by the activated neurotensin receptor when coupled to trimeric Gαißγ complex in the presence of PtdIns(4,5)P2 provide further evidence for a specific effect of PtdIns(4,5)P2 on coupling. We identify key residues on cognate Gα subunits through which PtdIns(4,5)P2 forms bridging interactions with basic residues on class A GPCRs. These modulating effects of lipids on receptors suggest consequences for understanding function, G-protein selectivity and drug targeting of class A GPCRs.

Influence of electronic polarization on the binding of anions to a chloride-pumping rhodopsin.

The functional properties of some biological ion channels and membrane transport proteins are proposed to exploit anion-hydrophobic interactions. Here, we investigate a chloride-pumping rhodopsin as an example of a membrane protein known to contain a defined anion binding site composed predominantly of hydrophobic residues. Using molecular dynamics simulations, we explore Cl- binding to this hydrophobic site and compare the dynamics arising when electronic polarization is neglected (CHARMM36 [c36] fixed-charge force field), included implicitly (via the prosECCo force field), or included explicitly (through the polarizable force field, AMOEBA). Free energy landscapes of Cl- moving out of the binding site and into bulk solution demonstrate that the inclusion of polarization results in stronger ion binding and a second metastable binding site in chloride-pumping rhodopsin. Simulations focused on this hydrophobic binding site also indicate longer binding durations and closer ion proximity when polarization is included. Furthermore, simulations reveal that Cl- within this binding site interacts with an adjacent loop to facilitate rebinding events that are not observed when polarization is neglected. These results demonstrate how the inclusion of polarization can influence the behavior of anions within protein binding sites and can yield results comparable with more accurate and computationally demanding methods.

Lipid-Dependent Regulation of Ion Channels and G Protein-Coupled Receptors: Insights from Structures and Simulations.

Ion channels and G protein-coupled receptors (GPCRs) are regulated by lipids in their membrane environment. Structural studies combined with biophysical and molecular simulation investigations reveal interaction sites for specific lipids on membrane protein structures. For K channels, PIP2 plays a key role in regulating Kv and Kir channels. Likewise, several recent cryo-EM structures of TRP channels have revealed bound lipids, including PIP2 and cholesterol. Among the pentameric ligand-gated ion channel family, structural and biophysical studies suggest the M4 TM helix may act as a lipid sensor, e.g., forming part of the binding sites for neurosteroids on the GABAA receptor. Structures of GPCRs have revealed multiple cholesterol sites, which may modulate both receptor dynamics and receptor oligomerization. PIP2 also interacts with GPCRs and may modulate their interactions with G proteins. Overall, it is evident that multiple lipid binding sites exist on channels and receptors that modulate their function allosterically and are potential druggable sites.

Defining how multiple lipid species interact with inward rectifier potassium (Kir2) channels.

Protein-lipid interactions are a key element of the function of many integral membrane proteins. These potential interactions should be considered alongside the complexity and diversity of membrane lipid composition. Inward rectifier potassium channel (Kir) Kir2.2 has multiple interactions with plasma membrane lipids: Phosphatidylinositol (4, 5)-bisphosphate (PIP2) activates the channel; a secondary anionic lipid site has been identified, which augments the activation by PIP2; and cholesterol inhibits the channel. Molecular dynamics simulations are used to characterize in molecular detail the protein-lipid interactions of Kir2.2 in a model of the complex plasma membrane. Kir2.2 has been simulated with multiple, functionally important lipid species. From our simulations we show that PIP2 interacts most tightly at the crystallographic interaction sites, outcompeting other lipid species at this site. Phosphatidylserine (PS) interacts at the previously identified secondary anionic lipid interaction site, in a PIP2 concentration-dependent manner. There is interplay between these anionic lipids: PS interactions are diminished when PIP2 is not present in the membrane, underlining the need to consider multiple lipid species when investigating protein-lipid interactions.

The MscS-like channel YnaI has a gating mechanism based on flexible pore helices.

The mechanosensitive channel of small conductance (MscS) is the prototype of an evolutionarily diversified large family that fine-tunes osmoregulation but is likely to fulfill additional functions. Escherichia coli has six osmoprotective paralogs with different numbers of transmembrane helices. These helices are important for gating and sensing in MscS but the role of the additional helices in the paralogs is not understood. The medium-sized channel YnaI was extracted and delivered in native nanodiscs in closed-like and open-like conformations using the copolymer diisobutylene/maleic acid (DIBMA) for structural studies. Here we show by electron cryomicroscopy that YnaI has an extended sensor paddle that during gating relocates relative to the pore concomitant with bending of a GGxGG motif in the pore helices. YnaI is the only one of the six paralogs that has this GGxGG motif allowing the sensor paddle to move outward. Access to the pore is through a vestibule on the cytosolic side that is fenestrated by side portals. In YnaI, these portals are obstructed by aromatic side chains but are still fully hydrated and thus support conductance. For comparison with large-sized channels, we determined the structure of YbiO, which showed larger portals and a wider pore with no GGxGG motif. Further in silico comparison of MscS, YnaI, and YbiO highlighted differences in the hydrophobicity and wettability of their pores and vestibule interiors. Thus, MscS-like channels of different sizes have a common core architecture but show different gating mechanisms and fine-tuned conductive properties.

Influence of effective polarization on ion and water interactions within a biomimetic nanopore.

Interactions between ions and water at hydrophobic interfaces within ion channels and nanopores are suggested to play a key role in the movement of ions across biological membranes. Previous molecular-dynamics simulations have shown that anion affinity for aqueous/hydrophobic interfaces can be markedly influenced by including polarization effects through an electronic continuum correction. Here, we designed a model biomimetic nanopore to imitate the polar pore openings and hydrophobic gating regions found in pentameric ligand-gated ion channels. Molecular-dynamics simulations were then performed using both a non-polarizable force field and the electronic-continuum-correction method to investigate the behavior of water, Na+, and Cl- ions confined within the hydrophobic region of the nanopore. Number-density distributions revealed preferential Cl- adsorption to the hydrophobic pore walls, with this interfacial layer largely devoid of Na+. Free-energy profiles for Na+ and Cl- permeating the pore also display an energy-barrier reduction associated with the localization of Cl- to this hydrophobic interface, and the hydration-number profiles reflect a corresponding reduction in the first hydration shell of Cl-. Crucially, these ion effects were only observed through inclusion of effective polarization, which therefore suggests that polarizability may be essential for an accurate description for the behavior of ions and water within hydrophobic nanoscale pores, especially those that conduct Cl-.

Membrane-binding mechanism of the EEA1 FYVE domain revealed by multi-scale molecular dynamics simulations.

Early Endosomal Antigen 1 (EEA1) is a key protein in endosomal trafficking and is implicated in both autoimmune and neurological diseases. The C-terminal FYVE domain of EEA1 binds endosomal membranes, which contain phosphatidylinositol-3-phosphate (PI(3)P). Although it is known that FYVE binds PI(3)P specifically, it has not previously been described of how FYVE attaches and binds to endosomal membranes. In this study, we employed both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations to determine how FYVE binds to PI(3)P-containing membranes. CG-MD showed that the dominant membrane binding mode resembles the crystal structure of EEA1 FYVE domain in complex with inositol-1,3-diphospate (PDB ID 1JOC). FYVE, which is a homodimer, binds the membrane via a hinge mechanism, where the C-terminus of one monomer first attaches to the membrane, followed by the C-terminus of the other monomer. The estimated total binding energy is ~70 kJ/mol, of which 50-60 kJ/mol stems from specific PI(3)P-interactions. By AT-MD, we could partition the binding mode into two types: (i) adhesion by electrostatic FYVE-PI(3)P interaction, and (ii) insertion of amphipathic loops. The AT simulations also demonstrated flexibility within the FYVE homodimer between the C-terminal heads and coiled-coil stem. This leads to a dynamic model whereby the 200 nm long coiled coil attached to the FYVE domain dimer can amplify local hinge-bending motions such that the Rab5-binding domain at the other end of the coiled coil can explore an area of 0.1 µm2 in the search for a second endosome with which to interact.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective.

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and ß-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

Structural basis of Smoothened regulation by its extracellular domains.

Developmental signals of the Hedgehog (Hh) and Wnt families are transduced across the membrane by Frizzledclass G-protein-coupled receptors (GPCRs) composed of both a heptahelical transmembrane domain (TMD) and an extracellular cysteine-rich domain (CRD). How the large extracellular domains of GPCRs regulate signalling by the TMD is unknown. We present crystal structures of the Hh signal transducer and oncoprotein Smoothened, a GPCR that contains two distinct ligand-binding sites: one in its TMD and one in the CRD. The CRD is stacked a top the TMD, separated by an intervening wedge-like linker domain. Structure-guided mutations show that the interface between the CRD, linker domain and TMD stabilizes the inactive state of Smoothened. Unexpectedly, we find a cholesterol molecule bound to Smoothened in the CRD binding site. Mutations predicted to prevent cholesterol binding impair the ability of Smoothened to transmit native Hh signals. Binding of a clinically used antagonist, vismodegib, to the TMD induces a conformational change that is propagated to the CRD, resulting in loss of cholesterol from the CRD-linker domain-TMD interface. Our results clarify the structural mechanism by which the activity of a GPCR is controlled by ligand-regulated interactions between its extracellular and transmembrane domains.

A heuristic derived from analysis of the ion channel structural proteome permits the rapid identification of hydrophobic gates.

Ion channel proteins control ionic flux across biological membranes through conformational changes in their transmembrane pores. An exponentially increasing number of channel structures captured in different conformational states are now being determined; however, these newly resolved structures are commonly classified as either open or closed based solely on the physical dimensions of their pore, and it is now known that more accurate annotation of their conductive state requires additional assessment of the effect of pore hydrophobicity. A narrow hydrophobic gate region may disfavor liquid-phase water, leading to local dewetting, which will form an energetic barrier to water and ion permeation without steric occlusion of the pore. Here we quantify the combined influence of radius and hydrophobicity on pore dewetting by applying molecular dynamics simulations and machine learning to nearly 200 ion channel structures. This allows us to propose a simple simulation-free heuristic model that rapidly and accurately predicts the presence of hydrophobic gates. This not only enables the functional annotation of new channel structures as soon as they are determined, but also may facilitate the design of novel nanopores controlled by hydrophobic gates.

The morphogen Sonic hedgehog inhibits its receptor Patched by a pincer grasp mechanism.

Hedgehog (HH) ligands, classical morphogens that pattern embryonic tissues in all animals, are covalently coupled to two lipids-a palmitoyl group at the N terminus and a cholesteroyl group at the C terminus. While the palmitoyl group binds and inactivates Patched 1 (PTCH1), the main receptor for HH ligands, the function of the cholesterol modification has remained mysterious. Using structural and biochemical studies, along with reassessment of previous cryo-electron microscopy structures, we find that the C-terminal cholesterol attached to Sonic hedgehog (Shh) binds the first extracellular domain of PTCH1 and promotes its inactivation, thus triggering HH signaling. Molecular dynamics simulations show that this interaction leads to the closure of a tunnel through PTCH1 that serves as the putative conduit for sterol transport. Thus, Shh inactivates PTCH1 by grasping its extracellular domain with two lipidic pincers, the N-terminal palmitate and the C-terminal cholesterol, which are both inserted into the PTCH1 protein core.

Atomistic mechanism of transmembrane helix association.

Transmembrane helix association is a fundamental step in the folding of helical membrane proteins. The prototypical example of this association is formation of the glycophorin dimer. While its structure and stability have been well-characterized experimentally, the detailed assembly mechanism is harder to obtain. Here, we use all-atom simulations within phospholipid membrane to study glycophorin association. We find that initial association results in the formation of a non-native intermediate, separated by a significant free energy barrier from the dimer with a native binding interface. We have used transition-path sampling to determine the association mechanism. We find that the mechanism of the initial bimolecular association to form the intermediate state can be mediated by many possible contacts, but seems to be particularly favoured by formation of non-native contacts between the C-termini of the two helices. On the other hand, the contacts which are key to determining progression from the intermediate to the native state are those which define the native binding interface, reminiscent of the role played by native contacts in determining folding of globular proteins. As a check on the simulations, we have computed association and dissociation rates from the transition-path sampling. We obtain results in reasonable accord with available experimental data, after correcting for differences in native state stability. Our results yield an atomistic description of the mechanism for a simple prototype of helical membrane protein folding.

Computational Modeling of Realistic Cell Membranes.

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Supramolecular assemblies underpin turnover of outer membrane proteins in bacteria.

Gram-negative bacteria inhabit a broad range of ecological niches. For Escherichia coli, this includes river water as well as humans and animals, where it can be both a commensal and a pathogen. Intricate regulatory mechanisms ensure that bacteria have the right complement of ß-barrel outer membrane proteins (OMPs) to enable adaptation to a particular habitat. Yet no mechanism is known for replacing OMPs in the outer membrane, an issue that is further confounded by the lack of an energy source and the high stability and abundance of OMPs. Here we uncover the process underpinning OMP turnover in E. coli and show it to be passive and binary in nature, in which old OMPs are displaced to the poles of growing cells as new OMPs take their place. Using fluorescent colicins as OMP-specific probes, in combination with ensemble and single-molecule fluorescence microscopy in vivo and in vitro, as well as molecular dynamics simulations, we established the mechanism for binary OMP partitioning. OMPs clustered to form ∼0.5-µm diameter islands, where their diffusion is restricted by promiscuous interactions with other OMPs. OMP islands were distributed throughout the cell and contained the Bam complex, which catalyses the insertion of OMPs in the outer membrane. However, OMP biogenesis occurred as a gradient that was highest at mid-cell but largely absent at cell poles. The cumulative effect is to push old OMP islands towards the poles of growing cells, leading to a binary distribution when cells divide. Hence, the outer membrane of a Gram-negative bacterium is a spatially and temporally organized structure, and this organization lies at the heart of how OMPs are turned over in the membrane.

Influence of water models on water movement through AQP1.

Water diffusion through membrane proteins is a key aspect of cellular function. Essential processes of cellular metabolism are driven by osmotic pressure, which depends on water channels. Membrane proteins such as aquaporins (AQPs) are responsible for enabling water permeation through the cell membrane. AQPs are highly selective, allowing only water and relatively small polar molecules to cross the membrane. Experimentally, estimation of water flux through membrane proteins is still a challenge, and hence, accurate simulations of water permeation are of particular importance. We present a numerical study of water diffusion through AQP1 comparing three water models: TIP3P, OPC, and TIP4P/2005. Bulk diffusion, diffusion permeability, and osmotic permeability are computed and compared among all models. The results show that there are significant differences between TIP3P (a particularly widespread model for simulations of biological systems) and the more recently developed TIP4P/2005 and OPC models. We demonstrate that OPC and TIP4P/2005 reproduce protein-water interactions and dynamics in very good agreement with experimental data. From this study, we find that the choice of the water model has a significant effect on the computed water dynamics as well as its molecular behavior within a biological nanopore.

The MemProtMD database: a resource for membrane-embedded protein structures and their lipid interactions.

Integral membrane proteins fulfil important roles in many crucial biological processes, including cell signalling, molecular transport and bioenergetic processes. Advancements in experimental techniques are revealing high resolution structures for an increasing number of membrane proteins. Yet, these structures are rarely resolved in complex with membrane lipids. In 2015, the MemProtMD pipeline was developed to allow the automated lipid bilayer assembly around new membrane protein structures, released from the Protein Data Bank (PDB). To make these data available to the scientific community, a web database (http://memprotmd.bioch.ox.ac.uk) has been developed. Simulations and the results of subsequent analysis can be viewed using a web browser, including interactive 3D visualizations of the assembled bilayer and 2D visualizations of lipid contact data and membrane protein topology. In addition, ensemble analyses are performed to detail conserved lipid interaction information across proteins, families and for the entire database of 3506 PDB entries. Proteins may be searched using keywords, PDB or Uniprot identifier, or browsed using classification systems, such as Pfam, Gene Ontology annotation, mpstruc or the Transporter Classification Database. All files required to run further molecular simulations of proteins in the database are provided.

The Glycosphingolipid GM3 Modulates Conformational Dynamics of the Glucagon Receptor.

The extracellular domain (ECD) of class B1 G-protein-coupled receptors (GPCRs) plays a central role in signal transduction and is uniquely positioned to sense both the extracellular and membrane environments. Although recent studies suggest a role for membrane lipids in the modulation of class A and class F GPCR signaling properties, little is known about the effect of lipids on class B1 receptors. In this study, we employed multiscale molecular dynamics simulations to access the dynamics of the glucagon receptor (GCGR) ECD in the presence of native-like membrane bilayers. Simulations showed that the ECD could move about a hinge region formed by residues Q122-E126 to adopt both closed and open conformations relative to the transmembrane domain. ECD movements were modulated by binding of the glycosphingolipid GM3. These large-scale fluctuations in ECD conformation may affect the ligand binding and receptor activation properties. We also identify a unique phosphatidylinositol (4,5)-bisphosphate (PIP2) interaction profile near intracellular loop (ICL) 2/TM3 at the G-protein-coupling interface, suggesting a mechanism of engaging G-proteins that may have a distinct dependence on PIP2 compared with class A GPCRs. Given the structural conservation of class B1 GPCRs, the modulatory effects of GM3 and PIP2 on GCGR may be conserved across these receptors, offering new insights into potential therapeutic targeting.