Coarse-grained simulations of proton-dependent conformational changes in lactose permease

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Substrate-dependent transport mechanism in AcrB of multidrug resistant bacteria

The multidrug resistance (MDR) system effectively expels antibiotics out of bacteria causing serious issues during bacterial infection. In addition to drug, indole, a common metabolic waste of bacteria, is expelled by MDR system of gram-negative bacteria for their survival. Experimental results suggest that AcrB, one of the key components of MDR system, undergoes large scale conformation changes during the pumping due to proton-motive process. However, due to extremely short time scale, it is difficult to observe (experimentally) those changes in the AcrB, which might facilitate the pumping process. Molecular simulations can shed light to understand the conformational changes for transport of indole in AcrB. Examination of conformational changes using all-atom simulation is, however, impractical. Here, we develop a hybrid coarse-grained force field to study the conformational changes of AcrB in presence of indole in the porter domain of monomer II. Using the coarse-grained force field, we investigated the conformational changes of AcrB for a number of model systems considering the effect of protonation in aspartic acid (Asp) residues Asp407 and Asp408 in the transmembrane domain of monomer II. Our results show that in the presence of indole, protonation of Asp408 or Asp407 residue causes conformational changes from binding state to extrusion state in monomer II, while remaining two monomers (I and III) approach access state in AcrB protein. We also observed that all three AcrB monomers prefer to go back to access state in the absence of indole. Steered molecular dynamics simulations were performed to demonstrate the feasibility of indole transport mechanism for protonated systems. Identification of indole transport pathway through AcrB can be very helpful in understanding the drug efflux mechanism used by the MDR bacteria.

Free Energy Landscape of the Complete Transport Cycle in a Key Bacterial Transporter

PepT_So is a proton-coupled bacterial symporter, from the major facilitator superfamily (MFS), which transports di-/tripeptide molecules. The recently obtained crystal structure of PepT_So provides an unprecedented opportunity to gain an understanding of functional insights of the substrate transport mechanism. Binding of the proton and peptide molecule induces conformational changes into occluded (OC) and outward-facing (OF) states, which we are able to characterize using molecular dynamics (MD) simulations. The structural knowledge of the OC and OF state is important to fully understand the major energy barrier associated with the transport cycle. In order to gain functional insight into the interstate dynamics, we performed extensive all atom MD simulations. The Markov state model was constructed to identify the free energy barriers between the states, and kinetic information on intermediate pathways was obtained using the transition pathway theory (TPT). TPT shows that the OF state is obtained by the movement of TM1 and TM7 at the extracellular side approximately 12-16 Å away from each other, and the inward movement of TM4 and TM10 at the intracellular halves to 3-4 Å characterizes the OC state. Helix distance distributions obtained from MD simulations were compared with experimental double electron-electron resonance spectroscopy and were found to be in excellent agreement with previous studies. We also predicted the optimal positions for placement of methane thiosulfonate spin label probes to capture the slowest protein dynamics. Our finding sheds light on the conformational cycle of this key membrane transporter and the functional relationships between the multiple intermediate states.

Uptake dynamics in the Lactose permease (LacY) membrane protein transporter

The sugar transporter Lactose permease (LacY) of Escherichia coli has become a prototype to understand the underlying molecular details of membrane transport. Crystal structures have trapped the protein in sugar-bound states facing the periplasm, but with narrow openings unable to accommodate sugar. Therefore, the molecular details of sugar uptake remain elusive. In this work, we have used extended simulations and metadynamics sampling to explore a putative sugar-uptake pathway and associated free energy landscape. We found an entrance at helix-pair 2 and 11, which involved lipid head groups and residues Gln 241 and Gln 359. Furthermore, the protein displayed high flexibility on the periplasmic side of Phe 27, which is located at the narrowest section of the pathway. Interactions to Phe 27 enabled passage into the binding site, which was associated with a 24 ± 4 kJ/mol binding free energy in excellent agreement with an independent binding free energy calculation and experimental data. Two free energy minima corresponding to the two possible binding poses of the lactose analog β-D-galactopyranosyl-1-thio-β-D-galactopyranoside (TDG) were aligned with the crystal structure-binding pocket. This work outlines the chemical environment of a putative periplasmic sugar pathway and paves way for understanding substrate affinity and specificity in LacY.

Coarse-grained simulations of conformational changes in the multidrug efflux transporter AcrB

The multidrug resistance (MDR) system actively pumps antibiotics out of cells causing serious health problems. During the pumping, AcrB (one of the key components of MDR) undergoes a series of large-scale and proton-motive conformational changes. Capturing the conformational changes through all-atom simulations is challenging. Here, we implement a hybrid coarse-grained force field to investigate the conformational changes of AcrB in the porter domain under different protonation states of Asp407/Asp408 in the trans-membrane domain. Our results show that protonation of Asp408 in monomer III (extrusion) stabilizes the asymmetric structure of AcrB; deprotonation of Asp408 induces clear opening of the entrance and closing of the exit leading to the transition from extrusion to access state. The structural changes in the porter domain of AcrB are strongly coupled with the proton translocation stoichiometry in the trans-membrane domain. Moreover, our simulations support the postulation that AcrB should adopt the symmetric resting state in a substrate-free situation.

Exploration of conformational changes in lactose permease upon sugar binding and proton transfer through coarse-grained simulations

Escherichia coli lactose permease (LacY) actively transports lactose and other galactosides across cell membranes through lactose/H⁺ symport process. Lactose/H⁺ symport is a highly complex process that involves sugar translocation, H⁺ transfer, and large-scale protein conformational changes. The complete picture of lactose/H⁺ symport is largely unclear due to the complexity and multiscale nature of the process. In this work, we develop the force field for sugar molecules compatible with PACE, a hybrid and coarse-grained force field that couples the united-atom protein models with the coarse-grained MARTINI water/lipid. After validation, we implement the new force field to investigate the binding of a β-d-galactopyranosyl-1-thio- β-d-galactopyranoside (TDG) molecule to a wild-type LacY. Results show that the local interactions between TDG and LacY at the binding pocket are consistent with the X-ray experiment. Transitions from inward-facing to outward-facing conformations upon TDG binding and protonation of Glu269 have been achieved from ∼5.5 µs simulations. Both the opening of the periplasmic side and closure of the cytoplasmic side of LacY are consistent with double electron-electron resonance and thiol cross-linking experiments. Our analysis suggests that the conformational changes of LacY are a cumulative consequence of interdomain H-bonds breaking at the periplasmic side, interdomain salt-bridge formation at the cytoplasmic side, and the TDG orientational changes during the transition.