Enzyme dynamics point to stepwise conformational selection in catalysis

Deciphering the Mystery in p300 Taz2-p53 TAD2 Recognition

Intrinsically disordered proteins (IDPs) engage in various fundamental biological activities, and their behavior is of particular importance for a better understanding of the verbose but well-organized signal transduction in cells. IDPs exhibit uniquely paradoxical features with low affinity but simultaneously high specificity in recognizing their binding targets. The transcription factor p53 plays a crucial role in cancer suppression, carrying out some of its biological functions using its disordered regions, such as N-terminal transactivation domain 2 (TAD2). Exploration of the binding and unbinding processes between proteins is challenging, and the inherently disordered properties of these regions further complicate the issue. Computer simulations are a powerful tool to complement the experiments to fill gaps to explore the binding/unbinding processes between proteins. Here, we investigated the binding mechanism between p300 Taz2 and p53 TAD2 through extensive molecular dynamics (MD) simulations using the physics-based UNited RESidue (UNRES) force field with additional Go̅-like potentials. Distance restraints extracted from the NMR-resolved structures were imposed on intermolecular residue pairs to accelerate binding simulations, in which Taz2 was immobilized in a native-like conformation and disordered TAD2 was fully free. Starting from six structures with TAD2 placed at different positions around Taz2, we observed a metastable intermediate state in which the middle helical segment of TAD2 is anchored in the binding pocket, highlighting the significance of the TAD2 helix in directing protein recognition. Physics-based binding simulations show that successful binding is achieved after a series of stages, including (1) protein collisions to initiate the formation of encounter complexes, (2) partial attachment of TAD2, and finally (3) full attachment of TAD2 to the correct binding pocket of Taz2. Furthermore, machine-learning-based PathDetect-SOM was used to identify two binding pathways, the encounter complexes, and the intermediate states.

Selection of a promiscuous minimalist cAMP phosphodiesterase from a library of de novo designed proteins

The ability of unevolved amino acid sequences to become biological catalysts was key to the emergence of life on Earth. However, billions of years of evolution separate complex modern enzymes from their simpler early ancestors. To probe how unevolved sequences can develop new functions, we use ultrahigh-throughput droplet microfluidics to screen for phosphoesterase activity amidst a library of more than one million sequences based on a de novo designed 4-helix bundle. Characterization of hits revealed that acquisition of function involved a large jump in sequence space enriching for truncations that removed >40% of the protein chain. Biophysical characterization of a catalytically active truncated protein revealed that it dimerizes into an α-helical structure, with the gain of function accompanied by increased structural dynamics. The identified phosphodiesterase is a manganese-dependent metalloenzyme that hydrolyses a range of phosphodiesters. It is most active towards cyclic AMP, with a rate acceleration of ~109 and a catalytic proficiency of >1014 M-1, comparable to larger enzymes shaped by billions of years of evolution.

A unified view on enzyme catalysis by cryo-EM study of a DNA topoisomerase

The theories for substrate recognition in enzyme catalysis have evolved from lock-key to induced fit, then conformational selection, and conformational selection followed by induced fit. However, the prevalence and consensus of these theories require further examination. Here we use cryogenic electron microscopy and African swine fever virus type 2 topoisomerase (AsfvTop2) to demonstrate substrate binding theories in a joint and ordered manner: catalytic selection by the enzyme, conformational selection by the substrates, then induced fit. The apo-AsfvTop2 pre-exists in six conformers that comply with the two-gate mechanism directing DNA passage and release in the Top2 catalytic cycle. The structures of AsfvTop2-DNA-inhibitor complexes show that substantial induced-fit changes occur locally from the closed apo-conformer that however is too far-fetched for the open apo-conformer. Furthermore, the ATPase domain of AsfvTop2 in the MgAMP-PNP-bound crystal structures coexist in reduced and oxidized forms involving a disulfide bond, which can regulate the AsfvTop2 function.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Ratcheting synthesis

Synthetic chemistry has traditionally relied on reactions between reactants of high chemical potential and transformations that proceed energetically downhill to either a global or local minimum (thermodynamic or kinetic control). Catalysts can be used to manipulate kinetic control, lowering activation energies to influence reaction outcomes. However, such chemistry is still constrained by the shape of one-dimensional reaction coordinates. Coupling synthesis to an orthogonal energy input can allow ratcheting of chemical reaction outcomes, reminiscent of the ways that molecular machines ratchet random thermal motion to bias conformational dynamics. This fundamentally distinct approach to synthesis allows multi-dimensional potential energy surfaces to be navigated, enabling reaction outcomes that cannot be achieved under conventional kinetic or thermodynamic control. In this Review, we discuss how ratcheted synthesis is ubiquitous throughout biology and consider how chemists might harness ratchet mechanisms to accelerate catalysis, drive chemical reactions uphill and programme complex reaction sequences.

Reactivity Tracking of an Enzyme Progress Coordinate

The reactivity of individual solvent-coupled protein configurations is used to track and resolve the progress coordinate for the core reaction sequence of substrate radical rearrangement and hydrogen atom transfer in the ethanolamine ammonia-lyase (EAL) enzyme from Salmonella enterica. The first-order decay of the substrate radical intermediate is the monitored reaction. Heterogeneous confinement from sucrose hydrates in the mesophase solvent surrounding the cryotrapped protein introduces distributed kinetics in the non-native decay of the substrate radical pair capture substate, which arise from an ensemble of configurational microstates. Reaction rates increase by >103-fold across the distribution to approach that for the native enabled substate for radical rearrangement, which reacts with monotonic kinetics. The native progress coordinate thus involves a collapse of the configuration space to generate optimized reactivity. Reactivity tracking reveals fundamental features of solvent-protein-reaction configurational coupling and leads to a model that refines the ensemble paradigm of enzyme catalysis for strongly adiabatic chemical steps.

Combining Protein Conformational Diversity and Phylogenetic Information Using CoDNaS and CoDNaS-Q

CoDNaS (http://ufq.unq.edu.ar/codnas/) and CoDNaS-Q (http://ufq.unq.edu.ar/codnasq) are repositories of proteins with different degrees of conformational diversity. Following the ensemble nature of the native state, conformational diversity represents the structural differences between the conformers in the ensemble. Each entry in CoDNaS and CoDNaS-Q contains a redundant collection of experimentally determined conformers obtained under different conditions. These conformers represent snapshots of the protein dynamism. While CoDNaS contains examples of conformational diversity at the tertiary level, a recent development, CoDNaS-Q, contains examples at the quaternary level. In the emerging age of accurate protein structure prediction by machine learning approaches, many questions remain open regarding the characterization of protein dynamism. In this context, most bioinformatics resources take advantage of distinct features derived from protein alignments, however, the complexity and heterogeneity of information makes it difficult to recover reliable biological signatures. Here we present five protocols to explore tertiary and quaternary conformational diversity at the individual protein level as well as for the characterization of the distribution of conformational diversity at the protein family level in a phylogenetic context. These protocols can provide curated protein families with experimentally known conformational diversity, facilitating the exploration of sequence determinants of protein dynamism. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Assessing conformational diversity with CoDNaS Alternate Protocol 1: Assessing conformational diversity at the quaternary level with CoDNaS-Q Basic Protocol 2: Exploring conformational diversity in a protein family Alternate Protocol 2: Exploring quaternary conformational diversity in a protein family Basic Protocol 3: Representing conformational diversity in a phylogenetic context.

Dynamical activation of function in metalloenzymes

Formulations of hydrogen tunneling in enzyme-catalysed C-H activation reactions indicate enthalpic barriers to reaction that are independent of chemical steps and dependent on the protein scaffold. A tool to identify catalytically relevant site-specific protein thermal networks has emerged from temperature-dependent hydrogen deuterium exchange (TDHDX). Focusing on mutant enzyme forms with altered activation energies for catalysis, TDHDX provides a comparative analysis of the impact of mutation on Ea for local protein unfolding. Identified thermal networks appear unrelated to protein scaffold conservation and track to the dictates of the catalysed reaction, including sites for metal binding. The positions of thermal networks provide a framework for further understanding of time-dependent, functionally relevant protein motions. Measurement of nanosecond Stokes shifts at the surface of the thermal network in soybean lipoxygenase yields activation energies that are identical to Ea values measured for kcat . This finding identifies a rapid (> nanosecond), long-range and cooperative structural reorganization as the thermal barrier to catalysis. A model for protein dynamics is put forward that integrates broadly distributed protein conformational sampling with protein embedded thermal networks.

Structure-Based Function and Regulation of NCX Variants: Updates and Challenges

The plasma-membrane homeostasis Na⁺/Ca²⁺ exchangers (NCXs) mediate Ca²⁺ extrusion/entry to dynamically shape Ca²⁺ signaling/in biological systems ranging from bacteria to humans. The NCX gene orthologs, isoforms, and their splice variants are expressed in a tissue-specific manner and exhibit nearly 10⁴-fold differences in the transport rates and regulatory specificities to match the cell-specific requirements. Selective pharmacological targeting of NCX variants could benefit many clinical applications, although this intervention remains challenging, mainly because a full-size structure of eukaryotic NCX is unavailable. The crystal structure of the archaeal NCX_Mj, in conjunction with biophysical, computational, and functional analyses, provided a breakthrough in resolving the ion transport mechanisms. However, NCX_Mj (whose size is nearly three times smaller than that of mammalian NCXs) cannot serve as a structure-dynamic model for imitating high transport rates and regulatory modules possessed by eukaryotic NCXs. The crystal structures of isolated regulatory domains (obtained from eukaryotic NCXs) and their biophysical analyses by SAXS, NMR, FRET, and HDX-MS approaches revealed structure-based variances of regulatory modules. Despite these achievements, it remains unclear how multi-domain interactions can decode and integrate diverse allosteric signals, thereby yielding distinct regulatory outcomes in a given ortholog/isoform/splice variant. This article summarizes the relevant issues from the perspective of future developments.

Temperature-dependent hydrogen deuterium exchange shows impact of analog binding on adenosine deaminase flexibility but not embedded thermal networks

The analysis of hydrogen deuterium exchange by mass spectrometry as a function of temperature and mutation has emerged as a generic and efficient tool for the spatial resolution of protein networks that are proposed to function in the thermal activation of catalysis. In this work, we extend temperature-dependent hydrogen deuterium exchange from apo-enzyme structures to protein-ligand complexes. Using adenosine deaminase as a prototype, we compared the impacts of a substrate analog (1-deaza-adenosine) and a very tight-binding inhibitor/transition state analog (pentostatin) at single and multiple temperatures. At a single temperature, we observed different hydrogen deuterium exchange-mass spectrometry properties for the two ligands, as expected from their 106-fold differences in strength of binding. By contrast, analogous patterns for temperature-dependent hydrogen deuterium exchange mass spectrometry emerge in the presence of both 1-deaza-adenosine and pentostatin, indicating similar impacts of either ligand on the enthalpic barriers for local protein unfolding. We extended temperature-dependent hydrogen deuterium exchange to a function-altering mutant of adenosine deaminase in the presence of pentostatin and revealed a protein thermal network that is highly similar to that previously reported for the apo-enzyme (Gao et al., 2020, JACS 142, 19936-19949). Finally, we discuss the differential impacts of pentostatin binding on overall protein flexibility versus site-specific thermal transfer pathways in the context of models for substrate-induced changes to a distributed protein conformational landscape that act in synergy with embedded protein thermal networks to achieve efficient catalysis.

Identification of novel functional mini-receptors by combinatorial screening of split-WW domains

β-Sheet motifs such as the WW domain are increasingly being explored as building blocks for synthetic biological applications. Since the sequence-structure relationships of β-sheet motifs are generally complex compared to the well-studied α-helical coiled coil (CC), other approaches such as combinatorial screening should be included to vary the function of the peptide. In this study, we present a combinatorial approach to identify novel functional mini-proteins based on the WW-domain scaffold, which takes advantage of the successful reconstitution of the fragmented WW domain of hPin1 (hPin1WW) by CC association. Fragmentation of hPin1WW was performed in both loop 1 (CC-hPin1WW-L1) and loop 2 (CC-hPin1WW-L2), and the respective fragments were linked to the strands of an antiparallel heterodimeric CC. Structural analysis by CD and NMR spectroscopy revealed structural reconstitution of the WW-domain scaffold only in CC-hPin1WW-L1, but not in CC-hPin1WW-L2. Furthermore, by using 1H-15N HSQC NMR, fluorescence and CD spectroscopy, we demonstrated that binding properties of fragmented hPin1WW in CC-hPin1WW-L1 were fully restored by CC association. To demonstrate the power of this approach as a combinatorial screening platform, we synthesized a four-by-six library of N- and C-terminal hPin1WW-CC peptide fragments that was screened for a WW domain that preferentially binds to ATP over cAMP, phophocholine, or IP6. Using this screening platform, we identified one WW domain, which specifically binds ATP, and a phosphorylcholine-specific WW-based mini-receptor, both having binding dissociation constants in the lower micromolar range.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

Chemical fuels for molecular machinery

Chemical reaction networks that transform out-of-equilibrium 'fuel' to 'waste' are the engines that power the biomolecular machinery of the cell. Inspired by such systems, autonomous artificial molecular machinery is being developed that functions by catalysing the decomposition of chemical fuels, exploiting kinetic asymmetry to harness energy released from the fuel-to-waste reaction to drive non-equilibrium structures and dynamics. Different aspects of chemical fuels profoundly influence their ability to power molecular machines. Here we consider the structure and properties of the fuels that biology has evolved and compare their features with those of the rudimentary synthetic chemical fuels that have so far been used to drive autonomous non-equilibrium molecular-level dynamics. We identify desirable, but context-specific, traits for chemical fuels together with challenges and opportunities for the design and invention of new chemical fuels to power synthetic molecular machinery and other dissipative nanoscale processes.

Leveraging intrinsic flexibility to engineer enhanced enzyme catalytic activity

Dynamic motions of enzymes occurring on a broad range of timescales play a pivotal role in all steps of the reaction pathway, including substrate binding, catalysis, and product release. However, it is unknown whether structural information related to conformational flexibility can be exploited for the directed evolution of enzymes with higher catalytic activity. Here, we show that mutagenesis of residues exclusively located at flexible regions distal to the active site of Homo sapiens kynureninase (HsKYNase) resulted in the isolation of a variant (BF-HsKYNase) in which the rate of the chemical step toward kynurenine was increased by 45-fold. Mechanistic pre–steady-state kinetic analysis of the wild type and the evolved enzyme shed light on the underlying effects of distal mutations (>10 Å from the active site) on the rate-limiting step of the catalytic cycle. Hydrogen-deuterium exchange coupled to mass spectrometry and molecular dynamics simulations revealed that the amino acid substitutions in BF-HsKYNase allosterically affect the flexibility of the pyridoxal-5′-phosphate (PLP) binding pocket, thereby impacting the rate of chemistry, presumably by altering the conformational ensemble and sampling states more favorable to the catalyzed reaction.

Autonomous fuelled directional rotation about a covalent single bond

Biology operates through autonomous chemically fuelled molecular machinery¹, including rotary motors such as adenosine triphosphate synthase² and the bacterial flagellar motor³. Chemists have long sought to create analogous molecular structures with chemically powered, directionally rotating, components^4-17. However, synthetic motor molecules capable of autonomous 360° directional rotation about a single bond have proved elusive, with previous designs lacking either autonomous fuelling^7,10,12 or directionality⁶. Here we show that 1-phenylpyrrole 2,2'-dicarboxylic acid^18,19 (1a) is a catalysis-driven^20,21 motor that can continuously transduce energy from a chemical fuel^9,20-27 to induce repetitive 360° directional rotation of the two aromatic rings around the covalent N-C bond that connects them. On treatment of 1a with a carbodiimide^21,25-27, intramolecular anhydride formation between the rings and the anhydride's hydrolysis both occur incessantly. Both reactions are kinetically gated^28-30 causing directional bias. Accordingly, catalysis of carbodiimide hydration by the motor molecule continuously drives net directional rotation around the N-C bond. The directionality is determined by the handedness of both an additive that accelerates anhydride hydrolysis and that of the fuel, and is easily reversed additive³¹. More than 97% of fuel molecules are consumed through the chemical engine cycle²⁴ with a directional bias of up to 71:29 with a chirality-matched fuel and additive. In other words, the motor makes a 'mistake' in direction every three to four turns. The 26-atom motor molecule's simplicity augurs well for its structural optimization and the development of derivatives that can be interfaced with other components for the performance of work and tasks^32-36.

Insights of conformational dynamics on catalytic activity in the computational stability design of Bacillus subtilis LipA

In protein engineering, the contributions of individual mutations to designed combinatorial mutants are unpredictable. Screening designed mutations that affect enzyme catalytic activity enables evolutions towards efficient activities. Here, Bacillus subtilis LipA (BSLA) was selected as a model protein for thermostabilization designs, and the circular dichroism measurements showed six combinatorial designs with improved stability (from 5.81 °C to 13.61 °C). Based on molecular dynamic simulations, the conformational dynamics of the mutants revealed that mutations alter the populations of conformational states and the increased ensembles of inactive conformations might lead to a reduction in activity. We further demonstrated that the mutations responsible for the reduced enzyme catalytic activity involved a short dynamic correlation path to disturbing the equilibrium conformation of active sites. By removing N82V, which had a close dynamic correlation to the active sites in mutant D3, the redesigned mutant RD3 had an increased activity of 57.6%. By combining computational simulation with experimental verification, this work established that essential sites to counteract the activity-stability trade-off in multipoint combinatorial mutants could be computationally predicted and thus provide a possible strategy by which to indirectly or directly guide protein design.

Resolution and characterization of contributions of select protein and coupled solvent configurational fluctuations to radical rearrangement catalysis in coenzyme B12-dependent ethanolamine ammonia-lyase

Coenzyme B12 (adenosylcobalamin) -dependent ethanolamine ammonia-lyase (EAL) is the signature enzyme in ethanolamine utilization metabolism associated with microbiome homeostasis and disease conditions in the human gut. The enzyme conducts a complex choreography of bond-making/bond-breaking steps that rearrange substrate to products through a radical mechanism, with themes common to other coenzyme B12-dependent and radical enzymes. The methods presented are targeted to test the hypothesis that particular, select protein and coupled solvent configurational fluctuations contribute to enzyme function. The general approach is to correlate enzyme function with an introduced perturbation that alters the properties (for example, degree of concertedness, or collectiveness) of protein and coupled solvent dynamics. Methods for sample preparation and low-temperature kinetic measurements by using temperature-step reaction initiation and time-resolved, full-spectrum electron paramagnetic resonance spectroscopy are detailed. A framework for interpretation of results obtained in ensemble systems under conditions of statistical equilibrium within the reacting, globally unstable state is presented. The temperature-dependence of the first-order rate constants for decay of the cryotrapped paramagnetic substrate radical state in EAL, through the chemical step of radical rearrangement, displays a piecewise-continuous Arrhenius dependence from 203 to 295K, punctuated by a kinetic bifurcation over 219-220K. The results reveal the obligatory contribution of a class of select collective protein and coupled solvent fluctuations to the interconversion of two resolved, sequential configurational substates, on the decay time scale. The select class of collective fluctuations also contributes to the chemical step. The methods and analysis are generally applicable to other coenzyme B12-dependent and related radical enzymes.

Myricetin Allosterically Inhibits the Dengue NS2B-NS3 Protease by Disrupting the Active and Locking the Inactive Conformations

The dengue NS2B-NS3 protease existing in equilibrium between the active and inactive forms is essential for virus replication, thus representing a key drug target. Here, myricetin, a plant flavonoid, was characterized to noncompetitively inhibit the dengue protease. Further NMR study identified the protease residues perturbed by binding to myricetin, which were utilized to construct the myricetin-protease complexes. Strikingly, in the active form, myricetin binds to a new allosteric site (AS2) far away from the active site pocket and the allosteric site (AS1) for binding curcumin, while in the inactive form, it binds to both AS1 and AS2. To decipher the mechanism for the allosteric inhibition by myricetin, we conducted molecular dynamics simulations on different forms of dengue NS2B-NS3 proteases. Unexpectedly, the binding of myricetin to AS2 is sufficient to disrupt the active conformation by displacing the characteristic NS2B C-terminal β-hairpin from the active site pocket. By contrast, the binding of myricetin to AS1 and AS2 results in locking the inactive conformation. Therefore, myricetin represents the first small molecule, which allosterically inhibits the dengue protease by both disrupting the active conformation and locking the inactive conformation. The results enforce the notion that a global allosteric network exists in the dengue NS2B-NS3 protease, which is susceptible to allosteric inhibition by small molecules such as myricetin and curcumin. As myricetin has been extensively used as a food additive, it might be directly utilized to fight the dengue infections and as a promising starting material for further design of potent allosteric inhibitors.

Lysine-Dependent Entropy Effects in the B. subtilis Lysine Riboswitch: Insights from Single-Molecule Thermodynamic Studies

Riboswitches play an important role in RNA-based sensing/gene regulation control for many bacteria. In particular, the accessibility of multiple conformational states at physiological temperatures allows riboswitches to selectively bind a cognate ligand in the aptamer domain, which triggers secondary structural changes in the expression platform, and thereby "switching" between on or off transcriptional or translational states for the downstream RNA. The present work exploits temperature-controlled, single-molecule total internal reflection fluorescence (TIRF) microscopy to study the thermodynamic landscape of such ligand binding/folding processes, specifically for the Bacillus subtilis lysine riboswitch. The results confirm that the riboswitch folds via an induced-fit (IF) mechanism, in which cognate lysine ligand first binds to the riboswitch before structural rearrangement takes place. The transition state to folding is found to be enthalpically favored (ΔH_fold^‡ < 0), yet with a free-energy barrier that is predominantly entropic (-TΔS_fold^‡ > 0), which results in folding (unfolding) rate constants strongly dependent (independent) of lysine concentration. Analysis of the single-molecule kinetic "trajectories" reveals this rate constant dependence of k_fold on lysine to be predominantly entropic in nature, with the additional lysine conferring preferential advantage to the folding process by the presence of ligands correctly oriented with respect to the riboswitch platform. By way of contrast, van't Hoff analysis reveals enthalpic contributions to the overall folding thermodynamics (ΔH⁰) to be surprisingly constant and robustly independent of lysine concentration. The results demonstrate the crucial role of hydrogen bonding between the ligand and riboswitch platform but with only a relatively modest fraction (45%) of the overall enthalpy change needed to access the transition state and initiate transcriptional switching.

GTP Cyclohydrolase I as a Potential Drug Target: New Insights into Its Allosteric Modulation via Normal Mode Analysis

Guanosine triphosphate (GTP) cyclohydrolase I (GCH1) catalyzes the conversion of GTP into dihydroneopterin triphosphate (DHNP). DHNP is the first intermediate of the folate de novo biosynthesis pathway in prokaryotic and lower eukaryotic microorganisms and the tetrahydrobiopterin (BH4) biosynthesis pathway in higher eukaryotes. The de novo folate biosynthesis provides essential cofactors for DNA replication, cell division, and synthesis of key amino acids in rapidly replicating pathogen cells, such as Plasmodium falciparum (P. falciparum), a causative agent of malaria. In eukaryotes, the product of the BH4 biosynthesis pathway is essential for the production of nitric oxide and several neurotransmitter precursors. An increased copy number of the malaria parasite P. falciparum GCH1 gene has been reported to influence antimalarial antifolate drug resistance evolution, whereas mutations in the human GCH1 are associated with neuropathic and inflammatory pain disorders. Thus, GCH1 stands as an important and attractive drug target for developing therapeutics. The GCH1 intrinsic dynamics that modulate its activity remains unclear, and key sites that exert allosteric effects across the structure are yet to be elucidated. This study employed the anisotropic network model to analyze the intrinsic motions of the GCH1 structure alone and in complex with its regulatory partner protein. We showed that the GCH1 tunnel-gating mechanism is regulated by a global shear motion and an outward expansion of the central five-helix bundle. We further identified hotspot residues within sites of structural significance for the GCH1 intrinsic allosteric modulation. The obtained results can provide a solid starting point to design novel antineuropathic treatments for humans and novel antimalarial drugs against the malaria parasite P. falciparum GCH1 enzyme.

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

Dimerization of PHGDH via the catalytic unit is essential for its enzymatic function

Human D-3-phosphoglycerate dehydrogenase (PHGDH), a key enzyme in de novo serine biosynthesis, is amplified in various cancers and serves as a potential target for anticancer drug development. To facilitate this process, more information is needed on the basic biochemistry of this enzyme. For example, PHGDH was found to form tetramers in solution and the structure of its catalytic unit (sPHGDH) was solved as a dimer. However, how the oligomeric states affect PHGDH enzyme activity remains elusive. We studied the dependence of PHGDH enzymatic activity on its oligomeric states. We found that sPHGDH forms a mixture of monomers and dimers in solution with a dimer dissociation constant of ∼0.58 μM, with the enzyme activity depending on the dimer content. We computationally identified hotspot residues at the sPHGDH dimer interface. Single-point mutants at these sites disrupt dimer formation and abolish enzyme activity. Molecular dynamics simulations showed that dimer formation facilitates substrate binding and maintains the correct conformation required for enzyme catalysis. We further showed that the full-length PHGDH exists as a dynamic mixture of monomers, dimers, and tetramers in solution with enzyme concentration-dependent activity. Mutations that can completely disrupt the sPHGDH dimer show different abilities to interrupt the full-length PHGDH tetramer. Among them, E108A and I121A can also disrupt the oligomeric structures of the full-length PHGDH and abolish its enzyme activity. Our study indicates that disrupting the oligomeric structure of PHGDH serves as a novel strategy for PHGDH drug design and the hotspot residues identified can guide the design process.

Conformational Ensemble of TteAdoCbl Riboswitch Provides Stable Structural Elements for Conformation Selection and Population Shift in Cobalamin Recognition

Cobalamin riboswitch is a cis-regulatory element widely found in the 5'-UTRs of the vitamin B12-associated genes in bacteria, resulting in modulation and production of a particular protein. Thermoanaerobacter tengcongensis (Tte) AdoCbl riboswitches are the largest of the known riboswitches with 210 nucleotides, partially due to its long peripheral P6-extension, which enable high affinity of AdoCbl. Two structural elements, T-loop/T-looplike motif and kissing loop are key to the global folding of the RNA. While the structure of the TteAdoCbl riboswitch complex is known, we still do not understand the structure and conformation before AdoCbl ligand recognition. In order to delineate the conformational changes and the stabilities of long-range interactions, we have performed extensive all-atom replica-exchange molecular dynamics simulations of the TteAdoCbl riboswitch with a total simulation time of 2296 ns. We found that both the T-loop/T-looplike motif and kissing loop are very stable with ligand binding. The gating conformation changes of P6-extension allow the ligand to bind to the preorganized kissing loop binding pocket. The T-loop/T-looplike motif has much more hydrogen bonds than observed in TteAdoCbl riboswitch complex crystal structure, indicating an allosteric response of the T-loop/T-looplike motif. Our study demonstrated that the conformational ensemble of TteAdoCbl riboswitch provides stable structural elements for conformation selection and population shift in cobalamin recognition.

Drug discovery for enzymes

Enzymes are essential, physiological catalysts involved in all processes of life, including metabolism, cellular signaling and motility, as well as cell growth and division. They are attractive drug targets because of the presence of defined substrate-binding pockets, which can be exploited as binding sites for pharmaceutical enzyme inhibitors. Understanding the reaction mechanisms of enzymes and the molecular mode of action of enzyme inhibitors is indispensable for the discovery and development of potent, efficacious, and safe novel drugs. The combination of classical concepts of enzymology with new experimental and data analysis methods opens new routes for drug discovery.

Conformational Dynamics of Deubiquitinase A and Functional Implications

Deubiquitinase A (DUBA) belongs to the ovarian tumor family of deubiquitinating enzymes and was initially identified as a negative regulator of type I interferons, whose overproduction has been linked to autoimmune diseases. The deubiquitinating activity of DUBA is positively regulated by phosphorylation at a single serine residue, S177, which results in minimal structural changes. We have previously shown that phosphorylation induces a two-state conformational equilibrium observed only in the active form of DUBA, highlighting the functional importance of DUBA dynamics. Here, we report the conformational dynamics of DUBA on the microsecond-to-millisecond time scales characterized by nuclear magnetic resonance relaxation dispersion experiments. We found that motions on these time scales are highly synchronized in the phosphorylated and nonphosphorylated DUBA. Despite the overall similarity of these two forms, different dynamic properties were observed in helix α1 and the neighboring regions, including residue S177, which likely contribute to the activation of DUBA by phosphorylation. Moreover, our data suggest that transient unfolding of helix α6 drives the global conformational process and that mutations can be introduced to modulate this process, which provides a basis for future studies to define the exact functional roles of motions in DUBA activation and substrate specificity.

Curcumin Allosterically Inhibits the Dengue NS2B-NS3 Protease by Disrupting Its Active Conformation

Flaviviruses including dengue virus and Zika virus encode a unique two-component NS2B-NS3 protease essential for maturation/infectivity, thus representing a key target for designing antiflavivirus drugs. Here, for the first time, by NMR and molecular docking, we reveal that curcumin allosterically inhibits the dengue protease by binding to a cavity with no overlap with the active site. Further molecular dynamics simulations decode that the binding of curcumin leads to unfolding/displacing the characteristic β-hairpin of the C-terminal NS2B and consequently disrupting the closed (active) conformation of the protease. Our study identified a cavity most likely conserved in all flaviviral NS2B-NS3 proteases, which could thus serve as a therapeutic target for the discovery/design of small-molecule allosteric inhibitors. Moreover, as curcumin has been used as a food additive for thousands of years in many counties, it can be directly utilized to fight the flaviviral infections and as a promising starting for further design of potent allosteric inhibitors.

Schistosomal Sulfotransferase Interaction with Oxamniquine Involves Hybrid Mechanism of Induced-fit and Conformational Selection

Background:

Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention.

Objective:

The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations.

Methods:

Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA.

Results:

The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol.

Conclusion:

The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity-enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.

Protein-Protein Binding as a Two-Step Mechanism: Preselection of Encounter Poses during the Binding of BPTI and Trypsin

Biomolecular recognition between proteins follows complex mechanisms, the understanding of which can substantially advance drug discovery efforts. Here, we track each step of the binding process in atomistic detail with molecular dynamics simulations using trypsin and its inhibitor bovine pancreatic trypsin inhibitor (BPTI) as a model system. We use umbrella sampling to cover a range of unbinding pathways. Starting from these simulations, we subsequently seed classical simulations at different stages of the process and combine them to a Markov state model. We clearly identify three kinetically separated states (an unbound state, an encounter state, and the final complex) and describe the mechanisms that dominate the binding process. From our model, we propose the following sequence of events. The initial formation of the encounter complex is driven by long-range interactions because opposite charges in trypsin and BPTI draw them together. The encounter complex features the prealigned binding partners with binding sites still partially surrounded by solvation shells. Further approaching leads to desolvation and increases the importance of van der Waals interactions. The native binding pose is adopted by maximizing short-range interactions. Thereby side-chain rearrangements ensure optimal shape complementarity. In particular, BPTI’s P1 residue adapts to the S1 pocket and prime site residues reorient to optimize interactions. After the paradigm of conformation selection, binding-competent conformations of BPTI and trypsin are already present in the apo ensembles and their probabilities increase during this proposed two-step association process. This detailed characterization of the molecular forces driving the binding process includes numerous aspects that have been discussed as central to the binding of trypsin and BPTI and protein complex formation in general. In this study, we combine all these aspects into one comprehensive model of protein recognition. We thereby contribute to enhance our general understanding of this fundamental mechanism, which is particularly critical as the development of biopharmaceuticals continuously gains significance.

Enzyme-Substrate-Cofactor Dynamical Networks Revealed by High-Resolution Field Cycling Relaxometry

The remarkable power and specificity of enzyme catalysis rely on the dynamic alignment of the enzyme, substrates, and cofactors, yet the role of dynamics has usually been approached from the perspective of the protein. We have been using an underappreciated NMR technique, subtesla high-resolution field cycling 31P NMR relaxometry, to investigate the dynamics of the enzyme-bound substrates and cofactor on guanosine-5'-monophosphate reductase (GMPR). GMPR forms two dead end, yet catalytically competent, complexes that mimic distinct steps in the catalytic cycle: E·IMP·NADP+ undergoes a partial hydride transfer reaction, while E·GMP·NADP+ undergoes a partial deamination reaction. A different cofactor conformation is required for each partial reaction. Here we report the effects of mutations designed to perturb cofactor conformation and ammonia binding with the goal of identifying the structural features that contribute to the distinct dynamic signatures of the hydride transfer and deamination complexes. These experiments suggest that Asp129 is a central cog in a dynamic network required for both hydride transfer and deamination. In contrast, Lys77 modulates the conformation and mobility of substrates and cofactors in a reaction-specific manner. Thr105 and Tyr318 are part of a deamination-specific dynamic network that includes the 2'-OH of GMP. These residues have comparatively little effect on the dynamic properties of the hydride transfer complex. These results further illustrate the potential of high-resolution field cycling NMR relaxometry for the investigation of ligand dynamics. In addition, exchange experiments indicate that NH3/NH4+ has a high affinity for the deamination complex but a low affinity for the hydride transfer complex, suggesting that the movement of ammonia may gate the cofactor conformational change. Collectively, these experiments reinforce the view that the enzyme, substrates, and cofactor are linked in intricate, reaction-specific, dynamic networks and demonstrate that distal portions of the substrates and cofactors are critical features in these networks.

Kramers’ Theory and the Dependence of Enzyme Dynamics on Trehalose-Mediated Viscosity

The disaccharide trehalose is accumulated in the cytoplasm of some organisms in response to harsh environmental conditions. Trehalose biosynthesis and accumulation are important for the survival of such organisms by protecting the structure and function of proteins and membranes. Trehalose affects the dynamics of proteins and water molecules in the bulk and the protein hydration shell. Enzyme catalysis and other processes dependent on protein dynamics are affected by the viscosity generated by trehalose, as described by the Kramers’ theory of rate reactions. Enzyme/protein stabilization by trehalose against thermal inactivation/unfolding is also explained by the viscosity mediated hindering of the thermally generated structural dynamics, as described by Kramers’ theory. The analysis of the relationship of viscosity–protein dynamics, and its effects on enzyme/protein function and other processes (thermal inactivation and unfolding/folding), is the focus of the present work regarding the disaccharide trehalose as the viscosity generating solute. Finally, trehalose is widely used (alone or in combination with other compounds) in the stabilization of enzymes in the laboratory and in biotechnological applications; hence, considering the effect of viscosity on catalysis and stability of enzymes may help to improve the results of trehalose in its diverse uses/applications.

Molecular dynamics simulation study of Plasmodium falciparum and Escherichia coli SufA: Exploration of conformational changes possibly involved in iron-sulfur cluster transfer

Iron-sulfur (Fe-S) clusters are one of the earliest known metal complexes in biological molecules. Suf system is one of the Fe-S biogenesis pathways. SufA belongs to the Suf pathway. It is an A-type carrier protein that transfers Fe-S clusters from the scaffold to target proteins. Structural studies were performed for the Suf pathway protein, SufA, in order to explore the conformational changes that probably aid in the transfer of Fe-S clusters to target proteins. Three-dimensional (3D) structure of Plasmodium falciparum (Pf) SufA homodimer was obtained by homology modeling using 3D structure of Escherichia coli (Ec) SufA as template. Molecular dynamics (MD) simulation of Pf SufA and Ec SufA homodimers followed by trajectory and pocket analyses were carried out. A co-ordinated displacement of the homodimeric chains in the interfacial region, resembling a swinging trapeze-like movement was observed. Potential involvement of this swinging trapeze-like movement of the residues belonging to the interfacial region has been proposed as a probable mechanism that assists in the transfer of Fe-S cluster from SufA to apo proteins. This was substantiated by protein-protein interaction studies in Pf SufA by performing molecular docking of 3D conformations of Pf SufA obtained from MD trajectory at every 1 ns interval with Pf ferredoxin.Communicated by Ramaswamy H. Sarma.

Building a macro-mixing dual-basin Gō model using the Multistate Bennett Acceptance Ratio

The dual-basin Gō-model is a structural-based coarsegrained model for simulating a conformational transition between two known structures of a protein. Two parameters are required to produce a dual-basin potential mixed using two single-basin potentials, although the determination of mixing parameters is usually not straightforward. Here, we have developed an efficient scheme to determine the mixing parameters using the Multistate Bennett Acceptance Ratio (MBAR) method after short simulations with a set of parameters. In the scheme, MBAR allows us to predict observables at various unsimulated conditions, which are useful to improve the mixing parameters in the next round of iterative simulations. The number of iterations that are necessary for obtaining the converged mixing parameters are significantly reduced in the scheme. We applied the scheme to two proteins, the glutamine binding protein and the ribose binding protein, for showing the effectiveness in the parameter determination. After obtaining the converged parameters, both proteins show frequent conformational transitions between open and closed states, providing the theoretical basis to investigate structure-dynamics-function relationships of the proteins.

Regional conformational flexibility couples substrate specificity and scissile phosphate diester selectivity in human flap endonuclease 1

Human flap endonuclease-1 (hFEN1) catalyzes the divalent metal ion-dependent removal of single-stranded DNA protrusions known as flaps during DNA replication and repair. Substrate selectivity involves passage of the 5'-terminus/flap through the arch and recognition of a single nucleotide 3'-flap by the α2-α3 loop. Using NMR spectroscopy, we show that the solution conformation of free and DNA-bound hFEN1 are consistent with crystal structures; however, parts of the arch region and α2-α3 loop are disordered without substrate. Disorder within the arch explains how 5'-flaps can pass under it. NMR and single-molecule FRET data show a shift in the conformational ensemble in the arch and loop region upon addition of DNA. Furthermore, the addition of divalent metal ions to the active site of the hFEN1-DNA substrate complex demonstrates that active site changes are propagated via DNA-mediated allostery to regions key to substrate differentiation. The hFEN1-DNA complex also shows evidence of millisecond timescale motions in the arch region that may be required for DNA to enter the active site. Thus, hFEN1 regional conformational flexibility spanning a range of dynamic timescales is crucial to reach the catalytically relevant ensemble.

Interaction vs Preorganization in Enzyme Catalysis. A Dispute That Calls for Resolution

This essay focuses on the debate between Warshel et al. (proponents of preorganization) and Menger and Nome (proponents of spatiotemporal effects) over the source of fast enzyme catalysis. The Warshel model proposes that the main function of enzymes is to push the solvent coordinate toward the transition state. Other physical-organic factors (e.g., desolvation, entropic effects, ground state destabilization, etc.) do not, ostensibly, contribute substantially to the rate. Indeed, physical organic chemistry in its entirety was claimed to be "irrelevant to an enzyme's active site". Preorganization had been applied by Warshel to his "flagship" enzyme, ketosteroid isomerase, but we discuss troubling issues with their ensuing analysis. For example, the concepts of "general acid" and "general base", known to play a role in this enzyme's mechanism, are ignored in the text. In contrast, the spatiotemporal theory postulates that enzyme-like rates (i.e., accelerations >10⁸) occur when two functionalities are held rigidly at contact distances less than ca. 3 Å. Numerous diverse organic systems are shown to bear this out experimentally. Many of these are intramolecular systems where distances between functionalities are known. Among them are fast intramolecular systems where strain is actually generated during the reaction, thereby excluding steric compression as a source of the observed enzyme-like rates. Finally, the account ends with structural data from four active sites of enzymes, obtained by others, all showing contact distances between substrate analogues and enzyme. To our knowledge, contact distances less than the diameter of water are found universally among enzymes, and it is to this fact that we attribute their extremely fast rates given the assumption that enzymes, whatever their particular mechanism, obey elementary chemical principles.

Oxidative damage to epigenetically methylated sites affects DNA stability, dynamics and enzymatic demethylation

DNA damage can affect various regulatory elements of the genome, with the consequences for DNA structure, dynamics, and interaction with proteins remaining largely unexplored. We used solution NMR spectroscopy, restrained and free molecular dynamics to obtain the structures and investigate dominant motions for a set of DNA duplexes containing CpG sites permuted with combinations of 5-methylcytosine (mC), the primary epigenetic base, and 8-oxoguanine (oxoG), an abundant DNA lesion. Guanine oxidation significantly changed the motion in both hemimethylated and fully methylated DNA, increased base pair breathing, induced BI→BII transition in the backbone 3′ to the oxoG and reduced the variability of shift and tilt helical parameters. UV melting experiments corroborated the NMR and molecular dynamics results, showing significant destabilization of all methylated contexts by oxoG. Notably, some dynamic and thermodynamic effects were not additive in the fully methylated oxidized CpG, indicating that the introduced modifications interact with each other. Finally, we show that the presence of oxoG biases the recognition of methylated CpG dinucleotides by ROS1, a plant enzyme involved in epigenetic DNA demethylation, in favor of the oxidized DNA strand. Thus, the conformational and dynamic effects of spurious DNA oxidation in the regulatory CpG dinucleotide can have far-reaching biological consequences.

Modulating Hinge Flexibility in the APP Transmembrane Domain Alters γ-Secretase Cleavage

Intramembrane cleavage of the β-amyloid precursor protein C99 substrate by γ-secretase is implicated in Alzheimer's disease pathogenesis. Biophysical data have suggested that the N-terminal part of the C99 transmembrane domain (TMD) is separated from the C-terminal cleavage domain by a di-glycine hinge. Because the flexibility of this hinge might be critical for γ-secretase cleavage, we mutated one of the glycine residues, G38, to a helix-stabilizing leucine and to a helix-distorting proline. Both mutants impaired γ-secretase cleavage and also altered its cleavage specificity. Circular dichroism, NMR, and backbone amide hydrogen/deuterium exchange measurements as well as molecular dynamics simulations showed that the mutations distinctly altered the intrinsic structural and dynamical properties of the substrate TMD. Although helix destabilization and/or unfolding was not observed at the initial ε-cleavage sites of C99, subtle changes in hinge flexibility were identified that substantially affected helix bending and twisting motions in the entire TMD. These resulted in altered orientation of the distal cleavage domain relative to the N-terminal TMD part. Our data suggest that both enhancing and reducing local helix flexibility of the di-glycine hinge may decrease the occurrence of enzyme-substrate complex conformations required for normal catalysis and that hinge mobility can thus be conducive for productive substrate-enzyme interactions.

Discovery of processive catalysis by an exo-hydrolase with a pocket-shaped active site

Substrates associate and products dissociate from enzyme catalytic sites rapidly, which hampers investigations of their trajectories. The high-resolution structure of the native Hordeum exo-hydrolase HvExoI isolated from seedlings reveals that non-covalently trapped glucose forms a stable enzyme-product complex. Here, we report that the alkyl β-d-glucoside and methyl 6-thio-β-gentiobioside substrate analogues perfused in crystalline HvExoI bind across the catalytic site after they displace glucose, while methyl 2-thio-β-sophoroside attaches nearby. Structural analyses and multi-scale molecular modelling of nanoscale reactant movements in HvExoI reveal that upon productive binding of incoming substrates, the glucose product modifies its binding patterns and evokes the formation of a transient lateral cavity, which serves as a conduit for glucose departure to allow for the next catalytic round. This path enables substrate-product assisted processive catalysis through multiple hydrolytic events without HvExoI losing contact with oligo- or polymeric substrates. We anticipate that such enzyme plasticity could be prevalent among exo-hydrolases.

Enzyme substrates and products often diffuse too rapidly to assess the catalytic implications of these movements. Here, the authors characterise the structural basis of product and substrate diffusion for an exo-hydrolase and discover a substrate-product assisted processive catalytic mechanism.

Molecular dynamics based improvement of the solubilizing self-cleavable tag Zbasic-ΔI-CM application in the preparation of recombinant proteins in Escherichia coli

Z_basic-ΔI-CM is a novel intein-based self-cleavable tag we developed to accelerate the soluble expression of recombinant proteins in Escherichia coli (E. coli). Previously we found that intein activity could be interfered by its flanking exteins, and thus reducing the production efficiency and final yield. In this work, we used CXC-chemokine 9 (CXCL9) as a model C-extein, which fusion with Z_basic-ΔI-CM showed high intein activity. When the fusion protein got soluble expression, CXCL9 was released immediately and purified directly from cell lysis supernatant. The results demonstrated that Z_basic-ΔI-CM tag had successfully mediated the efficient production of high-quality CXCL9 with reduced time and resources consumption in comparison with inclusion bodies expression. Molecular dynamics simulations suggested that the improved cleavage activity of Z_basic-ΔI-CM upon fusion with CXCL9 may be due to the higher dynamics of the first half loop and stabilization of the second half loop of intein. Our results proved that the self-cleavable Z_basic-ΔI-CM mediated soluble expression could be a feasible process for cytokines like CXCL9, thus of attractive potentials for production of therapeutic proteins using E. coli expression system.

Increased H-Bond Stability Relates to Altered ε-Cleavage Efficiency and Aβ Levels in the I45T Familial Alzheimer's Disease Mutant of APP

Cleavage of the amyloid precursor protein's (APP) transmembrane domain (TMD) by γ-secretase is a crucial step in the aetiology of Alzheimer's Disease (AD). Mutations in the APP TMD alter cleavage and lead to familial forms of AD (FAD). The majority of FAD mutations shift the preference of initial cleavage from ε49 to ε48, thus raising the AD-related Aβ42/Aβ40 ratio. The I45T mutation is among the few FAD mutations that do not alter ε-site preference, while it dramatically reduces the efficiency of ε-cleavage. Here, we investigate the impact of the I45T mutation on the backbone dynamics of the substrate TMD. Amide exchange experiments and molecular dynamics simulations in solvent and a lipid bilayer reveal an increased stability of amide hydrogen bonds at the ζ- and γ-cleavage sites. Stiffening of the H-bond network is caused by an additional H-bond between the T45 side chain and the TMD backbone, which alters dynamics within the cleavage domain. In particular, the increased H-bond stability inhibits an upward movement of the ε-sites in the I45T mutant. Thus, an altered presentation of ε-sites to the active site of γ-secretase as a consequence of restricted local flexibility provides a rationale for reduced ε-cleavage efficiency of the I45T mutant.

Studying the phosphoryl transfer mechanism of the E. coli phosphofructokinase-2: from X-ray structure to quantum mechanics/molecular mechanics simulations

Phosphofructokinases catalyze the ATP-dependent phosphorylation of fructose-6-phosphate and they are highly regulated.

Phosphofructokinases (Pfks) catalyze the ATP-dependent phosphorylation of fructose-6-phosphate (F6P) and they are regulated in a wide variety of organisms. Although numerous aspects of the kinetics and regulation have been characterized for Pfks, the knowledge about the mechanism of the phosphoryl transfer reaction and the transition state lags behind. In this work, we describe the X-ray crystal structure of the homodimeric Pfk-2 from E. coli, which contains products in one site and reactants in the other, as well as an additional ATP molecule in the inhibitory allosteric site adjacent to the reactants. This complex was previously predicted when studying the kinetic mechanism of ATP inhibition. After removing the allosteric ATP, molecular dynamic (MD) simulations revealed conformational changes related to domain packing, as well as stable interactions of Lys27 and Asp256 with donor (ATP) and acceptor (fructose-6-) groups, and of Asp166 with Mg²⁺. The phosphoryl transfer reaction mechanism catalyzed by Pfk-2 was investigated through Quantum Mechanics/Molecular Mechanics (QM/MM) simulations using a combination of the string method and a path-collective variable for the exploration of its free energy surface. The calculated activation free energies showed that a dissociative mechanism, occurring with a metaphosphate intermediate formation followed by a proton transfer to Asp256, is more favorable than an associative one. The structural analysis reveals the role of Asp256 acting as a catalytic base and Lys27 stabilizing the transition state of the dissociative mechanism.

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Elucidation of conformational diversity of druggable enzymes and classification of chemical modulators based on inhibitor-bound structures

Communicated by Ramaswamy H. Sarma.

Structurally- and dynamically-driven allostery of the chymotrypsin-like proteases of SARS, Dengue and Zika viruses

Coronavirus 3C-like and Flavivirus NS2B-NS3 proteases utilize the chymotrypsin fold to harbor their catalytic machineries but also contain additional domains/co-factors. Over the past decade, we aimed to decipher how the extra domains/co-factors mediate the catalytic machineries of SARS 3C-like, Dengue and Zika NS2B-NS3 proteases by characterizing their folding, structures, dynamics and inhibition with NMR, X-ray crystallography and MD simulations, and the results revealed: 1) the chymotrypsin fold of the SARS 3C-like protease can independently fold, while, by contrast, those of Dengue and Zika proteases lack the intrinsic capacity to fold without co-factors. 2) Mutations on the extra domain of SARS 3C-like protease can transform the active catalytic machinery into the inactive collapsed state by structurally-driven allostery. 3) Amazingly, even without detectable structural changes, mutations on the extra domain are sufficient to either inactivate or enhance the catalytic machinery of SARS 3C-like protease by dynamically-driven allostery. 4) Global networks of correlated motions have been identified: for SARS 3C-like protease, N214A inactivates the catalytic machinery by decoupling the network, while STI/A and STIF/A enhance by altering the patterns of the network. The global networks of Dengue and Zika proteases are coordinated by their NS2B-cofactors. 5) Natural products were identified to allosterically inhibit Zika and Dengue proteases through binding a pocket on the back of the active site. Therefore, by introducing extra domains/cofactors, nature develops diverse strategies to regulate the catalytic machinery embedded on the chymotrypsin fold through folding, structurally- and dynamically-driven allostery, all of which might be exploited to develop antiviral drugs.

Evolutionary repurposing of a sulfatase: A new Michaelis complex leads to efficient transition state charge offset

The recruitment and evolutionary optimization of promiscuous enzymes is key to the rapid adaptation of organisms to changing environments. Our understanding of the precise mechanisms underlying enzyme repurposing is, however, limited: What are the active-site features that enable the molecular recognition of multiple substrates with contrasting catalytic requirements? To gain insights into the molecular determinants of adaptation in promiscuous enzymes, we performed the laboratory evolution of an arylsulfatase to improve its initially weak phenylphosphonate hydrolase activity. The evolutionary trajectory led to a 100,000-fold enhancement of phenylphosphonate hydrolysis, while the native sulfate and promiscuous phosphate mono- and diester hydrolyses were only marginally affected (≤50-fold). Structural, kinetic, and in silico characterizations of the evolutionary intermediates revealed that two key mutations, T50A and M72V, locally reshaped the active site, improving access to the catalytic machinery for the phosphonate. Measured transition state (TS) charge changes along the trajectory suggest the creation of a new Michaelis complex (E•S, enzyme-substrate), with enhanced leaving group stabilization in the TS for the promiscuous phosphonate (β_leavinggroup from -1.08 to -0.42). Rather than altering the catalytic machinery, evolutionary repurposing was achieved by fine-tuning the molecular recognition of the phosphonate in the Michaelis complex, and by extension, also in the TS. This molecular scenario constitutes a mechanistic alternative to adaptation solely based on enzyme flexibility and conformational selection. Instead, rapid functional transitions between distinct chemical reactions rely on the high reactivity of permissive active-site architectures that allow multiple substrate binding modes.

Dissecting conformational changes in APP's transmembrane domain linked to ε-efficiency in familial Alzheimer's disease

The mechanism by which familial Alzheimer's disease (FAD) mutations within the transmembrane domain (TMD) of the Amyloid Precursor Protein (APP) affect ε-endoproteolysis is only poorly understood. Thereby, mutations in the cleavage domain reduce ε-efficiency of γ-secretase cleavage and some even shift entry into production lines. Since cleavage occurs within the TMD, a relationship between processing and TMD structure and dynamics seems obvious. Using molecular dynamic simulations, we dissect the dynamic features of wild-type and seven FAD-mutants into local and global components. Mutations consistently enhance hydrogen-bond fluctuations upstream of the ε-cleavage sites but maintain strong helicity there. Dynamic perturbation-response scanning reveals that FAD-mutants target backbone motions utilized in the bound state. Those motions, obscured by large-scale motions in the pre-bound state, provide (i) a dynamic mechanism underlying the proposed coupling between binding and ε-cleavage, (ii) key sites consistent with experimentally determined docking sites, and (iii) the distinction between mutants and wild-type.