Analysis of the Free Energy Landscapes for the Opening-Closing Dynamics of the Maltose Transporter ATPase MalK2 Using Enhanced-Sampling Molecular Dynamics Simulation

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

The Transporter-Mediated Cellular Uptake and Efflux of Pharmaceutical Drugs and Biotechnology Products: How and Why Phospholipid Bilayer Transport Is Negligible in Real Biomembranes

Over the years, my colleagues and I have come to realise that the likelihood of pharmaceutical drugs being able to diffuse through whatever unhindered phospholipid bilayer may exist in intact biological membranes in vivo is vanishingly low. This is because (i) most real biomembranes are mostly protein, not lipid, (ii) unlike purely lipid bilayers that can form transient aqueous channels, the high concentrations of proteins serve to stop such activity, (iii) natural evolution long ago selected against transport methods that just let any undesirable products enter a cell, (iv) transporters have now been identified for all kinds of molecules (even water) that were once thought not to require them, (v) many experiments show a massive variation in the uptake of drugs between different cells, tissues, and organisms, that cannot be explained if lipid bilayer transport is significant or if efflux were the only differentiator, and (vi) many experiments that manipulate the expression level of individual transporters as an independent variable demonstrate their role in drug and nutrient uptake (including in cytotoxicity or adverse drug reactions). This makes such transporters valuable both as a means of targeting drugs (not least anti-infectives) to selected cells or tissues and also as drug targets. The same considerations apply to the exploitation of substrate uptake and product efflux transporters in biotechnology. We are also beginning to recognise that transporters are more promiscuous, and antiporter activity is much more widespread, than had been realised, and that such processes are adaptive (i.e., were selected by natural evolution). The purpose of the present review is to summarise the above, and to rehearse and update readers on recent developments. These developments lead us to retain and indeed to strengthen our contention that for transmembrane pharmaceutical drug transport "phospholipid bilayer transport is negligible".

ATP hydrolysis and nucleotide exit enhance maltose translocation in the MalFGK2E importer

ATP binding cassette (ABC) transporters employ ATP hydrolysis to harness substrate translocation across membranes. The Escherichia coli MalFGK₂E maltose importer is an example of a type I ABC importer and a model system for this class of ABC transporters. The MalFGK₂E importer is responsible for the intake of malto-oligossacharides in E.coli. Despite being extensively studied, little is known about the effect of ATP hydrolysis and nucleotide exit on substrate transport. In this work, we studied this phenomenon using extensive molecular dynamics simulations (MD) along with potential of mean force calculations of maltose transport across the pore, in the pre-hydrolysis, post-hydrolysis and nucleotide-free states. We concluded that ATP hydrolysis and nucleotide exit trigger conformational changes that result in the decrease of energetic barriers to maltose translocation towards the cytoplasm, with a concomitant increase of the energy barrier in the periplasmic side of the pore, contributing for the irreversibility of the process. We also identified key residues that aid in positioning and orientation of maltose, as well as a novel binding pocket for maltose in MalG. Additionally, ATP hydrolysis leads to conformations similar to the nucleotide-free state. This study shows the contribution of ATP hydrolysis and nucleotide exit in the transport cycle, shedding light on ABC type I importer mechanisms.

Structural dynamics of ABC transporters: molecular simulation studies

The biological activities of living organisms involve various inputs and outputs. The ATP-driven substances (biomolecules) responsible for these kinds of activities through membrane (i.e. uptake and efflux of substrates) include ATP-binding cassette (ABC) transporters, some of which play important roles in multidrug resistance. The basic architecture of ABC transporters comprises transmembrane domains (TMDs) and nucleotide-binding domains (NBDs). The functional dynamics (substrate transport) of ABC transporters are realized by concerted motions, such as NBD dimerization, mechanical transmission via coupling helices (CHs), and the translocation of substrates through TMDs, which are induced by the binding and/or hydrolysis of ATP molecules and substrates. In this mini-review, we briefly discuss recent progresses in the structural dynamics as revealed by molecular simulation studies at all-atom (AA), coarse-grained (CG), and quantum mechanics/molecular mechanics (QM/MM) levels.

Development of an ATP force field for coarse-grained simulation of ATPases and its application to the maltose transporter

The biological functions of ATPases, such as myosin, kinesin, and ABC transporter, are due to large conformational motions driven by energy obtained from ATP. Elucidation of the mechanisms underlying these ATP-driven movements is one of the greatest challenges in computational chemistry. It has been shown that the MARTINI coarse-grained method is a promising tool for the investigation of large conformational motions in various proteins. However, this method has not yet been applied to ATPases because of the lack of a force field for the ATP molecule. Here, we developed force field parameters for the ATP molecule and conducted simulations using these parameters for the subunits (MalK₂ ) and the full-length structure (MalFGK₂ -E) of a maltose transporter. It was found for both targets that the dimerization of the nucleotide binding domains (NBDs) is induced upon ATP binding. Moreover, for the full-length transporter, the conformational transition from the pre-translocation state to the outward-facing state was observed and was accompanied by an initial transport motion of the substrate. It is expected that coarse-grained simulations utilizing the parameters for the ATP molecule developed here will serve as a powerful tool for investigating other ATPases as well. © 2019 Wiley Periodicals, Inc.

The mechanism of nucleotide-binding domain dimerization in the intact maltose transporter as studied by all-atom molecular dynamics simulations

The Escherichia coli maltose transporter MalFGK₂ -E belongs to the protein superfamily of ATP-binding cassette (ABC) transporters. This protein is composed of heterodimeric transmembrane domains (TMDs) MalF and MalG, and the homodimeric nucleotide-binding domains (NBDs) MalK₂ . In addition to the TMDs and NBDs, the periplasmic maltose binding protein MalE captures maltose and shuttle it to the transporter. In this study, we performed all-atom molecular dynamics (MD) simulations on the maltose transporter and found that both the binding of MalE to the periplasmic side of the TMDs and binding of ATP to the MalK₂ are necessary to facilitate the conformational change from the inward-facing state to the occluded state, in which MalK₂ is completely dimerized. MalE binding suppressed the fluctuation of the TMDs and MalF periplasmic region (MalF-P2), and thus prevented the incorrect arrangement of the MalF C-terminal (TM8) helix. Without MalE binding, the MalF TM8 helix showed a tendency to intrude into the substrate translocation pathway, hindering the closure of the MalK₂ . This observation is consistent with previous mutagenesis experimental results on MalF and provides a new point of view regarding the understanding of the substrate translocation mechanism of the maltose transporter.

Analysis of an ATP-induced conformational transition of ABC transporter MsbA using a coarse-grained model

Upon the binding of ATP molecules to nucleotide binding domains (NBDs), ATP-binding cassette (ABC) exporters undergo a conformational transition from an inward-facing (IF) to an outward-facing (OF) state. This molecular event is a typical example of chemo-mechanical coupling. However, the underlying mechanism remains unclear. In this study, we analyzed the IF→OF transition of a representative ABC exporter, MsbA, by solving the equation of motion under an elastic network model (ENM). ATP was represented as a single node in ENM or replaced by external forces. When two ATP nodes were added to the ENM of the IF state protein, the two NBDs dimerized; subsequently, the two transmembrane domains opened toward the extracellular side, resulting in the formation of the OF structure. Such a conformational transition was also reproduced by applying external forces, which caused the rotational motion of the NBDs instead of the addition of ATP nodes. The process of the conformational transition was analyzed in detail using cross-correlation maps for node-node interactions. More importantly, it was revealed that the ATP binding energy is converted into distortion energy of several transmembrane helices. These results are useful for understanding the chemo-mechanical coupling in ABC transporters.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Negative Stain Single-particle EM of the Maltose Transporter in Nanodiscs Reveals Asymmetric Closure of MalK2 and Catalytic Roles of ATP, MalE, and Maltose

The Escherichia coli MalE-MalFGK₂ complex is one of the best characterized members of the large and ubiquitous family of ATP-binding cassette (ABC) transporters. It is composed of a membrane-spanning heterodimer, MalF-MalG; a homodimeric ATPase, MalK₂; and a periplasmic maltose receptor, MalE. Opening and closure of MalK₂ is coupled to conformational changes in MalF-MalG and the alternate exposition of the substrate-binding site to either side of the membrane. To further define this alternate access mechanism and the impact of ATP, MalE, and maltose on the conformation of the transporter during the transport cycle, we have reconstituted MalFGK₂ in nanodiscs and analyzed its conformations under 10 different biochemical conditions using negative stain single-particle EM. EM map results (at 15-25 Å resolution) indicate that binding of ATP to MalK₂ promotes an asymmetric, semi-closed conformation in accordance with the low ATPase activity of MalFGK₂ In the presence of MalE, the MalK dimer becomes fully closed, gaining the ability to hydrolyze ATP. In the presence of ADP or maltose, MalE·MalFGK₂ remains essentially in a semi-closed symmetric conformation, indicating that release of these ligands is required for the return to the initial state. Taken together, this structural information provides a rationale for the stimulation of MalK ATPase activity by MalE as well as by maltose.

Coupling between ATP hydrolysis and protein conformational change in maltose transporter

As the intracellular part of maltose transporter, MalK dimer utilizes the energy of ATP hydrolysis to drive protein conformational change, which then facilitates substrate transport. Free energy evaluation of the complete conformational change before and after ATP hydrolysis is helpful to elucidate the mechanism of chemical-to-mechanical energy conversion in MalK dimer, but is lacking in previous studies. In this work, we used molecular dynamics simulations to investigate the structural transition of MalK dimer among closed, semi-open and open states. We observed spontaneous structural transition from closed to open state in the ADP-bound system and partial closure of MalK dimer from the semi-open state in the ATP-bound system. Subsequently, we calculated the reaction pathways connecting the closed and open states for the ATP- and ADP-bound systems and evaluated the free energy profiles along the paths. Our results suggested that the closed state is stable in the presence of ATP but is markedly destabilized when ATP is hydrolyzed to ADP, which thus explains the coupling between ATP hydrolysis and protein conformational change of MalK dimer in thermodynamics. Proteins 2017; 85:207-220. © 2016 Wiley Periodicals, Inc.

ATP Hydrolysis Mechanism in a Maltose Transporter Explored by QM/MM Metadynamics Simulation

Translocation of substrates across the cell membrane by adenosine 5'-triphosphate (ATP)-binding cassette (ABC) transporters depends on the energy provided by ATP hydrolysis within the nucleotide-binding domains (NBDs). However, the detailed mechanism remains unclear. In this study, we focused on maltose transporter NBDs (MalK₂) and performed a quantum mechanical/molecular mechanical (QM/MM) well-tempered metadynamics simulation to address this issue. We explored the free-energy profile along an assigned collective variable. As a result, it was determined that the activation free energy is approximately 10.5 kcal/mol, and the reaction released approximately 3.8 kcal/mol of free energy, indicating that the reaction of interest is a one-step exothermic reaction. The dissociation of the ATP γ-phosphate seems to be the rate-limiting step, which supports the so-called dissociative model. Moreover, Glu159, located in the Walker B motif, acts as a base to abstract the proton from the lytic water, but is not the catalytic base, which corresponds to an atypical general base catalysis model. We also observed two interesting proton transfers: transfer from the His192 ε-position nitrogen to the dissociated inorganic phosphate, Pi, and transfer from the Lys42 side chain to adenosine 5'-diphosphate β-phosphate. These proton transfers would stabilize the posthydrolysis state. Our study provides significant insight into the ATP hydrolysis mechanism in MalK₂ from a dynamical viewpoint, and this insight would be applicable to other ABC transporters.