The substrate import mechanism of the human serotonin transporter

Structural Models of Human Norepinephrine Transporter Ensemble Reveal the Allosteric Sites and Ligand-Binding Mechanism

The norepinephrine transporter (NET) plays a pivotal role in recycling norepinephrine (NE) from the synaptic cleft. However, the structures referring to the conformational heterogeneity of NET during the transport cycle remain poorly understood. Here, three structural models of NE bound to the orthosteric site of NET in outward-open (OOholo), outward-occluded (OCholo), and inward-open (IOholo) conformations were first obtained using the multistate structures of serotonin transporter as templates and further characterized through Gaussian-accelerated molecular dynamics and free energy reweighting. Analysis of the structures revealed eight potential allosteric sites on the functional-specific states of NET. One of the pharmacologically relevant pockets located at the extracellular vestibule was further verified by simulating the binding behaviors of a clinical trial drug χ-MrIA that is allosterically regulating NET. These structural and energetic insights into NET advanced our understanding of NE reuptake and paved the way for discovering novel molecules targeting the allosteric sites.

Free energy profile of the substrate-induced occlusion of the human serotonin transporter

The serotonin transporter (SERT) is a member of the Solute Carrier 6 (SLC6) family and is responsible for maintaining the appropriate level of serotonin in the brain. Dysfunction of SERT has been linked to several neuropsychiatric disorders, including depression, anxiety and obsessive-compulsive disorder. Therefore, an in-depth understanding of the mechanism on an atomistic level, coupled with a quantification of transporter dynamics and the associated free energies is required. Here, we constructed Markov state models (MSMs) from extensive unbiased molecular dynamics simulations to quantify the free energy profile of serotonin (5HT) triggered SERT occlusion and explored the driving forces of the mechanism of occlusion. Our results reveal that SERT occludes via multiple intermediate conformations and show that the motion of occlusion is energetically downhill for the 5HT-bound transporter. Force distribution analyses show that the interactions of 5HT with the bundle domain are crucial. During occlusion, attractive forces steadily increase and pull on the bundle domain, which leads to SERT occlusion. Some interactions become repulsive upon full occlusion, suggesting that SERT creates pressure on 5HT to promote its movement towards the cytosol.

APACE: AlphaFold2 and advanced computing as a service for accelerated discovery in biophysics

The prediction of protein 3D structure from amino acid sequence is a computational grand challenge in biophysics and plays a key role in robust protein structure prediction algorithms, from drug discovery to genome interpretation. The advent of AI models, such as AlphaFold, is revolutionizing applications that depend on robust protein structure prediction algorithms. To maximize the impact, and ease the usability, of these AI tools we introduce APACE, AlphaFold2 and advanced computing as a service, a computational framework that effectively handles this AI model and its TB-size database to conduct accelerated protein structure prediction analyses in modern supercomputing environments. We deployed APACE in the Delta and Polaris supercomputers and quantified its performance for accurate protein structure predictions using four exemplar proteins: 6AWO, 6OAN, 7MEZ, and 6D6U. Using up to 300 ensembles, distributed across 200 NVIDIA A100 GPUs, we found that APACE is up to two orders of magnitude faster than off-the-self AlphaFold2 implementations, reducing time-to-solution from weeks to minutes. This computational approach may be readily linked with robotics laboratories to automate and accelerate scientific discovery.

The Arabidopsis AtSWEET13 transporter discriminates sugars by selective facial and positional substrate recognition

Transporters are targeted by endogenous metabolites and exogenous molecules to reach cellular destinations, but it is generally not understood how different substrate classes exploit the same transporter's mechanism. Any disclosure of plasticity in transporter mechanism when treated with different substrates becomes critical for developing general selectivity principles in membrane transport catalysis. Using extensive molecular dynamics simulations with an enhanced sampling approach, we select the Arabidopsis sugar transporter AtSWEET13 as a model system to identify the basis for glucose versus sucrose molecular recognition and transport. Here we find that AtSWEET13 chemical selectivity originates from a conserved substrate facial selectivity demonstrated when committing alternate access, despite mono-/di-saccharides experiencing differing degrees of conformational and positional freedom throughout other stages of transport. However, substrate interactions with structural hallmarks associated with known functional annotations can help reinforce selective preferences in molecular transport.

Energetics of substrate transport in proton-dependent oligopeptide transporters

The PepT So transporter mediates the transport of peptides across biological membranes. Despite advancements in structural biology, including cryogenic electron microscopy structures resolving PepT So in different states, the molecular basis of peptide recognition and transport by PepT So is not fully elucidated. In this study, we employed molecular dynamics simulations, Markov State Models (MSMs), and Transition Path Theory (TPT) to investigate the transport mechanism of an alanine-alanine peptide (Ala-Ala) through the PepT So transporter. Our simulations revealed conformational changes and key intermediate states involved in peptide translocation. We observed that the presence of the Ala-Ala peptide substrate lowers the free energy barriers associated with transition to the inward-facing state. Furthermore, we elucidated the proton transport model and analyzed the pharmacophore features of intermediate states, providing insights for rational drug design. These findings highlight the significance of substrate binding in modulating the conformational dynamics of PepT So and identify critical residues that facilitate transport.

Identification of the potassium-binding site in serotonin transporter

Clearance of serotonin (5-hydroxytryptamine, 5-HT) from the synaptic cleft after neuronal signaling is mediated by serotonin transporter (SERT), which couples this process to the movement of a Na+ ion down its chemical gradient. After release of 5-HT and Na+ into the cytoplasm, the transporter faces a rate-limiting challenge of resetting its conformation to be primed again for 5-HT and Na+ binding. Early studies of vesicles containing native SERT revealed that K+ gradients can provide an additional driving force, via K+ antiport. Moreover, under appropriate conditions, a H+ ion can replace K+. Intracellular K+ accelerates the resetting step. Structural studies of SERT have identified two binding sites for Na+ ions, but the K+ site remains enigmatic. Here, we show that K+ antiport can drive substrate accumulation into vesicles containing SERT extracted from a heterologous expression system, allowing us to study the residues responsible for K+ binding. To identify candidate binding residues, we examine many cation binding configurations using molecular dynamics simulations, predicting that K+ binds to the so-called Na2 site. Site-directed mutagenesis of residues in this site can eliminate the ability of both K+ and H+ to drive 5-HT accumulation into vesicles and, in patch clamp recordings, prevent the acceleration of turnover rates and the formation of a channel-like state by K+ or H+. In conclusion, the Na2 site plays a pivotal role in orchestrating the sequential binding of Na+ and then K+ (or H+) ions to facilitate 5-HT uptake in SERT.

Ligand coupling mechanism of the human serotonin transporter differentiates substrates from inhibitors

The presynaptic serotonin transporter (SERT) clears extracellular serotonin following vesicular release to ensure temporal and spatial regulation of serotonergic signalling and neurotransmitter homeostasis. Prescription drugs used to treat neurobehavioral disorders, including depression, anxiety, and obsessive-compulsive disorder, trap SERT by blocking the transport cycle. In contrast, illicit drugs of abuse like amphetamines reverse SERT directionality, causing serotonin efflux. Both processes result in increased extracellular serotonin levels. By combining molecular dynamics simulations with biochemical experiments and using a homologous series of serotonin analogues, we uncovered the coupling mechanism between the substrate and the transporter, which triggers the uptake of serotonin. Free energy analysis showed that only scaffold-bound substrates could initiate SERT occlusion through attractive long-range electrostatic interactions acting on the bundle domain. The associated spatial requirements define substrate and inhibitor properties, enabling additional possibilities for rational drug design approaches.

Probing binding and occlusion of substrate in the human creatine transporter-1 by computation and mutagenesis

In chordates, energy buffering is achieved in part through phosphocreatine, which requires cellular uptake of creatine by the membrane-embedded creatine transporter (CRT1/SLC6A8). Mutations in human slc6a8 lead to creatine transporter deficiency syndrome, for which there is only limited treatment. Here, we used a combined homology modeling, molecular dynamics, and experimental approach to generate a structural model of CRT1. Our observations support the following conclusions: contrary to previous proposals, C144, a key residue in the substrate binding site, is not present in a charged state. Similarly, the side chain D458 must be present in a protonated form to maintain the structural integrity of CRT1. Finally, we identified that the interaction chain Y148-creatine-Na+ is essential to the process of occlusion, which occurs via a "hold-and-pull" mechanism. The model should be useful to study the impact of disease-associated point mutations on the folding of CRT1 and identify approaches which correct folding-deficient mutants.

Oligomerization of Monoamine Transporters

Transporters of the monoamine transporter (MAT) family regulate the uptake of important neurotransmitters like dopamine, serotonin, and norepinephrine. The MAT family functions using the electrochemical gradient of ions across the membrane and comprises three transporters, dopamine transporter (DAT), serotonin transporter (SERT), and norepinephrine transporter (NET). MAT transporters have been observed to exist in monomeric states to higher-order oligomeric states. Structural features, allosteric modulation, and lipid environment regulate the oligomerization of MAT transporters. NET and SERT oligomerization are regulated by levels of PIP2 present in the membrane. The kink present in TM12 in the MAT family is crucial for dimer interface formation. Allosteric modulation in the dimer interface hinders dimer formation. Oligomerization also influences the transporters' function, trafficking, and regulation. This chapter will focus on recent studies on monoamine transporters and discuss the factors affecting their oligomerization and its impact on their function.

Adaptive Sampling Methods for Molecular Dynamics in the Era of Machine Learning

Molecular dynamics (MD) simulations are fundamental computational tools for the study of proteins and their free energy landscapes. However, sampling protein conformational changes through MD simulations is challenging due to the relatively long time scales of these processes. Many enhanced sampling approaches have emerged to tackle this problem, including biased sampling and path-sampling methods. In this Perspective, we focus on adaptive sampling algorithms. These techniques differ from other approaches because the thermodynamic ensemble is preserved and the sampling is enhanced solely by restarting MD trajectories at particularly chosen seeds rather than introducing biasing forces. We begin our treatment with an overview of theoretically transparent methods, where we discuss principles and guidelines for adaptive sampling. Then, we present a brief summary of select methods that have been applied to realistic systems in the past. Finally, we discuss recent advances in adaptive sampling methodology powered by deep learning techniques, as well as their shortcomings.

Dynamic extracellular vestibule of human SERT: Unveiling druggable potential with high-affinity allosteric inhibitors

The serotonin transporter (SERT) tightly regulates synaptic serotonin levels and has been the primary target of antidepressants. Binding of inhibitors to the allosteric site of human SERT (hSERT) impedes the dissociation of antidepressants bound at the central site and may enhance the efficacy of such antidepressants to potentially reduce their dosage and side effects. Here, we report the identification of a series of high-affinity allosteric inhibitors of hSERT in a unique scaffold, with the lead compound, Lu AF88273 (3-(1-(2-(1H-indol-3-yl)ethyl)piperidin-4-yl)-6-chloro-1H-indole), having 2.1 nM allosteric potency in inhibiting imipramine dissociation. In addition, we find that Lu AF88273 also inhibits serotonin transport in a noncompetitive manner. The binding pose of Lu AF88273 in the allosteric site of hSERT is determined with extensive molecular dynamics simulations and rigorous absolute binding free energy perturbation (FEP) calculations, which show that a part of the compound occupies a dynamically formed small cavity. The predicted binding location and pose are validated by site-directed mutagenesis and can explain much of the structure-activity relationship of these inhibitors using the relative binding FEP calculations. Together, our findings provide a promising lead compound and the structural basis for the development of allosteric drugs targeting hSERT. Further, they demonstrate that the divergent allosteric sites of neurotransmitter transporters can be selectively targeted.

Activation mechanism of the human Smoothened receptor

Smoothened (SMO) is a membrane protein of the class F subfamily of G protein-coupled receptors (GPCRs) and maintains homeostasis of cellular differentiation. SMO undergoes conformational change during activation, transmitting the signal across the membrane, making it amenable to bind to its intracellular signaling partner. Receptor activation has been studied at length for class A receptors, but the mechanism of class F receptor activation remains unknown. Agonists and antagonists bound to SMO at sites in the transmembrane domain (TMD) and the cysteine-rich domain have been characterized, giving a static view of the various conformations SMO adopts. Although the structures of the inactive and active SMO outline the residue-level transitions, a kinetic view of the overall activation process remains unexplored for class F receptors. We describe SMO's activation process in atomistic detail by performing 300 μs of molecular dynamics simulations and combining it with Markov state model theory. A molecular switch, conserved across class F and analogous to the activation-mediating D-R-Y motif in class A receptors, is observed to break during activation. We also show that this transition occurs in a stage-wise movement of the transmembrane helices: TM6 first, followed by TM5. To see how modulators affect SMO activity, we simulated agonist and antagonist-bound SMO. We observed that agonist-bound SMO has an expanded hydrophobic tunnel in SMO's core TMD, whereas antagonist-bound SMO shrinks this tunnel, further supporting the hypothesis that cholesterol travels through a tunnel inside Smoothened to activate it. In summary, this study elucidates the distinct activation mechanism of class F GPCRs and shows that SMO's activation process rearranges the core TMD to open a hydrophobic conduit for cholesterol transport.

Reconciling membrane protein simulations with experimental DEER spectroscopy data

Spectroscopy experiments are crucial to study membrane proteins for which traditional structure determination methods still prove challenging. Double electron-electron resonance (DEER) spectroscopy experiments provide protein residue-pair distance distributions that are indicative of their conformational heterogeneity. Atomistic molecular dynamics (MD) simulations are another tool that have been proven to be vital to study the structural dynamics of membrane proteins such as to identify inward-open, occluded, and outward-open conformations of transporter membrane proteins, among other partially open or closed states of the protein. Yet, studies have reported that there is no direct consensus between the distributional data from DEER experiments and MD simulations, which has challenged validation of structures obtained from long-timescale simulations and using simulations to design experiments. Current coping strategies for comparisons rely on heuristics, such as mapping the nearest matching peaks between two ensembles or biased simulations. Here we examine the differences in residue-pair distance distributions arising due to the choice of membranes around the protein and covalent modification of a pair of residues to nitroxide spin labels in DEER experiments. Through comparing MD simulations of two proteins, PepT_So and LeuT-both of which have been characterized using DEER experiments previously-we show that the proteins' dynamics are similar despite the choice of the detergent micelle as a membrane mimetic in DEER experiments. On the other hand, covalently modified residues show slight local differences in their dynamics and a huge divergence when the oxygen atom pair distances between spin labeled residues are measured rather than protein backbone distances. Given the computational expense associated with pairwise MTSSL labeled MD simulations, we examine the use of biased simulations to explore the conformational dynamics of the spin labels only to reveal that such simulations alter the underlying protein dynamics. Our study identifies the main cause for the mismatch between DEER experiments and MD simulations and will accelerate the development of potential mitigation strategies to improve the match.

The Bright Side of Psychedelics: Latest Advances and Challenges in Neuropharmacology

The need to identify effective therapies for the treatment of psychiatric disorders is a particularly important issue in modern societies. In addition, difficulties in finding new drugs have led pharmacologists to review and re-evaluate some past molecules, including psychedelics. For several years there has been growing interest among psychotherapists in psilocybin or lysergic acid diethylamide for the treatment of obsessive-compulsive disorder, of depression, or of post-traumatic stress disorder, although results are not always clear and definitive. In fact, the mechanisms of action of psychedelics are not yet fully understood and some molecular aspects have yet to be well defined. Thus, this review aims to summarize the ethnobotanical uses of the best-known psychedelic plants and the pharmacological mechanisms of the main active ingredients they contain. Furthermore, an up-to-date overview of structural and computational studies performed to evaluate the affinity and binding modes to biologically relevant receptors of ibogaine, mescaline, N,N-dimethyltryptamine, psilocin, and lysergic acid diethylamide is presented. Finally, the most recent clinical studies evaluating the efficacy of psychedelic molecules in some psychiatric disorders are discussed and compared with drugs already used in therapy.

Multiagent Reinforcement Learning-Based Adaptive Sampling for Conformational Dynamics of Proteins

Machine learning is increasingly applied to improve the efficiency and accuracy of molecular dynamics (MD) simulations. Although the growth of distributed computer clusters has allowed researchers to obtain higher amounts of data, unbiased MD simulations have difficulty sampling rare states, even under massively parallel adaptive sampling schemes. To address this issue, several algorithms inspired by reinforcement learning (RL) have arisen to promote exploration of the slow collective variables (CVs) of complex systems. Nonetheless, most of these algorithms are not well-suited to leverage the information gained by simultaneously sampling a system from different initial states (e.g., a protein in different conformations associated with distinct functional states). To fill this gap, we propose two algorithms inspired by multiagent RL that extend the functionality of closely related techniques (REAP and TSLC) to situations where the sampling can be accelerated by learning from different regions of the energy landscape through coordinated agents. Essentially, the algorithms work by remembering which agent discovered each conformation and sharing this information with others at the action-space discretization step. A stakes function is introduced to modulate how different agents sense rewards from discovered states of the system. The consequences are three-fold: (i) agents learn to prioritize CVs using only relevant data, (ii) redundant exploration is reduced, and (iii) agents that obtain higher stakes are assigned more actions. We compare our algorithm with other adaptive sampling techniques (least counts, REAP, TSLC, and AdaptiveBandit) to show and rationalize the gain in performance.

Recent Advances in Machine Learning Variant Effect Prediction Tools for Protein Engineering

Proteins are Nature's molecular machinery and comprise diverse roles while consisting of chemically similar building blocks. In recent years, protein engineering and design have become important research areas, with many applications in the pharmaceutical, energy, and biocatalysis fields, among others-where the aim is to ultimately create a protein given desired structural and functional properties. It is often critical to model the relationship between a protein's sequence, folded structure, and biological function to assist in such protein engineering pursuits. However, significant challenges remain in concretely mapping an amino acid sequence to specific protein properties and biological activities. Mutations may enhance or diminish molecular protein function, and the epistatic interactions between mutations result in an inherently complex mapping between genetic modifications and protein function. Therefore, estimating the quantitative effects of mutations on protein function(s) remains a grand challenge of biology, bioinformatics, and many related fields and would rapidly accelerate protein engineering tasks when successful. Such estimation is often known as variant effect prediction (VEP). However, progress has been demonstrated in recent years with the development of machine learning (ML) methods in modeling the relationship between mutations and protein function. In this Review, recent advances in variant effect prediction (VEP) are discussed as tools for protein engineering, focusing on techniques incorporating gains from the broader ML community and challenges in estimating biomolecular functional differences. Primary developments highlighted include convolutional neural networks, graph neural networks, and natural language embeddings for protein sequences.

Structural Rearrangement of the Serotonin Transporter Intracellular Gate Induced by Thr276 Phosphorylation

The reuptake of the neurotransmitter serotonin from the synaptic cleft by the serotonin transporter, SERT, is essential for proper neurological signaling. Biochemical studies have shown that Thr276 of transmembrane helix 5 is a site of PKG-mediated SERT phosphorylation, which has been proposed to shift the SERT conformational equilibria to promote inward-facing states, thus enhancing 5-HT transport. Recent structural and simulation studies have provided insights into the conformation transitions during substrate transport but have not shed light on SERT regulation via post-translational modifications. Using molecular dynamics simulations and Markov state models, we investigate how Thr276 phosphorylation impacts the SERT mechanism and its role in enhancing transporter stability and function. Our simulations show that Thr276 phosphorylation alters the hydrogen-bonding network involving residues on transmembrane helix 5. This in turn decreases the free energy barriers for SERT to transition to the inward-facing state, thus facilitating 5-HT import. The results provide atomistic insights into in vivo SERT regulation and can be extended to other pharmacologically important transporters in the solute carrier family.