Large Conformational Changes Research Articles

Architecture and regulation of filamentous human cystathionine beta-synthase

Cystathionine beta-synthase (CBS) is an essential metabolic enzyme across all domains of life for the production of glutathione, cysteine, and hydrogen sulfide. Appended to the conserved catalytic domain of human CBS is a regulatory domain that modulates activity by S-adenosyl-L-methionine (SAM) and promotes oligomerisation. Here we show using cryo-electron microscopy that full-length human CBS in the basal and SAM-bound activated states polymerises as filaments mediated by a conserved regulatory domain loop. In the basal state, CBS regulatory domains sterically block the catalytic domain active site, resulting in a low-activity filament with three CBS dimers per turn. This steric block is removed when in the activated state, one SAM molecule binds to the regulatory domain, forming a high-activity filament with two CBS dimers per turn. These large conformational changes result in a central filament of SAM-stabilised regulatory domains at the core, decorated with highly flexible catalytic domains. Polymerisation stabilises CBS and reduces thermal denaturation. In PC-3 cells, we observed nutrient-responsive CBS filamentation that disassembles when methionine is depleted and reversed in the presence of SAM. Together our findings extend our understanding of CBS enzyme regulation, and open new avenues for investigating the pathogenic mechanism and therapeutic opportunities for CBS-associated disorders.

Using atomistic solution scattering modelling to elucidate the role of the Fc glycans in human IgG4.

Human immunoglobulin G (IgG) exists as four subclasses IgG1-4, each of which has two Fab subunits joined by two hinges to a Fc subunit. IgG4 has the shortest hinge with 12 residues. The Fc subunit has two glycan chains, but the importance of glycosylation is not fully understood in IgG4. Here, to evaluate the stability and structure of non-glycosylated IgG4, we performed a multidisciplinary structural study of glycosylated and deglycosylated human IgG4 A33 for comparison with our similar study of human IgG1 A33. After deglycosylation, IgG4 was found to be monomeric by analytical ultracentrifugation; its sedimentation coefficient of 6.52 S was reduced by 0.27 S in reflection of its lower mass. X-ray and neutron solution scattering showed that the overall Guinier radius of gyration RG and its cross-sectional values after deglycosylation were almost unchanged. In the P(r) distance distribution curves, the two M1 and M2 peaks that monitor the two most common distances within IgG4 were unchanged following deglycosylation. Further insight from Monte Carlo simulations for glycosylated and deglycosylated IgG4 came from 111,382 and 117,135 possible structures respectively. Their comparison to the X-ray and neutron scattering curves identified several hundred best-fit models for both forms of IgG4. Principal component analyses showed that glycosylated and deglycosylated IgG4 exhibited different conformations from each other. Within the constraint of unchanged RG and M1-M2 values, the glycosylated IgG4 models showed more restricted Fc conformations compared to deglycosylated IgG4, but no other changes. Kratky plots supported this interpretation of greater disorder upon deglycosylation, also observed in IgG1. Overall, these more variable Fc conformations may demonstrate a generalisable impact of deglycosylation on Fc structures, but with no large conformational changes in IgG4 unlike those seen in IgG1.

A model for cooperativity of kinesin-4 motors by communicating through the microtubule track

Kinesin-4 motor proteins at nanomolar concentrations can accumulate at microtubule plus ends, regulating the microtubule length. Recent experimental studies showed that at picomolar motor concentrations, under which nearly no motor crowding can occur, the run length of the kinesin-4 motor increases whereas the velocity decreases with the motor concentration, indicating the presence of a long-range communication between the motors through the microtubule track. However, the mechanism of the long-range communication is illusive. Here, based on the proposal that large conformational changes of the tubulin induced by the strong interaction with a kinesin head can propagate along the microtubule such that influence the tubulin-binding energy of other motors distant away, the dynamics of the kinesin-4 motor at the picomolar motor concentrations is studied theoretically, explaining quantitatively the experimental data. Moreover, the lifetime of the motor when multiple motors are jammed near the plus ends at high motor concentrations is studied theoretically.

Structural insights into thermophilic chaperonin complexes

Group I chaperonins are dual heptamer protein complexes that play significant roles in protein homeostasis. The structure and function of the Escherichia coli chaperonin are well characterized. However, the dynamic properties of chaperonins, such as large ATPase-dependent conformational changes by binding of lid-like co-chaperonin GroES, have made structural analyses challenging, and our understanding of these changes during the turnover of chaperonin complex formation is limited. In this study, we used single-particle cryogenic electron microscopy to investigate the structures of GroES-bound chaperonin complexes from the thermophilic hydrogen-oxidizing bacteria Hydrogenophilus thermoluteolus and Hydrogenobacter thermophilus in the presence of ATP and AMP-PNP. We captured the structure of an intermediate state chaperonin complex, designated as an asymmetric football-shaped complex, and performed analyses to decipher the dynamic structural variations. Our structural analyses of inter- and intra-subunit communications revealed a unique mechanism of complex formation through the binding of a second GroES to a bullet-shaped complex.

Multiscale Study on the Intramolecular C-S Bond Formation Catalyzed by P450 Monooxygenase CxnD Involved in the Biosynthesis of Chuangxinmycin: The Critical Roles of Noncrystal Water Molecule and Conformational Change.

Cytochrome P450 monooxygenase CxnD catalyzes intramolecular C-S bond formation in the biosynthesis of chuangxinmycin, which is representative of the synthesis of sulfur-containing natural heterocyclic compounds. The intramolecular cyclization usually requires the activation of two reaction sites and a large conformational change; thus, illuminating its detailed reaction mechanism remains challengeable. Here, the reaction pathway of CxnD-catalyzed C-S bond formation was clarified by a series of calculations, including Gaussian accelerated molecular dynamics simulations and quantum mechanical-molecular mechanical calculations. Our results revealed that the C-S formation follows a diradical coupling mechanism. CxnD first employs Cpd I to abstract the hydrogen atom from the imino group of the indole ring, and then, the resulted Cpd II further extracts another hydrogen atom from the thiol group of the side chain to afford a diradical intermediate, in which a noncrystal water molecule entering into the active site after the formation of Cpd I was proved to play an indispensable role. Moreover, the diradical intermediate cannot directly perform the coupling reaction. It should first undergo a series of conformational changes leading to the proximity of two reaction sites. It is the flexibility of the active site of the enzyme and the side chain of the substrate that makes the diradical coupling to be successful.

Myomerger Induces Membrane Hemifusion and Regulates Fusion Pore Expansion.

Membrane fusion is a crucial mechanism in a wide variety of important events in cell biology from viral infection to exocytosis. However, despite many efforts and much progress, cell-cell fusion has remained elusive to our understanding. Along the life of the fusion pore, large conformational changes take place from the initial lipid bilayer bending, passing through the hemifusion intermediates, and ending with the formation of the first nascent fusion pore. In this sense, computer simulations are an ideal technique for describing such complex lipid remodeling at the molecular level. In this work, we studied the role played by the muscle-specific membrane protein Myomerger during the formation of the fusion pore. We have conducted μs length atomistic and coarse-grained molecular dynamics, together with free-energy calculations using ad hoc collective variables. Our results show that Myomerger favors the hemifusion diaphragm-stalk transition, reduces the nucleation-expansion energy difference, and promotes the formation of nonenlarging fusion pores.

Machine Learning Methods as a Cost-Effective Alternative to Physics-Based Binding Free Energy Calculations.

The rank ordering of ligands remains one of the most attractive challenges in drug discovery. While physics-based in silico binding affinity methods dominate the field, they still have problems, which largely revolve around forcefield accuracy and sampling. Recent advances in machine learning have gained traction for protein-ligand binding affinity predictions in early drug discovery programs. In this article, we perform retrospective binding free energy evaluations for 172 compounds from our internal collection spread over four different protein targets and five congeneric ligand series. We compared multiple state-of-the-art free energy methods ranging from physics-based methods with different levels of complexity and conformational sampling to state-of-the-art machine-learning-based methods that were available to us. Overall, we found that physics-based methods behaved particularly well when the ligand perturbations were made in the solvation region, and they did not perform as well when accounting for large conformational changes in protein active sites. On the other end, machine-learning-based methods offer a good cost-effective alternative for binding free energy calculations, but the accuracy of their predictions is highly dependent on the experimental data available for training the model.

PH-dependent reaction triggering in PmHMGR crystals for time-resolved crystallography

Serial crystallography and time-resolved data collection can readily be employed to investigate the catalytic mechanism of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl (HMG)-coenzyme-A (CoA) reductase (PmHMGR) by changing the environmental conditions in the crystal and so manipulating the reaction rate. This enzyme uses a complex mechanism to convert mevalonate to HMG-CoA using the co-substrate CoA and cofactor NAD+. The multi-step reaction mechanism involves an exchange of bound NAD+ and large conformational changes by a 50-residue subdomain. The enzymatic reaction can be run in both forward and reverse directions in solution and is catalytically active in the crystal for multiple reaction steps. Initially, the enzyme was found to be inactive in the crystal starting with bound mevalonate, CoA, and NAD+. To observe the reaction from this direction, we examined the effects of crystallization buffer constituents and pH on enzyme turnover, discovering a strong inhibition in the crystallization buffer and a controllable increase in enzyme turnover as a function of pH. The inhibition is dependent on ionic concentration of the crystallization precipitant ammonium sulfate but independent of its ionic composition. Crystallographic studies show that the observed inhibition only affects the oxidation of mevalonate but not the subsequent reactions of the intermediate mevaldehyde. Calculations of the pKa values for the enzyme active site residues suggest that the effect of pH on turnover is due to the changing protonation state of His381. We have now exploited the changes in ionic inhibition in combination with the pH-dependent increase in turnover as a novel approach for triggering the PmHMGR reaction in crystals and capturing information about its intermediate states along the reaction pathway.

DynamicBind: predicting ligand-specific protein-ligand complex structure with a deep equivariant generative model

While significant advances have been made in predicting static protein structures, the inherent dynamics of proteins, modulated by ligands, are crucial for understanding protein function and facilitating drug discovery. Traditional docking methods, frequently used in studying protein-ligand interactions, typically treat proteins as rigid. While molecular dynamics simulations can propose appropriate protein conformations, they’re computationally demanding due to rare transitions between biologically relevant equilibrium states. In this study, we present DynamicBind, a deep learning method that employs equivariant geometric diffusion networks to construct a smooth energy landscape, promoting efficient transitions between different equilibrium states. DynamicBind accurately recovers ligand-specific conformations from unbound protein structures without the need for holo-structures or extensive sampling. Remarkably, it demonstrates state-of-the-art performance in docking and virtual screening benchmarks. Our experiments reveal that DynamicBind can accommodate a wide range of large protein conformational changes and identify cryptic pockets in unseen protein targets. As a result, DynamicBind shows potential in accelerating the development of small molecules for previously undruggable targets and expanding the horizons of computational drug discovery.

Tailoring the Dynamics of Excited‐State Charge Transfer through Conformational Engineering to Improve Second Near‐Infrared Fluorescence for High‐Resolution Osteosarcoma Imaging

The dynamics of excited‐state charge transfer (CT) in second near‐infrared (NIR‐II) fluorophore proceeds with conformational change that govern fluorescence. Unveiling the relationship between CT dynamics and conformational change in excited state is of great fundamental significance in optimizing NIR‐II fluorescence but remains less explored. This study unveils the pivotal role of large conformational change in promoting the evolution of emissive CT state into nonemissive CT state in NIR‐II fluorophore (CA‐BBT). Spectroscopic and computational results reveal that large conformational rotation/twisting leads to a pronounced increase in the CT character in excited state. This heightened CT character in CA‐BBT enables a rapid evolution of emissive CT state into a nonemissive CT state within 1.4 ps, as observed by ultrafast spectroscopy. Subsequently, this nonemissive state dominates nonradiative decay, resulting in very low NIR‐II fluorescence. Preventing such detrimental evolution by constraining the conformational dynamics of CA‐BBT results in a 10‐fold enhancement of NIR‐II fluorescence, enabling high‐resolution dynamic visualization of vasculature within osteosarcoma. This study offers a profound understanding of the relationship between CT dynamics and conformational changes in NIR‐II fluorescence, presenting innovative perspectives to develop bright NIR‐II fluorophores.

Nanomechanical footprint of SARS-CoV-2 variants in complex with a potent nanobody by molecular simulations.

Rational design of novel antibody therapeutics against viral infections such as coronavirus relies on surface complementarity and high affinity for their effectiveness. Here, we explore an additional property of protein complexes, the intrinsic mechanical stability, in SARS-CoV-2 variants when complexed with a potent antibody. In this study, we utilized a recent implementation of the GōMartini 3 approach to investigate large conformational changes in protein complexes with a focus on the mechanostability of the receptor-binding domain (RBD) from WT, Alpha, Delta, and XBB.1.5 variants in complex with the H11-H4 nanobody. The analysis revealed moderate differences in mechanical stability among these variants. Also, we identified crucial residues in both the RBD and certain protein segments in the nanobody that contribute to this property. By performing pulling simulations and monitoring the presence of specific native and non-native contacts across the protein complex interface, we provided mechanistic insights into the dissociation process. Force-displacement profiles indicate a tensile force clamp mechanism associated with the type of protein complex. Our computational approach not only highlights the key mechanostable interactions that are necessary to maintain overall stability, but it also paves the way for the rational design of potent antibodies that are mechanostable and effective against emergent SARS-CoV-2 variants.

Arthropod repellent interactions with olfactory receptors and ionotropic receptors analyzed by molecular modeling

The main insect chemoreceptors are olfactory receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs). The odorant binding sites of many insect ORs appear to be occluded and inaccessible from the surface of the receptor protein, based on the three-dimensional structure of OR5 from the jumping bristletail Machilis hrabei (MhraOR5) and a survey of a sample of vinegar fly (Drosophila melanogaster) OR structures obtained from artificial intellegence (A.I.) modeling. Molecular dynamics simulations revealed that the occluded site can become accessible through tunnels that transiently open and close. The present study extends this analysis to examine seventeen ORs and one GR docking with ligands that have known valence: nine that signal attraction and nine that signal aversion. All but one of the receptors displayed occluded ligand binding sites analogous to MhraOR5, and docking software predicted the known attractant and repellent ligands will bind to the occluded sites. Docking of the repellent DEET was examined, and more than half of the OR ligand sites were predicted to bind DEET, including receptors that signal aversion as well as those that signal attraction. However, DEET may not actually have access to all the attractant binding sites. The larger size and lower flexibility of repellent molecules may restrict their passage through the tunnel bottlenecks, which could act as filters to select access to the ligand binding sites. In contrast to ORs and GRs, the IR ligand binding site is in an extracellular domain known to undergo a large conformational change from an open to a closed state. A.I. models of two D. melanogaster IRs of known valence and two blacklegged tick (Ixodes scapularis) IRs having unknown ligands were computationally tested for attractant and repellent binding. The ligand-binding sites in the closed state appear inaccessible to the protein surface, so attractants and repellents must bind initially at an accessible site in the open state before triggering the conformational change. In some IRs, repellent binding sites were identified at exterior sites adjacent to the ligand-binding site. These may be allosteric sites that, when occupied by repellents, can stabilize the open state of an attractant IR, or stabilize the closed state of an IR in the absence of its activating ligand. The model of D. melanogaster IR64a suggests a possible molecular mechanism for the activation of this IR by H+. The amino acids involved in this proposed mechanism are conserved in IR64a from several Dipteran pest species and disease vectors, potentially offering a route to discovery of new repellents that act via the allosteric site.

CASP15 cryo-EM protein and RNA targets: Refinement and analysis using experimental maps.

CASP assessments primarily rely on comparing predicted coordinates with experimental reference structures. However, experimental structures by their nature are only models themselves-their construction involves a certain degree of subjectivity in interpreting density maps and translating them to atomic coordinates. Here, we directly utilized density maps to evaluate the predictions by employing a method for ranking the quality of protein chain predictions based on their fit into the experimental density. The fit-based ranking was found to correlate well with the CASP assessment scores. Overall, the evaluation against the density map indicated that the models are of high accuracy, and occasionally even better than the reference structure in some regions of the model. Local assessment of predicted side chains in a 1.52 Å resolution map showed that side-chains are sometimes poorly positioned. Additionally, the top 118 predictions associated with 9 protein target reference structures were selected for automated refinement, in addition to the top 40 predictions for 11 RNA targets. For both proteins and RNA, the refinement of CASP15 predictions resulted in structures that are close to the reference target structure. This refinement was successful despite large conformational changes often being required, showing that predictions from CASP-assessed methods could serve as a good starting point for building atomic models in cryo-EM maps for both proteins and RNA. Loop modeling continued to pose a challenge for predictors, and together with the lack of consensus amongst models in these regions suggests that modeling, in combination with model-fit to the density, holds the potential for identifying more flexible regions within the structure.

The Onset of Molecule-Spanning Dynamics in Heat Shock Protein Hsp90.

Protein dynamics have been investigated on a wide range of time scales. Nano- and picosecond dynamics have been assigned to local fluctuations, while slower dynamics have been attributed to larger conformational changes. However, it is largely unknown how fast (local) fluctuations can lead to slow global (allosteric) changes. Here, fast molecule-spanning dynamics on the 100to 200ns time scale in the heat shock protein 90 (Hsp90) are shown. Global real-space movements are assigned to dynamic modes on this time scale, which is possible by a combination of single-molecule fluorescence, quasi-elastic neutron scattering and all-atom molecular dynamics (MD) simulations. The time scale of these dynamic modes depends on the conformational state of the Hsp90 dimer. In addition, the dynamic modes are affected to various degrees by Sba1, a co-chaperone of Hsp90, depending on the location within Hsp90, which is in very good agreement with MD simulations. Altogether, this data is best described by fast molecule-spanning dynamics, which precede larger conformational changes in Hsp90 and might be the molecular basis for allostery. This integrative approach provides comprehensive insights into molecule-spanning dynamics on the nanosecond time scale for a multi-domain protein.

An allosteric redox switch involved in oxygen protection in a CO2 reductase.

Metal-dependent formate dehydrogenases reduce CO2 with high efficiency and selectivity, but are usually very oxygen sensitive. An exception is Desulfovibrio vulgaris W/Sec-FdhAB, which can be handled aerobically, but the basis for this oxygen tolerance was unknown. Here we show that FdhAB activity is controlled by a redox switch based on an allosteric disulfide bond. When this bond is closed, the enzyme is in an oxygen-tolerant resting state presenting almost no catalytic activity and very low formate affinity. Opening this bond triggers large conformational changes that propagate to the active site, resulting in high activity and high formate affinity, but also higher oxygen sensitivity. We present the structure of activated FdhAB and show that activity loss is associated with partial loss of the metal sulfido ligand. The redox switch mechanism is reversible in vivo and prevents enzyme reduction by physiological formate levels, conferring a fitness advantage during O2 exposure.

Conformation state-specific monobodies regulate the functions of flexible proteins through conformation trapping.

Synthetic binding proteins have emerged as modulators of protein functions through protein-protein interactions (PPIs). Because PPIs are influenced by the structural dynamics of targeted proteins, investigating whether the synthetic-binders-based strategy is applicable for proteins with large conformational changes is important. This study demonstrates the applicability of monobodies (fibronectin type-III domain-based synthetic binding proteins) in regulating the functions of proteins that undergo tens-of-angstroms-scale conformational changes, using an example of the A55C/C77S/V169C triple mutant (Adktm ; a phosphoryl transfer-catalyzing enzyme with a conformational change between OPEN/CLOSED forms). Phage display successfully developed monobodies that recognize the OPEN form (substrate-unbound form), but not the CLOSED form of Adktm . Two OPEN form-specific clones (OP-2 and OP-4) inhibited Adktm kinase activity. Epitope mapping with a yeast-surface display/flow cytometry indicated that OP-2 binds to the substrate-entry side of Adktm , whereas OP-4 binding occurs at another site. Small angle X-ray scattering coupled with size-exclusion chromatography (SEC-SAXS) indicated that OP-4 binds to the hinge side opposite to the substrate-binding site of Adktm , retaining the whole OPEN-form structure of Adktm . Titration of the OP-4-Adktm complex with Ap5 A, a transition-state analog of Adktm , showed that the conformational shift to the CLOSED form was suppressed although Adktm retained the OPEN-form (i.e., substrate-binding ready form). These results show that OP-4 captures and stabilizes the OPEN-form state, thereby affecting the hinge motion. These experimental results indicate that monobody-based modulators can regulate the functions of proteins that show tens-of-angstroms-scale conformational changes, by trapping specific conformational states generated during large conformational change process that is essential for function exertion.

Structural characterization of a novel cyclic 2,3-diphosphoglycerate synthetase involved in extremolyte production in the archaeon Methanothermus fervidus.

The enzyme cyclic di-phosphoglycerate synthetase that is involved in the production of the osmolyte cyclic 2,3-diphosphoglycerate has been studied both biochemically and structurally. Cyclic 2,3-diphosphoglycerate is found exclusively in the hyperthermophilic archaeal methanogens, such as Methanothermus fervidus, Methanopyrus kandleri, and Methanothermobacter thermoautotrophicus. Its presence increases the thermostability of archaeal proteins and protects the DNA against oxidative damage caused by hydroxyl radicals. The cyclic 2,3-diphosphoglycerate synthetase enzyme has been crystallized and its structure solved to 1.7 Å resolution by experimental phasing. It has also been crystallized in complex with its substrate 2,3 diphosphoglycerate and the co-factor ADP and this structure has been solved to 2.2 Å resolution. The enzyme structure has two domains, the core domain shares some structural similarity with other NTP-dependent enzymes. A significant proportion of the structure, including a 127 amino acid N-terminal domain, has no structural similarity to other known enzyme structures. The structure of the complex shows a large conformational change that occurs in the enzyme during catalytic turnover. The reaction involves the transfer of the γ-phosphate group from ATP to the substrate 2,3 -diphosphoglycerate and the subsequent SN2 attack to form a phosphoanhydride. This results in the production of the unusual extremolyte cyclic 2,3 -diphosphoglycerate which has important industrial applications.

Multistep conformational changes leading to the gate opening of light-driven sodium pump rhodopsin

Membrane transport proteins require a gating mechanism that opens and closes the substrate transport pathway to carry out unidirectional transport. The "gating" involves large conformational changes and is achieved via multistep reactions. However, these elementary steps have not been clarified for most transporters due to the difficulty of detecting the individual steps. Here, we propose these steps for the gate opening of the bacterial Na+ pump rhodopsin, which outwardly pumps Na+ upon illumination. We herein solved an asymmetric dimer structure of Na+ pump rhodopsin from the bacterium Indibacter alkaliphilus. In one protomer, the Arg108 sidechain is oriented toward the protein center and appears to block a Na+ release pathway to the extracellular (EC) medium. In the other protomer, however, this sidechain swings to the EC side and then opens the release pathway. Assuming that the latter protomer mimics the Na+-releasing intermediate, we examined the mechanism for the swing motion of the Arg108 sidechain. On the EC surface of the first protomer, there is a characteristic cluster consisting of Glu10, Glu159, and Arg242 residues connecting three helices. In contrast, this cluster is disrupted in the second protomer. Our experimental results suggested that this disruption is a key process. The cluster disruption induces the outward movement of the Glu159-Arg242 pair and simultaneously rotates the seventh transmembrane helix. This rotation resultantly opens a space for the swing motion of the Arg108 sidechain. Thus, cluster disruption might occur during the photoreaction and then trigger sequential conformation changes leading to the gate-open state.

Autoantibody binding and unique enzyme-substrate intermediate conformation of human transglutaminase 3

Transglutaminase 3 (TG3), the autoantigen of dermatitis herpetiformis (DH), is a calcium dependent enzyme that targets glutamine residues in polypeptides for either transamidation or deamidation modifications. To become catalytically active TG3 requires proteolytic cleavage between the core domain and two C-terminal β-barrels (C1C2). Here, we report four X-ray crystal structures representing inactive and active conformations of human TG3 in complex with a TG3-specific Fab fragment of a DH patient derived antibody. We demonstrate that cleaved TG3, upon binding of a substrate-mimicking inhibitor, undergoes a large conformational change as a β-sheet in the catalytic core domain moves and C1C2 detaches. The unique enzyme-substrate conformation of TG3 without C1C2 is recognized by DH autoantibodies. The findings support a model where B-cell receptors of TG3-specific B cells bind and internalize TG3-gluten enzyme-substrate complexes thereby facilitating gluten-antigen presentation, T-cell help and autoantibody production.

Structure-based discovery of pyrazole-benzothiadiazole hybrid as human HPPD inhibitors

4-Hydroxyphenylpyruvate dioxygenase (HPPD) has attracted increasing attention as a target for treating type I tyrosinemia and other diseases with defects in tyrosine catabolism. Only one commercial drug, 2-(2-nitro-4-trifluoromethylbenzoyl)-1, 3-cyclohexanedione (NTBC), clinically treat type I tyrosinemia, but show some severe side effects in clinical application. Here, we determined the structure of human HPPD-NTBC complex, and developed new pyrazole-benzothiadiazole 2,2-dioxide hybrids from the binding of NTBC. These compounds showed improved inhibition against human HPPD, among which compound a10 was the most active candidate. The Absorption Distribution Metabolism Excretion Toxicity (ADMET) predicted properties suggested that a10 had good druggability, and was with lower toxicity than NTBC. The structure comparison between inhibitor-bound and ligand-free form human HPPD showed a large conformational change of the C-terminal helix. Furthermore, the loop 1 and α7 helix were found adopting different conformations to assist the gating of the cavity, which explains the gating mechanism of human HPPD.