Balanced activation of IspG and IspH to eliminate MEP intermediate accumulation and improve isoprenoids production in Escherichia coli

The MEP pathway genes were modulated to investigate whether there were new rate-limiting steps and toxic intermediates in this pathway. Activating IspG led to significant decrease of cell growth and β-carotene production. It was found that ispG overexpression led to accumulation of intermediate HMBPP, which seriously interfered with synthesis machinery of nucleotide and protein in Escherichia coli. Activation of the downstream enzyme IspH could solve HMBPP accumulation problem and eliminate the negative effects of ispG overexpression. In addition, intermediate MECPP accumulated in the starting strain, while balanced activation of IspG and IspH could push the carbon flux away from MECPP and led to 73% and 77% increase of β-carotene and lycopene titer respectively. Our work for the first time identified HMBPP to be a cytotoxic intermediate in MEP pathway and demonstrated that balanced activation of IspG and IspH could eliminate accumulation of HMBPP and MECPP and improve isoprenoids production.     

Structural Evidence for Enhancement of Sequential Vitamin D3 Hydroxylation Activities by Directed Evolution of Cytochrome P450 Vitamin D3 Hydroxylase

Vitamin D₃ hydroxylase (Vdh) isolated from actinomycete Pseudonocardia autotrophica is a cytochrome P450 (CYP) responsible for the biocatalytic conversion of vitamin D₃ (VD₃) to 1α,25-dihydroxyvitamin D₃ (1α,25(OH)₂VD₃) by P. autotrophica. Although its biological function is unclear, Vdh is capable of catalyzing the two-step hydroxylation of VD₃, i.e. the conversion of VD₃ to 25-hydroxyvitamin D₃ (25(OH)VD₃) and then of 25(OH)VD₃ to 1α,25(OH)₂VD₃, a hormonal form of VD₃. Here we describe the crystal structures of wild-type Vdh (Vdh-WT) in the substrate-free form and of the highly active quadruple mutant (Vdh-K1) generated by directed evolution in the substrate-free, VD₃-bound, and 25(OH)VD₃-bound forms. Vdh-WT exhibits an open conformation with the distal heme pocket exposed to the solvent both in the presence and absence of a substrate, whereas Vdh-K1 exhibits a closed conformation in both the substrate-free and substrate-bound forms. The results suggest that the conformational equilibrium was largely shifted toward the closed conformation by four amino acid substitutions scattered throughout the molecule. The substrate-bound structure of Vdh-K1 accommodates both VD₃ and 25(OH)VD₃ but in an anti-parallel orientation. The occurrence of the two secosteroid binding modes accounts for the regioselective sequential VD₃ hydroxylation activities. Moreover, these structures determined before and after directed evolution, together with biochemical and spectroscopic data, provide insights into how directed evolution has worked for significant enhancement of both the VD₃ 25-hydroxylase and 25(OH)VD₃ 1α-hydroxylase activities.     

Nucleotide pocket thermodynamics measured by EPR reveal how energy partitioning relates myosin speed to efficiency

We have used spin-labeled ADP to investigate the dynamics of the nucleotide-binding pocket in a series of myosins, which have a range of velocities. Electron paramagnetic resonance spectroscopy reveals that the pocket is in equilibrium between open and closed conformations. In the absence of actin, the closed conformation is favored. When myosin binds actin, the open conformation becomes more favored, facilitating nucleotide release. We found that faster myosins favor a more closed pocket in the actomyosin•ADP state, with smaller values of ΔH⁰ and ΔS⁰, even though these myosins release ADP at a faster rate. A model involving a partitioning of free energy between work-generating steps prior to rate-limiting ADP release explains both the unexpected correlation between velocity and opening of the pocket and the observation that fast myosins are less efficient than slow myosins.     

Structural Assignment of Spectra by Characterization of Conformational Substates in Bound MbCO

Residue motions of the distal heme pocket and bound CO ligand of carbonmonoxy Myoglobin are studied using a combination of molecular dynamics simulations and quantum chemical methods. Using mixed quantum mechanics/molecular mechanics calculations together with sampling from molecular dynamics simulations (QM/MM(MD)), the experimentally observed spectroscopic A₀ and A₁ substates of the bound CO ligand are assigned to the open and closed conformation of His⁶⁴ and the His_ɛ⁶⁴ tautomer, respectively. Several previously proposed origins of the A₃ substate, including rotamers of the doubly protonated His⁶⁴H⁺ side chain, His⁶⁴H⁺ inside the distal pocket, and cooperative motions with Arg⁴⁵, are investigated with QM/MM(MD). However, the signatures of the calculated infrared spectra do not agree with the experimentally observed ones. For additional insight on this, extensive molecular dynamics simulations are used together with improved electrostatics for the bound ligand. A CO fluctuating charge model is developed to describe the ab initio dipole and quadrupole moments of the bound ligand. CO absorption spectra are then obtained directly from the dynamics simulations. Finally, the electrostatics of the heme pocket is examined in detail in an attempt to determine the structural origins of the observed spectroscopic A-states from MD simulations. However, contrary to related simulations for unbound CO in myoglobin, the shifts and splittings for carbonmonoxy Myoglobin are generally small and difficult to relate to structural change. This suggests that coupling of the CO motion to other degrees of freedom, such as the Fe-CO stretching and bending, is important to correctly describe the dynamics of bound CO in myoglobin.     

Porter domain opening and closing motions in the multi-drug efflux transporter AcrB

Acriflavine resistance protein B acts as the active transporter in the multi-drug efflux pump Acriflavine resistance proteins A / B - Tolerance to colicins protein in Escherichia coli. Within the same reaction cycle intermediate all Acriflavine resistance protein B X-ray structures display highly similar conformations of the substrate-recruiting and transporting porter domain. To assess if this structural homogeneity is an intrinsic feature of Acriflavine resistance protein B or stems from other causes we performed a series of six independent, unbiased 100 ns molecular dynamics simulations of membrane-embedded, asymmetric, substrate-free wild type Acriflavine resistance protein B in a 150 mM NaCl solution. We find the porter domain more flexible than previously assumed displaying clear opening and closing motions of the proximal binding pocket (L and T-state) and the exit of the drug transport channels (O-intermediate). Concurrently the hydrophobic binding pocket favors a closed conformation in all three protomers. Our findings suggest that the conformational homogeneity seen in the crystal structures is likely an effect of bound but structurally unresolved substrate. Our simulations further imply that each of the known three reaction cycle intermediates occurs in at least two variants, the Thr676 loop independently regulates porter domain access and likely plays a key role in substrate transport. On a 100 ns time scale we find no evidence supporting the proposed LLL resting state in the absence of substrate. If the proximal binding pocket dynamics have an inhibiting effect on Acriflavine resistance protein B pump activity lowering the life time of substrate-accessible conformations, the observed dynamics could provide a structural explanation for the Acriflavine resistance protein B activity-enhancing effect of the adaptor protein Acriflavine resistance protein A stabilizing PC1 and PC2 subdomain orientations.     

A closer look at the spectroscopic properties of possible reaction intermediates in wild-type and mutant (E)-4-hydroxy-3-methylbut-2-enyl diphosphate reductase

(E)-4-Hydroxy-3-methylbut-2-enyl diphosphate reductase (IspH or LytB) catalyzes the terminal step of the MEP/DOXP pathway where it converts (E)-4-hydroxy-3-methylbut-2-enyl diphosphate (HMBPP) into the two products, isopentenyl diphosphate and dimethylallyl diphosphate. The reaction involves the reductive elimination of the C4 hydroxyl group, using a total of two electrons. Here we show that the active form of IspH contains a [4Fe-4S] cluster and not the [3Fe-4S] form. Our studies show that the cluster is the direct electron source for the reaction and that a reaction intermediate is bound directly to the cluster. This active form has been trapped in a state, dubbed FeS(A), that was detected by electron paramagnetic resonance (EPR) spectroscopy when one-electron-reduced IspH was incubated with HMBPP. In addition, three mutants of IspH have been prepared and studied, His42, His124, and Glu126 (Aquifex aeolicus numbering), with particular attention paid to the effects on the cluster properties and possible reaction intermediates. None of the mutants significantly affected the properties of the [4Fe-4S](+) cluster, but different effects were observed when one-electron-reduced forms were incubated with HMBPP. Replacing His42 led to an increased K(M) value and a much lower catalytic efficiency, confirming the role of this residue in substrate binding. Replacing the His124 also resulted in a lower catalytic efficiency. In this case, however, the enzyme showed the loss of the [4Fe-4S](+) EPR signal upon addition of HMBPP without the subsequent formation of the FeS(A) signal. Instead, a radical-type signal was observed in some of the samples, indicating that this residue plays a role in the correct positioning of the substrate. The incorrect orientation in the mutant leads to the formation of substrate-based radicals instead of the cluster-bound intermediate complex FeS(A). Replacing the Glu126 also resulted in a lower catalytic efficiency, with yet a third type of EPR signal being detected upon incubation with HMBPP. (31)P and (2)H ENDOR measurements of the FeS(A) species incubated with regular and (2)H-C4-labeled HMBPP reveal that the substrate binds to the enzyme in the proximity of the active-site cluster with C4 adjacent to the site of linkage between the FeS cluster and HMBPP. Comparison of the spectroscopic properties of this intermediate to those of intermediates detected in (E)-4-hydroxy-3-methylbut-2-enyl diphosphate synthase and ferredoxin:thioredoxin reductase suggests that HMBPP binds to the FeS cluster via its hydroxyl group instead of a side-on binding as previously proposed for the species detected in the inactive Glu126 variant. Consequences for the IspH reaction mechanism are discussed.     

Reconstruction and Evaluation of the Synthetic Bacterial MEP Pathway in Saccharomyces cerevisiae

Isoprenoids, which are a large group of natural and chemical compounds with a variety of applications as e.g. fragrances, pharmaceuticals and potential biofuels, are produced via two different metabolic pathways, the mevalonate (MVA) pathway and the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway. Here, we attempted to replace the endogenous MVA pathway in Saccharomyces cerevisiae by a synthetic bacterial MEP pathway integrated into the genome to benefit from its superior properties in terms of energy consumption and productivity at defined growth conditions. It was shown that the growth of a MVA pathway deficient S. cerevisiae strain could not be restored by the heterologous MEP pathway even when accompanied by the co-expression of genes erpA, hISCA1 and CpIscA involved in the Fe-S trafficking routes leading to maturation of IspG and IspH and E. coli genes fldA and fpr encoding flavodoxin and flavodoxin reductase believed to be responsible for electron transfer to IspG and IspH.     

Comparative structural dynamics of Tyrosyl-tRNA synthetase complexed with different substrates explored by molecular dynamics

Several molecular dynamics simulations of S. aureus Tyrosyl-tRNA synthetase (TyrRS) in its free form and complexed with Tyr, ATP, tyrosyl adenylate and inhibitor respectively have been carried out to investigate the ligand-linked conformational stability changes associated with its catalytic cycle. The results show that unliganded S. aureus TyrRS samples a more relaxed conformational space than substrate-bound TyrRS. There are three high flexibility regions encompassing residues 114–118, 128–133, and 226–238 respectively. The region which includes the KMSKS motif (KFGKS in S. aureus TyrRS) shows the highest difference in fluctuations. Hydrogen bond network formed by Tyr, ATP, tyrosyl adenylate and inhibitor with S. aureus TyrRS is discussed. Our simulations suggest the induced-fit conformational changes of the KMSKS loop as follows: the KMSKS loop of substrate-free S. aureus TyrRS adopts an open conformation. The tyrosine binds in the pocket with the KMSKS loop balancing between semi-open and open forms. The ATP binding induces the KMSKS loop to the open form. After the Tyr-AMP is formed, the first two residues of KMSKS loop twists in an anticlockwise direction and drives the loop in a conformation similar to the closed one, while those of the last three GKS residues adopt a conformation between semi-open and open conformation. This conformational change may probably be involved in the initial tRNA binding. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structural and Functional Characterization of the Kindlin-1 Pleckstrin Homology Domain

Inside-out activation of integrins is mediated via the binding of talin and kindlin to integrin β-subunit cytoplasmic tails. The kindlin FERM domain is interrupted by a pleckstrin homology (PH) domain within its F2 subdomain. Here, we present data confirming the importance of the kindlin-1 PH domain for integrin activation and its x-ray crystal structure at a resolution of 2.1 Å revealing a C-terminal second α-helix integral to the domain but found only in the kindlin protein family. An isoform-specific salt bridge occludes the canonical phosphoinositide binding site, but molecular dynamics simulations display transient switching to an alternative open conformer. Molecular docking reveals that the opening of the pocket would enable potential ligands to bind within it. Although lipid overlay assays suggested the PH domain binds inositol monophosphates, surface plasmon resonance demonstrated weak affinities for inositol 3,4,5-triphosphate (Ins(3,4,5)P₃; K_D ∼100 μm) and no monophosphate binding. Removing the salt bridge by site-directed mutagenesis increases the PH domain affinity for Ins(3,4,5)P₃ as measured by surface plasmon resonance and enables it to bind PtdIns(3,5)P₂ on a dot-blot. Structural comparison with other PH domains suggests that the phosphate binding pocket in the kindlin-1 PH domain is more occluded than in kindlins-2 and -3 due to its salt bridge. In addition, the apparent affinity for Ins(3,4,5)P₃ is affected by the presence of PO₄ ions in the buffer. We suggest the physiological ligand of the kindlin-1 PH domain is most likely not an inositol phosphate but another phosphorylated species.     

Structure of the GcpE (IspG)-MEcPP complex from Thermus thermophilus

Isoprenoid precursor biosynthesis occurs through the mevalonate or the methylerythritol phosphate (MEP) pathway, used i.e., by humans and by many human pathogens, respectively. In the MEP pathway, 2-C-methyl-d-erythritol-2,4-cyclo-diphosphate (MEcPP) is converted to (E)-1-hydroxy-2-methyl-but-2-enyl-4-diphosphate (HMBPP) by the iron-sulfur cluster enzyme HMBPP synthase (GcpE). The presented X-ray structure of the GcpE-MEcPP complex from Thermus thermophilus at 1.55 Å resolution provides valuable information about the catalytic mechanism and for rational inhibitor design. MEcPP binding inside the TIM-barrel funnel induces a 60° rotation of the [4Fe-4S] cluster containing domain onto the TIM-barrel entrance. The apical iron of the [4Fe-4S] cluster ligates with the C3 oxygen atom of MEcPP.     

A specific interdomain interaction preserves the structural and binding properties of the ModA protein from the phytopathogen Xanthomonas citri domain interaction and transport in ModA

The periplasmic-binding proteins in ATP-binding cassette systems (ABC Transporters) are responsible for the capture and delivery of ligands to their specific transporters, triggering a series of ATP-driven conformational changes that leads to the transport of the ligand. Structurally consisting of two lobes, the proteins change conformation after interaction with the ligand. The structure of the molybdate-binding protein (ModA) from Xanthomonas citri, bound to molybdate, was previously solved by our group and an interdomain interaction, mediated by a salt bridge between K127 and D59, apparently supports the binding properties and keeps the domains closed. To determinate the importance of this interaction, we built two ModA mutants, K127S and D59A, and analysed their functional and structural properties. Based on a set of spectroscopic experiments, crystallisation trials, structure determination and molecular dynamics (MD) simulations, we showed that the salt bridge is essential to maintain the structure and binding properties. Additionally, the MD simulations revealed that this mutant adopted a more compact structure that packed down the ligand-binding pocket. From the closed bound to open structure, the positioning of the helices forming the dipole and the salt bridge are essential to induce an intermediate state.     

Comparative analysis of structural and dynamic properties of the loaded and unloaded hemophore HasA: functional implications

A heme-acquisition system present in several Gram-negative bacteria requires the secretion of hemophores. These extracellular carrier proteins capture heme and deliver it to specific outer membrane receptors. The Serratia marcescens HasA hemophore is a monodomain protein that binds heme with a very high affinity. Its α/β structure, as that of its binding pocket, has no common features with other iron- or heme-binding proteins. Heme is held by two loops L1 and L2 and coordinated to iron by an unusual ligand pair, H32/Y75. Two independent regions of the hemophore β-sheet are involved in HasA-HasR receptor interaction. Here, we report the 3-D NMR structure of apoHasA and the backbone dynamics of both loaded and unloaded hemophore. While the overall structure of HasA is very similar in the apo and holo forms, the hemophore presents a transition from an open to a closed form upon ligand binding, through a large movement, of up to 30 Å, of loop L1 bearing H32. Comparison of loaded and unloaded HasA dynamics on different time scales reveals striking flexibility changes in the binding pocket. We propose a mechanism by which these structural and dynamic features provide the dual function of heme binding and release to the HasR receptor.     

Influence of Cholesterol and Its Stereoisomers on Members of the Serotonin Receptor Family

Despite the ubiquity of cholesterol within the cell membrane, the mechanism by which it influences embedded proteins remains elusive. Numerous G-protein coupled receptors exhibit dramatic responses to membrane cholesterol with regard to the ligand-binding affinity and functional properties, including the 5-HT receptor family. Here, we use over 25 μs of unbiased atomistic molecular dynamics simulations to identify cholesterol interaction sites in the 5-HT_1B and 5-HT_2B receptors and evaluate their impact on receptor structure. Susceptibility to membrane cholesterol is shown to be subtype dependent and determined by the quality of interactions between the extracellular loops. Charged residues are essential for maintaining the arrangement of the extracellular surface in 5-HT_2B; in the absence of such interactions, the extracellular surface of the 5-HT_1B is malleable, populating a number of distinct conformations. Elevated cholesterol density near transmembrane helix 4 is considered to be conducive to the conformation of extracellular loop 2. Occupation of this site is also shown to be stereospecific, illustrated by differential behavior of nat-cholesterol isomers, ent- and epi-cholesterol. In simulations containing the endogenous agonist, serotonin, cholesterol binding at transmembrane helix 4 biases bound serotonin molecules toward an unexpected binding mode in the extended binding pocket. The results highlight the capability of membrane cholesterol to influence the mobility of the extracellular surface in the 5-HT₁ receptor family and manipulate the architecture of the extracellular ligand-binding pocket.     

Brownian dynamics theory for predicting internal and external blockages of tetraethylammonium in the KcsA potassium channel

The theory of Brownian dynamics is used to model permeation and the blocking of KcsA potassium channels by tetraethylammonium (TEA). A novel Brownian dynamics simulation algorithm is implemented that comprises two free energy profiles; one profile is seen by the potassium ions and the other by the TEA molecules whose shape is approximated by a sphere. Our simulations reveal that internally applied TEA blocks the passage of K⁺ ions by physically occluding the pore. A TEA molecule in the external reservoir encounters an attractive energy-well created by four tyrosine residues at position 82, in addition to all other attractive and repulsive forces impinging on it. Using Brownian dynamics, we investigate how deep the energy-well needs to be to reproduce the experimentally determined inhibitory constant k_i for the TEA blockade of KcsA or the mutant Shaker T449Y. The one-dimensional free energy profile obtained from molecular dynamics is first converted into a one-dimensional potential energy profile, and is then transformed into a three-dimensional free energy profile in Brownian dynamics by adding the short-range potential from the channel walls. When converted, the free energy profile calculated from molecular dynamics gives a well-depth of ∼10 kT. We systematically alter the depths of the profiles, and then use Brownian dynamics simulations to numerically determine the current versus TEA-concentration curves. We show that the sequence of binding and unbinding events of the TEA molecule to the binding pocket can be modeled by a first-order Markov process. The Brownian dynamics simulations also reveal that the probability of a TEA molecule binding to the binding pocket in KcsA potassium channels increases exponentially with TEA concentration and depends also on the applied potential and the K⁺ concentration in the simulation assembly.     

Electrostatic interactions determine entrance/release order of substrates in the catalytic cycle of adenylate kinase

Adenylate kinase is a monomeric phosphotransferase with important biological function in regulating concentration of adenosine triphosphate (ATP) in cells, by transferring the terminal phosphate group from ATP to adenosine monophosphate (AMP) and forming two adenosine diphosphate (ADP) molecules. During this reaction, the kinase may undergo a large conformational transition, forming different states with its substrates. Although many structures of the protein are available, atomic details of the whole process remain unclear. In this article, we use both conventional molecular dynamics (MD) simulation and an enhanced sampling technique called parallel cascade selection MD simulation to explore different conformational states of the Escherichia coli adenylate kinase. Based on the simulation results, we propose a possible entrance/release order of substrates during the catalytic cycle. The substrate-free protein prefers an open conformation, but changes to a closed state once ATP·Mg enters into its binding pocket first and then AMP does. After the reaction of ATP transferring the terminal phosphate group to AMP, ADP·Mg and ADP are released sequentially, and finally the whole catalyze cycle is completed. Detailed contact and distance analysis reveals that the entrance/release order of substrates may be largely controlled by electrostatic interactions between the protein and the substrates.     

In silico experiments of single-chain antibody fragment against drugs of abuse

Three sets of in silico experiments have been conducted to elucidate the binding mechanics of two drugs, (+)-methamphetamine (METH) and amphetamine (AMP) to the single-chain variable fragment (scFv) recently engineered from anti-METH monoclonal antibody mAb6H4 (IgG, κlight chain, K_d = 11 nM). The first set of in silico experiments are long time equilibration runs of scFv:drug complexes and of drug-free scFv both in the solution. They demonstrate how the solution structures of scFv deviate from its crystallographic form with or without drug molecules bound to it. They lead to the prediction that the Arrhenius activation barrier is nearly zero for transitions from the dissociated state to the bound state. The second set of in silico experiments are nonequilibrium dynamics of pulling the drug molecules out of the binding pocket of scFv and the equilibration runs for drugs to fall back into the binding pocket. They demonstrate that extra water molecules (in addition to the two crystallographic waters) exist inside the binding pocket, underneath the drug molecules. These extra waters must have been evaporated from the binding pockets during the crystallization process of the in vitro experiments of structural determination. The third set of in silico experiments are nonequilibrium steered molecular dynamics simulations to determine the absolute binding free energies of METH and AMP to scFv. The center of mass of a drug molecule (METH or AMP) is steered (pulled) towards (forward) and away from (reverse) the binding site, sampling forward and reverse pulling paths. Mechanic work is measured along the pulling paths. The work measurements are averaged through the Brownian dynamics fluctuation dissipation theorem to produce the free-energy profiles of the scFv:drug complexes as a function of the drug-scFv separation. These experiments lead to the theoretical prediction of absolute binding energies of METH and AMP that are in agreement with the in vitro experimental results.     

Analysis of the structure and dynamics of human serum albumin

A computational approach to designing a peptide-based ligand for the purification of human serum albumin (HSA) was undertaken using molecular docking and molecular dynamics (MD) simulation. A three-step procedure was performed to design a specific ligand for HSA. Based on the candidate pocket structure of HSA (warfarin binding site), a peptide library was built. These peptides were then docked into the pocket of HSA using the GOLD program. The GOLDscore values were used to determine the affinity of peptides for HSA. Consequently, the dipeptide Trp–Trp, which shows a high GOLDscore value, was selected and linked to a spacer arm of Lys[CO(CH₂)₅NH] on the surface of ECH-lysine sepharose 4 gel. For further evaluation, the Autodock Vina program was used to dock the linked compound into the pocket of HSA. The docking simulation was performed to obtain a first guess of the binding structure of the spacer–Trp–Trp–HSA complex and subsequently analyzed by MD simulations to assess the reliability of the docking results. These MD simulations indicated that the ligand–HSA complex remains stable, and water molecules can bridge between the ligand and the protein by hydrogen bonds. Finally, absorption spectroscopic studies were performed to illustrate the appropriateness of the binding affinity of the designed ligand toward HSA. These studies demonstrate that the designed dipeptide can bind preferentially to the warfarin binding site. Graphical Abstract

Three-step computational approach to the design of a dipeptide ligand for human serum albumin purification exploiting structure-based docking and molecular dynamics simulation     

Identification of new Hsp90 inhibitors by structure-based virtual screening

Structure-based virtual screening identified pyrimidine-2,4,6-trione and 4H-1,2,4-triazole-3-thiol as novel scaffolds of Hsp90 ATPase inhibitors. Their binding modes in the ATP-binding pocket of Hsp90 were analyzed using AutoDoc program combined with molecular dynamics (MD) simulations.     

MD simulations reveal alternate conformations of the oxyanion hole in the Zika virus NS2B/NS3 protease

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple 100-ns conventional MD simulations in both the presence and absence of inhibitors. The S1, S2, and S3 pockets are stable as well. However, in two of eight simulations, the A132-G133 peptide bond in the binding pocket of S1' spontaneously flips to form a 3₁₀-helix that corresponds to the inactive conformation of the oxyanion hole, and then maintains this conformation until the end of the 100-ns conventional MD simulations without inversion of the flip. This conformational change affects the S1' pocket in ZIKV NS2B/NS3 protease active site, which is important for small molecule binding. The simulation results provide evidence at the atomic level that the inactive conformation of the oxyanion hole is more favored energetically when no specific interactions are formed between substrate/inhibitor and oxyanion hole residues. Interestingly, however, transition between the active and inactive conformation of the oxyanion hole can be observed by boosting the valley potential in accelerated MD simulations. This supports a proposed induced-fit mechanism of ZIKV NS2B/NS3 protease from computational methods and provides useful direction to enhance inhibitor binding predictions in structure-based drug design.