Molecular mechanism of proton-coupled ligand translocation by the bacterial efflux pump EmrE

Alchemical Transformations and Beyond: Recent Advances and Real-World Applications of Free Energy Calculations in Drug Discovery

Computational methods constitute efficient strategies for screening and optimizing potential drug molecules. A critical factor in this process is the binding affinity between candidate molecules and targets, quantified as binding free energy. Among various estimation methods, alchemical transformation methods stand out for their theoretical rigor. Despite challenges in force field accuracy and sampling efficiency, advancements in algorithms, software, and hardware have increased the application of free energy perturbation (FEP) calculations in the pharmaceutical industry. Here, we review the practical applications of FEP in drug discovery projects since 2018, covering both ligand-centric and residue-centric transformations. We show that relative binding free energy calculations have steadily achieved chemical accuracy in real-world applications. In addition, we discuss alternative physics-based simulation methods and the incorporation of deep learning into free energy calculations.

Molecular stripping underpins derepression of a toxin-antitoxin system

Transcription factors control gene expression; among these, transcriptional repressors must liberate the promoter for derepression to occur. Toxin-antitoxin (TA) modules are bacterial elements that autoregulate their transcription by binding the promoter in a T:A ratio-dependent manner, known as conditional cooperativity. The molecular basis of how excess toxin triggers derepression has remained elusive, largely because monitoring the rearrangement of promoter-repressor complexes, which underpin derepression, is challenging. Here, we dissect the autoregulation of the Salmonella enterica tacAT3 module. Using a combination of assays targeting DNA binding and promoter activity, as well as structural characterization, we determine the essential TA and DNA elements required to control transcription, and we reconstitute a repression-to-derepression path. We demonstrate that excess toxin triggers molecular stripping of the repressor complex off the DNA through multiple allosteric changes causing DNA distortion and ultimately leading to derepression. Thus, our work provides important insight into the mechanisms underlying conditional cooperativity.

All-atoms molecular dynamics study to screen potent efflux pump inhibitors against KpnE protein of Klebsiella pneumoniae

The Small Multidrug Resistance efflux pump protein KpnE, plays a pivotal role in multi-drug resistance in Klebsiella pneumoniae. Despite well-documented study of its close homolog, EmrE, from Escherichia coli, the mechanism of drug binding to KpnE remains obscure due to the absence of a high-resolution experimental structure. Herein, we exclusively elucidate its structure-function mechanism and report some of the potent inhibitors through drug repurposing. We used molecular dynamics simulation to develop a dimeric structure of KpnE and explore its dynamics in lipid-mimetic bilayers. Our study identified both semi-open and open conformations of KpnE, highlighting its importance in transport process. Electrostatic surface potential map suggests a considerable degree of similarity between KpnE and EmrE at the binding cleft, mostly occupied by negatively charged residues. We identify key amino acids Glu14, Trp63 and Tyr44, indispensable for ligand recognition. Molecular docking and binding free energy calculations recognizes potential inhibitors like acarbose, rutin and labetalol. Further validations are needed to confirm the therapeutic role of these compounds. Altogether, our membrane dynamics study uncovers the crucial charged patches, lipid-binding sites and flexible loop that could potentiate substrate recognition, transport mechanism and pave the way for development of novel inhibitors against K. pneumoniae.Communicated by Ramaswamy H. Sarma.

Comprehending conformational changes in EmrE, multidrug transporter at different pH: insights from molecular dynamics simulations

EmrE is a small multidrug resistance (SMR) pump of antiparallel topology that confers resistance to a broad range of polyaromatic cations in Escherichia coli. Atomic-level understanding of conformational changes for the selectivity of substrate and transport of a diverse array of drugs through the smallest known efflux pumps is crucial to multi-drug resistance. Therefore, the present study aims to provide insights into conformational changes during the transport through EmrE transporter at different pH. Molecular dynamics simulations have been carried out on the complete structure of EmrE in the absence of substrate. Computational analyses such as secondary structure, principal component, dynamic cross-correlation matrix, and hydrogen bond calculations have been performed. Analysis of MD trajectories in this study revealed pH-dependent interactions that influenced the structural dynamics of EmrE. Notably, at high pH, Glu14 and Tyr60 in both monomers formed electrostatic interactions, while these interactions decreased significantly at a low pH. Interestingly, a kink at helix 3 (H3) and dual open conformation of EmrE at low pH were also observed in contrast to a closed state discerned towards the periplasmic side at high pH. Significant interactions between C-terminal residues and residues at the edge of H1 & Loop1 and H3 & Loop3 were identified, suggesting their role in stabilizing the closed conformation of EmrE at the periplasmic end under high pH conditions. The present study enhances our understanding of EmrE's conformational changes, shedding light on the pH-dependent mechanisms that are likely to impact its function in multi-drug resistance.Communicated by Ramaswamy H. Sarma.

Investigation of molecular details of a bacterial cationic amino acid transporter (GkApcT) during arginine transportation using molecular dynamics simulation and umbrella sampling techniques

CONTEXT

Cationic amino acid transporters (CATs) facilitate arginine transport across membranes and maintain its levels in various tissues and organs, but their overexpression has been associated with severe cancers. A recent study identified the alternating access mechanism and critical residues involved in arginine transportation in a cationic amino acid transporter from Geobacillus kaustophilus (GkApcT). Here, we used molecular dynamics (MD) simulation methods to investigate the transportation mechanism of arginine (Arg) through GkApcT. The results revealed that arginine strongly interacts with specific binding site residues (Thr43, Asp111, Glu115, Lys191, Phe231, Ile234, and Asp237). Based on the umbrella sampling, the main driving force for arginine transport is the polar interactions of the arginine with channel-lining residues. An in-depth description of the dissociation mechanism and binding energy analysis brings valuable insight into the interactions between arginine and transporter residues, facilitating the design of effective CAT inhibitors in cancer cells.

METHODS

The membrane-protein system was constructed by uploading the prokaryotic CAT (PDB ID: 6F34) to the CHARMM-GUI web server. Molecular dynamics simulations were done using the GROMACS package, version 5.1.4, with the CHARMM36 force field and TIP3P water model. The MM-PBSA approach was performed for determining the arginine binding free energy. Furthermore, the hotspot residues were identified through per-residue decomposition analysis. The characteristics of the channel such as bottleneck radius and channel length were analyzed using the CaverWeb 1.1 web server. The proton wire inside the transporter was investigated based on the classic Grotthuss mechanism. We also investigated the atomistic details of arginine transportation using the path-based free energy umbrella sampling technique (US).

Still rocking in the structural era: A molecular overview of the small multidrug resistance (SMR) transporter family

The small multidrug resistance (SMR) family is composed of widespread microbial membrane proteins that fulfill different transport functions. Four functional SMR subtypes have been identified, which variously transport the small, charged metabolite guanidinium, bulky hydrophobic drugs and antiseptics, polyamines, and glycolipids across the membrane bilayer. The transporters possess a minimalist architecture, with ∼100-residue subunits that require assembly into homodimers or heterodimers for transport. In part because of their simple construction, the SMRs are a tractable system for biochemical and biophysical analysis. Studies of SMR transporters over the last 25 years have yielded deep insights for diverse fields, including membrane protein topology and evolution, mechanisms of membrane transport, and bacterial multidrug resistance. Here, we review recent advances in understanding the structures and functions of SMR transporters. New molecular structures of SMRs representing two of the four functional subtypes reveal the conserved structural features that have permitted the emergence of disparate substrate transport functions in the SMR family and illuminate structural similarities with a distantly related membrane transporter family, SLC35/DMT.

Crystal structures of bacterial small multidrug resistance transporter EmrE in complex with structurally diverse substrates

Proteins from the bacterial small multidrug resistance (SMR) family are proton-coupled exporters of diverse antiseptics and antimicrobials, including polyaromatic cations and quaternary ammonium compounds. The transport mechanism of the Escherichia coli transporter, EmrE, has been studied extensively, but a lack of high-resolution structural information has impeded a structural description of its molecular mechanism. Here, we apply a novel approach, multipurpose crystallization chaperones, to solve several structures of EmrE, including a 2.9 Å structure at low pH without substrate. We report five additional structures in complex with structurally diverse transported substrates, including quaternary phosphonium, quaternary ammonium, and planar polyaromatic compounds. These structures show that binding site tryptophan and glutamate residues adopt different rotamers to conform to disparate structures without requiring major rearrangements of the backbone structure. Structural and functional comparison to Gdx-Clo, an SMR protein that transports a much narrower spectrum of substrates, suggests that in EmrE, a relatively sparse hydrogen bond network among binding site residues permits increased sidechain flexibility.

A General Picture of Cucurbit[8]uril Host-Guest Binding

Describing, understanding, and designing complex interaction networks within macromolecular systems remain challenging in modern chemical research. Host-guest systems, despite their relative simplicity in both the structural feature and interaction patterns, still pose problems in theoretical modeling. The barrel-shaped supramolecular container cucurbit[8]uril (CB8) shows promising functionalities in various areas, e.g., catalysis and molecular recognition. It can stably coordinate a series of structurally diverse guests with high affinities. In this work, we examine the binding of seven commonly abused drugs to the CB8 host, aiming at providing a general picture of CB8-guest binding. Extensive sampling of the configurational space of these host-guest systems is performed, and the binding pathway and interaction patterns of CB8-guest complexes are investigated. A thorough comparison of widely used fixed-charge models for drug-like molecules is presented. Iterative refitting of the atomic charges suggests significant conformation dependence of charge generation. The initial model generated at the original conformation could be inaccurate for new conformations explored during conformational search, and the newly fitted charge set improves the prediction-experiment correlation significantly. Our investigations of the configurational space of CB8-drug complexes suggest that the host-guest interactions are more complex than expected. Despite the structural simplicities of these molecules, the conformational fluctuations of the host and the guest molecules and orientations of functional groups lead to the existence of an ensemble of binding modes. The insights of the binding thermodynamics, performance of fixed-charge models, and binding patterns of the CB8-guest systems are useful for studying and elucidating the binding mechanism of other host-guest complexes.