Molecular mechanism of prestin electromotive signal amplification.

Hearing involves two fundamental processes: mechano-electrical transduction and signal amplification. Despite decades of studies, the molecular bases for both remain elusive. Here, we show how prestin, the electromotive molecule of outer hair cells (OHCs) that senses both voltage and membrane tension, mediates signal amplification by coupling conformational changes to alterations in membrane surface area. Cryoelectron microscopy (cryo-EM) structures of human prestin bound with chloride or salicylate at a common "anion site" adopt contracted or expanded states, respectively. Prestin is ensconced within a perimeter of well-ordered lipids, through which it induces dramatic deformation in the membrane and couples protein conformational changes to the bulk membrane. Together with computational studies, we illustrate how the anion site is allosterically coupled to changes in the transmembrane domain cross-sectional area and the surrounding membrane. These studies provide insight into OHC electromotility by providing a structure-based mechanism of the membrane motor prestin.

Structures of the TMC-1 complex illuminate mechanosensory transduction.

The initial step in the sensory transduction pathway underpinning hearing and balance in mammals involves the conversion of force into the gating of a mechanosensory transduction channel1. Despite the profound socioeconomic impacts of hearing disorders and the fundamental biological significance of understanding mechanosensory transduction, the composition, structure and mechanism of the mechanosensory transduction complex have remained poorly characterized. Here we report the single-particle cryo-electron microscopy structure of the native transmembrane channel-like protein 1 (TMC-1) mechanosensory transduction complex isolated from Caenorhabditis elegans. The two-fold symmetric complex is composed of two copies each of the pore-forming TMC-1 subunit, the calcium-binding protein CALM-1 and the transmembrane inner ear protein TMIE. CALM-1 makes extensive contacts with the cytoplasmic face of the TMC-1 subunits, whereas the single-pass TMIE subunits reside on the periphery of the complex, poised like the handles of an accordion. A subset of complexes additionally includes a single arrestin-like protein, arrestin domain protein (ARRD-6), bound to a CALM-1 subunit. Single-particle reconstructions and molecular dynamics simulations show how the mechanosensory transduction complex deforms the membrane bilayer and suggest crucial roles for lipid-protein interactions in the mechanism by which mechanical force is transduced to ion channel gating.

Cryo-EM structures of human ABCA7 provide insights into its phospholipid translocation mechanisms.

Phospholipid extrusion by ABC subfamily A (ABCA) exporters is central to cellular physiology, although the specifics of the underlying substrate interactions and transport mechanisms remain poorly resolved at the molecular level. Here we report cryo-EM structures of lipid-embedded human ABCA7 in an open state and in a nucleotide-bound, closed state at resolutions between 3.6 and 4.0 Å. The former reveals an ordered patch of bilayer lipids traversing the transmembrane domain (TMD), while the latter reveals a lipid-free, closed TMD with a small extracellular opening. These structures offer a structural framework for both substrate entry and exit from the ABCA7 TMD and highlight conserved rigid-body motions that underlie the associated conformational transitions. Combined with functional analysis and molecular dynamics (MD) simulations, our data also shed light on lipid partitioning into the ABCA7 TMD and localized membrane perturbations that underlie ABCA7 function and have broader implications for other ABCA family transporters.

Differential dynamics and direct interaction of bound ligands with lipids in multidrug transporter ABCG2.

ABCG2 is an ATP-binding cassette (ABC) transporter that extrudes a wide range of xenobiotics and drugs from the cell and contributes to multidrug resistance in cancer cells. Following our recent structural characterization of topotecan-bound ABCG2, here, we present cryo-EM structures of ABCG2 under turnover conditions in complex with a special modulator and slow substrate, tariquidar, in nanodiscs. The structures reveal that similar to topotecan, tariquidar induces two distinct ABCG2 conformations under turnover conditions (turnover-1 and turnover-2). µs-scale molecular dynamics simulations of drug-bound and apo ABCG2 in native-like lipid bilayers, in both topotecan- and tariquidar-bound states, characterize the ligand size as a major determinant of its binding stability. The simulations highlight direct lipid-drug interactions for the smaller topotecan, which exhibits a highly dynamic binding mode. In contrast, the larger tariquidar occupies most of the available volume in the binding pocket, thus leaving little space for lipids to enter the cavity and interact with it. Similarly, when simulating ABCG2 in the apo inward-open state, we also observe spontaneous penetration of phospholipids into the binding cavity. The captured phospholipid diffusion pathway into ABCG2 offers a putative general path to recruit any hydrophobic/amphiphilic substrates directly from the membrane. Our simulations also reveal that ABCG2 rejects cholesterol as a substrate, which is omnipresent in plasma membranes that contain ABCG2. At the same time, cholesterol is found to prohibit the penetration of phospholipids into ABCG2. These molecular findings have direct functional ramifications on ABCG2's function as a transporter.

Tumor-acquired somatic mutation affects conformation to abolish ABCG2-mediated drug resistance.

ABCG2 is an important ATP-binding cassette transporter impacting the absorption and distribution of over 200 chemical toxins and drugs. ABCG2 also reduces the cellular accumulation of diverse chemotherapeutic agents. Acquired somatic mutations in the phylogenetically conserved amino acids of ABCG2 might provide unique insights into its molecular mechanisms of transport. Here, we identify a tumor-derived somatic mutation (Q393K) that occurs in a highly conserved amino acid across mammalian species. This ABCG2 mutant seems incapable of providing ABCG2-mediated drug resistance. This was perplexing because it is localized properly and retained interaction with substrates and nucleotides. Using a conformationally sensitive antibody, we show that this mutant appears "locked" in a non-functional conformation. Structural modeling and molecular dynamics simulations based on ABCG2 cryo-EM structures suggested that the Q393K interacts with the E446 to create a strong salt bridge. The salt bridge is proposed to stabilize the inward-facing conformation, resulting in an impaired transporter that lacks the flexibility to readily change conformation, thereby disrupting the necessary communication between substrate binding and transport.

Membrane Mixer: A Toolkit for Efficient Shuffling of Lipids in Heterogeneous Biological Membranes.

Molecular dynamics (MD) simulations of biological membranes have achieved such levels of sophistication that are commonly used to predict unresolved structures and various properties of lipids and to substantiate experimental data. While achieving sufficient sampling of lipid dynamics remains a major challenge, a commonly used method to improve lipid sampling, e.g., in terms of specific interactions with membrane-associated proteins, is to randomize the initial arrangement of lipid constituents in multiple replicas of simulations, without changing the overall lipid composition of the membrane of interest. Here, we introduce a method that can rapidly generate multiple replicas of lipid bilayers with different spatial and conformational configurations for any given lipid composition. The underlying algorithm, which allows one to shuffle lipids at any desired level, relies on the application of an external potential, here referred to as the "carving potential", that removes clashes/entanglements before lipid positions are exchanged (shuffled), thereby minimizing the energy penalty due to abrupt lipid repositioning. The method is implemented as "Membrane Mixer Plugin (MMP) 1.0" in VMD, with a convenient graphical user interface that guides the user in setting various options and parameters. The plugin is fully automated and generates new membrane replicas more rapidly and conveniently than other analogous tools. The plugin and its capabilities introduced here can be extended to include additional features in future versions.

Complex self-propelled rings: a minimal model for cell motility.

Collective behavior of active matter is observed for self-propelled particles, such as vibrated disks and active Brownian particles, as well as for cytoskeletal filaments in motile cells. Here, a system of quasi two-dimensional penetrable self-propelled rods inside rigid rings is used to construct a complex self-propelled particle. The rods interact sterically with each other and with a stationary or mobile ring via a separation-shifted Lennard-Jones potential. They either have a sliding attachment to the inside of the ring at one of their ends, or can move freely within the ring confinement. We study the inner structure and dynamics of the mobile self-propelled rings. We find that these complex particles cannot only be characterized as active Brownian particles, but can also exhibit cell-like motility: random walks, persistent motion, circling, and run-and-circle motion.

Dissecting the Interaction Fingerprints and Binding Affinity of BYL719 Analogs Targeting PI3Kα.

Phosphatidylinositol-3-kinase Alpha (PI3Kα) is a lipid kinase which regulates signaling pathways involved in cell proliferation. Dysregulation of these pathways promotes several human cancers, pushing for the development of anticancer drugs to target PI3Kα. One such medicinal chemistry campaign at Novartis led to the discovery of BYL719 (Piqray, Alpelicib), a PI3Kα inhibitor approved by the FDA in 2019 for treatment of HR+/HER2-advanced breast cancer with a PIK3CA mutation. Structure-based drug design played a key role in compound design and optimization throughout the discovery process. However, further characterization of potency drivers via structural dynamics and energetic analyses can be advantageous for ensuing PI3Kα programs. Here, our goal is to employ various in-silico techniques, including molecular simulations and machine learning, to characterize 14 ligands from the BYL719 analogs and predict their binding affinities. The structural insights from molecular simulations suggest that although the ligand-hinge interaction is the primary driver of ligand stability at the pocket, the R group positioning at C2 or C6 of pyridine/pyrimidine also plays a major role. Binding affinities predicted via thermodynamic integration (TI) are in good agreement with previously reported IC50s. Yet, computationally demanding techniques such as TI might not always be the most efficient approach for affinity prediction, as in our case study, fast high-throughput techniques were capable of classifying compounds as active or inactive, and one docking approach showed accuracy comparable to TI.

Ionizable Amino Lipids Distribution and Effects on DSPC/Cholesterol Membranes: Implications for Lipid Nanoparticle Structure.

Lipid nanoparticles (LNPs) containing ionizable aminolipids are among the leading platforms for the successful delivery of nucleic-acid-based therapeutics, including messenger RNA (mRNA). The two recently FDA-approved COVID-19 vaccines developed by Moderna and Pfizer/BioNTech belong to this category. Ionizable aminolipids, cholesterol, and DSPC lipids are among the key components of such formulations, crucially modulating physicochemical properties of these formulations and, consequently, the potency of these therapeutics. Despite the importance of these components, the distribution of these molecules in LNPs containing mRNA is not clear. In this study, we used all-atom molecular dynamics (MD) simulations to investigate the distribution and effects of the Lipid-5 (apparent pKa of the lipid nanoparticle = 6.56), a rationally designed and previously reported ionizable aminolipid by Moderna, on lipid bilayers [Mol. Ther. 2018, 26, 1509-1519]. The simulations were conducted with half of the aminolipids charged and half neutral approximately to the expected ionization in the microenvironment of the LNP surface. In all five simulated systems in this work, the cholesterol content was kept constant, whereas the DSPC and Lipid-5 concentrations were changed systematically. We found that at higher concentrations of the ionizable aminolipids, the neutral aminolipids form a disordered aggregate in the membrane interior that preferentially includes cholesterol. The rules underlying the lipid redistribution could be used to rationally choose lipids to optimize the LNP function.

Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications.

Protein lipidations are vital co/post-translational modifications that tether lipid tails to specific protein amino acids, allowing them to anchor to biological membranes, switch their subcellular localization, and modulate association with other proteins. Such lipidations are thus crucial for multiple biological processes including signal transduction, protein trafficking, and membrane localization and are implicated in various diseases as well. Examples of lipid-anchored proteins include the Ras family of proteins that undergo farnesylation; actin and gelsolin that are myristoylated; phospholipase D that is palmitoylated; glycosylphosphatidylinositol-anchored proteins; and others. Here, we develop parameters for cysteine-targeting farnesylation, geranylgeranylation, and palmitoylation, as well as glycine-targeting myristoylation for the latest version of the Martini 3 coarse-grained force field. The parameters are developed using the CHARMM36m all-atom force field parameters as reference. The behavior of the coarse-grained models is consistent with that of the all-atom force field for all lipidations and reproduces key dynamical and structural features of lipid-anchored peptides, such as the solvent-accessible surface area, bilayer penetration depth, and representative conformations of the anchors. The parameters are also validated in simulations of the lipid-anchored peripheral membrane proteins Rheb and Arf1, after comparison with independent all-atom simulations. The parameters, along with mapping schemes for the popular martinize2 tool, are available for download at 10.5281/zenodo.7849262 and also as supporting information.

Lipid-mediated prestin organization in outer hair cell membranes and its implications in sound amplification.

Prestin is a high-density motor protein in the outer hair cells (OHCs), whose conformational response to acoustic signals alters the shape of the cell, thereby playing a major role in sound amplification by the cochlea. Despite recent structures, prestin's intimate interactions with the membrane, which are central to its function remained unresolved. Here, employing a large set (collectively, more than 0.5 ms) of coarse-grained molecular dynamics simulations, we demonstrate the impact of prestin's lipid-protein interactions on its organization at densities relevant to the OHCs and its effectiveness in reshaping OHCs. Prestin causes anisotropic membrane deformation, which mediates a preferential membrane organization of prestin where deformation patterns by neighboring copies are aligned constructively. The resulting reduced membrane rigidity is hypothesized to maximize the impact of prestin on OHC reshaping. These results demonstrate a clear case of protein-protein cooperative communication in membrane, purely mediated by interactions with lipids.

Proton-driven alternating access in a spinster lipid transporter.

Spinster (Spns) lipid transporters are critical for transporting sphingosine-1-phosphate (S1P) across cellular membranes. In humans, Spns2 functions as the main S1P transporter in endothelial cells, making it a potential drug target for modulating S1P signaling. Here, we employed an integrated approach in lipid membranes to identify unknown conformational states of a bacterial Spns from Hyphomonas neptunium (HnSpns) and to define its proton- and substrate-coupled conformational dynamics. Our systematic study reveals conserved residues critical for protonation steps and their regulation, and how sequential protonation of these proton switches coordinates the conformational transitions in the context of a noncanonical ligand-dependent alternating access. A conserved periplasmic salt bridge (Asp60TM2:Arg289TM7) keeps the transporter in a closed conformation, while proton-dependent conformational dynamics are significantly enhanced on the periplasmic side, providing a pathway for ligand exchange.

Backbone amides are determinants of Cl- selectivity in CLC ion channels.

Chloride homeostasis is regulated in all cellular compartments. CLC-type channels selectively transport Cl- across biological membranes. It is proposed that side-chains of pore-lining residues determine Cl- selectivity in CLC-type channels, but their spatial orientation and contributions to selectivity are not conserved. This suggests a possible role for mainchain amides in selectivity. We use nonsense suppression to insert α-hydroxy acids at pore-lining positions in two CLC-type channels, CLC-0 and bCLC-k, thus exchanging peptide-bond amides with ester-bond oxygens which are incapable of hydrogen-bonding. Backbone substitutions functionally degrade inter-anion discrimination in a site-specific manner. The presence of a pore-occupying glutamate side chain modulates these effects. Molecular dynamics simulations show backbone amides determine ion energetics within the bCLC-k pore and how insertion of an α-hydroxy acid alters selectivity. We propose that backbone-ion interactions are determinants of Cl- specificity in CLC channels in a mechanism reminiscent of that described for K+ channels.

Conformational changes in the nucleotide-binding domains of P-glycoprotein induced by ATP hydrolysis.

P-glycoprotein (Pgp) is a member of the ABC transporter superfamily with high physiological importance. Pgp nucleotide-binding domains (NBDs) drive the transport cycle through ATP binding and hydrolysis. We use molecular dynamics simulations to investigate the ATP hydrolysis-induced conformational changes in NBDs. Five systems, including all possible ATP/ADP combinations in the NBDs and the APO system, are simulated. ATP/ADP exchange induces conformational changes mostly within the conserved signature motif of the NBDs, resulting in relative orientational changes in the NBDs. Nucleotide removal leads to additional orientational changes in the NBDs, allowing their dissociation. Furthermore, we capture putative hydrolysis-competent configurations in which the conserved glutamate in the Walker-B motif acts as a catalytic base capturing a water molecule likely initiating ATP hydrolysis.