Testing geometrical discrimination within an enzyme active site: constrained hydrogen bonding in the ketosteroid isomerase oxyanion hole

Identifying the toxic mechanisms of emerging electronic contaminations liquid crystal monomers and the construction of a priority control list for graded control

Liquid crystal monomers (LCMs) are identified as emerging organic contaminations with largely unexplored health impacts. To elucidate their toxic mechanisms, support the establishment of environmental discharge and management standards, and promote effective LCMs control, this study constructs a database covering 20,545 potential targets of 1431 LCMs, highlighting 9 key toxic target proteins that disrupt the nervous system and metabolic functions. GO and KEGG pathway analysis suggests LCMs severely affect nervous system, linked to neurodegenerative diseases and mental health disorders, with toxicity variations driven by electronegativity and structural complexity of LCM terminal groups. To achieve tiered control of LCMs, construct toxicity risk control lists for 9 key toxic target proteins, suitable for the graded control of LCMs, management recommendations are provided based on toxicity levels. These lists were validated for reliability and offer reliable toxicity predictions for LCMs. SHAP analysis points to electronic properties, molecular shape, and structural characteristics of LCMs as primary health impact factors. As the first study integrating machine learning with computational toxicology to outline LCMs health impacts, it aims to enhance public understanding of LCM toxicity risks and support the development of environmental standards, effective management of LCM production and emissions, and reduction of public exposure risks.

Structural and functional effects of phosphopriming and scaffolding in the kinase GSK-3β

Glycogen synthase kinase 3β (GSK-3β) targets specific signaling pathways in response to distinct upstream signals. We used structural and functional studies to dissect how an upstream phosphorylation step primes the Wnt signaling component β-catenin for phosphorylation by GSK-3β and how scaffolding interactions contribute to this reaction. Our crystal structure of GSK-3β bound to a phosphoprimed β-catenin peptide confirmed the expected binding mode of the phosphoprimed residue adjacent to the catalytic site. An aspartate phosphomimic in the priming site of β-catenin adopted an indistinguishable structure but reacted approximately 1000-fold slower than the native phosphoprimed substrate. This result suggests that substrate positioning alone is not sufficient for catalysis and that native phosphopriming interactions are necessary. We also obtained a structure of GSK-3β with an extended peptide from the scaffold protein Axin that bound with greater affinity than that of previously crystallized Axin fragments. This structure neither revealed additional contacts that produce the higher affinity nor explained how substrate interactions in the GSK-3β active site are modulated by remote Axin binding. Together, our findings suggest that phosphopriming and scaffolding produce small conformational changes or allosteric effects, not captured in the crystal structures, that activate GSK-3β and facilitate β-catenin phosphorylation. These results highlight limitations in our ability to predict catalytic activity from structure and have potential implications for the role of natural phosphomimic mutations in kinase regulation and phosphosite evolution.

Insights into Modeling Approaches in Chemistry: Assessing Ligand-Protein Binding Thermodynamics Based on Rigid-Flexible Model Molecules

In the field of chemistry, model compounds find extensive use for investigating complex objects. One prime example of such object is the protein-ligand supramolecular interaction. Prediction the enthalpic and entropic contribution to the free energy associated with this process, as well as the structural and dynamic characteristics of protein-ligand complexes poses considerable challenges. This review exemplifies modeling approaches used to study protein-ligand binding (PLB) thermodynamics by employing pairs of conformationally constrained/flexible model molecules. Strategically designing the model molecules can reduce the number of variables that influence thermodynamic parameters. This enables scientists to gain deeper insights into the enthalpy and entropy of PLB, which is relevant for medicinal chemistry and drug design. The model studies reviewed here demonstrate that rigidifying ligands may induce compensating changes in the enthalpy and entropy of binding. Some "rules of thumb" have started to emerge on how to minimize entropy-enthalpy compensation and design efficient rigidified or flexible ligands.

Bridging the functional gap between reactivity and inhibition in dehaloperoxidase B from Amphitrite ornata: Mechanistic and structural studies with 2,4- and 2,6-dihalophenols

The multifunctional catalytic globin dehaloperoxidase (DHP) from the marine worm Amphitrite ornata was shown to catalyze the H₂O₂-dependent oxidation of 2,4- and 2,6-dihalophenols (DXP; X = F, Cl, Br). Product identification by LC-MS revealed multiple monomeric products with varying degrees of oxidation and/or dehalogenation, as well as oligomers with n up to 6. Mechanistic and ¹⁸O-labeling studies demonstrated sequential dihalophenol oxidation via peroxidase and peroxygenase activities. Binding studies established that 2,4-DXP (X = Cl, Br) have the highest affinities of any known DHP substrate. X-ray crystallography identified different binding positions for 2,4- and 2,6-DXP substrates in the hydrophobic distal pocket of DHP. Correlation between the number of halogens and the substrate binding orientation revealed a halogen-dependent binding motif for mono- (4-halophenol), di- (2,4- and 2,6-dihalophenol) and trihalophenols (2,4,6-trihalopenol). Taken together, the findings here on dihalophenol reactivity with DHP advance our understanding of how these compounds bridge the inhibitory and oxidative functions of their mono- and trihalophenol counterparts, respectively, and provide further insight into the protein structure-function paradigm relevant to multifunctional catalytic globins in comparison to their monofunctional analogs.

High throughput and quantitative enzymology in the genomic era

Accurate predictions from models based on physical principles are the ultimate metric of our biophysical understanding. Although there has been stunning progress toward structure prediction, quantitative prediction of enzyme function has remained challenging. Realizing this goal will require large numbers of quantitative measurements of rate and binding constants and the use of these ground-truth data sets to guide the development and testing of these quantitative models. Ground truth data more closely linked to the underlying physical forces are also desired. Here, we describe technological advances that enable both types of ground truth measurements. These advances allow classic models to be tested, provide novel mechanistic insights, and place us on the path toward a predictive understanding of enzyme structure and function.

Heme-binding enables allosteric modulation in an ancient TIM-barrel glycosidase

Glycosidases are phylogenetically widely distributed enzymes that are crucial for the cleavage of glycosidic bonds. Here, we present the exceptional properties of a putative ancestor of bacterial and eukaryotic family-1 glycosidases. The ancestral protein shares the TIM-barrel fold with its modern descendants but displays large regions with greatly enhanced conformational flexibility. Yet, the barrel core remains comparatively rigid and the ancestral glycosidase activity is stable, with an optimum temperature within the experimental range for thermophilic family-1 glycosidases. None of the ∼5500 reported crystallographic structures of ∼1400 modern glycosidases show a bound porphyrin. Remarkably, the ancestral glycosidase binds heme tightly and stoichiometrically at a well-defined buried site. Heme binding rigidifies this TIM-barrel and allosterically enhances catalysis. Our work demonstrates the capability of ancestral protein reconstructions to reveal valuable but unexpected biomolecular features when sampling distant sequence space. The potential of the ancestral glycosidase as a scaffold for custom catalysis and biosensor engineering is discussed.

Family 1 glycosidases (GH1) are present in the three domains of life and share classical TIM-barrel fold. Structural and biochemical analyses of a resurrected ancestral GH1 enzyme reveal heme binding, not known in its modern descendants. Heme rigidifies the TIM-barrel and allosterically enhances catalysis.

Structure-based reconstruction of a Mycobacterium hypothetical protein into an active Δ5-3-ketosteroid isomerase

Protein engineering based on structure homology holds the potential to engineer steroid-transforming enzymes on demand. Based on the genome sequencing analysis of industrial Mycobacterium strain HGMS2 to produce 4-androstene-3,17-dione (4-AD), three hypothetical proteins were predicted as putative Δ⁵-3-ketosteroid isomerases (KSIs) to catalyze an intramolecular proton transfer involving the transformation of 5-androstene-3,17-dione (5-AD) into 4-AD, which were defined as mKSI228, mKSI291 and mKSI753. Activity assays indicated that mKSI228 and mKSI291 exhibited weak activity, as low as 0.7% and 1.5%, respectively, of a well-studied and highly active KSI from Pseudomonas putida KSI (pKSI), while mKSI753 had no activity similar to Mycobacterium tuberculosis KSI (mtKSI). Although the 3D structures of the putative mKSIs were homologous to pKSI, their amino acid sequences were significantly different from those of pKSI and tKSI. Thus, by use of these two KSIs as homology models, we were able to convert the low-active mKSI291 into a high-active active KSI by site-directed mutagenesis. On the other hand, an X-ray crystallographic structure of mKSI291 identified a water molecule in its active site. This unique water molecule might function as a bridge to connect Ser-OH, Tyr57-OH and C3O of the intermediate form a hydrogen-bonding network that was responsible for its weak activity, compared with that of mtKSI. Our results not only demonstrated the use of a protein engineering approach to understanding KSI catalytic mechanism, but also provided an example for engineering the catalytic active sites and gaining a functional enzyme based on homologous structures.

Self-assembled organic nanotube promoted allylation of ketones in aqueous phase

A self-assembled organic nanotube was found to promote the allylation of ketones in the aqueous phase.

The Recognition of Unrelated Ligands by Identical Proteins

Unrelated ligands, often found in drug discovery campaigns, can bind to the same receptor, even with the same protein residues. To investigate how this might occur, and whether it might be typically possible to find unrelated ligands for the same drug target, we sought examples of topologically unrelated ligands that bound to the same protein in the same site. Seventy-six pairs of ligands, each bound to the same protein (152 complexes total), were considered, classified into three groups. In the first (31 pairs of complexes), unrelated ligands interacted largely with the same pocket residues through different functional groups. In the second group (39 pairs), the unrelated ligand in each pair engaged different residues, though still within the same pocket. The smallest group (6 pairs) contained ligands with different scaffolds but with shared functional groups interacting with the same residues. We found that there are multiple chemically unrelated but physically similar functional groups that can complement any given local protein pocket; when these functional group substitutions are combined within a single molecule, they lead to topologically unrelated ligands that can each well-complement a site. It may be that many active and orthosteric sites can recognize topologically unrelated ligands.

Efficiently Photocontrollable or Not? Biological Activity of Photoisomerizable Diarylethenes

Diarylethene derivatives, the biological activity of which can be reversibly changed by irradiation with light of different wavelengths, have shown promise as scientific tools and as candidates for photocontrollable drugs. However, examples demonstrating efficient photocontrol of their biological activity are still relatively rare. This concept article discusses the possible reasons for this situation and presents a critical analysis of existing data and hypotheses in this field, in order to extract the design principles enabling the construction of efficient photocontrollable diarylethene-based molecules. Papers addressing biologically relevant interactions between diarylethenes and biomolecules are analyzed; however, in most published cases, the efficiency of photocontrol in living systems remains to be demonstrated. We hope that this article will encourage further discussion of design principles, primarily among pharmacologists, synthetic and medicinal chemists.

Hydrogen Bonds: Simple after All?

Hydrogen bonds play integral roles in biological structure, function, and conformational dynamics and are fundamental to life as it has evolved on Earth. However, our understanding of these fundamental and ubiquitous interactions has seemed fractured and incomplete, and it has been difficult to extract generalities and principles about hydrogen bonds despite thousands of papers published on this topic, perhaps in part because of the expanse of this subject and the density of studies. Fortunately, recent hydrogen bond proposals, discussions, and debates have stimulated new tests and models and have led to a remarkably simple picture of the structure of hydrogen bonds. This knowledge also provides clarity concerning hydrogen bond energetics, limiting and simplifying the factors that need be considered. Herein we recount the advances that have led to this simpler view of hydrogen bond structure, dynamics, and energetics. A quantitative predictive model for hydrogen bond length can now be broadly and deeply applied to evaluate current proposals and to uncover structural features of proteins, their conformational restraints, and their correlated motions. In contrast, a quantitative energetic description of molecular recognition and catalysis by proteins remains an important ongoing challenge, although our improved understanding of hydrogen bonds may aid in testing predictions from current and future models. We close by codifying our current state of understanding into five "Rules for Hydrogen Bonding" that may provide a foundation for understanding and teaching about these vital interactions and for building toward a deeper understanding of hydrogen bond energetics.

Electric Fields and Fast Protein Dynamics in Enzymes

In recent years, there has been much discussion regarding the origin of enzymatic catalysis and whether including protein dynamics is necessary for understanding catalytic enhancement. An important contribution in this debate was made with the application of the vibrational Stark effect spectroscopy to measure electric fields in the active site. This provided a window on electric fields at the transition state in enzymatic reactions. We performed computational studies on two enzymes where we have shown that fast dynamics is part of the reaction mechanism and calculated the electric field near the bond-breaking event. We found that the fast motions that we had identified lead to an increase of the electric field, thus preparing an enzymatic configuration that is electrostatically favorable for the catalytic chemical step. We also studied the enzyme that has been the subject of Stark spectroscopy, ketosteroid isomerase, and found electric fields of a similar magnitude to the two previous examples.

19 F NMR isotropic chemical shift for efficient screening of fluorinated fragments which are racemates and/or display multiple conformers

Fluorine ligand-based NMR spectroscopy is now an established method for performing binding screening against a macromolecular target. Typically, the transverse relaxation rate of the fluorine signals is monitored in the absence and presence of the target. However, useful structural information can sometimes be obtained from the analysis of the fluorine isotropic chemical shift. This is particularly relevant for molecules that are racemates and/or display multiple conformers. The large difference in fluorine isotropic chemical shift between free and bound state deriving mainly from the breaking and/or making of intramolecular and/or intermolecular hydrogen bonds allows the detection of very weak affinity ligands. According to our experimental results, racemates should always be included in the generation of the fluorinated fragment libraries. The selection or the availability of only one of the enantiomers for the fluorinated screening library could result in missing relevant chemical scaffold motifs.

Electric Fields and Enzyme Catalysis

What happens inside an enzyme's active site to allow slow and difficult chemical reactions to occur so rapidly? This question has occupied biochemists' attention for a long time. Computer models of increasing sophistication have predicted an important role for electrostatic interactions in enzymatic reactions, yet this hypothesis has proved vexingly difficult to test experimentally. Recent experiments utilizing the vibrational Stark effect make it possible to measure the electric field a substrate molecule experiences when bound inside its enzyme's active site. These experiments have provided compelling evidence supporting a major electrostatic contribution to enzymatic catalysis. Here, we review these results and develop a simple model for electrostatic catalysis that enables us to incorporate disparate concepts introduced by many investigators to describe how enzymes work into a more unified framework stressing the importance of electric fields at the active site.

De novo endo-functionalized organic cages as cooperative multi-hydrogen-bond-donating catalysts

Two endo-functionalized organic cages as oxyanion hole mimics were achieved via dynamic covalent chemistry, which exhibit good size selectivity, catalytic activity and broad substrate scopes for Friedel-Crafts reactions. The modular character of the synthesis provides an easy way to modify the size, shape and inner function of the molecular cavity.

Kemp Eliminase Activity of Ketosteroid Isomerase

Kemp eliminases represent the most successful class of computationally designed enzymes, with rate accelerations of up to 10⁹-fold relative to the rate of the same reaction in aqueous solution. Nevertheless, several other systems such as micelles, catalytic antibodies, and cavitands are known to accelerate the Kemp elimination by several orders of magnitude. We found that the naturally occurring enzyme ketosteroid isomerase (KSI) also catalyzes the Kemp elimination. Surprisingly, mutations of D38, the residue that acts as a general base for its natural substrate, produced variants that catalyze the Kemp elimination up to 7000-fold better than wild-type KSI does, and some of these variants accelerate the Kemp elimination more than the computationally designed Kemp eliminases. Analysis of the D38N general base KSI variant suggests that a different active site carboxylate residue, D99, performs the proton abstraction. Docking simulations and analysis of inhibition by active site binders suggest that the Kemp elimination takes place in the active site of KSI and that KSI uses the same catalytic strategies of the computationally designed enzymes. In agreement with prior observations, our results strengthen the conclusion that significant rate accelerations of the Kemp elimination can be achieved with very few, nonspecific interactions with the substrate if a suitable catalytic base is present in a hydrophobic environment. Computational design can fulfill these requirements, and the design of more complex and precise environments represents the next level of challenges for protein design.

Handicap-Recover Evolution Leads to a Chemically Versatile, Nucleophile-Permissive Protease

Mutation of the tobacco etch virus (TEV) protease nucleophile from cysteine to serine causes an approximately ∼104 -fold loss in activity. Ten rounds of directed evolution of the mutant, TEVSer , overcame the detrimental effects of nucleophile exchange to recover near-wild-type activity in the mutant TEVSer X. Rather than respecialising TEV to the new nucleophile, all the enzymes along the evolutionary trajectory also retained the ability to use the original cysteine nucleophile. Therefore the adaptive evolution of TEVSer is paralleled by a neutral trajectory for TEVCys , in which mutations that increase serine nucleophile reactivity hardly affect the reactivity of cysteine. This apparent nucleophile permissiveness explains how nucleophile switches can occur in the phylogeny of the chymotrypsin-like protease PA superfamily. Despite the changed key component of their chemical mechanisms, the evolved variants TEVSer X and TEVCys X have similar activities; this could potentially facilitate escape from adaptive conflict to enable active-site evolution.

Differential Distortion of Purine Substrates by Human and Plasmodium falciparum Hypoxanthine-Guanine Phosphoribosyltransferase to Catalyse the Formation of Mononucleotides

Plasmodium falciparum (Pf) hypoxanthine-guanine phosphoribosyltransferase (HGPRT) is a potential therapeutic target. Compared to structurally homologous human enzymes, it has expanded substrate specificity. In this study, 9-deazapurines are used as in situ probes of the active sites of human and Pf HGPRTs. Through the use of these probes it is found that non-covalent interactions stabilise the pre-transition state of the HGPRT-catalysed reaction. Vibrational spectra reveal that the bound substrates are extensively distorted, the carbonyl bond of nucleobase moiety is weakened and the substrate is destabilised along the reaction coordinate. Raman shifts of the human and Pf enzymes are used to quantify the differing degrees of hydrogen bonding in the homologues. A decreased Raman cross-section in enzyme-bound 9-deazaguanine (9DAG) shows that the phenylalanine residue (Phe186 in human and Phe197 in Pf) of HGPRT stacks with the nucleobase. Differential loss of the Raman cross-section suggests that the active site is more compact in human HGPRT as compared to the Pf enzyme, and is more so in the phosphoribosyl pyrophosphate (PRPP) complex 9DAG-PRPP-HGPRT than in 9-deazahypoxanthine (9DAH)-PRPP-HGPRT.

Determination of Hydrogen Bond Structure in Water versus Aprotic Environments To Test the Relationship Between Length and Stability

Hydrogen bonds profoundly influence the architecture and activity of biological macromolecules. Deep appreciation of hydrogen bond contributions to biomolecular function thus requires a detailed understanding of hydrogen bond structure and energetics and the relationship between these properties. Hydrogen bond formation energies (ΔGf) are enormously more favorable in aprotic solvents than in water, and two classes of contributing factors have been proposed to explain this energetic difference, focusing respectively on the isolated and hydrogen-bonded species: (I) water stabilizes the dissociated donor and acceptor groups much better than aprotic solvents, thereby reducing the driving force for hydrogen bond formation; and (II) water lengthens hydrogen bonds compared to aprotic environments, thereby decreasing the potential energy within the hydrogen bond. Each model has been proposed to provide a dominant contribution to ΔGf, but incisive tests that distinguish the importance of these contributions are lacking. Here we directly test the structural basis of model II. Neutron crystallography, NMR spectroscopy, and quantum mechanical calculations demonstrate that O-H···O hydrogen bonds in crystals, chloroform, acetone, and water have nearly identical lengths and very similar potential energy surfaces despite ΔGf differences >8 kcal/mol across these solvents. These results rule out a substantial contribution from solvent-dependent differences in hydrogen bond structure and potential energy after association (model II) and thus support the conclusion that differences in hydrogen bond ΔGf are predominantly determined by solvent interactions with the dissociated groups (model I). These findings advance our understanding of universal hydrogen-bonding interactions and have important implications for biology and engineering.

Modulation of inherent dynamical tendencies of the bisabolyl cation via preorganization in epi-isozizaene synthase

The relative importance of preorganization, selective transition state stabilization and inherent reactivity are assessed through quantum chemical and docking calculations for a sesquiterpene synthase (epi-isozizaene synthase, EIZS). Inherent reactivity of the bisabolyl cation, both static and dynamic, appears to determine the pathway to product, although preorganization and selective binding of the final transition state structure in the multi-step carbocation cascade that forms epi-isozizaene appear to play important roles.

Quantum delocalization of protons in the hydrogen-bond network of an enzyme active site

Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.

Probing the Origins of Catalytic Discrimination between Phosphate and Sulfate Monoester Hydrolysis: Comparative Analysis of Alkaline Phosphatase and Protein Tyrosine Phosphatases

Catalytic promiscuity, the ability of enzymes to catalyze multiple reactions, provides an opportunity to gain a deeper understanding of the origins of catalysis and substrate specificity. Alkaline phosphatase (AP) catalyzes both phosphate and sulfate monoester hydrolysis reactions with a ∼10¹⁰-fold preference for phosphate monoester hydrolysis, despite the similarity between these reactions. The preponderance of formal positive charge in the AP active site, particularly from three divalent metal ions, was proposed to be responsible for this preference by providing stronger electrostatic interactions with the more negatively charged phosphoryl group versus the sulfuryl group. To test whether positively charged metal ions are required to achieve a high preference for the phosphate monoester hydrolysis reaction, the catalytic preference of three protein tyrosine phosphatases (PTPs), which do not contain metal ions, were measured. Their preferences ranged from 5 × 10⁶ to 7 × 10⁷, lower than that for AP but still substantial, indicating that metal ions and a high preponderance of formal positive charge within the active site are not required to achieve a strong catalytic preference for phosphate monoester over sulfate monoester hydrolysis. The observed ionic strength dependences of k_cat/K_M values for phosphate and sulfate monoester hydrolysis are steeper for the more highly charged phosphate ester with both AP and the PTP Stp1, following the dependence expected based on the charge difference of these two substrates. However, the dependences for AP were not greater than those of Stp1 and were rather shallow for both enzymes. These results suggest that overall electrostatics from formal positive charge within the active site is not the major driving force in distinguishing between these reactions and that substantial discrimination can be attained without metal ions. Thus, local properties of the active site, presumably including multiple positioned dipolar hydrogen bond donors within the active site, dominate in defining this reaction specificity.

Fluorine as a hydrogen-bond acceptor: experimental evidence and computational calculations

Hydrogen-bonding interactions play an important role in many chemical and biological systems. Fluorine acting as a hydrogen-bond acceptor in intermolecular and intramolecular interactions has been the subject of many controversial discussions and there are different opinions about it. Recently, we have proposed a correlation between the propensity of fluorine to be involved in hydrogen bonds and its (19)F NMR chemical shift. We now provide additional experimental and computational evidence for this correlation. The strength of hydrogen-bond complexes involving the fluorine moieties CH2F, CHF2, and CF3 was measured and characterized in simple systems by using established and novel NMR methods and compared to the known hydrogen-bond complex formed between acetophenone and p-fluorophenol. Implications of these results for (19)F NMR screening are analyzed in detail. Computed values of the molecular electrostatic potential at the different fluorine atoms and the analysis of the electron density topology at bond critical points correlate well with the NMR results.

The Hydractinia echinata test-system. III: Structure-toxicity relationship study of some azo-, azo-anilide, and diazonium salt derivatives

Structure-toxicity relationships for a series of 75 azo and azo-anilide dyes and five diazonium salts were developed using Hydractinia echinata (H. echinata) as model species. In addition, based on these relationships, predictions for 58 other azo-dyes were made. The experimental results showed that the measured effectiveness Mlog(1/MRC₅₀) does not depend on the number of azo groups or the ones corresponding to metobolites, but it is influenced by the number of anilide groups, as well as by the substituents’ positions within molecules. The conformational analysis pointed out the intramolecular hydrogen bonds, especially the simple tautomerization of quinoidic (ST_OH) or aminoidic (ST_NH2) type. The effectiveness is strongly influenced by the “push-pull” electronic effect, specific to two hydroxy or amino groups separated by an azo moiety (double alternate tautomery, (DAT), to the –COOH or –SO₃H groups which are located in ortho or para position with respect to the azo group. The levels of the lipophylic/hydrophilic, electronic and steric equilibriums, pointed out by the Mlog(1/MRC₅₀) values, enabled the calculation of their average values Clog(1/MRC₅₀) (“Köln model”), characteristic to one derivative class (class isotoxicity). The azo group reduction and the hydrolysis of the amido/peptidic group are two concurrent enzymatic reactions, which occur with different reaction rates and mechanisms. The products of the partial biodegradation are aromatic amines. No additive or synergic effects are noticed among them.

Heterogeneous preferential solvation of water and trifluoroethanol in homologous lysozymes

Cytoplasmic osmolytes can significantly alter the thermodynamic and kinetic properties of proteins relative to those under dilute solution conditions. Spectroscopic experiments of lysozymes in cosolvents indicate that such changes may arise from the heterogeneous, site-specific hydrophobic interactions between protein surface residues and individual solvent molecules. In pursuit of an accurate and predictive model for explaining biomolecular interactions, we study the averaged structural characteristics of mixed solvents with homologous lysozyme solutes using all-atom molecular dynamics. By observing the time-averaged densities of different aqueous solutions of trifluoroethanol, we deduce trends in the heterogeneous solvent interactions over each protein’s surface, and investigate how the homology of protein structure does not necessarily translate to similarities in solvent structure and composition—even when observing identical side chains.

Using unnatural amino acids to probe the energetics of oxyanion hole hydrogen bonds in the ketosteroid isomerase active site

Hydrogen bonds are ubiquitous in enzyme active sites, providing binding interactions and stabilizing charge rearrangements on substrate groups over the course of a reaction. But understanding the origin and magnitude of their catalytic contributions relative to hydrogen bonds made in aqueous solution remains difficult, in part because of complexities encountered in energetic interpretation of traditional site-directed mutagenesis experiments. It has been proposed for ketosteroid isomerase and other enzymes that active site hydrogen bonding groups provide energetic stabilization via "short, strong" or "low-barrier" hydrogen bonds that are formed due to matching of their pKa or proton affinity to that of the transition state. It has also been proposed that the ketosteroid isomerase and other enzyme active sites provide electrostatic environments that result in larger energetic responses (i.e., greater "sensitivity") to ground-state to transition-state charge rearrangement, relative to aqueous solution, thereby providing catalysis relative to the corresponding reaction in water. To test these models, we substituted tyrosine with fluorotyrosines (F-Tyr's) in the ketosteroid isomerase (KSI) oxyanion hole to systematically vary the proton affinity of an active site hydrogen bond donor while minimizing steric or structural effects. We found that a 40-fold increase in intrinsic F-Tyr acidity caused no significant change in activity for reactions with three different substrates. F-Tyr substitution did not change the solvent or primary kinetic isotope effect for proton abstraction, consistent with no change in mechanism arising from these substitutions. The observed shallow dependence of activity on the pKa of the substituted Tyr residues suggests that the KSI oxyanion hole does not provide catalysis by forming an energetically exceptional pKa-matched hydrogen bond. In addition, the shallow dependence provides no indication of an active site electrostatic environment that greatly enhances the energetic response to charge accumulation, consistent with prior experimental results.

Acceleration of an aromatic Claisen rearrangement via a designed spiroligozyme catalyst that mimics the ketosteroid isomerase catalytic dyad

A series of hydrogen-bonding catalysts have been designed for the aromatic Claisen rearrangement of a 1,1-dimethylallyl coumarin. These catalysts were designed as mimics of the two-point hydrogen-bonding interaction present in ketosteroid isomerase that has been proposed to stabilize a developing negative charge on the ether oxygen in the migration of the double bond.1 Two hydrogen bond donating groups, a phenol alcohol and a carboxylic acid, were grafted onto a conformationally restrained spirocyclic scaffold, and together they enhance the rate of the Claisen rearrangement by a factor of 58 over the background reaction. Theoretical calculations correctly predict the most active catalyst and suggest that both preorganization and favorable interactions with the transition state of the reaction are responsible for the observed rate enhancement.

Quantitative dissection of hydrogen bond-mediated proton transfer in the ketosteroid isomerase active site

Hydrogen bond networks are key elements of protein structure and function but have been challenging to study within the complex protein environment. We have carried out in-depth interrogations of the proton transfer equilibrium within a hydrogen bond network formed to bound phenols in the active site of ketosteroid isomerase. We systematically varied the proton affinity of the phenol using differing electron-withdrawing substituents and incorporated site-specific NMR and IR probes to quantitatively map the proton and charge rearrangements within the network that accompany incremental increases in phenol proton affinity. The observed ionization changes were accurately described by a simple equilibrium proton transfer model that strongly suggests the intrinsic proton affinity of one of the Tyr residues in the network, Tyr16, does not remain constant but rather systematically increases due to weakening of the phenol-Tyr16 anion hydrogen bond with increasing phenol proton affinity. Using vibrational Stark spectroscopy, we quantified the electrostatic field changes within the surrounding active site that accompany these rearrangements within the network. We were able to model these changes accurately using continuum electrostatic calculations, suggesting a high degree of conformational restriction within the protein matrix. Our study affords direct insight into the physical and energetic properties of a hydrogen bond network within a protein interior and provides an example of a highly controlled system with minimal conformational rearrangements in which the observed physical changes can be accurately modeled by theoretical calculations.

Fundamental challenges in mechanistic enzymology: progress toward understanding the rate enhancements of enzymes

Enzymes are remarkable catalysts that lie at the heart of biology, accelerating chemical reactions to an astounding extent with extraordinary specificity. Enormous progress in understanding the chemical basis of enzymatic transformations and the basic mechanisms underlying rate enhancements over the past decades is apparent. Nevertheless, it has been difficult to achieve a quantitative understanding of how the underlying mechanisms account for the energetics of catalysis, because of the complexity of enzyme systems and the absence of underlying energetic additivity. We review case studies from our own work that illustrate the power of precisely defined and clearly articulated questions when dealing with such complex and multifaceted systems, and we also use this approach to evaluate our current ability to design enzymes. We close by highlighting a series of questions that help frame some of what remains to be understood, and we encourage the reader to define additional questions and directions that will deepen and broaden our understanding of enzymes and their catalysis.

QM/MM modelling of ketosteroid isomerase reactivity indicates that active site closure is integral to catalysis

Ketosteroid isomerase (Δ⁵-3-keto steroid isomerase or steroid Δ-isomerase) is a highly efficient enzyme at the centre of current debates on enzyme catalysis. We have modelled the reaction mechanism of the isomerization of 3-oxo-Δ⁵-steroids into their Δ⁴-conjugated isomers using high-level combined quantum mechanics/molecular mechanics (QM/MM) methods, and semi-empirical QM/MM molecular dynamics simulations. Energy profiles were obtained at various levels of QM theory (AM1, B3LYP and SCS-MP2). The high-level QM/MM profile is consistent with experimental data. QM/MM dynamics simulations indicate that active site closure and desolvation of the catalytic Asp38 occur before or during formation of dienolate intermediates. These changes have a significant effect on the reaction barrier. A low barrier to reaction is found only when the active site is closed, poising it for catalysis. This conformational change is thus integral to the whole process. The effects on the barrier are apparently largely due to changes in solvation. The combination of high-level QM/MM energy profiles and QM/MM dynamics simulation shows that the reaction involves active site closure, desolvation of the catalytic base, efficient isomerization and re-opening of the active site. These changes highlight the transition between the ligand binding/releasing form and the catalytic form of the enzyme. The results demonstrate that electrostatic interactions (as a consequence of pre-organization of the active site) are crucial for stabilization during the chemical reaction step, but closure of the active site is essential for efficient catalysis to occur.

Flexibility and reactivity in promiscuous enzymes

Best of both worlds: The interplay of active site reactivity and the dynamic character of proteins allows enzymes to be promiscuous and--sometimes--remarkably efficient at the same time. This review analyses the roles structural flexibility and chemical reactivity play in the catalytic mechanism of selected enzymes.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

Psychrophilic enzymes: from folding to function and biotechnology

Psychrophiles thriving permanently at near-zero temperatures synthesize cold-active enzymes to sustain their cell cycle. Genome sequences, proteomic, and transcriptomic studies suggest various adaptive features to maintain adequate translation and proper protein folding under cold conditions. Most psychrophilic enzymes optimize a high activity at low temperature at the expense of substrate affinity, therefore reducing the free energy barrier of the transition state. Furthermore, a weak temperature dependence of activity ensures moderate reduction of the catalytic activity in the cold. In these naturally evolved enzymes, the optimization to low temperature activity is reached via destabilization of the structures bearing the active site or by destabilization of the whole molecule. This involves a reduction in the number and strength of all types of weak interactions or the disappearance of stability factors, resulting in improved dynamics of active site residues in the cold. These enzymes are already used in many biotechnological applications requiring high activity at mild temperatures or fast heat-inactivation rate. Several open questions in the field are also highlighted.

The swipe card model of odorant recognition

Just how we discriminate between the different odours we encounter is not completely understood yet. While obviously a matter involving biology, the core issue isa matter for physics: what microscopic interactions enable the receptors in our noses-small protein switches—to distinguish scent molecules? We survey what is and is not known about the physical processes that take place when we smell things, highlighting the difficulties in developing a full understanding of the mechanics of odorant recognition. The main current theories, discussed here, fall into two major groups. One class emphasises the scent molecule's shape, and is described informally as a "lock and key" mechanism. But there is another category, which we focus on and which we call "swipe card" theories:the molecular shape must be good enough, but the information that identifies the smell involves other factors. One clearly-defined "swipe card" mechanism that we discuss here is Turin's theory, in which inelastic electron tunnelling is used to discern olfactant vibration frequencies. This theory is explicitly quantal, since it requires the molecular vibrations to take in or give out energy only in discrete quanta. These ideas lead to obvious experimental tests and challenges. We describe the current theory in a form that takes into account molecular shape as well as olfactant vibrations. It emerges that this theory can explain many observations hard to reconcile in other ways. There are still some important gaps in a comprehensive physics-based description of the central steps in odorant recognition. We also discuss how far these ideas carry over to analogous processes involving other small biomolecules, like hormones, steroids and neurotransmitters. We conclude with a discussion of possible quantum behaviours in biology more generally, the case of olfaction being just one example. This paper is presented in honour of Prof. Marshall Stoneham who passed away unexpectedly during its writing.

Three hydrogen bond donor catalysts: oxyanion hole mimics and transition state analogues

Enzymes and their mimics use hydrogen bonds to catalyze chemical transformations. Small-molecule transition state analogues of oxyanion holes have been characterized by computations, gas-phase IR and photoelectron spectroscopy, and determination of their binding constants in acetonitrile. A new class of hydrogen bond catalysts is proposed (donors that can contribute three hydrogen bonds to a single functional group) and demonstrated in a Friedel-Crafts reaction. The employed catalyst was observed to react 100 times faster than its rotamer that can employ only two hydrogen bonds. The former compound also binds anions more tightly and was found to have a thermodynamic advantage.

Accurate simulation and detection of coevolution signals in multiple sequence alignments

BACKGROUND

While the conserved positions of a multiple sequence alignment (MSA) are clearly of interest, non-conserved positions can also be important because, for example, destabilizing effects at one position can be compensated by stabilizing effects at another position. Different methods have been developed to recognize the evolutionary relationship between amino acid sites, and to disentangle functional/structural dependencies from historical/phylogenetic ones.

METHODOLOGY/PRINCIPAL FINDINGS

We have used two complementary approaches to test the efficacy of these methods. In the first approach, we have used a new program, MSAvolve, for the in silico evolution of MSAs, which records a detailed history of all covarying positions, and builds a global coevolution matrix as the accumulated sum of individual matrices for the positions forced to co-vary, the recombinant coevolution, and the stochastic coevolution. We have simulated over 1600 MSAs for 8 protein families, which reflect sequences of different sizes and proteins with widely different functions. The calculated coevolution matrices were compared with the coevolution matrices obtained for the same evolved MSAs with different coevolution detection methods. In a second approach we have evaluated the capacity of the different methods to predict close contacts in the representative X-ray structures of an additional 150 protein families using only experimental MSAs.

CONCLUSIONS/SIGNIFICANCE

Methods based on the identification of global correlations between pairs were found to be generally superior to methods based only on local correlations in their capacity to identify coevolving residues using either simulated or experimental MSAs. However, the significant variability in the performance of different methods with different proteins suggests that the simulation of MSAs that replicate the statistical properties of the experimental MSA can be a valuable tool to identify the coevolution detection method that is most effective in each case.

Computational protein engineering: bridging the gap between rational design and laboratory evolution

Enzymes are tremendously proficient catalysts, which can be used as extracellular catalysts for a whole host of processes, from chemical synthesis to the generation of novel biofuels. For them to be more amenable to the needs of biotechnology, however, it is often necessary to be able to manipulate their physico-chemical properties in an efficient and streamlined manner, and, ideally, to be able to train them to catalyze completely new reactions. Recent years have seen an explosion of interest in different approaches to achieve this, both in the laboratory, and in silico. There remains, however, a gap between current approaches to computational enzyme design, which have primarily focused on the early stages of the design process, and laboratory evolution, which is an extremely powerful tool for enzyme redesign, but will always be limited by the vastness of sequence space combined with the low frequency for desirable mutations. This review discusses different approaches towards computational enzyme design and demonstrates how combining newly developed screening approaches that can rapidly predict potential mutation “hotspots” with approaches that can quantitatively and reliably dissect the catalytic step can bridge the gap that currently exists between computational enzyme design and laboratory evolution studies.