Membrane mediated toppling mechanism of the folate energy coupling factor transporter

Assessing the Martini 3 protein model: A review of its path and potential

Coarse-grained (CG) protein models have become indispensable tools for studying many biological protein details, from conformational dynamics to the organization of protein macro-complexes, and even the interaction of proteins with other molecules. The Martini force field is one of the most widely used CG models for bio-molecular simulations, partly because of the enormous success of its protein model. With the recent release of a new and improved version of the Martini force field - Martini 3 - a new iteration of its protein model was also made available. The Martini 3 protein force field is an evolution of its Martini 2 counterpart, aimed at improving many of the shortcomings that had been previously identified. In this mini-review, we first provide a general overview of the model and then focus on the successful advances made in the short time since its release, many of which would not have been possible before. Furthermore, we discuss reported limitations, potential directions for model improvement and comment on what the likely future development and application avenues are.

Hit optimization by dynamic combinatorial chemistry on Streptococcus pneumoniae energy-coupling factor transporter ECF-PanT

Herein, we present the first application of target-directed dynamic combinatorial chemistry (tdDCC) to the whole complex of the highly dynamic transmembrane, energy-coupling factor (ECF) transporter ECF-PanT in Streptococcus pneumoniae. In addition, we successfully employed the tdDCC technique as a hit-identification and -optimization strategy that led to the identification of optimized ECF inhibitors with improved activity. We characterized the best compounds regarding cytotoxicity and performed computational modeling studies on the crystal structure of ECF-PanT to rationalize their binding mode. Notably, docking studies showed that the acylhydrazone linker is able to maintain the crucial interactions.

Identification of inhibitors targeting the energy-coupling factor (ECF) transporters

The energy-coupling factor (ECF) transporters are a family of transmembrane proteins involved in the uptake of vitamins in a wide range of bacteria. Inhibition of the activity of these proteins could reduce the viability of pathogens that depend on vitamin uptake. The central role of vitamin transport in the metabolism of bacteria and absence from humans make the ECF transporters an attractive target for inhibition with selective chemical probes. Here, we report on the identification of a promising class of inhibitors of the ECF transporters. We used coarse-grained molecular dynamics simulations on Lactobacillus delbrueckii ECF-FolT2 and ECF-PanT to profile the binding mode and mechanism of inhibition of this novel chemotype. The results corroborate the postulated mechanism of transport and pave the way for further drug-discovery efforts.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Expulsion mechanism of the substrate-translocating subunit in ECF transporters

Energy-coupling factor (ECF)-type transporters mediate the uptake of micronutrients in many bacteria. They consist of a substrate-translocating subunit (S-component) and an ATP-hydrolysing motor (ECF module) Previous data indicate that the S-component topples within the membrane to alternately expose the binding site to either side of the membrane. In many ECF transporters, the substrate-free S-component can be expelled from the ECF module. Here we study this enigmatic expulsion step by cryogenic electron microscopy and reveal that ATP induces a concave-to-convex shape change of two long helices in the motor, thereby destroying the S-component's docking site and allowing for its dissociation. We show that adaptation of the membrane morphology to the conformational state of the motor may favour expulsion of the substrate-free S-component when ATP is bound and docking of the substrate-loaded S-component after hydrolysis. Our work provides a picture of bilayer-assisted chemo-mechanical coupling in the transport cycle of ECF transporters.

Minimal Out-of-Equilibrium Metabolism for Synthetic Cells: A Membrane Perspective

Life-like systems need to maintain a basal metabolism, which includes importing a variety of building blocks required for macromolecule synthesis, exporting dead-end products, and recycling cofactors and metabolic intermediates, while maintaining steady internal physical and chemical conditions (physicochemical homeostasis). A compartment, such as a unilamellar vesicle, functionalized with membrane-embedded transport proteins and metabolic enzymes encapsulated in the lumen meets these requirements. Here, we identify four modules designed for a minimal metabolism in a synthetic cell with a lipid bilayer boundary: energy provision and conversion, physicochemical homeostasis, metabolite transport, and membrane expansion. We review design strategies that can be used to fulfill these functions with a focus on the lipid and membrane protein composition of a cell. We compare our bottom-up design with the equivalent essential modules of JCVI-syn3a, a top-down genome-minimized living cell with a size comparable to that of large unilamellar vesicles. Finally, we discuss the bottlenecks related to the insertion of a complex mixture of membrane proteins into lipid bilayers and provide a semiquantitative estimate of the relative surface area and lipid-to-protein mass ratios (i.e., the minimal number of membrane proteins) that are required for the construction of a synthetic cell.

Different membrane order measurement techniques are not mutually consistent

"Membrane order" is a term commonly used to describe the elastic and mechanical properties of the lipid bilayer, though its exact meaning is somewhat context- and method dependent. These mechanical properties of the membrane control many cellular functions and are measured using various biophysical techniques. Here, we ask if the results obtained from various techniques are mutually consistent. Such consistency cannot be assumed a priori because these techniques probe different spatial locations and different spatial and temporal scales. We evaluate the change of membrane order induced by serotonin using nine different techniques in lipid bilayers of three different compositions. Serotonin is an important neurotransmitter present at 100s of mM concentrations in neurotransmitter vesicles, and therefore its interaction with the lipid bilayer is biologically relevant. Our measurement tools include fluorescence of lipophilic dyes (Nile Red, Laurdan, TMA-DPH, DPH), whose properties are a function of membrane order; atomic force spectroscopy, which provides a measure of the force required to indent the lipid bilayer; 2H solid-state NMR spectroscopy, which measures the molecular order of the lipid acyl chain segments; fluorescence correlation spectroscopy, which provides a measure of the diffusivity of the probe in the membrane; and Raman spectroscopy, where spectral intensity ratios are affected by acyl chain order. We find that different measures often do not correlate with each other and sometimes even yield conflicting results. We conclude that no probe provides a general measure of membrane order and that any inference based on the change of membrane order measured by a particular probe may be unreliable.

Structural Determinants and Functional Significance of Dimerization for Osmosensing Transporter ProP in Escherichia coli

Osmosensing transporter ProP forestalls cellular dehydration by detecting environments with high osmotic pressure and mediating the accumulation of organic osmolytes by bacterial cells. It is composed of 12 transmembrane helices with cytoplasmic N- and C-termini. In Escherichia coli, dimers form when the C-terminal domains of ProP molecules form homodimeric, antiparallel, α-helical coiled coils. No dominant negative effect was detected when inactive and active ProP molecules formed heterodimers in vivo. Purification of ProP in detergent dodecylmaltoside yielded monomers, which were functional after reconstitution in proteoliposomes. With other evidence, this suggests that ProP monomers function independently whether in the monomeric or dimeric state. Amino acid replacements that disrupted or reversed the coiled coil did not prevent in vivo dimerization of ProP detected with a bacterial two-hybrid system. Maleimide labeling detected no osmolality-dependent variation in the reactivities of cysteine residues introduced to transmembrane helix (TM) XII. In contrast, coarse-grained molecular dynamic simulations detected deformation of the lipid around TMs III and VI, on the lipid-exposed protein surface opposite to TM XII. This suggests that the dimer interface of ProP includes the surfaces of TMs III and VI, not of TM XII as previously suggested by crosslinking data. Homology modeling suggested that coiled-coil formation and dimerization via such an interface are not mutually exclusive. In previous work, alterations to the C-terminal coiled coil blocked co-localization of ProP with phospholipid cardiolipin at E. coli cell poles. Thus, dimerization may contribute to ProP targeting, adjust its lipid environment, and hence indirectly modify its osmotic stress response.

ABCA1 is an extracellular phospholipid translocase

Production of high density lipoprotein (HDL) requires ATP-binding cassette transporter A1 (ABCA1) to drive phospholipid (PL) from the plasma membrane into extracellular apolipoprotein A-I. Here, we use simulations to show that domains of ABCA1 within the plasma membrane remove PL from the membrane’s outer leaflet. In our simulations, after the lipid diffuses into the interior of ABCA1’s outward-open cavity, PL extracted by the gateway passes through a ring-shaped domain, the annulus orifice, which forms the base of an elongated hydrophobic tunnel in the transporter’s extracellular domain. Engineered mutations in the gateway and annulus strongly inhibit lipid export by ABCA1 without affecting cell-surface expression levels. Our finding that ABCA1 extracts lipid from the outer face of the plasma membrane and forces it through its gateway and annulus into an elongated hydrophobic tunnel contrasts with the alternating access model, which proposes that ABCA1 flops PL substrate from the inner leaflet to the outer leaflet of the membrane. Consistent with our model, ABCA1 lacks the charged amino acid residues in the transmembrane domain found in the floppase members of the ABC transporter family.

ATP-binding cassette transporter A1 (ABCA1) drives phospholipid (PL) from the plasma membrane into extracellular apolipoprotein A-I, for the production of high density lipoprotein (HDL). Here, the authors use simulations to assess the mechanism of ABCA1 function and show that ABCA1 extracts lipid from the outer face of the plasma membrane.

Membrane transport of cobalamin

A wide variety of organisms encode cobalamin-dependent enzymes catalyzing essential metabolic reactions, but the cofactor cobalamin (vitamin B12) is only synthesized by a subset of bacteria and archaea. The biosynthesis of cobalamin is complex and energetically costly, making cobalamin variants and precursors metabolically valuable. Auxotrophs for these molecules have evolved uptake mechanisms to compensate for the lack of a synthesis pathway. Bacterial transport of cobalamin involves the passage over one or two lipidic membranes in Gram-positive and -negative bacteria, respectively. In higher eukaryotes, a complex system of carriers, receptors and transporters facilitates the delivery of the essential molecule to the tissues. Biochemical and genetic approaches have identified different transporter families involved in cobalamin transport. The majority of the characterized cobalamin transporters are active transport systems that belong to the ATP-binding cassette (ABC) superfamily of transporters. In this chapter, we describe the different cobalamin transport systems characterized to date that are present in bacteria and humans, as well as yet-to-be-identified transporters.

Insights into the bilayer-mediated toppling mechanism of a folate-specific ECF transporter by cryo-EM

Energy-coupling factor (ECF)-type transporters are small, asymmetric membrane protein complexes (∼115 kDa) that consist of a membrane-embedded, substrate-binding protein (S component) and a tripartite ATP-hydrolyzing module (ECF module). They import micronutrients into bacterial cells and have been proposed to use a highly unusual transport mechanism, in which the substrate is dragged across the membrane by a toppling motion of the S component. However, it remains unclear how the lipid bilayer could accommodate such a movement. Here, we used cryogenic electron microscopy at 200 kV to determine structures of a folate-specific ECF transporter in lipid nanodiscs and detergent micelles at 2.7- and 3.4-Å resolution, respectively. The structures reveal an irregularly shaped bilayer environment around the membrane-embedded complex and suggest that toppling of the S component is facilitated by protein-induced membrane deformations. In this way, structural remodeling of the lipid bilayer environment is exploited to guide the transport process.

The Role of the Membrane in Transporter Folding and Activity

The synthesis, folding, and function of membrane transport proteins are critical factors for defining cellular physiology. Since the stability of these proteins evolved amidst the lipid bilayer, it is no surprise that we are finding that many of these membrane proteins demonstrate coupling of their structure or activity in some way to the membrane. More and more transporter structures are being determined with some information about the surrounding membrane, and computational modeling is providing further molecular details about these solvation structures. Thus, the field is moving towards identifying which molecular mechanisms - lipid interactions, membrane perturbations, differential solvation, and bulk membrane effects - are involved in linking membrane energetics to transporter stability and function. In this review, we present an overview of these mechanisms and the growing evidence that the lipid bilayer is a major determinant of the fold, form, and function of membrane transport proteins in membranes.

Coacervate formation studied by explicit solvent coarse-grain molecular dynamics with the Martini model

Complex coacervates are liquid-liquid phase separated systems, typically containing oppositely charged polyelectrolytes. They are widely studied for their functional properties as well as their potential involvement in cellular compartmentalization as biomolecular condensates. Diffusion and partitioning of solutes into a coacervate phase are important to address because their highly dynamic nature is one of their most important functional characteristics in real-world systems, but are difficult to study experimentally or even theoretically without an explicit representation of every molecule in the system. Here, we present an explicit-solvent, molecular dynamics coarse-grain model of complex coacervates, based on the Martini 3.0 force field. We demonstrate the accuracy of the model by reproducing the salt dependent coacervation of poly-lysine and poly-glutamate systems, and show the potential of the model by simulating the partitioning of ions and small nucleotides between the condensate and surrounding solvent phase. Our model paves the way for simulating coacervates and biomolecular condensates in a wide range of conditions, with near-atomic resolution.

Evaluating the Efficiency of the Martini Force Field to Study Protein Dimerization in Aqueous and Membrane Environments

Protein-protein complex assembly is one of the major drivers of biological response. Understanding the mechanisms of protein oligomerization/dimerization would allow one to elucidate how these complexes participate in biological activities and could ultimately lead to new approaches in designing novel therapeutic agents. However, determining the exact association pathways and structures of such complexes remains a challenge. Here, we use parallel tempering metadynamics simulations in the well-tempered ensemble to evaluate the performance of Martini 2.2P and Martini open-beta 3 (Martini 3) force fields in reproducing the structure and energetics of the dimerization process of membrane proteins and proteins in an aqueous solution in reasonable accuracy and throughput. We find that Martini 2.2P systematically overestimates the free energy of association by estimating large barriers in distinct areas, which likely leads to overaggregation when multiple monomers are present. In comparison, the less viscous Martini 3 results in a systematic underestimation of the free energy of association for proteins in solution, while it performs well in describing the association of membrane proteins. In all cases, the near-native dimer complexes are identified as minima in the free energy surface albeit not always as the lowest minima. In the case of Martini 3, we find that the spurious supramolecular protein aggregation present in Martini 2.2P multimer simulations is alleviated and thus this force field may be more suitable for the study of protein oligomerization. We propose that the use of enhanced sampling simulations with a refined coarse-grained force field and appropriately defined collective variables is a robust approach for studying the protein dimerization process, although one should be cautious of the ranking of energy minima.

The Groovy TMEM16 Family: Molecular Mechanisms of Lipid Scrambling and Ion Conduction

The TMEM16 family of membrane proteins displays a remarkable functional dichotomy - while some family members function as Ca²⁺-activated anion channels, the majority of characterized TMEM16 homologs are Ca²⁺-activated lipid scramblases, which catalyze the exchange of phospholipids between the two membrane leaflets. Furthermore, some TMEM16 scramblases can also function as channels. Due to their involvement in important physiological processes, the family has been actively studied ever since their molecular identity was unraveled. In this review, we will summarize the recent advances in the field and how they influenced our view of TMEM16 family function and evolution. Structural, functional and computational studies reveal how relatively small rearrangements in the permeation pathway are responsible for the observed functional duality: while TMEM16 scramblases can adopt both ion- and lipid conductive conformations, TMEM16 channels can only populate the former. Recent data further provides the molecular details of a stepwise activation mechanism, which is initiated by Ca²⁺ binding and modulated by various cellular factors, including lipids. TMEM16 function and the surrounding membrane properties are inextricably intertwined, with the protein inducing bilayer deformations associated with scrambling, while the surrounding lipids modulate TMEM16 conformation and activity.

Elevator-type mechanisms of membrane transport

Membrane transporters are integral membrane proteins that mediate the passage of solutes across lipid bilayers. These proteins undergo conformational transitions between outward- and inward-facing states, which lead to alternating access of the substrate-binding site to the aqueous environment on either side of the membrane. Dozens of different transporter families have evolved, providing a wide variety of structural solutions to achieve alternating access. A sub-set of structurally diverse transporters operate by mechanisms that are collectively named 'elevator-type'. These transporters have one common characteristic: they contain a distinct protein domain that slides across the membrane as a rigid body, and in doing so it 'drags" the transported substrate along. Analysis of the global conformational changes that take place in membrane transporters using elevator-type mechanisms reveals that elevator-type movements can be achieved in more than one way. Molecular dynamics simulations and experimental data help to understand how lipid bilayer properties may affect elevator movements and vice versa.

In vitro reconstitution of dynamically interacting integral membrane subunits of energy-coupling factor transporters

Energy-coupling factor (ECF) transporters mediate import of micronutrients in prokaryotes. They consist of an integral membrane S-component (that binds substrate) and ECF module (that powers transport by ATP hydrolysis). It has been proposed that different S-components compete for docking onto the same ECF module, but a minimal liposome-reconstituted system, required to substantiate this idea, is lacking. Here, we co-reconstituted ECF transporters for folate (ECF-FolT2) and pantothenate (ECF-PanT) into proteoliposomes, and assayed for crosstalk during active transport. The kinetics of transport showed that exchange of S-components is part of the transport mechanism. Competition experiments suggest much slower substrate association with FolT2 than with PanT. Comparison of a crystal structure of ECF-PanT with previously determined structures of ECF-FolT2 revealed larger conformational changes upon binding of folate than pantothenate, which could explain the kinetic differences. Our work shows that a minimal in vitro system with two reconstituted transporters recapitulates intricate kinetics behaviour observed in vivo.

What monomeric nucleotide binding domains can teach us about dimeric ABC proteins

The classic conceptualization of ATP binding cassette (ABC) transporter function is an ATP-dependent conformational change coupled to transport of a substrate across a biological membrane via the transmembrane domains (TMDs). The binding of two ATP molecules within the transporter's two nucleotide binding domains (NBDs) induces their dimerization. Despite retaining the ability to bind nucleotides, isolated NBDs frequently fail to dimerize. ABC proteins without a TMD, for example ABCE and ABCF, have NBDs tethered via elaborate linkers, further supporting that NBD dimerization does not readily occur for isolated NBDs. Intriguingly, even in full-length transporters, the NBD-dimerized, outward-facing state is not as frequently observed as might be expected. This leads to questions regarding what drives NBD interaction and the role of the TMDs or linkers. Understanding the NBD-nucleotide interaction and the subsequent NBD dimerization is thus pivotal for understanding ABC transporter activity in general. Here, we hope to provide new insights into ABC protein function by discussing the perplexing issue of (missing) NBD dimerization in isolation and in the context of full-length ABC proteins.