Evaluations of the Absolute and Relative Free Energies for Antidepressant Binding to the Amino Acid Membrane Transporter LeuT with Free Energy Simulations

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

The Pore-Lipid Interface: Role of Amino-Acid Determinants of Lipophilic Access by Ivabradine to the hERG1 Pore Domain

Abnormal cardiac electrical activity is a common side effect caused by unintended block of the promiscuous drug target human ether-à-go-go-related gene (hERG1), the pore-forming domain of the delayed rectifier K⁺ channel in the heart. hERG1 block leads to a prolongation of the QT interval, a phase of the cardiac cycle that underlies myocyte repolarization detectable on the electrocardiogram. Even newly released drugs such as heart-rate lowering agent ivabradine block the rapid delayed rectifier current I_Kr, prolong action potential duration, and induce potentially lethal arrhythmia known as torsades de pointes. In this study, we describe a critical drug-binding pocket located at the lateral pore surface facing the cellular membrane. Mutations of the conserved M651 residue alter ivabradine-induced block but not by the common hERG1 blocker dofetilide. As revealed by molecular dynamics simulations, binding of ivabradine to a lipophilic pore access site is coupled to a state-dependent reorientation of aromatic residues F557 and F656 in the S5 and S6 helices. We show that the M651 mutation impedes state-dependent dynamics of F557 and F656 aromatic cassettes at the protein-lipid interface, which has a potential to disrupt drug-induced block of the channel. This fundamentally new mechanism coupling the channel dynamics and small-molecule access from the membrane into the hERG1 intracavitary site provides a simple rationale for the well established state-dependence of drug blockade. SIGNIFICANCE STATEMENT: The drug interference with the function of the cardiac hERG channels represents one of the major sources of drug-induced heart disturbances. We found a novel and a critical drug-binding pocket adjacent to a lipid-facing surface of the hERG1 channel, which furthers our molecular understanding of drug-induced QT syndrome.

Characterization of the Cation Binding Sites in the NCKX2 Na+/Ca2+-K+ Exchanger

NCKX1-5 are proteins involved in K⁺-dependent Na⁺/Ca²⁺ exchange in various signal tissues. Here we present a homology model of NCKX2 based on the crystal structure of the NCX_Mj transporter found in Methanoccocus jannaschii. Molecular dynamics simulations were performed on the resultant wild-type NCKX2 model and two mutants (D548N and D575N) loaded with either four Na⁺ ions or one Ca²⁺ ion and one K⁺ ion, in line with the experimentally observed transport stoichiometry. The selectivity of the active site in wild-type NCKX2 for Na⁺, K⁺, and Li⁺ and the electrostatic interactions of the positive Na⁺ ions in the negatively charged active site of wild-type NCKX2 and the two mutants were evaluated from free energy perturbation calculations. For validation of the homology model, our computational results were compared to available experimental data obtained from numerous prior functional studies. The NCKX2 homology model is in good agreement with the discussed experimental data and provides valuable insights into the structure of the active site, which is lined with acidic and polar residues. The binding of the potassium and calcium ions is accomplished via Asp 575 and 548, respectively. Mutation of these residues to Asn alters the functionality of NCKX2 because of the elimination of the favorable carboxylate-cation interactions. The knowledge obtained from the NCKX2 model can be transferred to other isoforms of the NCKX family: newly discovered pathological mutations in NCKX4 and NCKX5 affect residues that are involved in ion binding and/or transport according to our homology model.

Effect of External Electric Field on Substrate Transport of a Secondary Active Transporter

Substrate transport across a membrane accomplished by a secondary active transporter (SAT) is essential to the normal physiological function of living cells. In the present research, a series of all-atom molecular dynamics (MD) simulations under different electric field (EF) strengths was performed to investigate the effect of an external EF on the substrate transport of an SAT. The results show that EF both affects the interaction between substrate and related protein's residues by changing their conformations and tunes the timeline of the transport event, which collectively reduces the height of energy barrier for substrate transport and results in the appearance of two intermediate conformations under the existence of an external EF. Our work spotlights the crucial influence of external EFs on the substrate transport of SATs and could provide a more penetrating understanding of the substrate transport mechanism of SATs.

Improved QM/MM Linear-Interaction Energy Model for Substrate Recognition in Zinc-Containing Metalloenzymes

One of the essential challenges in the description of receptor-drug interactions in the presence of various polyvalent cations (such as zinc, magnesium, or iron) is the accurate assessment of the electronic effects due to cofactor binding. The effects can range from partial electronic polarization of the proximal atoms in a receptor and bound substrate to long-range effects related to partial charge transfer and electronic delocalization effects between the cofactor and the drug. Here, we examine the role of the explicit account for electronic effects for a panel of small-molecule inhibitors binding to the zinc-aminopeptidase PfA-M1, an essential target for antimalarial drug development. Our study on PfA-M1:inhibitor interactions at the QM level reveals that the partial charge and proton transfer due to bound zinc ion are important mechanisms in the inhibitors' recognition and catalysis. The combination of classical MD simulations with a posteriori QM/MM corrections with novel DFTB parameters for the zinc cation and the linear-interaction energy (LIE) approach offers by far the most accurate estimates for the PfA-M1:inhibitor binding affinities, opening the door for future inhibitor design.

Accurate calculation of the absolute free energy of binding for drug molecules

Accurate prediction of binding affinities has been a central goal of computational chemistry for decades, yet remains elusive. Despite good progress, the required accuracy for use in a drug-discovery context has not been consistently achieved for drug-like molecules. Here, we perform absolute free energy calculations based on a thermodynamic cycle for a set of diverse inhibitors binding to bromodomain-containing protein 4 (BRD4) and demonstrate that a mean absolute error of 0.6 kcal mol-1 can be achieved. We also show a similar level of accuracy (1.0 kcal mol-1) can be achieved in pseudo prospective approach. Bromodomains are epigenetic mark readers that recognize acetylation motifs and regulate gene transcription, and are currently being investigated as therapeutic targets for cancer and inflammation. The unprecedented accuracy offers the exciting prospect that the binding free energy of drug-like compounds can be predicted for pharmacologically relevant targets.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

ATP transport through VDAC and the VDAC-tubulin complex probed by equilibrium and nonequilibrium MD simulations

Voltage-dependent anion channel (VDAC), the major channel of the mitochondrial outer membrane, serves as a principal pathway for ATP, ADP, and other respiratory substrates across this membrane. Using umbrella-sampling simulations, we established the thermodynamic and kinetic components governing ATP transport across the VDAC1 channel. We found that there are several low-affinity binding sites for ATP along the translocation pathway and that the main barrier for ATP transport is located around the center of the channel and is formed predominantly by residues in the N-terminus. The binding affinity of ATP to an open channel was found to be in the millimolar to micromolar range. However, we show that this weak binding increases the ATP translocation probability by about 10-fold compared with the VDAC pore in which attractive interactions were artificially removed. Recently, it was found that free dimeric tubulin induces a highly efficient, reversible blockage of VDAC reconstituted into planar lipid membranes. It was proposed that by blocking VDAC permeability for ATP/ADP and other mitochondrial respiratory substrates tubulin controls mitochondrial respiration. Using the Rosetta protein-protein docking algorithm, we established a tentative structure of the VDAC-tubulin complex. An extensive set of equilibrium and nonequilibrium (under applied electric field) molecular dynamics (MD) simulations was used to establish the conductance of the open and blocked channel. It was found that the presence of the unstructured C-terminal tail of tubulin in the VDAC pore decreases its conductance by more than 40% and switches its selectivity from anionic to cationic. The subsequent 1D potential of mean force (PMF) computations for the VDAC-tubulin complex show that the state renders ATP transport virtually impossible. A number of residues pivotal for tubulin binding to the channel were identified that help to clarify the molecular details of VDAC-tubulin interaction and to provide new insight into the mechanism of the control of mitochondria respiration by VDAC.

Insights from molecular dynamics: the binding site of cocaine in the dopamine transporter and permeation pathways of substrates in the leucine and dopamine transporters

The dopamine transporter (DAT) facilitates the regulation of synaptic neurotransmitter levels. As a target for therapeutic and illicit psycho-stimulant drugs like antidepressants and cocaine, DAT has been studied intensively. Despite a wealth of mutational and physiological data regarding DAT, the structure remains unsolved and details of the transport mechanism, binding sites and conformational changes remain debated. A bacterial homolog of DAT, the leucine transporter (LeuT(Aa)) has been used as a template and framework for modeling and understanding DAT. Free energy profiles obtained from Multi-Configuration Thermodynamic Integration simulations allowed us to correctly identify the primary and secondary binding pockets of LeuT(Aa). A comparison of free energy profiles for dopamine and cocaine in DAT suggests that the binding site of cocaine is located in a secondary pocket, not the primary substrate site. Two recurring primary pathways for intracellular substrate release from the primary pocket are identified in both transporters using the Random Acceleration Molecular Dynamics method. One pathway appears to follow transmembranes (TMs) 1a and 6b while the other pathway follows along TMs 6b and 8. Interestingly, we observe that a single sodium ion is co-transported with leucine during both simulation types.

Atomistic models of ion and solute transport by the sodium-dependent secondary active transporters

The recent determination of high-resolution crystal structures of several transporters offers unprecedented insights into the structural mechanisms behind secondary transport. These proteins utilize the facilitated diffusion of the ions down their electrochemical gradients to transport the substrate against its concentration gradient. The structural studies revealed striking similarities in the structural organization of ion and solute binding sites and a well-conserved inverted-repeat topology between proteins from several gene families. In this paper we will overview recent atomistic simulations applied to study the mechanisms of selective binding of ion and substrate in LeuT, Glt, vSGLT and hSERT as well as its consequences for the transporter conformational dynamics. This article is part of a Special Issue entitled: Membrane protein structure and function.

Alchemical free energy methods for drug discovery: progress and challenges

Improved rational drug design methods are needed to lower the cost and increase the success rate of drug discovery and development. Alchemical binding free energy calculations, one potential tool for rational design, have progressed rapidly over the past decade, but still fall short of providing robust tools for pharmaceutical engineering. Recent studies, especially on model receptor systems, have clarified many of the challenges that must be overcome for robust predictions of binding affinity to be useful in rational design. In this review, inspired by a recent joint academic/industry meeting organized by the authors, we discuss these challenges and suggest a number of promising approaches for overcoming them.

Evidence for a third sodium-binding site in glutamate transporters suggests an ion/substrate coupling model

Excitatory amino acid transporters (EAATs) remove glutamate from synapses. They maintain an efficient synaptic transmission and prevent glutamate from reaching neurotoxic levels. Glutamate transporters couple the uptake of one glutamate to the cotransport of three sodium ions and one proton and the countertransport of one potassium ion. The molecular mechanism for this coupled uptake of glutamate and its co- and counter-transported ions is not known. In a crystal structure of the bacterial glutamate transporter homolog, GltPh, only two cations are bound to the transporter, and there is no indication of the location of the third sodium site. In experiments using voltage clamp fluorometry and simulations based on molecular dynamics combined with grand canonical Monte Carlo and free energy simulations performed on different isoforms of GltPh as well on a homology model of EAAT3, we sought to locate the third sodium-binding site in EAAT3. Both experiments and computer simulations suggest that T370 and N451 (T314 and N401 in GltPh) form part of the third sodium-binding site. Interestingly, the sodium bound at T370 forms part of the binding site for the amino acid substrate, perhaps explaining both the strict coupling of sodium transport to uptake of glutamate and the ion selectivity of the affinity for the transported amino acid in EAATs.