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

An ATPase-Mimicking MXene nanozyme pharmacologically breaks the ironclad defense system for ferroptosis cancer therapy

Anticancer nanomedicines used for ferroptosis therapy generally rely on the direct delivery of Fenton catalysts to drive lipid peroxidation in cancer cells. However, the therapeutic efficacy is limited by the ferroptosis resistance caused by the intracellular anti-ferroptotic signals. Herein, we report the intrinsic ATPase-mimicking activity of a vanadium carbide MXene nanozyme (PVCMs) to pharmacologically modulate the nuclear factor erythroid 2-related factor 2 (Nrf2) program, which is the master anti-ferroptotic mediator in the ironclad defense system in triple-negative breast cancer (TNBC) cells. The PVCMs perform high ATPase-like activity that can effectively and selectively catalyze the dephosphorylation of ATP to generate ADP. Through a cascade mechanism initiated by falling energy status, PVCMs can powerfully hinder the Nrf2 program to selectively drive ferroptosis in TNBC cells in response to PVCMs-induced glutathione depletion. This study provides a paradigm for the use of pharmacologically active nanozymes to moderate specific cellular signals and elicit desirable pharmacological activities for therapeutic applications.

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.

Human ABCB1 with an ABCB11-like degenerate nucleotide binding site maintains transport activity by avoiding nucleotide occlusion

Several ABC exporters carry a degenerate nucleotide binding site (NBS) that is unable to hydrolyze ATP at a rate sufficient for sustaining transport activity. A hallmark of a degenerate NBS is the lack of the catalytic glutamate in the Walker B motif in the nucleotide binding domain (NBD). The multidrug resistance transporter ABCB1 (P-glycoprotein) has two canonical NBSs, and mutation of the catalytic glutamate E556 in NBS1 renders ABCB1 transport-incompetent. In contrast, the closely related bile salt export pump ABCB11 (BSEP), which shares 49% sequence identity with ABCB1, naturally contains a methionine in place of the catalytic glutamate. The NBD-NBD interfaces of ABCB1 and ABCB11 differ only in four residues, all within NBS1. Mutation of the catalytic glutamate in ABCB1 results in the occlusion of ATP in NBS1, leading to the arrest of the transport cycle. Here we show that despite the catalytic glutamate mutation (E556M), ABCB1 regains its ATP-dependent transport activity, when three additional diverging residues are also replaced. Molecular dynamics simulations revealed that the rescue of ATPase activity is due to the modified geometry of NBS1, resulting in a weaker interaction with ATP, which allows the quadruple mutant to evade the conformationally locked pre-hydrolytic state to proceed to ATP-driven transport. In summary, we show that ABCB1 can be transformed into an active transporter with only one functional catalytic site by preventing the formation of the ATP-locked pre-hydrolytic state in the non-canonical site.

ABC transporters are one of the largest membrane protein superfamilies, present in all organisms from archaea to humans. They transport a wide range of molecules including amino acids, sugars, vitamins, nucleotides, peptides, lipids, metabolites, antibiotics, and xenobiotics. ABC transporters energize substrate transport by hydrolyzing ATP in two symmetrically arranged nucleotide binding sites (NBSs). The human multidrug resistance transporter ABCB1 has two active NBSs, and it is generally believed that integrity and cooperation of both sites are needed for transport. Several human ABC transporters, such as the bile salt transporter ABCB11, have one degenerate NBS, which has significantly reduced ATPase activity. Interestingly, unilateral mutations affecting one of the two NBSs completely abolish the function of symmetrical ABC transporters. Here we engineered an ABCB1 variant with a degenerate, ABCB11-like NBS1, which can nevertheless transport substrates. Our results indicate that ABCB1 can mediate active transport with a single active site, questioning the validity of models assuming strictly alternating catalysis.

Enhanced sampling and free energy calculations for protein simulations

Molecular dynamics simulation is a powerful computational technique to study biomolecular systems, which complements experiments by providing insights into the structural dynamics relevant to biological functions at atomic scale. It can also be used to calculate the free energy landscapes of the conformational transitions to better understand the functions of the biomolecules. However, the sampling of biomolecular configurations is limited by the free energy barriers that need to be overcome, leading to considerable gaps between the timescales reached by MD simulation and those governing biological processes. To address this issue, many enhanced sampling methodologies have been developed to increase the sampling efficiency of molecular dynamics simulations and free energy calculations. Usually, enhanced sampling algorithms can be classified into methods based on collective variables (CV-based) and approaches which do not require predefined CVs (CV-free). In this chapter, the theoretical basis of free energy estimation is briefly reviewed first, followed by the reviews of the most common CV-based and CV-free methods including the presentation of some examples and recent developments. Finally, the combination of different enhanced sampling methods is discussed.

F508del disturbs the dynamics of the nucleotide binding domains of CFTR before and after ATP hydrolysis

The cystic fibrosis transmembrane conductance regulator (CFTR) channel is an ion channel responsible for chloride transport in epithelia and it belongs to the class of ABC transporters. The deletion of phenylalanine 508 (F508del) in CFTR is the most common mutation responsible for cystic fibrosis. Little is known about the effect of the mutation in the isolated nucleotide binding domains (NBDs), on dimer dynamics, ATP hydrolysis and even on nucleotide binding. Using molecular dynamics simulations of the human CFTR NBD dimer, we showed that F508del increases, in the prehydrolysis state, the inter-motif distance in both ATP binding sites (ABP) when ATP is bound. Additionally, a decrease in the number of catalytically competent conformations was observed in the presence of F508del. We used the subtraction technique to study the first 300 ps after ATP hydrolysis in the catalytic competent site and found that the F508del dimer evidences lower conformational changes than the wild type. Using longer simulation times, the magnitude of the conformational changes in both forms increases. Nonetheless, the F508del dimer shows lower C-α RMS values in comparison to the wild-type, on the F508del loop, on the residues surrounding the catalytic site and the portion of NBD2 adjacent to ABP1. These results provide evidence that F508del interferes with the NBD dynamics before and after ATP hydrolysis. These findings shed a new light on the effect of F508del on NBD dynamics and reveal a novel mechanism for the influence of F508del on CFTR.

A novel ring-shaped reaction pathway with interconvertible intermediates in chitinase A as revealed by QM/MM simulation combined with a one-dimensional projection technique

Efficient sampling achieved by the use of a one-dimensional projection technique reveals the catalytic mechanism of chitinase A from Serratia marcescens.

Molecular Mechanism of ATP Hydrolysis in an ABC Transporter

Hydrolysis of nucleoside triphosphate (NTP) plays a key role for the function of many biomolecular systems. However, the chemistry of the catalytic reaction in terms of an atomic-level understanding of the structural, dynamic, and free energy changes associated with it often remains unknown. Here, we report the molecular mechanism of adenosine triphosphate (ATP) hydrolysis in the ATP-binding cassette (ABC) transporter BtuCD-F. Free energy profiles obtained from hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations show that the hydrolysis reaction proceeds in a stepwise manner. First, nucleophilic attack of an activated lytic water molecule at the ATP γ-phosphate yields ADP + HPO₄ ^2- as intermediate product. A conserved glutamate that is located very close to the γ-phosphate transiently accepts a proton and thus acts as catalytic base. In the second step, the proton is transferred back from the catalytic base to the γ-phosphate, yielding ADP + H₂PO₄ ^-. These two chemical reaction steps are followed by rearrangements of the hydrogen bond network and the coordination of the Mg²⁺ ion. The rate constant estimated from the computed free energy barriers is in very good agreement with experiments. The overall free energy change of the reaction is close to zero, suggesting that phosphate bond cleavage itself does not provide a power stroke for conformational changes. Instead, ATP binding is essential for tight dimerization of the nucleotide-binding domains and the transition of the transmembrane domains from inward- to outward-facing, whereas ATP hydrolysis resets the conformational cycle. The mechanism is likely relevant for all ABC transporters and might have implications also for other NTPases, as many residues involved in nucleotide binding and hydrolysis are strictly conserved.

Mapping Free Energy Pathways for ATP Hydrolysis in the E. coli ABC Transporter HlyB by the String Method

HlyB functions as an adenosine triphosphate (ATP)-binding cassette (ABC) transporter that enables bacteria to secrete toxins at the expense of ATP hydrolysis. Our previous work, based on potential energy profiles from combined quantum mechanical and molecular mechanical (QM/MM) calculations, has suggested that the highly conserved H-loop His residue H662 in the nucleotide binding domain (NBD) of E. coli HlyB may catalyze the hydrolysis of ATP through proton relay. To further test this hypothesis when entropic contributions are taken into account, we obtained QM/MM minimum free energy paths (MFEPs) for the HlyB reaction, making use of the string method in collective variables. The free energy profiles along the MFEPs confirm the direct participation of H662 in catalysis. The MFEP simulations of HlyB also reveal an intimate coupling between the chemical steps and a local protein conformational change involving the signature-loop residue S607, which may serve a catalytic role similar to an Arg-finger motif in many ATPases and GTPases in stabilizing the phosphoryl-transfer transition state.

Allostery in the dengue virus NS3 helicase: Insights into the NTPase cycle from molecular simulations

The C-terminus domain of non-structural 3 (NS3) protein of the Flaviviridae viruses (e.g. HCV, dengue, West Nile, Zika) is a nucleotide triphosphatase (NTPase) -dependent superfamily 2 (SF2) helicase that unwinds double-stranded RNA while translocating along the nucleic polymer. Due to these functions, NS3 is an important target for antiviral development yet the biophysics of this enzyme are poorly understood. Microsecond-long molecular dynamic simulations of the dengue NS3 helicase domain are reported from which allosteric effects of RNA and NTPase substrates are observed. The presence of a bound single-stranded RNA catalytically enhances the phosphate hydrolysis reaction by affecting the dynamics and positioning of waters within the hydrolysis active site. Coupled with results from the simulations, electronic structure calculations of the reaction are used to quantify this enhancement to be a 150-fold increase, in qualitative agreement with the experimental enhancement factor of 10–100. Additionally, protein-RNA interactions exhibit NTPase substrate-induced allostery, where the presence of a nucleotide (e.g. ATP or ADP) structurally perturbs residues in direct contact with the phosphodiester backbone of the RNA. Residue-residue network analyses highlight pathways of short ranged interactions that connect the two active sites. These analyses identify motif V as a highly connected region of protein structure through which energy released from either active site is hypothesized to move, thereby inducing the observed allosteric effects. These results lay the foundation for the design of novel allosteric inhibitors of NS3.

Non-structural protein 3 (NS3) is a Flaviviridae (e.g. Hepatitis C, dengue, and Zika viruses) helicase that unwinds double stranded RNA while translocating along the nucleic polymer during viral genome replication. As a member of superfamily 2 (SF2) helicases, NS3 utilizes the free energy of nucleotide triphosphate (NTP) binding, hydrolysis, and product unbinding to perform its functions. While much is known about SF2 helicases, the pathways and mechanisms through which free energy is transduced between the NTP hydrolysis active site and RNA binding cleft remains elusive. Here we present a multiscale computational study to characterize the allosteric effects induced by the RNA and NTPase substrates (ATP, ADP, and P_i) as well as the pathways of short-range, residue-residue interactions that connect the two active sites. Results from this body of molecular dynamics simulations and electronic structure calculations are highlighted in context to the NTPase enzymatic cycle, allowing for development of testable hypotheses for validation of these simulations. Our insights, therefore, provide novel details about the biophysics of NS3 and guide the next generation of experimental studies.

Comparison of mechanistic transport cycle models of ABC exporters

ABC (ATP binding cassette) transporters, ubiquitous in all kingdoms of life, carry out essential substrate transport reactions across cell membranes. Their transmembrane domains bind and translocate substrates and are connected to a pair of nucleotide binding domains, which bind and hydrolyze ATP to energize import or export of substrates. Over four decades of investigations into ABC transporters have revealed numerous details from atomic-level structural insights to their functional and physiological roles. Despite all these advances, a comprehensive understanding of the mechanistic principles of ABC transporter function remains elusive. The human multidrug resistance transporter ABCB1, also referred to as P-glycoprotein (P-gp), is one of the most intensively studied ABC exporters. Using ABCB1 as the reference point, we aim to compare the dominating mechanistic models of substrate transport and ATP hydrolysis for ABC exporters and to highlight the experimental and computational evidence in their support. In particular, we point out in silico studies that enhance and complement available biochemical data. “This article is part of a Special Issue entitled: Beyond the Structure Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.”

Importance of an Orchestrate Participation of each Individual Residue Present at a Catalytic Site

GTP hydrolysis is indispensable to keep a living cell healthy. Nature has evolved so many enzymes to enhance the slow GTP hydrolysis. Rab GTPases are evolved to regulate vesicle trafficking. GTPase activating proteins (GAPs) accelerates their intrinsic slow GTP hydrolysis in order to maintain the sustainability between cellular events. Any malfunction/interference in this hydrolysis disrupts normal cellular events and causes severe diseases. In this study, GTP hydrolysis mechanism of Rab33B catalyzed by TBC-domain GAP protein Gyp1p has been decoded using extensive ab initio QM/MM metadynamics simulations. An organized coupled movement of individual residues present at the catalytic site is found to be the key factor for this reaction. An unorganized coupled movement leads the hydrolysis through very high energy pathways. This also reveals that the chemical transformations occurring at a catalytic site are residue specific.

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.