Microsecond-timescale Molecular Dynamics Research Articles

Coupling structural, chemical composition and stress fluctuations with relaxation dynamics in metallic glasses

Understanding the atomistic mechanisms of relaxation dynamics in metallic glasses remains a longstanding challenge. Here, using microsecond time scale molecular dynamics simulations, three main relaxation stages in metallic glasses are identified. At about 0.6Tg, chemical composition plays a dominant role in the relaxation, manifested by stress accumulation and only minimal variations in structure. As Tg approaches, the confluence of structural heterogeneity and chemical composition leads to the decoupling of relaxation mechanisms. At this temperature, the relaxation results in a structure with a lower energy state and a lower level of stress. In the supercooled liquid regime, an extensive increase in the number of closed-packed icosahedral clusters is responsible for accelerated structural relaxation while their packing frustration leads to the accumulation of intrinsic residual stresses in the glass. The atomistic origin of these dynamic relaxation modes is discussed in terms of structural, chemical composition and stress variations.

Long-timescale atomistic simulations uncover loss-of-function mechanisms of uncharacterized Angiogenin mutants associated with ALS

Amyotrophic Lateral Sclerosis (ALS) is a devastating neurodegenerative disease characterized by progressive degeneration of motor neurons, resulting in respiratory failure and mortality within 3–5 years. Mutations in the Angiogenin (ANG) cause loss of ribonucleolytic and nuclear translocation activities, contributing to ALS pathogenesis. This study focused on investigating two uncharacterized ANG mutations, T11S and R122H, newly identified in the Project Mine consortium. Using extensive computational analysis, including structural modeling and microsecond-timescale molecular dynamics (MD) simulations, we observed conformational changes in the catalytic residue His114 of ANG induced by T11S and R122H mutations. These alterations impaired ribonucleolytic activity, as inferred through molecular docking and binding free energy calculations. Gibbs free energy landscape and residue-residue interaction network analysis further supported our findings, revealing the energetic states and allosteric pathway from the mutated site to His114. Additionally, we assessed the binding of NCI-65828, an inhibitor of ribonucleolytic activity of ANG, and found reduced effectiveness in binding to T11S and R122H mutants when His114 assumed a non-native conformation. This highlights the crucial role of His114 and its association with ALS. Elucidating the relationship between physical structure and functional dynamics of frequently mutated ANG mutants is essential for understanding ALS pathogenesis and developing more effective therapeutic interventions.

Computation-guided engineering of distal mutations in an artificial enzyme.

Artificial enzymes are valuable biocatalysts able to perform new-to-nature transformations with the precision and (enantio-)selectivity of natural enzymes. Although they are highly engineered biocatalysts, they often cannot reach catalytic rates akin those of their natural counterparts, slowing down their application in real-world industrial processes. Typically, their designs only optimise the chemistry inside the active site, while overlooking the role of protein dynamics on catalysis. In this work, we show how the catalytic performance of an already engineered artificial enzyme can be further improved by distal mutations that affect the conformational equilibrium of the protein. To this end, we subjected a specialised artificial enzyme based on the lactococcal multidrug resistance regulator (LmrR) to an innovative algorithm that quickly inspects the whole protein sequence space for hotpots which affect the protein dynamics. From an initial predicted selection of 73 variants, two variants with mutations distant by more than 11 Å from the catalytic pAF residue showed increased catalytic activity towards the new-to-nature hydrazone formation reaction. Their recombination displayed a 66% higher turnover number and 14 °C higher thermostability. Microsecond time scale molecular dynamics simulations evidenced a shift in the distribution of productive enzyme conformations, which are the result of a cascade of interactions initiated by the introduced mutations.

Exploring the Dynamics of Holo-Shikimate Kinase through Molecular Mechanics

Understanding the connection between local and global dynamics can provide valuable insights into enzymatic function and may contribute to the development of novel strategies for enzyme modulation. In this work, we investigated the dynamics at both the global and local (active site) levels of Shikimate Kinase (SK) through microsecond time-scale molecular dynamics (MD) simulations of the holoenzyme in the product state. Our focus was on the wild-type (WT) enzyme and two mutants (R116A and R116K) which are known for their reduced catalytic activity. Through exploring the dynamics of these variants, we gained insights into the role of residue R116 and its contribution to overall SK dynamics. We argue that the connection between local and global dynamics can be attributed to local frustration near the mutated residue which perturbs the global protein dynamics.

Dispersion-corrected DFT calculations and umbrella sampling simulations to investigate stability of Chrysin-cyclodextrin inclusion complexes

The study of inclusion complexes of Chrysin (ChR) with three forms of cyclodextrins (CDs) α-, β-, and γ-CD was accomplished to examine the stability of ChR inside the central cavities of CDs. The aim of study was to identify the most suitable form of CD to improve the hydro-solubility of poorly soluble ChR bioactive molecule. Microsecond timescale molecular dynamics (MD) simulations were performed on four inclusion complexes (α-CD/ChR, β-CD/ChR, and two conformations of γ-CD/ChR) to examine the dynamics of ChR inside the cavity of CDs. The first conformation of γ-CD/ChR inclusion complex (γ-CD1/ChR) was identified to possess the highest affinity between host and guest molecule on the basis of binding energy calculated by employing Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) and umbrella sampling simulations. To further strengthen the claims of classical and biased MD studies, Our own N-layered Integrated molecular Orbital and Molecular mechanics (ONIOM) (wB97XD/6–311+g(d,p):pm7) calculations were performed on the selected inclusion complexes. The ONIOM based complexation energy reaffirmed that ChR had highest affinity for the γ-CD1 host molecule. Further, the non-covalent interaction analysis was conducted using Multiwfn software on QM-optimized inclusion complexes with wB97XD/6–311+G(d,p) model chemistry, revealing non-covalent interactions between ChR and CDs. This atomic level information helped us to gain better insights into critical atoms of ChR and CD that participated in intermolecular interactions and identify γ-CD as a suitable host molecule for improving the hydro-solubulity of ChR. The structural insights would help to derive new derivatives of γ-CD with better host capacity.

Parameterization of the miniPEG-Modified γPNA Backbone: Toward Induced γPNA Duplex Dissociation.

γ-Modified peptide nucleic acids (γPNAs) serve as potential therapeutic agents against genetic diseases. Miniature poly(ethylene glycol) (miniPEG) has been reported to increase solubility and binding affinity toward genetic targets, yet details of γPNA structure and dynamics are not understood. Within our work, we parameterized missing torsional and electrostatic terms for the miniPEG substituent on the γ-carbon atom of the γPNA backbone in the CHARMM force field. Microsecond timescale molecular dynamics simulations were carried out on six miniPEG-modified γPNA duplexes from NMR structures (PDB ID: 2KVJ). Three NMR models for the γPNA duplex (PDB ID: 2KVJ) were simulated as a reference for structural and dynamic changes captured for the miniPEG-modified γPNA duplex. Principal component analysis performed on the γPNA backbone atoms identified a single isotropic conformational substate (CS) for the NMR simulations, whereas four anisotropic CSs were identified for the ensemble of miniPEG-modified γPNA simulations. The NMR structures were found to have a 23° helical bend toward the major groove, consistent with our simulated CS structure of 19.0°. However, a significant difference between simulated methyl- and miniPEG-modified γPNAs involved the opportunistic invasion of miniPEG through the minor and major groves. Specifically, hydrogen bond fractional analysis showed that the invasion was particularly prone to affect the second G-C base pair, reducing the Watson-Crick base pair hydrogen bond by 60% over the six simulations, whereas the A-T base pairs decreased by only 20%. Ultimately, the invasion led to base stack reshuffling, where the well-ordered base stacking was reduced to segmented nucleobase stacking interactions. Our 6 μs timescale simulations indicate that duplex dissociation suggests the onset toward γPNA single strands, consistent with the experimental observation of decreased aggregation. To complement the insight of miniPEG-modified γPNA structure and dynamics, the new miniPEG force field parameters allow for further exploration of such modified γPNA single strands as potential therapeutic agents against genetic diseases.

In Silico Evaluation of the Thr58-Associated Conserved Water with KRAS Switch-II Pocket Binders.

The KRAS switch-II pocket (SII-P) has proven to be one of the most successful tools for targeting KRAS with small molecules to date. This has been demonstrated with several KRAS(G12C)-targeting covalent inhibitors, already resulting in two FDA-approved drugs. Several earlier-stage compounds have also been reported to engage KRAS SII-P with other position 12 mutants, including G12D, G12S, and G12R. A highly conserved water molecule exists in the KRAS SII-P, linking Thr58 of switch-II and Gly10 of β1 sheet. This conserved water is also present in the cocrystal structures of most of the disclosed small-molecule inhibitors but is only displaced by a handful of SII-P binders. Here, we evaluated the conserved water molecule energetics by the WaterMap for the SII-P binders with publicly disclosed structures and studied the water behavior in the presence of selected inhibitors by microsecond timescale molecular dynamics (MD) simulations using two water models (total simulation time of 120 μs). Our data revealed the high-energy nature of this hydration site when coexisting with an SII-P binder and that there is a preference for a single isolated hydration site in this location within the most advanced compounds. Furthermore, water displacement was only achieved with a few disclosed compounds and was suboptimal, as for instance a cyanomethyl group as a water displacer appears to introduce repulsion with the native conformation of Thr58. These results suggested that this conserved water should be considered more central when designing new inhibitors, especially in the design of noncovalent inhibitors targeting the SII-P.

Front Cover: When Two Become One: Conformational Changes in FXR/RXR Heterodimers Bound to Steroidal Antagonists (ChemMedChem 4/2023)

The Front Cover features the farnesoid X receptor (FXR) which is a nuclear receptor with a role in regulating bile acid synthesis and cholesterol homeostasis. This work studies the ligand-mediated FXR activation and heterodimerization using microsecond timescale molecular dynamics simulations. Recently reported FXR antagonists and well stablished agonists were investigated. Changes in the heterodimerization interface induced by antagonists prevent recruitment of co-regulatory elements. These results provide insights into their conformational changes and show what happens when two nuclear receptors become one system. More information can be found in the Research Article by Alejandro Díaz-Holguín, Azam Rashidian, Thales Kronenberger et al.

The atomistic face of the human MHC-I peptide-loading complex

Antigens exposed at the cell surface by major histocompatibility complex class I (MHC-I) proteins enable self/non-self recognition by cytotoxic T cells, protecting the organism against viral infections and cancer-causing mutations. To perform their role, MHC-I must first be loaded with an antigenic peptide in the peptide-loading complex (PLC), a large multi-protein assembly whose atomic-level structure and dynamics are still poorly understood. Using all-atom molecular dynamics (MD) simulations, we studied key elements of the human PLC, how they stabilise MHC-I and catalyse antigen selection, and how they assemble to form the PLC. By combining microsecond-timescale MD simulations with a recent 6-Å-resolution cryo-EM structure of the PLC, we obtained an atomistic model of the complete complex, in explicit solvent and in a membrane environment (1.5 million atoms). This model offers unprecedented insights into the structure and dynamics of the human antigen-loading machinery. Our simulations explain how tapasin, a central component of the PLC, acts as both an MHC-I chaperone and a catalyst that accelerates the off-rate of low-affinity peptides to facilitate antigen triage (peptide editing). We also show how tapasin recruits the transporter associated with antigen processing (TAP) into the PLC via transmembrane interactions. Finally, truncating antigens or removing them from the MHC-I binding groove gives a spatially resolved map of MHC-I plasticity which reveals how peptide loading status affects key structural regions. Taken together, our MD simulations explain experimental kinetics and mutagenesis data, and represent the first in-depth, atomic-level study of the mechanisms underlying the biological function of the PLC, an important step towards a better understanding of adaptive immunity.

Formation of lipid raft nanodomains in homogeneous ternary lipid mixture of POPC/DPSM/cholesterol: Theoretical insights

Lipid rafts, in biological membranes, are cholesterol-rich nanodomains that regulate many protein activities and cellular processes. Understanding the formation of the lipid-raft nanodomains helps us elucidate many complex interactions in the cell. In this study, the formation of lipid-raft nanodomains in a ternary palmitoyl-oleoyl-phosphatidylcholine/stearoyl-sphingomyelin/cholesterol (POPC/DPSM/Chol) lipid mixture, the most realistic surrogate model for biological membranes, has been successfully observed for the first time in-silico using microsecond timescale molecular dynamics simulations. The model reveals the formation of cholesterol-induced nanodomains with raft-like characteristics and their underlying mechanism: the cholesterol molecules segregate themselves into cholesterol nanodomains and then enrich the cholesterol-rich domain with sphingomyelin molecules to form a lipid raft thanks to the weak bonding of cholesterol with sphingomyelin. Besides, it is found that the increase in cholesterol concentration enhances the biophysical properties (e.g., bilayer thickness, area per lipid headgroup, and order parameter) of the lipid raft nanodomains. Such findings suggest that the POPC/DPSM/Chol bilayer is a suitable model to fundamentally extend the nanodomain evolution to investigate their lifetime and kinetics as well as to study protein-lipid interaction, protein-protein interaction, and selection of therapeutic molecules in the presence of lipid rafts.

A ricin-based peptide BRIP from Hordeum vulgare inhibits Mpro of SARS-CoV-2

COVID-19 pandemic caused by SARS-CoV-2 led to the research aiming to find the inhibitors of this virus. Towards this world problem, an attempt was made to identify SARS-CoV-2 main protease (Mpro) inhibitory peptides from ricin domains. The ricin-based peptide from barley (BRIP) was able to inhibit Mpro in vitro with an IC50 of 0.52 nM. Its low and no cytotoxicity upto 50 µM suggested its therapeutic potential against SARS-CoV-2. The most favorable binding site on Mpro was identified by molecular docking and steered molecular dynamics (MD) simulations. The Mpro-BRIP interactions were further investigated by evaluating the trajectories for microsecond timescale MD simulations. The structural parameters of Mpro-BRIP complex were stable, and the presence of oppositely charged surfaces on the binding interface of BRIP and Mpro complex further contributed to the overall stability of the protein-peptide complex. Among the components of thermodynamic binding free energy, Van der Waals and electrostatic contributions were most favorable for complex formation. Our findings provide novel insight into the area of inhibitor development against COVID-19.

Mechanistic behavior and subtle key events during DNA clamp opening and closing in T4 bacteriophage

Clamp loaders ensure processive DNA replication by loading the toroidal shaped sliding clamps onto the DNA. The sliding clamps serve as a platform for the attachment of polymerases and several other proteins associated with the regulation of various cellular processes. Clamp loaders are fascinating as nanomachines that engage in protein-protein and protein-DNA interactions. The loading mechanism of the clamp around dsDNA at the atomic level has not yet been fully explored. We performed microsecond timescale molecular dynamics simulations to reveal the dynamics of two different intermediate complexes involved in loading of the clamps around DNA. We conducted various time-dependent MD-driven analyses including the highly robust Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) calculations to observe changes in the structural elements of the clamp loader-clamp-DNA complexes in open and closed states. Our studies revealed the structural consequences of ATP hydrolysis events at different subunits of the clamp loader. This study would help in a better understanding of the clamp loading mechanism and would allow tackling various complications that might arise due to irregularities in this process.

Structural insights into the activation of blood coagulation factor XI zymogen by thrombin: A computational molecular dynamics study

Activation of human blood coagulation factor XI zymogen to factor XIa plays a significant role in the upstream coagulation pathway, in which factor XIa activates factor IX zymogen. The mechanistic details of the proteolytic activation of factor XI by the activating enzyme thrombin are not well-understood at atomic level. In this study, we employed a combination of molecular docking and microsecond time-scale molecular dynamics simulations to identify the key regions of interaction between fXI and thrombin. The activating complex between the substrate and enzyme was modeled to represent the initial acylation step of the serine-protease hydrolysis mechanism. The proposed solution structural complex, fIX:fIIa, obtained from 3 microseconds of MD refinement, suggests that the activation of factor XI is mediated by thrombin's anion binding exosite-II interactions with A3 and A4 domains. We predict that the two positively charged arginine residues (Arg409 and Arg413) in the exosite-2 region, the β- and γ-insertion loops of thrombin play an important structural role in the initial activating complex between fXI and thrombin.

Structural origin of thermal shrinkage in soda-lime silicate glass below the glass transition temperature: A theoretical investigation by microsecond timescale molecular dynamics simulations.

Microscopic dynamical features in the relaxation of glass structures are one of the most important unsolved problems in condensed matter physics. Although the structural relaxation processes in the vicinity of glass transition temperature are phenomenologically expressed by the Kohlrausch-Williams-Watts function and the relaxation time can be successfully interpreted by Adam-Gibbs theory and/or Narayanaswamy's model, the atomic rearrangement, which is the origin of the volume change, and its driving force have not been elucidated. Using the microsecond time-scale molecular dynamics simulations, this study provides insights to quantitatively determine the origin of the thermal shrinkage below Tg in a soda-lime silicate glass. We found that during annealing below Tg, Na ions penetrate into the six-membered silicate rings, which remedies the acute O-O-O angles of the energetically unstable rings. The ring structure change makes the space to possess the cation inside the rings, but the ring volume is eventually reduced, which results in thermal shrinkage of the soda-lime silica glass. In conclusion, the dynamical structural relaxation due to the cation displacement evokes the overall volume relaxation at low temperature in the glassy material.

Addressing a Trapped High-Energy Water: Design and Synthesis of Highly Potent Pyrimidoindole-Based Glycogen Synthase Kinase-3β Inhibitors.

In small molecule binding, water is not a passive bystander but rather takes an active role in the binding site, which may be decisive for the potency of the inhibitor. Here, by addressing a high-energy water, we improved the IC50 value of our co-crystallized glycogen synthase kinase-3β (GSK-3β) inhibitor by nearly two orders of magnitude. Surprisingly, our results demonstrate that this high-energy water was not displaced by our potent inhibitor (S)-3-(3-((7-ethynyl-9H-pyrimido[4,5-b]indol-4-yl)(methyl)amino)piperidin-1-yl)propanenitrile ((S)-15, IC50 value of 6 nM). Instead, only a subtle shift in the location of this water molecule resulted in a dramatic decrease in the energy of this high-energy hydration site, as shown by the WaterMap analysis combined with microsecond timescale molecular dynamics simulations. (S)-15 demonstrated both a favorable kinome selectivity profile and target engagement in a cellular environment and reduced GSK-3 autophosphorylation in neuronal SH-SY5Y cells. Overall, our findings highlight that even a slight adjustment in the location of a high-energy water can be decisive for ligand binding.

Modeling structural interconversion in Alzheimers’ amyloid beta peptide with classical and intrinsically disordered protein force fields

A comprehensive understanding of the aggregation mechanism in amyloid beta 42 (Aβ42) peptide is imperative for developing therapeutic drugs to prevent or treat Alzheimer’s disease. Because of the high flexibility and lack of native tertiary structures of Aβ42, molecular dynamics (MD) simulations may help elucidate the peptide’s dynamics with atomic details and collectively improve ensembles not seen in experiments. We applied microsecond-timescale MD simulations to investigate the dynamics and conformational changes of Aβ42 by using a newly developed Amber force field (ff14IDPSFF). We compared the ff14IDPSFF and the regular ff14SB force field by examining the conformational changes of two distinct Aβ42 monomers in explicit solvent. Conformational ensembles obtained by simulations depend on the force field and initial structure, Aβ42α-helix or Aβ42β-strand. The ff14IDPSFF sampled a high ratio of disordered structures and diverse β-strand secondary structures; in contrast, ff14SB favored helicity during the Aβ42α-helix simulations. The conformations obtained from Aβ42β-strand simulations maintained a balanced content in the disordered and helical structures when simulated by ff14SB, but the conformers clearly favored disordered and β-sheet structures simulated by ff14IDPSFF. The results obtained with ff14IDPSFF qualitatively reproduced the NMR chemical shifts well. In-depth peptide and cluster analysis revealed some characteristic features that may be linked to early onset of the fibril-like structure. The C-terminal region (mainly M35-V40) featured in-registered anti-parallel β-strand (β-hairpin) conformations with tested systems. Our work should expand the knowledge of force field and structure dependency in MD simulations and reveals the underlying structural mechanism–function relationship in Aβ42 peptides. Communicated by Ramaswamy H. Sarma

Thermodynamic Dissection of the Intercalation Binding Process of Doxorubicin to dsDNA with Implications of Ionic and Solvent Effects

Doxorubicin (DOX) is a cancer drug that binds to dsDNA through intercalation. A comprehensive microsecond timescale molecular dynamics study is performed for DOX with 16 tetradecamer dsDNA sequences in explicit aqueous solvent, in order to investigate the intercalation process at both binding stages (conformational change and insertion binding stages). The molecular mechanics generalized Born surface area (MM-GBSA) method is adapted to quantify and break down the binding free energy (BFE) into its thermodynamic components, for a variety of different solution conditions as well as different DNA sequences. Our results show that the van der Waals interaction provides the largest contribution to the BFE at each stage of binding. The sequence selectivity depends mainly on the base pairs located downstream from the DOX intercalation site, with a preference for (AT)2 or (TA)2 driven by the favorable electrostatic and/or van der Waals interactions. Invoking the quartet sequence model proved to be most successful to predict the sequence selectivity. Our findings also indicate that the aqueous bathing solution (i.e., water and ions) opposes the formation of the DOX-DNA complex at every binding stage, thus implying that the complexation process preferably occurs at low ionic strength and is crucially dependent on solvent effects.

E487K-Induced Disorder in Functionally Relevant Dynamics of Mitochondrial Aldehyde Dehydrogenase 2

Mitochondrial aldehyde dehydrogenase 2 (ALDH2), which is a homotetramer assembled by two equivalent dimers, is an important enzyme that metabolizes ethanol-derived acetaldehyde to acetate in a coenzyme-dependent manner. The highly reactive acetaldehyde exhibits a toxic effect, indicating that the proper functioning of ALDH2 is essential to counteract aldehyde-associated diseases. It is known that the catalytic activity of ALDH2 is drastically impaired by a frequently observed mutation, E487K, in a dominant fashion. However, the molecular basis of the inactivation mechanism is elusive because of the complex nature of the dynamic behavior. Here, we performed microsecond-timescale molecular dynamics simulations of the proteins complexed with coenzymes. The E487K mutation elevated the conformational heterogeneity of the dimer interfaces, which are relatively distal from the substituted residue. Dynamic network analyses showed that Glu487 and the dimer interface were dynamically communicated, and the dynamic community further spanned throughout all of the subunits in the wild-type; however, this network was completely rearranged by the E487K mutation. The perturbation of the dynamic properties led to alterations of the global conformational motions and destabilization of the coenzyme binding required for receiving a proton from the catalytic nucleophile. The insights into the dynamic behavior of the dominant negative mutant in this work will provide clues to restore its function.

Microsecond-timescale simulations suggest 5-HT–mediated preactivation of the 5-HT3A serotonin receptor

Aided by efforts to improve their speed and efficiency, molecular dynamics (MD) simulations provide an increasingly powerful tool to study the structure-function relationship of pentameric ligand-gated ion channels (pLGICs). However, accurate reporting of the channel state and observation of allosteric regulation by agonist binding with MD remains difficult due to the timescales necessary to equilibrate pLGICs from their artificial and crystalized conformation to a more native, membrane-bound conformation in silico. Here, we perform multiple all-atom MD simulations of the homomeric 5-hydroxytryptamine 3A (5-HT3A) serotonin receptor for 15 to 20 μs to demonstrate that such timescales are critical to observe the equilibration of a pLGIC from its crystalized conformation to a membrane-bound conformation. These timescales, which are an order of magnitude longer than any previous simulation of 5-HT3A, allow us to observe the dynamic binding and unbinding of 5-hydroxytryptamine (5-HT) (i.e., serotonin) to the binding pocket located on the extracellular domain (ECD) and allosteric regulation of the transmembrane domain (TMD) from synergistic 5-HT binding. While these timescales are not long enough to observe complete activation of 5-HT3A, the allosteric regulation of ion gating elements by 5-HT binding is indicative of a preactive state, which provides insight into molecular mechanisms that regulate channel activation from a resting state. This mechanistic insight, enabled by microsecond-timescale MD simulations, will allow a careful examination of the regulation of pLGICs at a molecular level, expanding our understanding of their function and elucidating key structural motifs that can be targeted for therapeutic regulation.

Molecular Mechanism of Acetate Transport through the Acetate Channel SatP.

Acetate is a central metabolite that plays a key role in almost all organisms, and acetate channels are often essential for their survival. Recently solved structures of the acetate channel Succinate-Acetate Permease (SatP) provide an atomic view of its closed state. However, the open state of the channel, the key residue conformational changes that trigger the channel to open, and the free energy barrier of acetate transportation remain elusive. To address these questions, we performed microsecond time scale molecular dynamics (MD) simulations and umbrella sampling. Several acetate passing events were observed in the MD trajectories with the application of an external electric field. Further analyses reveal the molecular mechanism of the channel opening, which results from the repacking of key residues, such as Gln50 and Phe17, as well as the subsequent outward movement of all transmembrane helices. Our simulations show that acetate is always surrounded by several water molecules when passing through the channel. Furthermore, a high energy barrier of 15 kcal/mol was observed from the free energy profile generated by umbrella sampling on the closed state of the channel. Our study deepens the understanding of the molecular mechanism of acetate transport through the channel SatP and is expected to facilitate the drug discovery on this target.