Generating Concentration Gradients across Membranes for Molecular Dynamics Simulations of Periodic Systems.

The plasma membrane forms the boundary between a living entity and its environment and acts as a barrier to permeation and flow of substances. Several computational means of calculating permeability have been implemented for molecular dynamics (MD) simulations-based approaches. Except for double bilayer systems, most permeability studies have been performed under equilibrium conditions, in large part due to the challenges associated with creating concentration gradients in simulations utilizing periodic boundary conditions. To enhance the scientific understanding of permeation and complement the existing computational means of characterizing membrane permeability, we developed a non-equilibrium method that enables the generation and maintenance of steady-state gradients in MD simulations. We utilize PBCs advantageously by imposing a directional bias to the motion of permeants so that their crossing of the boundary replenishes the gradient, like a previous study on ions. Under these conditions, a net flow of permeants across membranes may be observed to determine bulk permeability by a direct application of J=PΔc. In the present study, we explore the results of its application to an exemplary O2 and POPC bilayer system, demonstrating accurate and precise permeability measurements. In addition, we illustrate the impact of permeant concentration and the choice of thermostat on the permeability. Moreover, we demonstrate that energetics of permeation can be closely examined by the dissipation of the gradient across the membrane to gain nuanced insights into the thermodynamics of permeability.

Beyond nothingness in the formation and functional relevance of voids in polymer films.

Voids-the nothingness-broadly exist within nanomaterials and impact properties ranging from catalysis to mechanical response. However, understanding nanovoids is challenging due to lack of imaging methods with the needed penetration depth and spatial resolution. Here, we integrate electron tomography, morphometry, graph theory and coarse-grained molecular dynamics simulation to study the formation of interconnected nanovoids in polymer films and their impacts on permeance and nanomechanical behaviour. Using polyamide membranes for molecular separation as a representative system, three-dimensional electron tomography at nanometre resolution reveals nanovoid formation from coalescence of oligomers, supported by coarse-grained molecular dynamics simulations. Void analysis provides otherwise inaccessible inputs for accurate fittings of methanol permeance for polyamide membranes. Three-dimensional structural graphs accounting for the tortuous nanovoids within, measure higher apparent moduli with polyamide membranes of higher graph rigidity. Our study elucidates the significance of nanovoids beyond the nothingness, impacting the synthesis‒morphology‒function relationships of complex nanomaterials.

Atomistic Characterization of Beta-2-Glycoprotein I Domain V Interaction with Anionic Membranes.

Background: Interaction of beta-2-glycoprotein I ( ß 2 GPI) with anionic membranes is crucial in antiphospholipid syndrome (APS), implicating the role of it's membrane bind-ing domain, Domain V (DV). The mechanism of DV binding to anionic lipids is not fully understood. Objectives: This study aims to elucidate the mechanism by which DV of ß 2 GPI binds to anionic membranes. Methods: We utilized molecular dynamics (MD) simulations to investigate the struc-tural basis of anionic lipid recognition by DV. To corroborate the membrane-binding mode identified in the HMMM simulations, we conducted additional simulations using a full mem-brane model. Results: The study identified critical regions in DV, namely the lysine-rich loop and the hydrophobic loop, essential for membrane association via electrostatic and hydrophobic interactions, respectively. A novel lysine pair contributing to membrane binding was also discovered, providing new insights into ß 2 GPI's membrane interaction. Simulations revealed two distinct binding modes of DV to the membrane, with mode 1 characterized by the insertion of the hydrophobic loop into the lipid bilayer, suggesting a dominant mechanism for membrane association. This interaction is pivotal for the pathogenesis of APS, as it facilitates the recognition of ß 2 GPI by antiphospholipid antibodies. Conclusion: The study advances our understanding of the molecular interactions be-tween ß 2 GPI's DV and anionic membranes, crucial for APS pathogenesis. It highlights the importance of specific regions in DV for membrane binding and reveals a predominant bind-ing mode. These findings have significant implications for APS diagnostics and therapeutics, offering a deeper insight into the molecular basis of the syndrome.

IHMCIF: An Extension of the PDBx/mmCIF Data Standard for Integrative Structure Determination Methods.

IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB.

WSB1/2 target chromatin-bound lysine-methylated RelA for proteasomal degradation and NF-κB termination.

Proteasome-mediated degradation of chromatin-bound NF-κB is critical in terminating the transcription of pro-inflammatory genes and can be triggered by Set9-mediated lysine methylation of the RelA subunit. However, the E3 ligase targeting methylated RelA remains unknown. Here, we find that two structurally similar substrate-recognizing components of Cullin-RING E3 ligases, WSB1 and WSB2, can recognize chromatin-bound methylated RelA for polyubiquitination and proteasomal degradation. We showed that WSB1/2 negatively regulated a subset of NF-κB target genes via associating with chromatin where they targeted methylated RelA for ubiquitination, facilitating the termination of NF-κB-dependent transcription. WSB1/2 specifically interacted with methylated lysines (K) 314 and 315 of RelA via their N-terminal WD-40 repeat (WDR) domains, thereby promoting ubiquitination of RelA. Computational modeling further revealed that a conserved aspartic acid (D) at position 158 within the WDR domain of WSB2 coordinates K314/K315 of RelA, with a higher affinity when either of the lysines is methylated. Mutation of D158 abolished WSB2's ability to bind to and promote ubiquitination of methylated RelA. Together, our study identifies a novel function and the underlying mechanism for WSB1/2 in degrading chromatin-bound methylated RelA and preventing sustained NF-κB activation, providing potential new targets for therapeutic intervention of NF-κB-mediated inflammatory diseases.

Structure of alpha-synuclein fibrils derived from human Lewy body dementia tissue.

The defining feature of Parkinson disease (PD) and Lewy body dementia (LBD) is the accumulation of alpha-synuclein (Asyn) fibrils in Lewy bodies and Lewy neurites. Here we develop and validate a method to amplify Asyn fibrils extracted from LBD postmortem tissue samples and use solid state nuclear magnetic resonance (SSNMR) studies to determine atomic resolution structure. Amplified LBD Asyn fibrils comprise a mixture of single protofilament and two protofilament fibrils with very low twist. The protofilament fold is highly similar to the fold determined by a recent cryo-electron microscopy study for a minority population of twisted single protofilament fibrils extracted from LBD tissue. These results expand the structural characterization of LBD Asyn fibrils and approaches for studying disease mechanisms, imaging agents and therapeutics targeting Asyn.

Elucidating the Mechanism of Metabolism of Cannabichromene by Human Cytochrome P450s.

Cannabichromene (CBC) is a nonpsychoactive phytocannabinoid well-known for its wide-ranging health advantages. However, there is limited knowledge regarding its human metabolism following CBC consumption. This research aimed to explore the metabolic pathways of CBC by various human liver cytochrome P450 (CYP) enzymes and support the outcomes using in vivo data from mice. The results unveiled two principal CBC metabolites generated by CYPs: 8'-hydroxy-CBC and 6',7'-epoxy-CBC, along with a minor quantity of 1″-hydroxy-CBC. Notably, among the examined CYPs, CYP2C9 demonstrated the highest efficiency in producing these metabolites. Moreover, through a molecular dynamics simulation spanning 1 µs, it was observed that CBC attains stability at the active site of CYP2J2 by forming hydrogen bonds with I487 and N379, facilitated by water molecules, which specifically promotes the hydroxy metabolite's formation. Additionally, the presence of cytochrome P450 reductase (CPR) amplified CBC's binding affinity to CYPs, particularly with CYP2C8 and CYP3A4. Furthermore, the metabolites derived from CBC reduced cytokine levels, such as IL6 and NO, by approximately 50% in microglia cells. This investigation offers valuable insights into the biotransformation of CBC, underscoring the physiological importance and the potential significance of these metabolites.

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.

Integrative analysis of pathogenic variants in glucose-6-phosphatase based on an AlphaFold2 model.

Mediating the terminal reaction of gluconeogenesis and glycogenolysis, the integral membrane protein glucose-6-phosphate catalytic subunit 1 (G6PC1) regulates hepatic glucose production by catalyzing hydrolysis of glucose-6-phosphate (G6P) within the lumen of the endoplasmic reticulum. Consistent with its vital contribution to glucose homeostasis, inactivating mutations in G6PC1 causes glycogen storage disease (GSD) type 1a characterized by hepatomegaly and severe hypoglycemia. Despite its physiological importance, the structural basis of G6P binding to G6PC1 and the molecular disruptions induced by missense mutations within the active site that give rise to GSD type 1a are unknown. In this study, we determine the atomic interactions governing G6P binding as well as explore the perturbations imposed by disease-linked missense variants by subjecting an AlphaFold2 G6PC1 structural model to molecular dynamics simulations and in silico predictions of thermodynamic stability validated with robust in vitro and in situ biochemical assays. We identify a collection of side chains, including conserved residues from the signature phosphatidic acid phosphatase motif, that contribute to a hydrogen bonding and van der Waals network stabilizing G6P in the active site. The introduction of GSD type 1a mutations modified the thermodynamic landscape, altered side chain packing and substrate-binding interactions, and induced trapping of catalytic intermediates. Our results, which corroborate the high quality of the AF2 model as a guide for experimental design and to interpret outcomes, not only confirm the active-site structural organization but also identify previously unobserved mechanistic contributions of catalytic and noncatalytic side chains.

A generative artificial intelligence framework based on a molecular diffusion model for the design of metal-organic frameworks for carbon capture.

Metal-organic frameworks (MOFs) exhibit great promise for CO2 capture. However, finding the best performing materials poses computational and experimental grand challenges in view of the vast chemical space of potential building blocks. Here, we introduce GHP-MOFassemble, a generative artificial intelligence (AI), high performance framework for the rational and accelerated design of MOFs with high CO2 adsorption capacity and synthesizable linkers. GHP-MOFassemble generates novel linkers, assembled with one of three pre-selected metal nodes (Cu paddlewheel, Zn paddlewheel, Zn tetramer) into MOFs in a primitive cubic topology. GHP-MOFassemble screens and validates AI-generated MOFs for uniqueness, synthesizability, structural validity, uses molecular dynamics simulations to study their stability and chemical consistency, and crystal graph neural networks and Grand Canonical Monte Carlo simulations to quantify their CO2 adsorption capacities. We present the top six AI-generated MOFs with CO2 capacities greater than 2m mol g-1, i.e., higher than 96.9% of structures in the hypothetical MOF dataset.

The carboxy terminus causes interfacial assembly of oleate hydratase on a membrane bilayer.

The soluble flavoprotein oleate hydratase (OhyA) hydrates the 9-cis double bond of unsaturated fatty acids. OhyA substrates are embedded in membrane bilayers; OhyA must remove the fatty acid from the bilayer and enclose it in the active site. Here, we show that the positively charged helix-turn-helix motif in the carboxy terminus (CTD) is responsible for interacting with the negatively charged phosphatidylglycerol (PG) bilayer. Super-resolution microscopy of Staphylococcus aureus cells expressing green fluorescent protein fused to OhyA or the CTD sequence shows subcellular localization along the cellular boundary, indicating OhyA is membrane-associated and the CTD sequence is sufficient for membrane recruitment. Using cryo-electron microscopy, we solved the OhyA dimer structure and conducted 3D variability analysis of the reconstructions to assess CTD flexibility. Our surface plasmon resonance experiments corroborated that OhyA binds the PG bilayer with nanomolar affinity and we found the CTD sequence has intrinsic PG binding properties. We determined that the nuclear magnetic resonance structure of a peptide containing the CTD sequence resembles the OhyA crystal structure. We observed intermolecular NOE from PG liposome protons next to the phosphate group to the CTD peptide. The addition of paramagnetic MnCl2 indicated the CTD peptide binds the PG surface but does not insert into the bilayer. Molecular dynamics simulations, supported by site-directed mutagenesis experiments, identify key residues in the helix-turn-helix that drive membrane association. The data show that the OhyA CTD binds the phosphate layer of the PG surface to obtain bilayer-embedded unsaturated fatty acids.

A computational spatial whole-Cell model for hepatitis B viral infection and drug interactions.

Despite a vaccine, hepatitis B virus (HBV) remains a world-wide source of infections and deaths. We develop a whole-cell computational platform combining spatial and kinetic models describing the infection cycle of HBV in a hepatocyte host. We simulate key parts of the infection cycle with this whole-cell platform for 10 min of biological time, to predict infection progression, map out virus-host and virus-drug interactions. We find that starting from an established infection, decreasing the copy number of the viral envelope proteins shifts the dominant infection pathway from capsid secretion to re-importing the capsids into the nucleus, resulting in more nuclear-localized viral covalently closed circular DNA (cccDNA) and boosting transcription. This scenario can mimic the consequence of drugs designed to manipulate viral gene expression. Mutating capsid proteins facilitates capsid destabilization and disassembly at nuclear pore complexes, resulting in an increase in cccDNA copy number. However, excessive destabilization leads to premature cytoplasmic disassembly and does not increase the cccDNA counts. Finally, our simulations can predict the best drug dosage and its administration timing to reduce the cccDNA counts. Our adaptable computational platform can be parameterized to study other viruses and identify the most central viral pathways that can be targeted by drugs.

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes.

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

Topological Learning Approach to Characterizing Biological Membranes.

Biological membranes play key roles in cellular compartmentalization, structure, and its signaling pathways. At varying temperatures, individual membrane lipids sample from different configurations, a process that frequently leads to higher-order phase behavior and phenomena. Here we present a persistent homology-based method for quantifying the structural features of individual and bulk lipids, providing local and contextual information on lipid tail organization. Our method leverages the mathematical machinery of algebraic topology and machine learning to infer temperature-dependent structural information of lipids from static coordinates. To train our model, we generated multiple molecular dynamics trajectories of DPPC membranes at varying temperatures. A fingerprint was then constructed for each set of lipid coordinates by a persistent homology filtration, in which interactions spheres were grown around the lipid atoms while tracking their intersections. The sphere filtration formed a simplicial complex that captures enduring key topological features of the configuration landscape, using homology, yielding persistence data. Following fingerprint extraction for physiologically relevant temperatures, the persistence data were used to train an attention-based neural network for assignment of effective temperature values to selected membrane regions. Our persistence homology-based method captures the local structural effects, via effective temperature, of lipids adjacent to other membrane constituents, e.g. sterols and proteins. This topological learning approach can predict lipid effective temperatures from static coordinates across multiple spatial resolutions. The tool, called MembTDA, can be accessed at https://github.com/hyunp2/Memb-TDA.

Myr-Arf1 conformational flexibility at the membrane surface sheds light on the interactions with ArfGAP ASAP1.

ADP-ribosylation factor 1 (Arf1) interacts with multiple cellular partners and membranes to regulate intracellular traffic, organelle structure and actin dynamics. Defining the dynamic conformational landscape of Arf1 in its active form, when bound to the membrane, is of high functional relevance and key to understanding how Arf1 can alter diverse cellular processes. Through concerted application of nuclear magnetic resonance (NMR), neutron reflectometry (NR) and molecular dynamics (MD) simulations, we show that, while Arf1 is anchored to the membrane through its N-terminal myristoylated amphipathic helix, the G domain explores a large conformational space, existing in a dynamic equilibrium between membrane-associated and membrane-distal conformations. These configurational dynamics expose different interfaces for interaction with effectors. Interaction with the Pleckstrin homology domain of ASAP1, an Arf-GTPase activating protein (ArfGAP), restricts motions of the G domain to lock it in what seems to be a conformation exposing functionally relevant regions.

Asymmetric conformations and lipid interactions shape the ATP-coupled cycle of a heterodimeric ABC transporter.

Here we used cryo-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy (DEER), and molecular dynamics (MD) simulations, to capture and characterize ATP- and substrate-bound inward-facing (IF) and occluded (OC) conformational states of the heterodimeric ATP binding cassette (ABC) multidrug exporter BmrCD in lipid nanodiscs. Supported by DEER analysis, the structures reveal that ATP-powered isomerization entails changes in the relative symmetry of the BmrC and BmrD subunits that propagates from the transmembrane domain to the nucleotide binding domain. The structures uncover asymmetric substrate and Mg2+ binding which we hypothesize are required for triggering ATP hydrolysis preferentially in one of the nucleotide-binding sites. MD simulations demonstrate that multiple lipid molecules differentially bind the IF versus the OC conformation thus establishing that lipid interactions modulate BmrCD energy landscape. Our findings are framed in a model that highlights the role of asymmetric conformations in the ATP-coupled transport with general implications to the mechanism of ABC transporters.

Asymmetric conformations and lipid interactions shape the ATP-coupled cycle of a heterodimeric ABC transporter.

To illuminate the structural origin of catalytic asymmetry of heterodimeric ABC transporters and how it shapes the energetics of their conformational cycles, we used cryo-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy (DEER), and molecular dynamics (MD) simulations, to capture and characterize conformational states of the heterodimeric ABC multidrug exporter BmrCD in lipid nanodiscs. In addition to multiple ATP- and substrate-bound inward-facing (IF) conformations, we obtained the structure of an occluded (OC) conformation wherein the unique extracellular domain (ECD) twists to partially open the extracellular gate. In conjunction with DEER analysis of the populations of these conformations, the structures reveal that ATP-powered isomerization entails changes in the relative symmetry of the BmrC and BmrD subunits that propagates from the transmembrane domain (TMD) to the nucleotide binding domain (NBD). The structures uncover asymmetric substrate and Mg 2+ binding which we hypothesize are required for triggering ATP hydrolysis preferentially in one of the nucleotide-binding sites. MD simulations demonstrated that multiple lipid molecules, identified from the cryo-EM density maps, differentially bind the IF versus the OC conformation thus modulating their relative stability. In addition to establishing how lipid interactions with BmrCD modulate the energy landscape, our findings are framed in a distinct transport model that highlights the role of asymmetric conformations in the ATP-coupled cycle with implications to the mechanism of ABC transporters in general.

Structure Reveals Homology in Elevator Transporters.

The elevator transport mechanism is one of the handful of canonical mechanisms by which transporters shuttle their substrates across the semi-permeable membranes that surround cells and organelles. Studies of molecular function are naturally guided by evolutionary context, but until now this context has been limited for elevator transporters because established evolutionary classification methods have organized them into several apparently unrelated families. Through comprehensive examination of the pertinent structures available in the Protein Data Bank, we show that 62 elevator transporters from 18 families share a conserved architecture in their transport domains consisting of 10 helices connected in 8 topologies. Through quantitative analysis of the structural similarity, structural complexity, and topologically-corrected sequence similarity among the transport domains, we provide compelling evidence that these elevator transporters are all homologous. Using our analysis, we have constructed a phylogenetic tree to enable quantification and visualization of the evolutionary relationships among elevator transporters and their families. We also report several examples of functional features that are shared by elevator transporters from different families. Our findings shed new light on the elevator transport mechanism and allow us to understand it in a far deeper and more nuanced manner.

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations.

Modeling and simulation of small molecules such as drugs and biological cofactors have been both a major focus of computational chemistry for decades and a growing need among computational biophysicists who seek to investigate the interaction of different types of ligands with biomolecules. Of particular interest in this regard are quantum mechanical (QM) calculations that are used to more accurately describe such small molecules, which can be of heterogeneous structures and chemistry, either in purely QM calculations or in hybrid QM/molecular mechanics (MM) simulations. QM programs are also used to develop MM force field parameters for small molecules to be used along with established force fields for biomolecules in classical simulations. With this growing need in mind, here we report a set of software tools developed and closely integrated within the broadly used molecular visualization/analysis program, VMD, that allow the user to construct, modify, and parametrize small molecules and prepare them for QM, hybrid QM/MM, or classical simulations. The tools also provide interactive analysis and visualization capabilities in an easy-to-use and integrated environment. In this paper, we briefly report on these tools and their major features and capabilities, along with examples of how they can facilitate molecular research in computational biophysics that might be otherwise prohibitively complex.

Structures and membrane interactions of native serotonin transporter in complexes with psychostimulants.

The serotonin transporter (SERT) is a member of the SLC6 neurotransmitter transporter family that mediates serotonin reuptake at presynaptic nerve terminals. SERT is the target of both therapeutic antidepressant drugs and psychostimulant substances such as cocaine and methamphetamines, which are small molecules that perturb normal serotonergic transmission by interfering with serotonin transport. Despite decades of studies, important functional aspects of SERT such as the oligomerization state of native SERT and its interactions with potential proteins remain unresolved. Here, we develop methods to isolate SERT from porcine brain (pSERT) using a mild, nonionic detergent, utilize fluorescence-detection size-exclusion chromatography to investigate its oligomerization state and interactions with other proteins, and employ single-particle cryo-electron microscopy to elucidate the structures of pSERT in complexes with methamphetamine or cocaine, providing structural insights into psychostimulant recognition and accompanying pSERT conformations. Methamphetamine and cocaine both bind to the central site, stabilizing the transporter in an outward open conformation. We also identify densities attributable to multiple cholesterol or cholesteryl hemisuccinate (CHS) molecules, as well as to a detergent molecule bound to the pSERT allosteric site. Under our conditions of isolation, we find that pSERT is best described as a monomeric entity, isolated without interacting proteins, and is ensconced by multiple cholesterol or CHS molecules.