Search

The force-from-lipids hypothesis of cellular mechanosensation posits that membrane channels open and close in response to changes in the physical state of the lipid bilayer, induced for example by lateral tension. Here, we investigate the molecular basis for this transduction mechanism by studying the mechanosensitive ion channel MscS from Escherichia coli and its eukaryotic homolog MSL1 from Arabidopsis thaliana. First, we use single-particle cryo-electron microscopy to determine the structure of a novel open conformation of wild-type MscS, stabilized in a thinned lipid nanodisc. Compared with the closed state, the structure shows a reconfiguration of helices TM1, TM2, and TM3a, and widening of the central pore. Based on these structures, we examined how the morphology of the membrane is altered upon gating, using molecular dynamics simulations. The simulations reveal that closed-state MscS causes drastic protrusions in the inner leaflet of the lipid bilayer, both in the absence and presence of lateral tension, and for different lipid compositions. These deformations arise to provide adequate solvation to hydrophobic crevices under the TM1-TM2 hairpin, and clearly reflect a high-energy conformation for the membrane, particularly under tension. Strikingly, these protrusions are largely eradicated upon channel opening. An analogous computational study of open and closed MSL1 recapitulates these findings. The gating equilibrium of MscS channels thus appears to be dictated by opposing conformational preferences, namely those of the lipid membrane and of the protein structure. We propose a membrane deformation model of mechanosensation, which posits that tension shifts the gating equilibrium towards the conductive state not because it alters the mode in which channel and lipids interact, but because it increases the energetic cost of the morphological perturbations in the membrane required by the closed state., Competing Interests: YP, BR, NB, EP No competing interests declared, JF Senior editor, eLife

Transmembrane protein 175 (TMEM175) is an evolutionarily distinct lysosomal cation channel whose mutation is associated with the development of Parkinson's disease. Here, we present a cryoelectron microscopy structure and molecular simulations of TMEM175 bound to 4-aminopyridine (4-AP), the only known small-molecule inhibitor of TMEM175 and a broad K + channel inhibitor, as well as a drug approved by the Food and Drug Administration against multiple sclerosis. The structure shows that 4-AP, whose mode of action had not been previously visualized, binds near the center of the ion conduction pathway, in the open state of the channel. Molecular dynamics simulations reveal that this binding site is near the middle of the transmembrane potential gradient, providing a rationale for the voltage-dependent dissociation of 4-AP from TMEM175. Interestingly, bound 4-AP rapidly switches between three predominant binding poses, stabilized by alternate interaction patterns dictated by the twofold symmetry of the channel. Despite this highly dynamic binding mode, bound 4-AP prevents not only ion permeation but also water flow. Together, these studies provide a framework for the rational design of novel small-molecule inhibitors of TMEM175 that might reveal the role of this channel in human lysosomal physiology both in health and disease.

Structures of the human lysosomal K + channel transmembrane protein 175 (TMEM175) in open and closed states revealed a novel architecture lacking the canonical K + selectivity filter motif present in previously known K + channel structures. A hydrophobic constriction composed of four isoleucine residues was resolved in the pore and proposed to serve as the gate in the closed state, and to confer ion selectivity in the open state. Here, we achieve higher-resolution structures of the open and closed states and employ molecular dynamics simulations to analyze the conducting properties of the putative open state, demonstrating that it is permeable to K + and, to a lesser degree, also Na + . Both cations must dehydrate significantly to penetrate the narrow hydrophobic constriction, but ion flow is assisted by a favorable electrostatic field generated by the protein that spans the length of the pore. The balance of these opposing energetic factors explains why permeation is feasible, and why TMEM175 is selective for K + over Na + , despite the absence of the canonical selectivity filter. Accordingly, mutagenesis experiments reveal an exquisite sensitivity of the channel to perturbations that mitigate the constriction. Together, these data reveal a novel mechanism for selective permeation of ions by TMEM175 that is unlike that of other K + channels., Competing Interests: SO, WZ, JL, HC, MK, RH The authors declare no competing interests, FM, JF No competing interests declared

Voltage-activated potassium (Kv) channels open upon membrane depolarization and proceed to spontaneously inactivate. Inactivation controls neuronal firing rates and serves as a form of short-term memory and is implicated in various human neurological disorders. Here, we use high-resolution cryo-electron microscopy and computer simulations to determine one of the molecular mechanisms underlying this physiologically crucial process. Structures of the activated Shaker Kv channel and of its W434F mutant in lipid bilayers demonstrate that C-type inactivation entails the dilation of the ion selectivity filter and the repositioning of neighboring residues known to be functionally critical. Microsecond-scale molecular dynamics trajectories confirm that these changes inhibit rapid ion permeation through the channel. This long-sought breakthrough establishes how eukaryotic K + channels self-regulate their functional state through the plasticity of their selectivity filters.

S-acylation, also known as palmitoylation, is the most abundant form of protein lipidation in humans. This reversible posttranslational modification, which targets thousands of proteins, is catalyzed by 23 members of the DHHC family of integral membrane enzymes. DHHC enzymes use fatty acyl-CoA as the ubiquitous fatty acyl donor and become autoacylated at a catalytic cysteine; this intermediate subsequently transfers the fatty acyl group to a cysteine in the target protein. Protein S-acylation intersects with almost all areas of human physiology, and several DHHC enzymes are considered as possible therapeutic targets against diseases such as cancer. These efforts would greatly benefit from a detailed understanding of the molecular basis for this crucial enzymatic reaction. Here, we combine X-ray crystallography with all-atom molecular dynamics simulations to elucidate the structure of the precatalytic complex of human DHHC20 in complex with palmitoyl CoA. The resulting structure reveals that the fatty acyl chain inserts into a hydrophobic pocket within the transmembrane spanning region of the protein, whereas the CoA headgroup is recognized by the cytosolic domain through polar and ionic interactions. Biochemical experiments corroborate the predictions from our structural model. We show, using both computational and experimental analyses, that palmitoyl CoA acts as a bivalent ligand where the interaction of the DHHC enzyme with both the fatty acyl chain and the CoA headgroup is important for catalytic chemistry to proceed. This bivalency explains how, in the presence of high concentrations of free CoA under physiological conditions, DHHC enzymes can efficiently use palmitoyl CoA as a substrate for autoacylation., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

Excitatory amino acid transporters (EAAT) play a key role in glutamatergic synaptic communication. Driven by transmembrane cation gradients, these transporters catalyze the reuptake of glutamate from the synaptic cleft once this neurotransmitter has been utilized for signaling. Two decades ago, pioneering studies in the Kanner lab identified a conserved methionine within the transmembrane domain as key for substrate turnover rate and specificity; later structural work, particularly for the prokaryotic homologs Glt Ph and Glt Tk , revealed that this methionine is involved in the coordination of one of the three Na + ions that are co-transported with the substrate. Albeit extremely atypical, the existence of this interaction is consistent with biophysical analyses of Glt Ph showing that mutations of this methionine diminish the binding cooperativity between substrates and Na + . It has been unclear, however, whether this intriguing methionine influences the thermodynamics of the transport reaction, i.e., its substrate:ion stoichiometry, or whether it simply fosters a specific kinetics in the binding reaction, which, while influential for the turnover rate, do not fundamentally explain the ion-coupling mechanism of this class of transporters. Here, studies of Glt Tk using experimental and computational methods independently arrive at the conclusion that the latter hypothesis is the most plausible, and lay the groundwork for future efforts to uncover the underlying mechanism., (© 2021. The Author(s).)

A methodology is proposed for the calculation of multidimensional free-energy landscapes of molecular systems, based on analysis of multiple molecular dynamics trajectories wherein adaptive biases have been applied to enhance the sampling of different collective variables. In this approach, which we refer to as the Force-Correction Analysis Method (FCAM), local averages of the total and biasing forces are evaluated post hoc, and the latter are subtracted from the former to obtain unbiased estimates of the mean force across collective-variable space. Multidimensional free-energy surfaces and minimum free-energy pathways are then derived by integrating the mean-force landscape with a kinetic Monte Carlo algorithm. To evaluate the proposed method, a series of numerical tests and comparisons with existing approaches were carried out for small molecules, peptides, and proteins, based on all-atom trajectories generated with standard, concurrent, and replica-exchange metadynamics in collective-variable spaces ranging from one to six dimensional. The tests confirm the correctness of the FCAM formulation and demonstrate that calculated mean forces and free energies converge rapidly and accurately, outperforming other methods used to unbias this kind of simulation data.

The nonpolarizable CHARMM force field is one of the most widely used energy functions for all-atom biomolecular simulations. Chloride is the only halide ion included in the latest version, CHARMM36m, and is used widely in simulation studies, often as an electrolyte ion but also as the biological substrate of transport proteins and enzymes. Here, we find that existing parameters systematically underestimate the interaction of Cl - with proteins and lipids. Accordingly, when examined in solution, little to no Cl - association can be observed with most components of the protein, including backbone, polar side chains and aromatic rings. The strength of the interaction with cationic side chains and with alkali ions is also incongruent with experimental measurements, specifically osmotic coefficients of concentrated solutions. Consistent with these findings, a 4-μs trajectory of the Cl - -specific transport protein CLC-ec1 shows irreversible Cl - dissociation from the so-called S cen binding site, even in a 150 mM NaCl buffer. To correct for these deficiencies, we formulate a series of pair-specific Lennard-Jones parameters that override those resulting from the conventional Lorentz-Berthelot combination rules. These parameters, referred to as NBFIX, are systematically calibrated against available experimental data as well as ab initio geometry optimizations and energy evaluations, for a wide set of binary and ternary Cl - complexes with protein and lipid analogs and alkali cations. Analogously, we also formulate parameter sets for the other three biological halide ions, namely, fluoride, bromide, and iodide. The resulting parameters are used to calculate the potential of mean force defining the interaction of each anion and each of the protein and lipid analogues in bulk water, revealing association free energies in the range of -0.3 to -3.3 kcal/mol, with the F - complexes being the least stable. The NBFIX corrections also preserve the Cl - occupancy of CLC-ec1 in a second 4-μs trajectory. We posit that these optimized molecular-mechanics models provide a more realistic foundation for all-atom simulation studies of processes entailing changes in hydration, recognition, or transport of halide anions.

Over two-thirds of integral membrane proteins of known structure assemble into oligomers. Yet, the forces that drive the association of these proteins remain to be delineated, as the lipid bilayer is a solvent environment that is both structurally and chemically complex. In this study, we reveal how the lipid solvent defines the dimerization equilibrium of the CLC-ec1 Cl - /H + antiporter. Integrating experimental and computational approaches, we show that monomers associate to avoid a thinned-membrane defect formed by hydrophobic mismatch at their exposed dimerization interfaces. In this defect, lipids are strongly tilted and less densely packed than in the bulk, with a larger degree of entanglement between opposing leaflets and greater water penetration into the bilayer interior. Dimerization restores the membrane to a near-native state and therefore, appears to be driven by the larger free-energy cost of lipid solvation of the dissociated protomers. Supporting this theory, we demonstrate that addition of short-chain lipids strongly shifts the dimerization equilibrium toward the monomeric state, and show that the cause of this effect is that these lipids preferentially solvate the defect. Importantly, we show that this shift requires only minimal quantities of short-chain lipids, with no measurable impact on either the macroscopic physical state of the membrane or the protein's biological function. Based on these observations, we posit that free-energy differentials for local lipid solvation define membrane-protein association equilibria. With this, we argue that preferential lipid solvation is a plausible cellular mechanism for lipid regulation of oligomerization processes, as it can occur at low concentrations and does not require global changes in membrane properties., Competing Interests: RC, NB, EK, ST, KG, AG, TÖ, LH, AF, VK No competing interests declared, JF Senior editor, eLife, JR Reviewing editor, eLife

In both prokaryotes and eukaryotes, multidrug and toxic-compound extrusion (MATE) transporters catalyze the efflux of a broad range of cytotoxic compounds, including human-made antibiotics and anticancer drugs. MATEs are secondary-active antiporters, i.e. their drug-efflux activity is coupled to, and powered by, the uptake of ions down a pre-existing transmembrane electrochemical gradient. Key aspects of this mechanism, however, remain to be delineated, such as its ion specificity and stoichiometry. We previously revealed the existence of a Na+-binding site in a MATE transporter from Pyroccocus furiosus (PfMATE) and hypothesized that this site might be broadly conserved among prokaryotic MATEs. Here, we evaluate this hypothesis by analyzing VcmN and ClbM, which along with PfMATE are the only three prokaryotic MATEs whose molecular structures have been determined at resolutions better than 3 Å. Analysis of available crystallographic data and molecular dynamics simulations indeed reveal an occupied Na+-binding site in the N-terminal lobe of both structures, analogous to that identified in PfMATE. We likewise find this site to be strongly selective against K+, suggesting it is mechanistically significant. Consistent with these computational results, DEER spectroscopy measurements for multiple doubly-spin-labeled VcmN constructs demonstrate Na+-dependent changes in protein conformation. The existence of this binding site in three MATE orthologs implicates Na+ in the ion-coupled drug-efflux mechanisms of this class of transporters. These results also imply that observations of H+-dependent activity stem either from a site elsewhere in the structure, or from H+ displacing Na+ under certain laboratory conditions, as has been noted for other Na+-driven transport systems., (Published under license by The American Society for Biochemistry and Molecular Biology, Inc.)

Hydrogen-deuterium exchange (HDX) is a comprehensive yet detailed probe of protein structure and dynamics and, coupled to mass spectrometry, has become a powerful tool for investigating an increasingly large array of systems. Computer simulations are often used to help rationalize experimental observations of exchange, but interpretations have frequently been limited to simple, subjective correlations between microscopic dynamical fluctuations and the observed macroscopic exchange behavior. With this in mind, we previously developed the HDX ensemble reweighting approach and associated software, HDXer, to aid the objective interpretation of HDX data using molecular simulations. HDXer has two main functions; first, to compute H-D exchange rates that describe each structure in a candidate ensemble of protein structures, for example from molecular simulations, and second, to objectively reweight the conformational populations present in a candidate ensemble to conform to experimental exchange data. In this article, we first describe the HDXer approach, theory, and implementation. We then guide users through a suite of tutorials that demonstrate the practical aspects of preparing experimental data, computing HDX levels from molecular simulations, and performing ensemble reweighting analyses. Finally we provide a practical discussion of the capabilities and limitations of the HDXer methods including recommendations for a user's own analyses. Overall, this article is intended to provide an up-to-date, pedagogical counterpart to the software, which is freely available at https://github.com/Lucy-Forrest-Lab/HDXer., Competing Interests: 8 Potentially Conflicting Interests The authors declare no potentially conflicting interests.

ATG9, the only transmembrane protein in the core macroautophagy/autophagy machinery, is a key player in the early stages of autophagosome formation. Yet, the lack of a high-resolution structure of ATG9 was a major impediment in understanding its three-dimensional organization and function. We recently solved a high-resolution cryoEM structure of the ubiquitously expressed human ATG9A isoform. The structure revealed that ATG9A is a domain-swapped homotrimer with a unique fold, and has an internal network of branched cavities. In cellulo analyses demonstrated the functional importance of the cavity-lining residues. These cavities could serve as conduits for transport of hydrophilic moieties, such as lipid headgroups, across the bilayer. Finally, structure-guided molecular dynamics predicted that ATG9A has membrane-bending properties, which is consistent with its localization to highly curved membranes.

Multidrug efflux pumps present a challenge to the treatment of bacterial infections, making it vitally important to understand their mechanism of action. Here, we investigate the nature of substrate binding within Lactococcus lactis LmrP, a prototypical multidrug transporter of the major facilitator superfamily. We determined the crystal structure of LmrP in a ligand-bound outward-open state and observed an embedded lipid in the binding cavity of LmrP, an observation supported by native mass spectrometry analyses. Molecular dynamics simulations suggest that the anionic lipid stabilizes the observed ligand-bound structure. Mutants engineered to disrupt binding of the embedded lipid display reduced transport of some, but not all, antibiotic substrates. Our results suggest that a lipid within the binding cavity could provide a malleable hydrophobic component that allows adaptation to the presence of different substrates, helping to explain the broad specificity of this protein and possibly other multidrug transporters.

S-acylation, whereby a fatty acid chain is covalently linked to a cysteine residue by a thioester linkage, is the most prevalent kind of lipid modification of proteins. Thousands of proteins are targets of this post-translational modification, which is catalyzed by a family of eukaryotic integral membrane enzymes known as DHHC protein acyltransferases (DHHC-PATs). Our knowledge of the repertoire of S-acylated proteins has been rapidly expanding owing to development of the chemoproteomic techniques. There has also been an increasing number of reports in the literature documenting the importance of S-acylation in human physiology and disease. Recently, the first atomic structures of two different DHHC-PATs were determined using X-ray crystallography. This review will focus on the insights gained into the molecular mechanism of DHHC-PATs from these structures and highlight representative data from the biochemical literature that they help explain., (Published by Elsevier Ltd.)

Bedaquiline (BDQ, Sirturo) has been approved to treat multidrug resistant forms of Mycobacterium tuberculosis. Prior studies suggested that BDQ was a selective inhibitor of the ATP synthase from M. tuberculosis. However, Sirturo treatment leads to an increased risk of cardiac arrhythmias and death, raising the concern that this adverse effect results from inhibition at a secondary site. Here we show that BDQ is a potent inhibitor of the yeast and human mitochondrial ATP synthases. Single-particle cryo-EM reveals that the site of BDQ inhibition partially overlaps with that of the inhibitor oligomycin. Molecular dynamics simulations indicate that the binding mode of BDQ to this site is similar to that previously seen for a mycobacterial enzyme, explaining the observed lack of selectivity. We propose that derivatives of BDQ ought to be made to increase its specificity toward the mycobacterial enzyme and thereby reduce the side effects for patients that are treated with Sirturo.

Hydrogen fluoride (HF) is the most polar diatomic molecule and one of the simplest molecules capable of hydrogen-bonding. HF deviates from ideality both in the gas phase and in solution and is thus of great interest from a fundamental standpoint. Pure and aqueous HF solutions are broadly used in chemical and industrial processes, despite their high toxicity. HF is a stable species also in some biological conditions, because it does not readily dissociate in water unlike other hydrogen halides; yet, little is known about how HF interacts with biomolecules. Here, we set out to develop a molecular-mechanics model to enable computer simulations of HF in chemical and biological applications. This model is based on a comprehensive high-level ab initio quantum chemical investigation of the structure and energetics of the HF monomer and dimer; (HF) n clusters, for n = 3-7; various clusters of HF and H 2 O; and complexes of HF with analogs of all 20 amino acids and of several commonly occurring lipids, both neutral and ionized. This systematic analysis explains the unique properties of this molecule: for example, that interacting HF molecules favor nonlinear geometries despite being diatomic and that HF is a strong H-bond donor but a poor acceptor. The ab initio data also enables us to calibrate a three-site molecular-mechanics model, with which we investigate the structure and thermodynamic properties of gaseous, liquid, and supercritical HF in a wide range of temperatures and pressures; the solvation structure of HF in water and of H 2 O in liquid HF; and the free diffusion of HF across a lipid bilayer, a key process underlying the high cytotoxicity of HF. Despite its inherent simplifications, the model presented significantly improves upon previous efforts to capture the properties of pure and aqueous HF fluids by molecular-mechanics methods and to our knowledge constitutes the first parameter set calibrated for biomolecular simulations.

Autophagy is a catabolic process involving capture of cytoplasmic materials into double-membraned autophagosomes that subsequently fuse with lysosomes for degradation of the materials by lysosomal hydrolases. One of the least understood components of the autophagy machinery is the transmembrane protein ATG9. Here, we report a cryoelectron microscopy structure of the human ATG9A isoform at 2.9-Å resolution. The structure reveals a fold with a homotrimeric domain-swapped architecture, multiple membrane spans, and a network of branched cavities, consistent with ATG9A being a membrane transporter. Mutational analyses support a role for the cavities in the function of ATG9A. In addition, structure-guided molecular simulations predict that ATG9A causes membrane bending, explaining the localization of this protein to small vesicles and highly curved edges of growing autophagosomes., Competing Interests: Declaration of Interests The authors declare no competing interest., (Published by Elsevier Inc.)

Hydrogen-deuterium exchange combined with mass spectrometry (HDX-MS) is a widely applied biophysical technique that probes the structure and dynamics of biomolecules without the need for site-directed modifications or bio-orthogonal labels. The mechanistic interpretation of HDX data, however, is often qualitative and subjective, owing to a lack of quantitative methods to rigorously translate observed deuteration levels into atomistic structural information. To help address this problem, we have developed a methodology to generate structural ensembles that faithfully reproduce HDX-MS measurements. In this approach, an ensemble of protein conformations is first generated, typically using molecular dynamics simulations. A maximum-entropy bias is then applied post hoc to the resulting ensemble such that averaged peptide-deuteration levels, as predicted by an empirical model, agree with target values within a given level of uncertainty. We evaluate this approach, referred to as HDX ensemble reweighting (HDXer), for artificial target data reflecting the two major conformational states of a binding protein. We demonstrate that the information provided by HDX-MS experiments and by the model of exchange are sufficient to recover correctly weighted structural ensembles from simulations, even when the relevant conformations are rarely observed. Degrading the information content of the target data-e.g., by reducing sequence coverage, by averaging exchange levels over longer peptide segments, or by incorporating different sources of uncertainty-reduces the structural accuracy of the reweighted ensemble but still allows for useful insights into the distinctive structural features reflected by the target data. Finally, we describe a quantitative metric to rank candidate structural ensembles according to their correspondence with target data and illustrate the use of HDXer to describe changes in the conformational ensemble of the membrane protein LeuT. In summary, HDXer is designed to facilitate objective structural interpretations of HDX-MS data and to inform experimental approaches and further developments of theoretical exchange models., (Published by Elsevier Inc.)

Cysteine palmitoylation, a form of S-acylation, is a key posttranslational modification in cellular signaling. This type of reversible lipidation occurs in both plasma and organellar membranes, and is catalyzed by a family of integral membrane proteins known as DHHC acyltransferases. The first step in the S-acylation process is the recognition of free acyl coenzyme A (acyl-CoA) from the lipid bilayer. The DHHC enzyme then becomes autoacylated at a site defined by a conserved Asp-His-His-Cys motif. This reaction entails ionization of the catalytic Cys. Intriguingly, in known DHHC structures, this catalytic Cys appears to be exposed to the hydrophobic interior of the lipid membrane, which would be highly unfavorable for a negatively charged nucleophile, thus hindering autoacylation. Here, we use biochemical and computational methods to reconcile these seemingly contradictory facts. First, we experimentally demonstrate that human DHHC20 is active when reconstituted in POPC nanodiscs. Microsecond-long all-atom molecular dynamics simulations are then calculated for human DHHC20 and for different acyl-CoA forms, also in a POPC membrane. Strikingly, we observe that human DHHC20 induces a drastic deformation in the membrane, particularly on the cytoplasmic side, where autoacylation occurs. As a result, the catalytic Cys becomes hydrated and optimally positioned to encounter the cleavage site in acyl-CoA. In summary, we hypothesize that DHHC enzymes locally reshape the membrane to foster a morphology that is specifically adapted for acyl-CoA recognition and autoacylation., (Published by Elsevier Inc.)

We report a methodology to calculate the free energy of a shape transformation in a lipid membrane directly from a molecular dynamics simulation. The bilayer need not be homogeneous or symmetric and can be atomically detailed or coarse grained. The method is based on a collective variable that quantifies the similarity between the membrane and a set of predefined density distributions. Enhanced sampling of this "Multi-Map" variable re-shapes the bilayer and permits the derivation of the corresponding potential of mean force. Calculated energies thus reflect the dynamic interplay of atoms and molecules, rather than postulated effects. Evaluation of deformations of different shape, amplitude, and range demonstrates that the macroscopic bending modulus assumed by the Helfrich-Canham model is increasingly unsuitable below the 100-Å scale. In this range of major biological significance, direct free-energy calculations reveal a much greater plasticity. We also quantify the stiffening effect of cholesterol on bilayers of different composition and compare with experiments. Lastly, we illustrate how this approach facilitates analysis of other solvent reorganization processes, such as hydrophobic hydration. Published 2019. This article is a U.S. Government work and is in the public domain in the USA., (Published 2019. This article is a U.S. Government work and is in the public domain in the USA.)

That channels and transporters can influence the membrane morphology is increasingly recognized. Less appreciated is that the extent and free-energy cost of these deformations likely varies among different functional states of a protein, and thus, that they might contribute significantly to defining its mechanism. We consider the trimeric Na + -aspartate symporter Glt Ph , a homolog of an important class of neurotransmitter transporters, whose mechanism entails one of the most drastic structural changes known. Molecular simulations indicate that when the protomers become inward-facing, they cause deep, long-ranged, and yet mutually-independent membrane deformations. Using a novel simulation methodology, we estimate that the free-energy cost of this membrane perturbation is in the order of 6-7 kcal/mol per protomer. Compensating free-energy contributions within the protein or its environment must thus stabilize this inward-facing conformation for the transporter to function. We discuss these striking results in the context of existing experimental observations for this and other transporters., Competing Interests: WZ, GF, CA, HK, HP No competing interests declared, LF Reviewing editor, eLife, JF Senior editor, eLife

P2X3 receptor channels expressed in sensory neurons are activated by extracellular ATP and serve important roles in nociception and sensory hypersensitization, making them attractive therapeutic targets. Although several P2X3 structures are known, it is unclear how physiologically abundant Ca 2+ -ATP and Mg 2+ -ATP activate the receptor, or how divalent cations regulate channel function. We used structural, computational and functional approaches to show that a crucial acidic chamber near the nucleotide-binding pocket in human P2X3 receptors accommodates divalent ions in two distinct modes in the absence and presence of nucleotide. The unusual engagement between the receptor, divalent ion and the γ-phosphate of ATP enables channel activation by ATP-divalent complex, cooperatively stabilizes the nucleotide on the receptor to slow ATP unbinding and recovery from desensitization, a key mechanism for limiting channel activity. These findings reveal how P2X3 receptors recognize and are activated by divalent-bound ATP, aiding future physiological investigations and drug development., Competing Interests: ML, YW, RB, FM, SS, MH No competing interests declared, JF, KS Reviewing editor, eLife

Many neurotoxins inflict pain by targeting receptors expressed on nociceptors, such as the polymodal cationic channel TRPV1. The tarantula double-knot toxin (DkTx) is a peptide with an atypical bivalent structure, providing it with the unique capability to lock TRPV1 in its open state and evoke an irreversible channel activation. Here, we describe a distinct gating mechanism of DkTx-evoked TRPV1 activation. Interestingly, DkTx evokes significantly smaller TRPV1 macroscopic currents than capsaicin, with a significantly lower unitary conductance. Accordingly, while capsaicin evokes aversive behaviors in TRPV1-transgenic Caenorhabditis elegans , DkTx fails to evoke such response at physiological concentrations. To determine the structural feature(s) responsible for this phenomenon, we engineered and evaluated a series of mutated toxins and TRPV1 channels. We found that elongating the DkTx linker, which connects its two knots, increases channel conductance compared with currents elicited by the native toxin. Importantly, deletion of the TRPV1 pore turret, a stretch of amino acids protruding out of the channel's outer pore region, is sufficient to produce both full conductance and aversive behaviors in response to DkTx. Interestingly, this deletion decreases the capsaicin-evoked channel activation. Taken together with structure modeling analysis, our results demonstrate that the TRPV1 pore turret restricts DkTx-mediated pore opening, probably through steric hindrance, limiting the current size and mitigating the evoked downstream physiological response. Overall, our findings reveal that DkTx and capsaicin elicit distinct TRPV1 gating mechanisms and subsequent pain responses. Our results also indicate that the TRPV1 pore turret regulates the mechanisms of channel gating and permeation., Competing Interests: The authors declare no conflict of interest.

Given its ability to measure multicomponent distance distributions between electron-spin probes, double electron-electron resonance (DEER) spectroscopy has become a leading technique to assess the structural dynamics of biomolecules. However, methodologies to evaluate the statistical error of these distributions are not standard, often hampering a rigorous interpretation of the experimental results. Distance distributions are often determined from the experimental DEER data through a mathematical method known as Tikhonov regularization, but this approach makes rigorous error estimates difficult. Here, we build upon an alternative, model-based approach in which the distance probability distribution is represented as a sum of Gaussian components, and use propagation of errors to calculate an associated confidence band. Our approach considers all sources of uncertainty, including the experimental noise, the uncertainty in the fitted background signal, and the limited time span of the data collection. The resulting confidence band reveals the most and least reliable features of the probability distribution, thereby informing the structural interpretation of DEER experiments. To facilitate this interpretation, we also generalize the molecular simulation method known as ensemble-biased metadynamics (EBMetaD). This method, originally designed to generate maximal-entropy structural ensembles consistent with one or more probability distributions, now also accounts for the uncertainty in those target distributions exactly as dictated by their confidence bands. After careful benchmarks, we demonstrate the proposed techniques using DEER results from spin-labeled T4 lysozyme., (Copyright © 2018 Biophysical Society. All rights reserved.)

Enzymes in the respiratory chain are increasingly seen as potential targets against multi-drug resistance of human pathogens and cancerous cells. However, a detailed understanding of the mechanism and specificity determinants of known inhibitors is still lacking. Oligomycin, for example, has been known to be an inhibitor of the membrane motor of the mitochondrial ATP synthase for over five decades, and yet little is known about its mode of action at the molecular level. In a recent breakthrough, a crystal structure of the S. cerevisiae c-subunit ring with bound oligomycin revealed the inhibitor docked on the outer face of the proton-binding sites, deep into the transmembrane region. However, the structure of the complex was obtained in an organic solvent rather than detergent or a lipid bilayer, and therefore it has been unclear whether this mode of recognition is physiologically relevant. Here, we use molecular dynamics simulations to address this question and gain insights into the mechanism of oligomycin inhibition. Our findings lead us to propose that oligomycin naturally partitions into the lipid/water interface, and that in this environment the inhibitor can indeed bind to any of the c-ring proton-carrying sites that are exposed to the membrane, thereby becoming an integral component of the proton-coordinating network. As the c-ring rotates within the membrane, driven either by downhill proton permeation or ATP hydrolysis, one of the protonated, oligomycin-bound sites eventually reaches the subunit-a interface and halts the rotary mechanism of the enzyme., (Copyright © 2018. Published by Elsevier B.V.)

Multidrug and toxic-compound extrusion (MATE) proteins comprise an important but largely uncharacterized family of secondary-active transporters. In both eukaryotes and prokaryotes, these transporters protect the cell by catalyzing the efflux of a broad range of cytotoxic compounds, including human-made antibiotics and anticancer drugs. MATEs are thus potential pharmacological targets against drug-resistant pathogenic bacteria and tumor cells. The activity of MATEs is powered by transmembrane electrochemical ion gradients, but their molecular mechanism and ion specificity are not understood, in part because high-quality structural information is limited. Here, we use computational methods to study PfMATE, from Pyrococcus furiosus , whose structure is the best resolved to date. Analysis of available crystallographic data and additional molecular dynamics simulations unequivocally reveal an occupied Na + -binding site in the N-lobe of this transporter, which had not been previously recognized. We find this site to be selective against K + and broadly conserved among prokaryotic MATEs, including homologs known to be Na + -dependent such as NorM-VC, VmrA, and ClbM, for which the location of the Na + site had been debated. We note, however, that the chemical makeup of the proposed Na + site indicates it is weakly specific against H + , explaining why MATEs featuring this Na + -binding motif may be solely driven by H + in laboratory conditions. We further posit that the concurrent coupling to H + and Na + gradients observed for some Na + -driven MATEs owes to a second H + -binding site, within the C-lobe. In summary, our study provides insights into the structural basis for the complex ion dependency of MATE transporters., Competing Interests: The authors declare no conflict of interest.

Mitochondrial adenosine triphosphate (ATP) synthase comprises a membrane embedded F o motor that rotates to drive ATP synthesis in the F 1 subunit. We used single-particle cryo-electron microscopy (cryo-EM) to obtain structures of the full complex in a lipid bilayer in the absence or presence of the inhibitor oligomycin at 3.6- and 3.8-angstrom resolution, respectively. To limit conformational heterogeneity, we locked the rotor in a single conformation by fusing the F6 subunit of the stator with the δ subunit of the rotor. Assembly of the enzyme with the F6-δ fusion caused a twisting of the rotor and a 9° rotation of the F o c 10 -ring in the direction of ATP synthesis, relative to the structure of isolated F o Our cryo-EM structures show how F 1 and F o are coupled, give insight into the proton translocation pathway, and show how oligomycin blocks ATP synthesis., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Adenosine triphosphate (ATP) synthases populate the inner membranes of mitochondria, where they produce the majority of the ATP required by the cell. From yeast to vertebrates, cryoelectron tomograms of these membranes have consistently revealed a very precise organization of these enzymes. Rather than being scattered throughout the membrane, the ATP synthases form dimers, and these dimers are organized into rows that extend for hundreds of nanometers. The rows are only observed in the membrane invaginations known as cristae, specifically along their sharply curved edges. Although the presence of these macromolecular structures has been irrefutably linked to the proper development of cristae morphology, it has been unclear what drives the formation of the rows and why they are specifically localized in the cristae. In this study, we present a quantitative molecular-simulation analysis that strongly suggests that the dimers of ATP synthases organize into rows spontaneously, driven by a long-range attractive force that arises from the relief of the overall elastic strain of the membrane. The strain is caused by the V-like shape of the dimers, unique among membrane protein complexes, which induces a strong deformation in the surrounding membrane. The process of row formation is therefore not a result of direct protein-protein interactions or a specific lipid composition of the membrane. We further hypothesize that, once assembled, the ATP synthase dimer rows prime the inner mitochondrial membrane to develop folds and invaginations by causing macroscopic membrane ridges that ultimately become the edges of cristae. In this way, mitochondrial ATP synthases would contribute to the generation of a morphology that maximizes the surface area of the inner membrane, and thus ATP production. Finally, we outline key experiments that would be required to verify or refute this hypothesis., (This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.)

Intracellular Ca 2+ signals control a wide array of cellular processes. These signals require spatial and temporal regulation of the intracellular Ca 2+ concentration, which is achieved in part by a class of ubiquitous membrane proteins known as sodium-calcium exchangers (NCXs). NCXs are secondary-active antiporters that power the translocation of Ca 2+ across the cell membrane by coupling it to the flux of Na + in the opposite direction, down an electrochemical gradient. Na + and Ca 2+ are translocated in separate steps of the antiport cycle, each of which is thought to entail a mechanism whereby ion-binding sites within the protein become alternately exposed to either side of the membrane. The prokaryotic exchanger NCX_Mj, the only member of this family with known structure, has been proposed to be a good functional and structural model of mammalian NCXs; yet our understanding of the functional properties of this protein remains incomplete. Here, we study purified NCX_Mj reconstituted into liposomes under well-controlled experimental conditions and demonstrate that this homologue indeed shares key functional features of the NCX family. Transport assays and reversal-potential measurements enable us to delineate the essential characteristics of this antiporter and establish that its ion-exchange stoichiometry is 3Na + :1Ca 2+ Together with previous studies, this work confirms that NCX_Mj is a valid model system to investigate the mechanism of ion recognition and membrane transport in sodium-calcium exchangers., (This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.)

The membrane transporter anion exchanger 1 (AE1), or band 3, is a key component in the processes of carbon-dioxide transport in the blood and urinary acidification in the renal collecting duct. In both erythrocytes and the basolateral membrane of the collecting-duct α-intercalated cells, the role of AE1 is to catalyze a one-for-one exchange of chloride for bicarbonate. After decades of biochemical and functional studies, the structure of the transmembrane region of AE1, which catalyzes the anion-exchange reaction, has finally been determined. Each protomer of the AE1 dimer comprises two repeats with inverted transmembrane topologies, but the structures of these repeats differ. This asymmetry causes the putative substrate-binding site to be exposed only to the extracellular space, consistent with the expectation that anion exchange occurs via an alternating-access mechanism. Here, we hypothesize that the unknown, inward-facing conformation results from inversion of this asymmetry, and we propose a model of this state constructed using repeat-swap homology modeling. By comparing this inward-facing model with the outward-facing experimental structure, we predict that the mechanism of AE1 involves an elevator-like motion of the substrate-binding domain relative to the nearly stationary dimerization domain and to the membrane plane. This hypothesis is in qualitative agreement with a wide range of biochemical and functional data, which we review in detail, and suggests new avenues of experimentation., (© 2017 Ficici et al.)

Recent years have witnessed a renewed interest in the ATP synthase as a drug target against human pathogens. Indeed, clinical, biochemical, and structural data indicate that hydrophobic inhibitors targeting the membrane-embedded proton-binding sites of the c-subunit ring could serve as last-resort antibiotics against multidrug resistant strains. However, because inhibition of the mitochondrial ATP synthase in humans is lethal, it is essential that these inhibitors be not only potent but also highly selective for the bacterial enzyme. To this end, a detailed understanding of the structure of this protein target is arguably instrumental. Here, we use computational methods to predict the atomic structures of the proton-binding sites in two prototypical c-rings: that of the ATP synthase from Saccharomyces cerevisiae, which is a model system for mitochondrial enzymes, and that from Escherichia coli, which can be pathogenic for humans. Our study reveals the structure of these binding sites loaded with protons and in the context of the membrane, that is, in the state that would mediate the recognition of a potential inhibitor. Both structures reflect a mode of proton coordination unlike those previously observed in other c-ring structures, whether experimental or modeled.

Pathological metabolic conditions such as ischemia induce the rupture of the mitochondrial envelope and the release of pro-apoptotic proteins, leading to cell death. At the onset of this process, the inner mitochondrial membrane becomes depolarized and permeable to osmolytes, proposedly due to the opening of a non-selective protein channel of unknown molecular identity. A recent study purports that this channel, referred to as Mitochondrial Permeability Transition Pore (MPTP), is formed within the c-subunit ring of the ATP synthase, upon its dissociation from the catalytic domain of the enzyme. Here, we examine this claim for two c-rings of different lumen width, through calculations of their ion conductance and selectivity based on all-atom molecular dynamics simulations. We also quantify the likelihood that the lumen of these c-rings is in a hydrated, potentially conducting state rather than empty or blocked by lipid molecules. These calculations demonstrate that the structure and biophysical properties of a correctly assembled c-ring are inconsistent with those attributed to the MPTP.

The atomic mechanisms of isomerization of ATP-Mg(2+) in solution are characterized using the recently developed String Method with Optimal Molecular Alignment (SOMA) and molecular-dynamics simulations. Bias-Exchange Metadynamics simulations are first performed to identify the primary conformers of the ATP-Mg(2+) complex and their connectivity. SOMA is then used to elucidate the minimum free-energy path (MFEP) for each transition, in a 48-dimensional space. Analysis of the per-atom contributions to the global free-energy profiles reveals that the mechanism of these transitions is controlled by the Mg(2+) ion and its coordinating oxygen atoms in the triphosphate moiety, as well as by the ion-hydration shell. Metadynamics simulations in path collective variables based on the MFEP demonstrate these isomerizations proceed across a narrow channel of configurational space, thus validating the premise underlying SOMA. This study provides a roadmap for the examination of conformational changes in biomolecules, based on complementary enhanced-sampling techniques with different strengths. © 2015 Wiley Periodicals, Inc., (© 2015 Wiley Periodicals, Inc.)

Numerous membrane transporters and enzymes couple their mechanisms to the permeation of Na(+) or H(+), thereby harnessing the energy stored in the form of transmembrane electrochemical potential gradients to sustain their activities. The molecular and environmental factors that control and modulate the ion specificity of most of these systems are, however, poorly understood. Here, we use isothermal titration calorimetry to determine the Na(+)/H(+) selectivity of the ion-driven membrane rotor of an F-type ATP synthase. Consistent with earlier theoretical predictions, we find that this rotor is significantly H(+) selective, although not sufficiently to be functionally coupled to H(+), owing to the large excess of Na(+) in physiological settings. The functional Na(+) specificity of this ATP synthase thus results from two opposing factors, namely its inherent chemical selectivity and the relative availability of the coupling ion. Further theoretical studies of this membrane rotor, and of two others with a much stronger and a slightly weaker H(+) selectivity, indicate that, although the inherent selectivity of their ion-binding sites is largely set by the balance of polar and hydrophobic groups flanking a conserved carboxylic side chain, subtle variations in their structure and conformational dynamics, for a similar chemical makeup, can also have a significant contribution. We propose that the principle of ion selectivity outlined here may provide a rationale for the differentiation of Na(+)- and H(+)-coupled systems in other families of membrane transporters and enzymes.

Na(+)/Ca(2+) exchangers (NCXs) are ubiquitous membrane transporters with a key role in Ca(2+) homeostasis and signaling. NCXs mediate the bidirectional translocation of either Na(+) or Ca(2+), and thus can catalyze uphill Ca(2+) transport driven by a Na(+) gradient, or vice versa. In a major breakthrough, a prokaryotic NCX homolog (NCX_Mj) was recently isolated and its crystal structure determined at atomic resolution. The structure revealed an intriguing architecture consisting of two inverted-topology repeats, each comprising five transmembrane helices. These repeats adopt asymmetric conformations, yielding an outward-facing occluded state. The crystal structure also revealed four putative ion-binding sites, but the occupancy and specificity thereof could not be conclusively established. Here, we use molecular-dynamics simulations and free-energy calculations to identify the ion configuration that best corresponds to the crystallographic data and that is also thermodynamically optimal. In this most probable configuration, three Na(+) ions occupy the so-called Sext, SCa, and Sint sites, whereas the Smid site is occupied by one water molecule and one H(+), which protonates an adjacent aspartate side chain (D240). Experimental measurements of Na(+)/Ca(2+) and Ca(2+)/Ca(2+) exchange by wild-type and mutagenized NCX_Mj confirm that transport of both Na(+) and Ca(2+) requires protonation of D240, and that this side chain does not coordinate either ion at Smid. These results imply that the ion exchange stoichiometry of NCX_Mj is 3:1 and that translocation of Na(+) across the membrane is electrogenic, whereas transport of Ca(2+) is not. Altogether, these findings provide the basis for further experimental and computational studies of the conformational mechanism of this exchanger.

We used electron cryotomography of mitochondrial membranes from wild-type and mutant Saccharomyces cerevisiae to investigate the structure and organization of ATP synthase dimers in situ. Subtomogram averaging of the dimers to 3.7 nm resolution revealed a V-shaped structure of twofold symmetry, with an angle of 86° between monomers. The central and peripheral stalks are well resolved. The monomers interact within the membrane at the base of the peripheral stalks. In wild-type mitochondria ATP synthase dimers are found in rows along the highly curved cristae ridges, and appear to be crucial for membrane morphology. Strains deficient in the dimer-specific subunits e and g or the first transmembrane helix of subunit 4 lack both dimers and lamellar cristae. Instead, cristae are either absent or balloon-shaped, with ATP synthase monomers distributed randomly in the membrane. Computer simulations indicate that isolated dimers induce a plastic deformation in the lipid bilayer, which is partially relieved by their side-by-side association. We propose that the assembly of ATP synthase dimer rows is driven by the reduction in the membrane elastic energy, rather than by direct protein contacts, and that the dimer rows enable the formation of highly curved ridges in mitochondrial cristae.

ATP synthase membrane rotors consist of a ring of c-subunits whose stoichiometry is constant for a given species but variable across different ones. We investigated the importance of c/c-subunit contacts by site-directed mutagenesis of a conserved stretch of glycines (GxGxGxGxG) in a bacterial c(11) ring. Structural and biochemical studies show a direct, specific influence on the c-subunit stoichiometry, revealing c(<11), c(12), c(13), c(14), and c(>14) rings. Molecular dynamics simulations rationalize this effect in terms of the energetics and geometry of the c-subunit interfaces. Quantitative data from a spectroscopic interaction study demonstrate that the complex assembly is independent of the c-ring size. Real-time ATP synthesis experiments in proteoliposomes show the mutant enzyme, harboring the larger c(12) instead of c(11), is functional at lower ion motive force. The high degree of compliance in the architecture of the ATP synthase rotor offers a rationale for the natural diversity of c-ring stoichiometries, which likely reflect adaptations to specific bioenergetic demands. These results provide the basis for bioengineering ATP synthases with customized ion-to-ATP ratios, by sequence modifications.

ATP synthases are the primary source of ATP in all living cells. To catalyze ATP synthesis, these membrane-associated complexes use a rotary mechanism powered by the transmembrane diffusion of ions down a concentration gradient. ATP synthases are assumed to be driven either by H(+) or Na(+), reflecting distinct structural motifs in their membrane domains, and distinct metabolisms of the host organisms. Here, we study the methanogenic archaeon Methanosarcina acetivorans using assays of ATP hydrolysis and ion transport in inverted membrane vesicles, and experimentally demonstrate that the rotary mechanism of its ATP synthase is coupled to the concurrent translocation of both H(+) and Na(+) across the membrane under physiological conditions. Using free-energy molecular simulations, we explain this unprecedented observation in terms of the ion selectivity of the binding sites in the membrane rotor, which appears to have been tuned via amino acid substitutions so that ATP synthesis in M. acetivorans can be driven by the H(+) and Na(+) gradients resulting from methanogenesis. We propose that this promiscuity is a molecular mechanism of adaptation to life at the thermodynamic limit.

Numerous membrane importers rely on accessory water-soluble proteins to capture their substrates. These substrate-binding proteins (SBP) have a strong affinity for their ligands; yet, substrate release onto the low-affinity membrane transporter must occur for uptake to proceed. It is generally accepted that release is facilitated by the association of SBP and transporter, upon which the SBP adopts a conformation similar to the unliganded state, whose affinity is sufficiently reduced. Despite the appeal of this mechanism, however, direct supporting evidence is lacking. Here, we use experimental and theoretical methods to demonstrate that an allosteric mechanism of enhanced substrate release is indeed plausible. First, we report the atomic-resolution structure of apo TeaA, the SBP of the Na(+)-coupled ectoine TRAP transporter TeaBC from Halomonas elongata DSM2581(T), and compare it with the substrate-bound structure previously reported. Conformational free-energy landscape calculations based upon molecular dynamics simulations are then used to dissect the mechanism that couples ectoine binding to structural change in TeaA. These insights allow us to design a triple mutation that biases TeaA toward apo-like conformations without directly perturbing the binding cleft, thus mimicking the influence of the membrane transporter. Calorimetric measurements demonstrate that the ectoine affinity of the conformationally biased triple mutant is 100-fold weaker than that of the wild type. By contrast, a control mutant predicted to be conformationally unbiased displays wild-type affinity. This work thus demonstrates that substrate release from SBPs onto their membrane transporters can be facilitated by the latter through a mechanism of allosteric modulation of the former.

As new atomic structures of membrane proteins are resolved, they reveal increasingly complex transmembrane topologies, and highly irregular surfaces with crevices and pores. In many cases, specific interactions formed with the lipid membrane are functionally crucial, as is the overall lipid composition. Compounded with increasing protein size, these characteristics pose a challenge for the construction of simulation models of membrane proteins in lipid environments; clearly, that these models are sufficiently realistic bears upon the reliability of simulation-based studies of these systems. Here, we introduce GRIFFIN, which uses a versatile framework to automate and improve a widely-used membrane-embedding protocol. Initially, GRIFFIN carves out lipid and water molecules from a volume equivalent to that of the protein, so as to conserve the system density. In the subsequent optimization phase GRIFFIN adds an implicit grid-based protein force-field to a molecular dynamics simulation of the pre-carved membrane. In this force-field, atoms inside the implicit protein volume experience an outward force that will expel them from that volume, whereas those outside are subject to electrostatic and van-der-Waals interactions with the implicit protein. At each step of the simulation, these forces are updated by GRIFFIN and combined with the intermolecular forces of the explicit lipid-water system. This procedure enables the construction of realistic and reproducible starting configurations of the protein-membrane interface within a reasonable timeframe and with minimal intervention. GRIFFIN is a standalone tool designed to work alongside any existing molecular dynamics package, such as NAMD or GROMACS.

CD1 molecules function to present lipid-based antigens to T cells. Here we present the crystal structure of CD1c at 2.5 Å resolution, in complex with the pathogenic Mycobacterium tuberculosis antigen mannosyl-β1-phosphomycoketide (MPM). CD1c accommodated MPM's methylated alkyl chain exclusively in the A' pocket, aided by a unique exit portal underneath the α1 helix. Most striking was an open F' pocket architecture lacking the closed cavity structure of other CD1 molecules, reminiscent of peptide binding grooves of classical major histocompatibility complex molecules. This feature, combined with tryptophan-fluorescence quenching during loading of a dodecameric lipopeptide antigen, provides a compelling model by which both the lipid and peptide moieties of the lipopeptide are involved in CD1c presentation of lipopeptides.

CD1 proteins mediate the trafficking and presentation of a diverse range of lipid antigens to T-cell receptors, and thus play a key role in our adaptive immune system. Crystal structures of several CD1 isoforms reveal a highly conserved tertiary structure, but also great variability in the anatomy of their binding pockets, reflecting their distinct ligand specificity. The structure of one important member of the family, CD1c, remains unknown. CD1c is of great interest as it can present an unusual and potent lipid antigen, mannosyl-beta(1)-phosphomycoketide (MPM) from Mycobacterium tuberculosis, the causative agent of tuberculosis. CD1c has also been reported to present acetylated 12-amino-acid-long peptides (lipo-12), an observation with broad immunological implications but difficult to rationalize on structural grounds. To gain insights into the structural basis for the ligand specificity of CD1c, we have generated an atomic model of its binding domain using a detailed position-specific multiple-template homology modeling approach. This model reveals structural features unique to this isoform, particularly with regard to the so-called pocket F', which provide a compelling rationale for the ability of CD1c to bind not only branched alkyl chains such as in MPM, but also long lipopeptides comparable to those presented by MHC proteins. A model of CD1c with bound MPM was constructed and analyzed through molecular dynamics simulations, showing marked structural stability in the time-scale of 100 ns. A model of CD1c in complex with lipo-12 is also presented.

Regulation of signaling pathways in the cell often involves multidomain allosteric enzymes that are able to adopt alternate active or inactive conformations in response to specific stimuli. It is therefore of great interest to elucidate the energetic and structural determinants that govern the conformational plasticity of these proteins. In this study, free-energy computations have been used to address this fundamental question, focusing on one important family of signaling enzymes, the Src tyrosine kinases. Inactivation of these enzymes depends on the formation of an assembly comprising a tandem of SH3 and SH2 modules alongside a catalytic domain. Activation results from the release of the SH3 and SH2 domains, which are then believed to be structurally uncoupled by virtue of a flexible peptide link. In contrast to this view, this analysis shows that inactivation depends critically on the intrinsic propensity of the SH3-SH2 tandem to adopt conformations that are conducive to the assembled inactive state, even when no interactions with the rest of the kinase are possible. This funneling of the available conformational space is encoded within the SH3-SH2 connector, which appears to have evolved to modulate the flexibility of the tandem in solution. To further substantiate this notion, we show how constitutively activating mutations in the SH3-SH2 connector shift the assembly equilibrium toward the disassembled, active state. Based on a similar analysis of several constructs of the kinase complex, we propose that assembly is characterized by the progressive optimization of the protein's conformational energy, with little or no energetic frustration.

Molecular dynamics (MD) simulations are used to explore the dynamics of a membrane protein in its crystal environment. A 50-ns-duration simulation (at a temperature of 300 K) is performed for the crystallographic unit cell of the bacterial outer membrane protein OmpA. The unit cell contains four protein molecules, plus detergent molecules and water. An excellent correlation between simulated and experimental values of crystallographic B factors is observed. Effectively, 0.2 micros of protein trajectories are obtained, allowing a critical assessment of simulation quality. Some deficiency in conformational sampling is demonstrated, but averaging over multiple trajectories improves this limitation. The previously undescribed structure and dynamics of detergent molecules in a unit cell are reported here, providing insight into the interactions important in the formation and stabilization of the crystalline environment at room temperature. In particular, we show that at room temperature the detergent molecules form a dynamic, extended micellar structure spreading over adjacent OmpA monomers within the crystal.

In this paper we describe a method for setting up an atomistic simulation of a membrane protein in a hydrated lipid bilayer and report the effect of differing electrostatic parameters on the drift in the protein structure during the subsequent simulation. The method aims to generate a suitable cavity in the interior of a lipid bilayer, using the solvent-accessible surface of the protein as a template, during the course of a short steered molecular dynamics simulation of a solvated lipid membrane. This is achieved by a two-stage process: firstly, lipid molecules whose headgroups are inside a cylindrical volume equivalent to that defined by the protein surface are removed; then the protein-lipid interface is optimized by applying repulsive forces perpendicular to the protein surface, and of gradually increased magnitude, to the remaining lipid atoms inside the volume occupied by the protein surface until it is emptied. The protein itself may then be inserted. Using the bacterial membrane proteins KcsA and FhuA as test cases, we show how the method achieves the formation of a suitable cavity in the interior of a dimyristoylphosphatidylcholine lipid bilayer without perturbing the configuration of the non-interfacial regions of the previously equilibrated lipid bilayer, even in cases of membrane proteins with irregular geometrical shapes. In addition, we compare subsequent simulations in which the long-range electrostatic interactions are treated via either a cut-off or a particle-mesh Ewald method. The results show that the drift from the initial structure is less in the latter case, especially for KcsA and for the non-core secondary structural elements (i.e. surface loops) of both proteins.