Continuum Electrostatic Calculations Research Articles

Delineating redox cooperativity in water-soluble and membrane multiheme cytochromes through protein design.

Nature has evolved diverse electron transport proteins and multiprotein assemblies essential to the generation and transduction of biological energy. However, substantially modifying or adapting these proteins for user-defined applications or to gain fundamental mechanistic insight can be hindered by their inherent complexity. De novo protein design offers an attractive route to stripping away this confounding complexity, enabling us to probe the fundamental workings of these bioenergetic proteins and systems, while providing robust, modular platforms for constructing completely artificial electron-conducting circuitry. Here, we use a set of de novo designed mono-heme and di-heme soluble and membrane proteins to delineate the contributions of electrostatic micro-environments and dielectric properties of the surrounding protein medium on the inter-heme redox cooperativity that we have previously reported. Experimentally, we find that the two heme sites in both the water-soluble and membrane constructs have broadly equivalent redox potentials in isolation, in agreement with Poisson-Boltzmann Continuum Electrostatics calculations. BioDC, a Python program for the estimation of electron transfer energetics and kinetics within multiheme cytochromes, also predicts equivalent heme sites, and reports that burial within the low dielectric environment of the membrane strengthens heme-heme electrostatic coupling. We conclude that redox cooperativity in our diheme cytochromes is largely driven by heme electrostatic coupling and confirm that this effect is greatly strengthened by burial in the membrane. These results demonstrate that while our de novo proteins present minimalist, new-to-nature constructs, they enable the dissection and microscopic examination of processes fundamental to the function of vital, yet complex, bioenergetic assemblies.

An expandable, modular de novo protein platform for precision redox engineering

The electron-conducting circuitry of life represents an as-yet untapped resource of exquisite, nanoscale biomolecular engineering. Here, we report the characterization and structure of a de novo diheme "maquette" protein, 4D2, which we subsequently use to create an expanded, modular platform for heme protein design. A well-folded monoheme variant was created by computational redesign, which was then utilized for the experimental validation of continuum electrostatic redox potential calculations. This demonstrates how fundamental biophysical properties can be predicted and fine-tuned. 4D2 was then extended into a tetraheme helical bundle, representing a 7 nm molecular wire. Despite a molecular weight of only 24 kDa, electron cryomicroscopy illustrated a remarkable level of detail, indicating the positioning of the secondary structure and the heme cofactors. This robust, expressible, highly thermostable and readily designable modular platform presents a valuable resource for redox protein design and the future construction of artificial electron-conducting circuitry.

Insights into receptor structure and dynamics at the surface of living cells

Evaluating protein structures in living cells remains a challenge. Here, we investigate Interleukin-4 receptor alpha (IL-4Rα) into which the non-canonical amino acid bicyclo[6.1.0]nonyne-lysine (BCNK) is incorporated by genetic code expansion. Bioorthogonal click labeling is performed with tetrazine-conjugated dyes. To quantify the reaction yield in situ, we develop brightness-calibrated ratiometric imaging, a protocol where fluorescent signals in confocal multi-color images are ascribed to local concentrations. Screening receptor mutants bearing BCNK in the extracellular domain uncovered site-specific variations of both click efficiency and Interleukin-4 binding affinity, indicating subtle well-defined structural perturbations. Molecular dynamics and continuum electrostatics calculations suggest solvent polarization to determine site-specific variations of BCNK reactivity. Strikingly, signatures of differential click efficiency, measured for IL-4Rα in ligand-bound and free form, mirror sub-angstrom deformations of the protein backbone at corresponding locations. Thus, click efficiency by itself represents a remarkably informative readout linked to protein structure and dynamics in the native plasma membrane.

Quantum Mechanical Study of Oxygen Ligands Protonation for the Stable States of the Laccase Active Site.

Laccases are enzymes catalyzing the oxidation of a wide range of organic and inorganic substrates accompanied by molecular oxygen reduction to water. Recently, oxygen reduction by laccases has been studied by single-crystal serial X-ray crystallography with increasing absorption doses at subatomic resolution. There were two determined structures corresponding to the reduced and oxidized stable states of the laccase active site. However, the protonation of the oxygen ligands involved cannot be determined even at subatomic resolution. In the present work, the protonation of oxygen ligands in the active site of laccase for the two stable states determined in the X-ray study was explored using quantum mechanical and continuum-electrostatics calculations. This is important for understanding the reaction of the oxygen reduction mechanism in laccases. The high precision of X-ray data at subatomic resolutions allowed us to optimize the quantum mechanical calculations.

Computational Assessment of Protein-Protein Binding Specificity within a Family of Synaptic Surface Receptors.

Atomic-level information is essential to explain the formation of specific protein complexes in terms of structure and dynamics. The set of Dpr and DIP proteins, which play a key role in the neuromorphogenesis in the nervous system of Drosophila melanogaster, offer a rich paradigm to learn about protein-protein recognition. Many members of the DIP subfamily cross-react with several members of the Dpr family and vice versa. While there exists a total of 231 possible Dpr-DIP heterodimer complexes from the 21 Dpr and 11 DIP proteins, only 57 "cognate" pairs have been detected by surface plasmon resonance (SPR) experiments, suggesting that the remaining 174 pairs have low or unreliable binding affinity. Our goal is to assess the performance of computational approaches to characterize the global set of interactions between Dpr and DIP proteins and identify the specificity of binding between each DIP with their corresponding Dpr binding partners. In addition, we aim to characterize how mutations influence the specificity of the binding interaction. In this work, a wide range of knowledge-based and physics-based approaches are utilized, including mutual information, linear discriminant analysis, homology modeling, molecular dynamics simulations, Poisson-Boltzmann continuum electrostatics calculations, and alchemical free energy perturbation to decipher the origin of binding specificity of the Dpr-DIP complexes examined. Ultimately, the results show that those two broad strategies are complementary, with different strengths and limitations. Biological inter-relations are more clearly revealed through knowledge-based approaches combining evolutionary and structural features, the molecular determinants controlling binding specificity can be predicted accurately with physics-based approaches based on atomic models.

Correcting electrostatic artifacts due to net-charge changes in the calculation of ligand binding free energies.

Alchemically derived free energies are artifacted when the perturbed moiety has a nonzero net charge. The source of the artifacts lies in the effective treatment of the electrostatic interactions within and between the perturbed atoms and remaining (partial) charges in the simulated system. To treat the electrostatic interactions effectively, lattice-summation (LS) methods or cutoff schemes in combination with a reaction-field contribution are usually employed. Both methods render the charging component of the calculated free energies sensitive to essential parameters of the system like the cutoff radius or the box side lengths. Here, we discuss the results of three previously published studies of ligand binding. These studies presented estimates of binding free energies that were artifacted due to the charged nature of the ligands. We show that the size of the artifacts can be efficiently calculated and raw simulation data can be corrected. We compare the corrected results with experimental estimates and nonartifacted estimates from path-sampling methods. Although the employed correction scheme involves computationally demanding continuum-electrostatics calculations, we show that the correction estimate can be deduced from a small sample of configurations rather than from the entire ensemble. This observation makes the calculations of correction terms feasible for complex biological systems. To show the general applicability of the proposed procedure, we also present results where the correction scheme was used to correct independent free energies obtained from simulations employing a cutoff scheme or LS electrostatics. In this work, we give practical guidelines on how to apply the appropriate corrections easily.

Computationally Modeling Electrostatic Binding Energetics in a Crowded, Dynamic Environment: Physical Insights from a Peptide-DNA System.

The cell is a crowded place, and it may be crucial at times to account for the local environment when studying determinants of molecular recognition. In this work, we use continuum electrostatics calculations on snapshots extracted from molecular dynamics simulations to understand how various aspects of a crowded environment affect electrostatic binding energies between the antimicrobial peptide buforin II and DNA. By comparing multiple models for representing crowding, sequentially introducing layers of model complexity for maximum control, we explore how electrostatic binding energetics depend on crowder physical properties, the sampling of the binding partners and crowder molecules, and the treatment of bulk solvent. We show that physical characteristics can combine to create an interplay of competing effects in this highly charged system. For example, increased ionic strength screening due to crowding partially cancels out the reduced solvent screening due to water depletion. We also quantify the effect of crowders' charge distributions on binding energetics. While we focus on electrostatic effects of crowding on binding, we begin to consider nonpolar components as well, and we implement a thermodynamic cycle accounting for both bound and unbound states to show the necessity of adequate crowder sampling in future studies. The insights developed here provide a rich starting point for experiments to further explore these competing effects and, ultimately, to rationally modulate molecular recognition in the complex cellular environment.

The H channel is not a proton transfer path in yeast cytochrome c oxidase

Cytochrome c oxidases (CcOs) in the respiratory chains of mitochondria and bacteria are primary consumers of molecular oxygen, converting it to water with the concomitant pumping of protons across the membrane to establish a proton electrochemical gradient. Despite a relatively well understood proton pumping mechanism of bacterial CcOs, the role of the H channel in mitochondrial forms of CcO remains debated. Here, we used site-directed mutagenesis to modify a central residue of the lower span of the H channel, Q413, in the genetically tractable yeast Saccharomyces cerevisiae. Exchange of Q413 to several different amino acids showed no effect on rates and efficiencies of respiratory cell growth, and redox potential measurements indicated minimal electrostatic interaction between the 413 locus and the nearest redox active component heme a. These findings clearly exclude a primary role of this section of the H channel in proton pumping in yeast CcO. In agreement with the experimental data, atomistic molecular dynamics simulations and continuum electrostatic calculations on wildtype and mutant yeast CcOs highlight potential bottlenecks in proton transfer through this route. Our data highlight the preference for neutral residues in the 413 locus, precluding sufficient hydration for formation of a proton conducting wire.

Redox-coupled quinone dynamics in the respiratory complex I

Complex I couples the free energy released from quinone (Q) reduction to pump protons across the biological membrane in the respiratory chains of mitochondria and many bacteria. The Q reduction site is separated by a large distance from the proton-pumping membrane domain. To address the molecular mechanism of this long-range proton-electron coupling, we perform here full atomistic molecular dynamics simulations, free energy calculations, and continuum electrostatics calculations on complex I from Thermus thermophilus We show that the dynamics of Q is redox-state-dependent, and that quinol, QH2, moves out of its reduction site and into a site in the Q tunnel that is occupied by a Q analog in a crystal structure of Yarrowia lipolytica We also identify a second Q-binding site near the opening of the Q tunnel in the membrane domain, where the Q headgroup forms strong interactions with a cluster of aromatic and charged residues, while the Q tail resides in the lipid membrane. We estimate the effective diffusion coefficient of Q in the tunnel, and in turn the characteristic time for Q to reach the active site and for QH2 to escape to the membrane. Our simulations show that Q moves along the Q tunnel in a redox-state-dependent manner, with distinct binding sites formed by conserved residue clusters. The motion of Q to these binding sites is proposed to be coupled to the proton-pumping machinery in complex I.

Protease‐protease interactions as a microenvironment‐dependent regulatory mechanism

Proteases are central to many physiological functions, including protein turnover, immunesystem response, and tissue remodeling. One family of proteases that have been of particular interest in recent years, are the cysteine cathepsins, which were originally characterized for their role in lysosomal protein turnover. Recent findings have shown that unchecked cathepsin activity contributes to the progression of diseases ranging from cancer to osteoporosis; however, all of the clinical trials for drugs targeting cathepsins, to date, have failed. Therapeutically targeting cathepsins is nontrivial due to the existence of a proteolytic network, defined by protease‐protease interactions that regulate concentrations of active enzymes through protease activation/degradation. Perturbation of the proteolytic network through the inhibition of cathepsins may lead to unexpected downstream off‐target effects. Further complexity arises when considering that the favorability of the interactions of the proteolytic network is affected by microenvironmental cues, such as pH. As a result, the direction and significance of connections within the proteolytic network are dependent on where the cathepsin is localized inside/outside the cell. Of particular interest to this work, are recently characterized cathepsin:cathepsin interactions that lead to cathepsin degradation, coined as “cannibalism”, which was first identified for the degradation of cathepsin K (catK) by cathepsin S (catS). We have utilized a hybrid computational and experimental approach to elucidate the pH‐dependence of catS:catK interactions and to quantify the effects of pH on catS on catK cannibalism. We performed computational protein docking, pKa predictions, and continuum electrostatics calculations to identify potential interactions that contribute to catS on catK cannibalism, as well as, to elucidate pH‐dependent effects. Our predictions indicate that interactions predicted for catS on catK cannibalism increase in favorability when increasing pH from 4 to 7.4. In vitro experiments using recombinant enzymes show that catK is degraded faster at pH 7.4 than at pH 4, suggesting that the strong predicted catS:catK interactions contribute to cannibalism. CatS, interestingly degrades catK slowest at pH 6, indicating that there is a preference for catK to be active at pH 6 despite cathepsins typically preferring more acidic environments, such as the lysosome. Our results illustrate a role for cathepsin cannibalism as a microenvironment‐dependent regulatory mechanism of cathepsin activity, which should be further investigated in the context of the proteolytic network and the development of cathepsin targeted therapeutics. Finally, the proposed hybrid approach lays the groundwork for investigating the role of microenvironmental pH in regulating other interactions in the proteolytic network inside and outside the cathepsin family of proteases.This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

In silico investigation of the molecular effects caused by R123H variant in secretory phospholipase A2-IIA associated with ARDS

Phospholipase A2-IIA catalyzes the hydrolysis of the sn-2 ester of glycerophospholipids. A rare c.428G > A (p.Arg143His) variant in PLA2G2A gene was found in two infants affected by acute respiratory distress syndrome (ARDS) by whole coding region and exon/intron boundaries sequencing. To obtain insights into the possible molecular effects of the rare R123H mutation in secretory PLA2-IIA (sPLA2-IIA), molecular modelling, molecular dynamics (MD) using principal component analysis (PCA) and continuum electrostatic calculations were conducted on the crystal structure of the wild type protein and on a generated model structure of the R123H mutant. Analysis of MD trajectories indicate that the overall stability of the protein is not affected by this mutation but nevertheless the catalytically crucial H-bond between Tyr51 and Asp91 as well as main electrostatic interactions in the region close to the mutation site are altered. PCA results indicate that the R123H replacement alter the internal molecular motions of the enzyme and that collective motions are increased. Electrostatic surface potential studies suggest that after mutation the interfacial binding to anionic phospholipid membranes and anionic proteins may be changed. The strengthening of electrostatic interactions may be propagated into the active site region thus potentially affecting the substrate recognition and enzymatic activity. Our findings provide the basis for further investigation and advances our understanding of the effects of mutations on sPLA2 structure and function.

Macromolecular Crowding Effects on Energetic Residue Contributions to Peptide-Nucleic Acid Interactions

Antimicrobial peptides (AMPs) are small, cationic peptides that are naturally found in a variety of organisms. AMPs have therapeutic potential as antibiotics due to their antimicrobial activity against a wide range of bacteria and other pathogens. Some AMPs function via translocation, which involves the undisruptive crossing of bacterial cell membranes followed by the targeting of an intracellular component. For example, buforin II (BF2) is believed to act by targeting intracellular nucleic acids. In an effort to better understand BF2 activity, previous studies have examined BF2-nucleic acid binding under dilute buffer conditions. These studies are expanded here to account for the macromolecular crowding environment present under physiological cellular conditions. In particular, we are interested in trying to determine the most effective computational models for considering these crowding effects. To this end, the binding of BF2 to DNA under crowded cell conditions is studied using a variety of computational crowding models. Crowding effects are accounted for through a representation of crowding as a lowered outer dielectric constant, as spherical crowders, or as explicit protein crowders. These models are applied to snapshots extracted from molecular dynamics (MD) simulations on which continuum electrostatics calculations are performed to examine any crowding effects on particular residue contributions to BF2-DNA binding. The use of MD enables better modeling of the dynamics and potential arrangements of crowders in solution. Preliminary analyses show that particular cationic residues exhibit an increased contribution to binding independent of the particular crowding model used. However, the magnitude of this increased contribution to binding appears to be dependent on the particular crowding model. These observations will help guide the models used by future work aimed at designing more potent peptides with enhanced nucleic acid binding.

Quantitative comparison of ABC membrane protein type I exporter structures in a standardized way

SummaryAn increasing number of ABC membrane protein structures are determined by cryo-electron microscopy and X-ray crystallography, consequently identifying differences between their conformations has become an arising issue. Therefore, we propose to define standardized measures for ABC Type I exporter structure characterization. We set conformational vectors, conftors, which describe the relative orientation of domains and can highlight structural differences. In addition, continuum electrostatics calculations were performed to characterize the energetics of membrane insertion illuminating functionally crucial regions. In summary, the proposed metrics contribute to deeper understanding of ABC membrane proteins' structural features, structure validation, and analysis of movements observed in a molecular dynamics trajectory. Moreover, the concept of standardized metrics can be applied not only to ABC membrane protein structures (http://conftors.hegelab.org).

Functional tuning of the catalytic residue pKa in a de novo designed esterase.

AlleyCatE is a de novo designed esterase that can be allosterically regulated by calcium ions. This artificial enzyme has been shown to hydrolyze p-nitrophenyl acetate (pNPA) and 4-nitrophenyl-(2-phenyl)-propanoate (pNPP) with high catalytic efficiency. AlleyCatE was created by introducing a single-histidine residue (His144 ) into a hydrophobic pocket of calmodulin. In this work, we explore the determinants of catalytic properties of AlleyCatE. We obtained the pKa value of the catalytic histidine using experimental measurements by NMR and pH rate profile and compared these values to those predicted from electrostatics pKa calculations (from both empirical and continuum electrostatics calculations). Surprisingly, the pKa value of the catalytic histidine inside the hydrophobic pocket of calmodulin is elevated as compared to the model compound pKa value of this residue in water. We determined that a short-range favorable interaction with Glu127 contributes to the elevated pKa of His144 . We have rationally modulated local electrostatic potential in AlleyCatE to decrease the pKa of its active nucleophile, His144 , by 0.7 units. As a direct result of the decrease in the His144 pKa value, catalytic efficiency of the enzyme increased by 45% at pH 6. This work shows that a series of simple NMR experiments that can be performed using low field spectrometers, combined with straightforward computational analysis, provide rapid and accurate guidance to rationally improve catalytic efficiency of histidine-promoted catalysis. Proteins 2017; 85:1656-1665. © 2017 Wiley Periodicals, Inc.

Insights into proton translocation in cbb3 oxidase from MD simulations

Heme-copper oxidases are membrane protein complexes that catalyse the final step of the aerobic respiration, namely the reduction of oxygen to water. The energy released during catalysis is coupled to the active translocation of protons across the membrane, which contributes to the establishment of an electrochemical gradient that is used for ATP synthesis. The distinctive C-type (or cbb3) cytochrome c oxidases, which are mostly present in proteobacteria, exhibit a number of unique structural and functional features, including high catalytic activity at low oxygen concentrations. At the moment, the functioning mechanism of C-type oxidases, in particular the proton transfer/pumping mechanism presumably via a single proton channel, is still poorly understood. In this work we used all-atom molecular dynamics simulations and continuum electrostatics calculations to obtain atomic-level insights into the hydration and dynamics of a cbb3 oxidase. We provide the details of the water dynamics and proton transfer pathways for both the “chemical” and “pumped” protons, and show that formation of protonic connections is strongly affected by the protonation state of key residues, namely H243, E323 and H337.

Assessment of Hydration Thermodynamics at Protein Interfaces with Grid Cell Theory.

Molecular dynamics simulations have been analyzed with the Grid Cell Theory (GCT) method to spatially resolve the binding enthalpies and entropies of water molecules at the interface of 17 structurally diverse proteins. Correlations between computed energetics and structural descriptors have been sought to facilitate the development of simple models of protein hydration. Little correlation was found between GCT-computed binding enthalpies and continuum electrostatics calculations. A simple count of contacts with functional groups in charged amino acids correlates well with enhanced water stabilization, but the stability of water near hydrophobic and polar residues depends markedly on its coordination environment. The positions of X-ray-resolved water molecules correlate with computed high-density hydration sites, but many unresolved waters are significantly stabilized at the protein surfaces. A defining characteristic of ligand-binding pockets compared to nonbinding pockets was a greater solvent-accessible volume, but average water thermodynamic properties were not distinctive from other interfacial regions. Interfacial water molecules are frequently stabilized by enthalpy and destabilized entropy with respect to bulk, but counter-examples occasionally occur. Overall detailed inspection of the local coordinating environment appears necessary to gauge the thermodynamic stability of water in protein structures.

Structural Determinants of Improved Fluorescence in a Family of Bacteriophytochrome-Based Infrared Fluorescent Proteins: Insights from Continuum Electrostatic Calculations and Molecular Dynamics Simulations.

Using X-ray crystallography, continuum electrostatic calculations, and molecular dynamics simulations, we have studied the structure, protonation behavior, and dynamics of the biliverdin chromophore and its molecular environment in a series of genetically engineered infrared fluorescent proteins (IFPs) based on the chromophore-binding domain of the Deinococcus radiodurans bacteriophytochrome. Our study suggests that the experimentally observed enhancement of fluorescent properties results from the improved rigidity and planarity of the biliverdin chromophore, in particular of the first two pyrrole rings neighboring the covalent linkage to the protein. We propose that the increases in the levels of both motion and bending of the chromophore out of planarity favor the decrease in fluorescence. The chromophore-binding pocket in some of the studied proteins, in particular the weakly fluorescent parent protein, is shown to be readily accessible to water molecules from the solvent. These waters entering the chromophore region form hydrogen bond networks that affect the otherwise planar conformation of the first three rings of the chromophore. On the basis of our simulations, the enhancement of fluorescence in IFPs can be achieved either by reducing the mobility of water molecules in the vicinity of the chromophore or by limiting the interactions of the nearby protein residues with the chromophore. Finally, simulations performed at both low and neutral pH values highlight differences in the dynamics of the chromophore and shed light on the mechanism of fluorescence loss at low pH.

Role of the Protonation State on the Structure and Dynamics of Albumin.

Human serum albumin undergoes reversible conformational transitions associated with ligand binding or pH changes. Among them, the neutral to base (N → B) transition occurring between pH 7 and pH 9 seems to be relevant for its function as a carrier. Unfortunately, a detailed atomic model for the B-form is still lacking, and several open questions remain concerning the charge distribution of the N-form. In this work, we report comparable molecular models for the N and B conformations that are built using continuum electrostatic calculations of pKa values and extended molecular dynamics (MD) simulations. Our computational models, which are critically assessed in terms of the available experimental observations relative to the N → B transition, reveal interesting similarities and differences between the N- and B-forms of HSA and highlight the importance of setting proper charge configurations in MD simulations.

The Hinge Region Strengthens the Nonspecific Interaction between Lac-Repressor and DNA: A Computer Simulation Study.

LacI is commonly used as a model to study the protein-DNA interaction and gene regulation. The headpiece of the lac-repressor (LacI) protein is an ideal system for investigation of nonspecific binding of the whole LacI protein to DNA. The hinge region of the headpiece has been known to play a key role in the specific binding of LacI to DNA, whereas its role in nonspecific binding process has not been elucidated. Here, we report the results of explicit solvent molecular dynamics simulation and continuum electrostatic calculations suggesting that the hinge region strengthens the nonspecific interaction, accounting for up to 50% of the micro-dissociation free energy of LacI from DNA. Consequently, the rate of microscopic dissociation of LacI from DNA is reduced by 2~3 orders of magnitude in the absence of the hinge region. We find the hinge region makes an important contribution to the electrostatic energy, the salt dependence of electrostatic energy, and the number of salt ions excluded from binding of the LacI-DNA complex.

Mechanisms for Two-Step Proton Transfer Reactions in the Outward-Facing Form of MATE Transporter

Bacterial pathogens or cancer cells can acquire multidrug resistance, which causes serious clinical problems. In cells with multidrug resistance, various drugs or antibiotics are extruded across the cell membrane by multidrug transporters. The multidrug and toxic compound extrusion (MATE) transporter is one of the five families of multidrug transporters. MATE from Pyrococcus furiosus uses H+ to transport a substrate from the cytoplasm to the outside of a cell. Crystal structures of MATE from P. furiosus provide essential information on the relevant H+-binding sites (D41 and D184). Hybrid quantum mechanical/molecular mechanical simulations and continuum electrostatic calculations on the crystal structures predict that D41 is protonated in one structure (Straight) and, both D41 and D184 protonated in another (Bent). All-atom molecular dynamics simulations suggest a dynamic equilibrium between the protonation states of the two aspartic acids and that the protonation state affects hydration in the substrate binding cavity and lipid intrusion in the cleft between the N- and C-lobes. This hypothesis is examined in more detail by quantum mechanical/molecular mechanical calculations on snapshots taken from the molecular dynamics trajectories. We find the possibility of two proton transfer (PT) reactions in Straight: the 1st PT takes place between side-chains D41 and D184 through a transient formation of low-barrier hydrogen bonds and the 2nd through another H+ from the headgroup of a lipid that intrudes into the cleft resulting in a doubly protonated (both D41 and D184) state. The 1st PT affects the local hydrogen bond network and hydration in the N-lobe cavity, which would impinge on the substrate-binding affinity. The 2nd PT would drive the conformational change from Straight to Bent. This model may be applicable to several prokaryotic H+-coupled MATE multidrug transporters with the relevant aspartic acids.