CD39, NTPDase 1, is attached to the plasma membrane by two transmembrane domains. Why?

Since the identification of CD39 and other members of the e-NTPDase (ecto-nucleoside triphosphate diphosphohydrolase) family as the primary enzymes responsible for cell surface nucleotide hydrolysis, one of their most intriguing features has been their unusual topology. The active site lies in the large extracellular region, but instead of being anchored in the membrane by a single transmembrane domain or lipid link like other ectoenzymes, CD39 has two transmembrane domains, one at each end. In this review we discuss evidence that the structure and dynamics of the transmembrane helices are intricately connected to enzymatic function. Removal of either or both transmembrane domains or disruption of their native state by detergent solubilization reduces activity by 90%, indicating that native function requires both transmembrane domains to be present and in the membrane. Enzymatic and mutational analysis of the native and truncated forms has shown that the active site can exist in distinct functional states characterized by different total activities, substrate specificities, hydrolysis mechanisms, and intermediate ADP release during ATP hydrolysis, depending on the state of the transmembrane domains. Disulfide crosslinking of cysteines introduced within the transmembrane helices revealed that they interact within and between molecules, in particular near the extracellular domain, and that activity depends on their organization. Both helices exhibit a high degree of rotational mobility, and the ability to undergo dynamic motions is required for activity and regulated by substrate binding. Recent reports suggest that membrane composition can regulate NTPDase activity. We propose that mechanical bilayer properties, potentially elasticity, might regulate CD39 by altering the balance between stability and mobility of its transmembrane domains.     

Dynamic motions of CD39 transmembrane domains regulate and are regulated by the enzymatic active site

The two transmembrane domains flanking the active site of CD39 regulate its activity, but little is known about the structural and dynamic features underlying their importance. Here we use a disulfide crosslinking strategy to examine transmembrane helix interactions and dynamics and to correlate these features with activity and substrate binding. We find strong intrasubunit TM1-TM2 interactions, as well as TM1-TM1' and TM2-TM2' interactions between dimer subunits, near the extracellular side of the membrane but only weak interactions near the cytoplasmic end. The specific helix faces that constitute each interface are highly flexible, indicating a significant degree of rotational mobility within the packed structure. Analysis of activity after locking the helices in various orientations via disulfide bonds suggests that not only the arrangement but also the ability of the helices to move relative to each other is crucial for enzyme function. Helix mobility is in turn modulated by substrate binding. These results suggest that rather than playing a static structural role to support an optimal active site conformation, the transmembrane domains undergo dynamic motions that underlie their functional relationship with the active site.     

Structural analysis of a rhomboid family intramembrane protease reveals a gating mechanism for substrate entry

Intramembrane proteolysis regulates diverse biological processes. Cleavage of substrate peptide bonds within the membrane bilayer is catalyzed by integral membrane proteases. Here we report the crystal structure of the transmembrane core domain of GlpG, a rhomboid-family intramembrane serine protease from Escherichia coli. The protein contains six transmembrane helices, with the catalytic Ser201 located at the N terminus of helix alpha4 approximately 10 A below the membrane surface. Access to water molecules is provided by a central cavity that opens to the extracellular region and converges on Ser201. One of the two GlpG molecules in the asymmetric unit has an open conformation at the active site, with the transmembrane helix alpha5 bent away from the rest of the molecule. Structural analysis suggests that substrate entry to the active site is probably gated by the movement of helix alpha5.     

Sequence-dependent backbone dynamics of a viral fusogen transmembrane helix

The transmembrane domains of membrane fusogenic proteins are known to contribute to lipid bilayer mixing as indicated by mutational studies and functional reconstitution of peptide mimics. Here, we demonstrate that mutations of a GxxxG motif or of Ile residues, that were previously shown to compromise the fusogenicity of the Vesicular Stomatitis virus G-protein transmembrane helix, reduce its backbone dynamics as determined by deuterium/hydrogen-exchange kinetics. Thus, the backbone dynamics of these helices may be linked to their fusogenicity which is consistent with the known over-representation of Gly and Ile in viral fusogen transmembrane helices. The transmembrane domains of membrane fusogenic proteins are known to contribute to lipid bilayer mixing. Our present results demonstrate that mutations of certain residues, that were previously shown to compromise the fusogenicity of the Vesicular Stomatitis virus G-protein transmembrane helix, reduce its backbone dynamics. Thus, the data suggest a relationship between sequence, backbone dynamics, and fusogenicity of transmembrane segments of viral fusogenic proteins.     

Side-Chain to Main-Chain Hydrogen Bonding Controls the Intrinsic Backbone Dynamics of the Amyloid Precursor Protein Transmembrane Helix

Many transmembrane helices contain serine and/or threonine residues whose side chains form intrahelical H-bonds with upstream carbonyl oxygens. Here, we investigated the impact of threonine side-chain/main-chain backbonding on the backbone dynamics of the amyloid precursor protein transmembrane helix. This helix consists of a N-terminal dimerization region and a C-terminal cleavage region, which is processed by γ-secretase to a series of products. Threonine mutations within this transmembrane helix are known to alter the cleavage pattern, which can lead to early-onset Alzheimer’s disease. Circular dichroism spectroscopy and amide exchange experiments of synthetic transmembrane domain peptides reveal that mutating threonine enhances the flexibility of this helix. Molecular dynamics simulations show that the mutations reduce intrahelical amide H-bonding and H-bond lifetimes. In addition, the removal of side-chain/main-chain backbonding distorts the helix, which alters bending and rotation at a diglycine hinge connecting the dimerization and cleavage regions. We propose that the backbone dynamics of the substrate profoundly affects the way by which the substrate is presented to the catalytic site within the enzyme. Changing this conformational flexibility may thus change the pattern of proteolytic processing.     

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.     

Transmembrane helices of membrane proteins may flex to satisfy hydrophobic mismatch

A novel mechanism for membrane modulation of transmembrane protein structure, and consequently function, is suggested in which mismatch between the hydrophobic surface of the protein and the hydrophobic interior of the lipid bilayer induces a flexing or bending of a transmembrane segment of the protein. Studies on model hydrophobic transmembrane peptides predict that helices tilt to submerge the hydrophobic surface within the lipid bilayer to satisfy the hydrophobic effect if the helix length exceeds the bilayer width. The hydrophobic surface of transmembrane helix 1 (TM1) of lactose permease, LacY, is accessible to the bilayer, and too long to be accommodated in the hydrophobic portion of a typical lipid bilayer if oriented perpendicular to the membrane surface. Hence, nuclear magnetic resonance (NMR) data and molecular dynamics simulations show that TM1 from LacY may flex as well as tilt to satisfy the hydrophobic mismatch with the bilayer. In an analogous study of the hydrophobic mismatch of TM7 of bovine rhodopsin, similar flexing of the transmembrane segment near the conserved NPxxY sequence is observed. As a control, NMR data on TM5 of lacY, which is much shorter than TM1, show that TM5 is likely to tilt, but not flex, consistent with the close match between the extent of hydrophobic surface of the peptide and the hydrophobic thickness of the bilayer. These data suggest mechanisms by which the lipid bilayer in which the protein is embedded modulates conformation, and thus function, of integral membrane proteins through interactions with the hydrophobic transmembrane helices.     

Membrane phospholipid bilayer as a determinant of monoacylglycerol lipase kinetic profile and conformational repertoire

The membrane‐associated serine hydrolase, monoacylglycerol lipase (MGL), is a well‐recognized therapeutic target that regulates endocannabinoid signaling. Crystallographic studies, while providing structural information about static MGL states, offer no direct experimental insight into the impact of MGL's membrane association upon its structure–function landscape. We report application of phospholipid bilayer nanodiscs as biomembrane models with which to evaluate the effect of a membrane system on the catalytic properties and conformational dynamics of human MGL (hMGL). Anionic and charge‐neutral phospholipid bilayer nanodiscs enhanced hMGL's kinetic properties [apparent maximum velocity (V_max) and substrate affinity (K_m)]. Hydrogen exchange mass spectrometry (HX MS) was used as a conformational analysis method to profile experimentally the extent of hMGL–nanodisc interaction and its impact upon hMGL structure. We provide evidence that significant regions of hMGL lid‐domain helix α4 and neighboring helix α6 interact with the nanodisc phospholipid bilayer, anchoring hMGL in a more open conformation to facilitate ligand access to the enzyme's substrate‐binding channel. Covalent modification of membrane‐associated hMGL by the irreversible carbamate inhibitor, AM6580, shielded the active site region, but did not increase solvent exposure of the lid domain, suggesting that the inactive, carbamylated enzyme remains intact and membrane associated. Molecular dynamics simulations generated conformational models congruent with the open, membrane‐associated topology of active and inhibited, covalently‐modified hMGL. Our data indicate that hMGL interaction with a phospholipid membrane bilayer induces regional changes in the enzyme's conformation that favor its recruiting lipophilic substrate/inhibitor from membrane stores to the active site via the lid, resulting in enhanced hMGL catalytic activity and substrate affinity.     

Side-Chain to Main-Chain Hydrogen Bonding Controls the Intrinsic Backbone Dynamics of the Amyloid Precursor Protein Transmembrane Helix

Many transmembrane helices contain serine and/or threonine residues whose side chains form intrahelical H-bonds with upstream carbonyl oxygens. Here, we investigated the impact of threonine side-chain/main-chain backbonding on the backbone dynamics of the amyloid precursor protein transmembrane helix. This helix consists of a N-terminal dimerization region and a C-terminal cleavage region, which is processed by γ-secretase to a series of products. Threonine mutations within this transmembrane helix are known to alter the cleavage pattern, which can lead to early-onset Alzheimer’s disease. Circular dichroism spectroscopy and amide exchange experiments of synthetic transmembrane domain peptides reveal that mutating threonine enhances the flexibility of this helix. Molecular dynamics simulations show that the mutations reduce intrahelical amide H-bonding and H-bond lifetimes. In addition, the removal of side-chain/main-chain backbonding distorts the helix, which alters bending and rotation at a diglycine hinge connecting the dimerization and cleavage regions. We propose that the backbone dynamics of the substrate profoundly affects the way by which the substrate is presented to the catalytic site within the enzyme. Changing this conformational flexibility may thus change the pattern of proteolytic processing.     

Transmembrane helix orientation and dynamics: insights from ensemble dynamics with solid-state NMR observables

As the major component of membrane proteins, transmembrane helices embedded in anisotropic bilayer environments adopt preferential orientations that are characteristic or related to their functional states. Recent developments in solid-state nuclear magnetic resonance (SSNMR) spectroscopy have made it possible to measure NMR observables that can be used to determine such orientations in a native bilayer environment. A quasistatic single conformer model is frequently used to interpret the SSNMR observables, but important motional information can be missing or misinterpreted in the model. In this work, we have investigated the orientation of the single-pass transmembrane domain of viral protein ”u“ (VpuTM) from HIV-1 by determining an ensemble of structures using multiple conformer models based on the SSNMR ensemble dynamics technique. The resulting structure ensemble shows significantly larger orientational fluctuations while the ensemble-averaged orientation is compatible with the orientation based on the quasistatic model. This observation is further corroborated by comparison with the VpuTM orientation from comparative molecular dynamics simulations in explicit bilayer membranes. SSNMR ensemble dynamics not only reveals the importance of transmembrane helix dynamics in interpretation of SSNMR observables, but also provides a means to simultaneously extract both transmembrane helix orientation and dynamics information from the SSNMR measurements.     

Molecular dynamics simulations of the ErbB-2 transmembrane domain within an explicit membrane environment: comparison with vacuum simulations

Two 500-ps molecular dynamics simulations performed on the single transmembrane domain of the ErbB-2 tyrosine kinase receptor immersed in a fully solvated dilauroylphosphatidyl-ethanolamine bilayer (DLPE) are compared to vacuum simulations. One membrane simulation shows that the initial alpha helix undergoes a local pi helix conversion in the peptide part embedded in the membrane core similar to that found in simulation vacuum. Lipid/water/peptide interaction analysis shows that in the helix core, the intramolecular peptide interactions are largely dominant compared to the interactions with water and lipids whereas the helix extremities are much more sensitive to these interactions at the membrane interfaces. Our results suggest that simulations in a lipid environment are required to understand the dynamics of transmembrane helices, but can be reasonably supplemented by in vacuo simulations to explore rapidly its conformational space and to describe the internal deformation of the hydrophobic core.     

Conformational dynamics of M2 helices in KirBac channels: helix flexibility in relation to gating via molecular dynamics simulations

KirBac1.1 and 3.1 are bacterial homologues of mammalian inward rectifier K channels. We have performed extended molecular dynamics simulations (five simulations, each of >20 ns duration) of the transmembrane domain of KirBac in two membrane environments, a palmitoyl oleoyl phosphatidylcholine bilayer and an octane slab. Analysis of these simulations has focused on the conformational dynamics of the pore-lining M2 helices, which form the cytoplasmic hydrophobic gate of the channel. Principal components analysis reveals bending of M2, with a molecular hinge at the conserved glycine (Gly134 in KirBac1.1, Gly120 in KirBac3.1). More detailed analysis reveals a dimer-of-dimers type motion. The first two eigenvectors describing the motions of M2 correspond to helix kink and swivel motions. The conformational flexibility of M2 seen in these simulations correlates with differences in M2 conformation between that seen in the X-ray structures of closed channels (KcsA and KirBac) in which the helix is undistorted, and in open channels (e.g. MthK) in which the M2 helix is kinked. Thus, the simulations, albeit on a time scale substantially shorter than that required for channel gating, suggest a gating model in which the intrinsic flexibility of M2 about a molecular hinge is coupled to conformational transitions of an intracellular 'gatekeeper' domain, the latter changing conformation in response to ligand binding.     

Structure, Dynamics, and Ion Conductance of the Phospholamban Pentamer

A 52-residue membrane protein, phospholamban (PLN) is an inhibitor of an adenosine-5′-triphosphate-driven calcium pump, the Ca²⁺-ATPase. Although the inhibition of Ca²⁺-ATPase involves PLN monomers, in a lipid bilayer membrane, PLN monomers form stable pentamers of unknown biological function. The recent NMR structure of a PLN pentamer depicts cytoplasmic helices extending normal to the bilayer in what is known as the bellflower conformation. The structure shows transmembrane helices forming a hydrophobic pore 4 Å in diameter, which is reminiscent of earlier reports of possible ion conductance through PLN pentamers. However, recent FRET measurements suggested an alternative structure for the PLN pentamer, known as the pinwheel model, which features a narrower transmembrane pore and cytoplasmic helices that lie against the bilayer. Here, we report on structural dynamics and conductance properties of the PLN pentamers from all-atom (AA) and coarse-grained (CG) molecular dynamics simulations. Our AA simulations of the bellflower model demonstrate that in a lipid bilayer membrane or a detergent micelle, the cytoplasmic helices undergo large structural fluctuations, whereas the transmembrane pore shrinks and becomes asymmetric. Similar asymmetry of the transmembrane region was observed in the AA simulations of the pinwheel model; the cytoplasmic helices remained in contact with the bilayer. Using the CG approach, structural dynamics of both models were investigated on a microsecond timescale. The cytoplasmic helices of the CG bellflower model were observed to fall against the bilayer, whereas in the CG pinwheel model the conformation of the cytoplasmic helices remained stable. Using steered molecular dynamics simulations, we investigated the feasibility of ion conductance through the pore of the bellflower model. The resulting approximate potentials of mean force indicate that the PLN pentamer is unlikely to function as an ion channel.     

Conformational Changes in Talin on Binding to Anionic Phospholipid Membranes Facilitate Signaling by Integrin Transmembrane Helices

Integrins are heterodimeric (αβ) cell surface receptors that are activated to a high affinity state by the formation of a complex involving the α/β integrin transmembrane helix dimer, the head domain of talin (a cytoplasmic protein that links integrins to actin), and the membrane. The talin head domain contains four sub-domains (F0, F1, F2 and F3) with a long cationic loop inserted in the F1 domain. Here, we model the binding and interactions of the complete talin head domain with a phospholipid bilayer, using multiscale molecular dynamics simulations. The role of the inserted F1 loop, which is missing from the crystal structure of the talin head, PDB:3IVF, is explored. The results show that the talin head domain binds to the membrane predominantly via cationic regions on the F2 and F3 subdomains and the F1 loop. Upon binding, the intact talin head adopts a novel V-shaped conformation which optimizes its interactions with the membrane. Simulations of the complex of talin with the integrin α/β TM helix dimer in a membrane, show how this complex promotes a rearrangement, and eventual dissociation of, the integrin α and β transmembrane helices. A model for the talin-mediated integrin activation is proposed which describes how the mutual interplay of interactions between transmembrane helices, the cytoplasmic talin protein, and the lipid bilayer promotes integrin inside-out activation.     

Differential backbone dynamics of companion helices in the extended helical coiled‐coil domain of a bacterial chemoreceptor

Cytoplasmic domains of transmembrane bacterial chemoreceptors are largely extended four‐helix coiled coils. Previous observations suggested the domain was structurally dynamic. We probed directly backbone dynamics of this domain of the transmembrane chemoreceptor Tar from Escherichia coli using site‐directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy. Spin labels were positioned on solvent‐exposed helical faces because EPR spectra for such positions reflect primarily polypeptide backbone movements. We acquired spectra for spin‐labeled, intact receptor homodimers solubilized in detergent or inserted into native E. coli lipid bilayers in Nanodiscs, characterizing 16 positions distributed throughout the cytoplasmic domain and on both helices of its helical hairpins, one amino terminal to the membrane‐distal tight turn (N‐helix), and the other carboxyl terminal (C‐helix). Detergent solubilization increased backbone dynamics for much of the domain, suggesting that loss of receptor activities upon solubilization reflects wide‐spread destabilization. For receptors in either condition, we observed an unanticipated difference between the N‐ and C‐helices. For bilayer‐inserted receptors, EPR spectra from sites in the membrane‐distal protein‐interaction region and throughout the C‐helix were typical of well‐structured helices. In contrast, for approximately two‐thirds of the N‐helix, from its origin as the AS‐2 helix of the membrane‐proximal HAMP domain to the beginning of the membrane‐distal protein‐interaction region, spectra had a significantly mobile component, estimated by spectral deconvolution to average approximately 15%. Differential helical dynamics suggests a four‐helix bundle organization with a pair of core scaffold helices and two more dynamic partner helices. This newly observed feature of chemoreceptor structure could be involved in receptor function.     

Nanomechanical properties of MscL α helices: A steered molecular dynamics study

Gating of mechanosensitive (MS) channels is driven by a hierarchical cascade of movements and deformations of transmembrane helices in response to bilayer tension. Determining the intrinsic mechanical properties of the individual transmembrane helices is therefore central to understanding the intricacies of the gating mechanism of MS channels. We used a constant-force steered molecular dynamics (SMD) approach to perform unidirectional pulling tests on all the helices of MscL in M. tuberculosis and E. coli homologs. Using this method, we could overcome the issues encountered with the commonly used constant-velocity SMD simulations, such as low mechanical stability of the helix during stretching and high dependency of the elastic properties on the pulling rate. We estimated Young's moduli of the α-helices of MscL to vary between 0.2 and 12.5 GPa with TM2 helix being the stiffest. We also studied the effect of water on the properties of the pore-lining TM1 helix. In the absence of water, this helix exhibited a much stiffer response. By monitoring the number of hydrogen bonds, it appears that water acts like a ‘lubricant’ (softener) during TM1 helix elongation. These data shed light on another physical aspect underlying hydrophobic gating of MS channels, in particular MscL.     

Orientational preferences of neighboring helices can drive ER insertion of a marginally hydrophobic transmembrane helix

α-helical integral membrane proteins critically depend on the correct insertion of their transmembrane α helices into the lipid bilayer for proper folding, yet a surprisingly large fraction of the transmembrane α helices in multispanning integral membrane proteins are not sufficiently hydrophobic to insert into the target membrane by themselves. How can such marginally hydrophobic segments nevertheless form transmembrane helices in the folded structure? Here, we show that a transmembrane helix with a strong orientational preference (N(cyt)-C(lum) or N(lum)-C(cyt)) can both increase and decrease the hydrophobicity threshold for membrane insertion of a neighboring, marginally hydrophobic helix. This effect helps explain the "missing hydrophobicity" in polytopic membrane proteins.     

High cationic charge and bilayer interface-binding helices in a regulatory lipid glycosyltransferase

In the prokaryote Acholeplasma laidlawii, membrane bilayer properties are sensed and regulated by two interface glycosyltransferases (GTs), synthesizing major nonbilayer- (alMGS GT) and bilayer-prone glucolipids. These enzymes are of similar structure, as many soluble GTs, but are sensitive to lipid charge and curvature stress properties. Multivariate and bioinformatic sequence analyses show that such interface enzymes, in relation to soluble ones of similar fold, are characterized by high cationic charge, certain distances between small and cationic amino acids, and by amphipathic helices. Varying surface contents of Lys/Arg pairs and Trp indicate different membrane-binding subclasses. A predicted potential (cationic) binding helix from alMGS was structurally verified by solution NMR and CD. The helix conformation was induced by a zwitterionic as well as anionic lipid environment, and the peptide was confined to the bilayer interface. Bilayer affinity of the peptide, analyzed by surface plasmon resonance, was higher than that for soluble membrane-seeking proteins/peptides and rose with anionic lipid content. Interface intercalation was supported by phase equilibria in membrane lipid mixtures, analyzed by 31P NMR and DSC. An analogous, potentially binding helix has a similar location in the structurally determined Escherichia coli cell wall precursor GT MurG. These two helices have little sequence conservation in alMGS and MurG homologues but maintain their amphipathic character. The evolutionary modification of the alMGS binding helix and its location close to the acceptor substrate site imply a functional importance in enzyme catalysis, potentially providing a mechanism by which glycolipid synthesis will be sensitive to membrane surface charge and intrinsic curvature strain.     

Coarse grain lipid-protein molecular interactions and diffusion with MsbA flippase

Coarse-grained (CG) modeling has proven effective for simulating lipid bilayer dynamics on scales of biological interest. Modeling the dynamics of flexible membrane proteins within the bilayer, on the other hand, poses a considerable challenge due to the complexity of the folding or conformational landscape. In the present work, the multiscale coarse-graining method is applied to atomistic peptide-lipid "soup" simulations to develop a general set of CG protein-lipid interaction potentials. The reduced model was constructed to be compatible with recent solvent-free CG models developed for protein-protein folding and lipid-lipid model bilayer interactions. The utility of the force field was demonstrated by molecular dynamics simulation of the MsbA ABC transporter in a mixed DOPC/DOPE bilayer. An elastic network was parameterized to restrain the MsbA dimer in its open, closed and hydrolysis intermediate conformations and its impact on domain flexibility was examined. Conformational stability enabled long-time dynamics simulation of MsbA freely diffusing in a 25 nm membrane patch. Three-dimensional density analysis revealed that a shell of weakly bound "annular lipids" solvate the membrane accessible surface of MsbA and its internal substrate-binding chamber. The annular lipid binding modes, along with local perturbations in head group structure, are a function of the orientation of grooves formed between transmembrane helices and may influence the alternating access mechanism of substrate entry and translocation.