Use of anion-aromatic interactions to position the general base in the ketosteroid isomerase active site

Design, synthesis, docking studies and bioactivity evaluation of 1,2,3-triazole eugenol derivatives

Aim: The design, synthesis, docking studies and evaluation of the in vitro antifungal and cytotoxic properties of eugenol (EUG) containing 1,2,3-triazole derivatives are reported. Most of the derivatives have not been reported.Materials & methods: The EUG derivatives were synthesized, molecular docked and tested for their antifungal activity.Results: The compounds showed potent antifungal activity against Trichophyton rubrum, associated with dermatophytosis. Compounds 2a and 2i exhibited promising results, with 2a being four-times more potent than EUG. The binding mode prediction was similar to itraconazole in the lanosterol-14-α-demethylase wild-type and G73E mutant binding sites. Additionally, the pharmacokinetic profile prediction suggests good gastrointestinal absorption and potential oral administration.Conclusion: Compound 2a is a promising antifungal agent against dermatophytosis caused by T. rubrum.

Tailoring Carbon Tails of Ligands on Au52(SR)32 Nanoclusters Enhances the Near-Infrared Photoluminescence Quantum Yield from 3.8 to 18.3

One of the important factors that determine the photoluminescence (PL) properties of gold nanoclusters pertain to the surface. In this study, four Au52(SR)32 nanoclusters that feature a series of aromatic thiolate ligands (-SR) with different bulkiness at the para-position are synthesized and investigated. The near-infrared (NIR) photoluminescence (peaks at 900-940 nm) quantum yield (QY) is largely enhanced with a decrease in the ligand's para-bulkiness. Specifically, the Au52(SR)32 capped with the least bulky p-methylbenzenethiolate (p-MBT) exhibits the highest PLQY (18.3% at room temperature in non-degassed dichloromethane), while Au52 with the bulkiest tert-butylbenzenethiolate (TBBT) only gives 3.8%. The large enhancement of QY with fewer methyl groups on the ligands implies a nonradiative decay via the multiphonon process mediated by C-H bonds. Furthermore, single-crystal X-ray diffraction (SCXRD) comparison of Au52(p-MBT)32 and Au52(TBBT)32 reveals that fewer methyl groups at the para-position lead to a stronger interligand π···π stacking on the Au52 core, thus restricting ligand vibrations and rotations. The emission nature is identified to be phosphorescence and thermally activated delayed fluorescence (TADF) based on the PL lifetime, 3O2 quenching, and temperature-dependent PL and absorption studies. The 1O2 generation efficiencies for the four Au52(SR)32 NCs follow the same trend as the observed PL performance. Overall, the highly NIR-luminescent Au52(p-MBT)32 nanocluster and the revealed mechanisms are expected to find future applications.

Molecular basis for selective uptake and elimination of organic anions in the kidney by OAT1

In mammals, the kidney plays an essential role in maintaining blood homeostasis through the selective uptake, retention or elimination of toxins, drugs and metabolites. Organic anion transporters (OATs) are responsible for the recognition of metabolites and toxins in the nephron and their eventual urinary excretion. Inhibition of OATs is used therapeutically to improve drug efficacy and reduce nephrotoxicity. The founding member of the renal organic anion transporter family, OAT1 (also known as SLC22A6), uses the export of α-ketoglutarate (α-KG), a key intermediate in the Krebs cycle, to drive selective transport and is allosterically regulated by intracellular chloride. However, the mechanisms linking metabolite cycling, drug transport and intracellular chloride remain obscure. Here, we present cryogenic-electron microscopy structures of OAT1 bound to α-KG, the antiviral tenofovir and clinical inhibitor probenecid, used in the treatment of Gout. Complementary in vivo cellular assays explain the molecular basis for α-KG driven drug elimination and the allosteric regulation of organic anion transport in the kidney by chloride.

Expanding the β-substitution reactions of serine synthase through mutagenesis of aromatic active site residues

The Gram-negative bacterium, Fusobacterium nucleatum, possesses a fold II type pyridoxal 5'-phosphate-dependent enzyme that catalyzes the reversible β-replacement of l-cysteine and l-serine, generating H2S and H2O, respectively. This enzyme, termed serine synthase (FN1055), contains an active site Asp232 that serves as a general base in the activation of a water molecule for nucleophilic attack of the ⍺-aminoacrylate intermediate. A network of hydrophobic residues surrounding Asp232 are key to catalysis as they increase the basicity of the side chain. However, these residues severely restrict the range of nucleophilic substrates that can react with the ⍺-aminoacrylate, making the enzyme an ineffective biocatalyst for noncanonical amino acid biosynthesis. Herein, we systematically substituted four aromatic active residues (Trp99, Phe125, Phe148 and Phe234) to an alanine to determine their catalytic importance in serine/cysteine synthase reactions and if their substitution could broaden the scope of nucleophiles that could react with the ⍺-aminoacrylate intermediate. All four single site mutants W99A, F125A, F148A, and F234A could form the ⍺-aminoacrylate intermediate upon reaction with either l-cysteine or l-serine; however, the rate constant associated with the elimination of the β-hydroxyl group from l-serine was 150 to 200-fold lower in the F125A and F148A variants. Substitution of Phe125 and Phe148, situated ∼3-4 Å from the general base, also abolished the serine synthase reaction due to their inability to activate a water molecule for nucleophilic attack of the ⍺-aminoacrylate. Overall, the mutational studies indicate that the clustering of aromatic residues disproportionately benefits the serine synthase reaction as they increase the binding affinity for l-cysteine, decrease the binding of the product, l-serine, and promote the activation of a water molecule. Notably, the aminoacrylate species present in F125A and F148A was able to react with thiophenol, signifying that serine synthase has biocatalytic potential in the synthesis of noncanonical amino acids.

Synthesis and evaluation of Fmoc-amino esters and amides bearing a substrate like quaternary ammonium group as selective butyrylcholinesterase inhibitors

The depletion of the neurotransmitter acetylcholine has been suggested to contribute to the reduced cognitive function observed in individuals suffering from neurodegenerative diseases such as Alzheimer's Disease (AD). For the two major cholinesterases, butyrylcholinesterase (BChE) and acetylcholinesterase (AChE), increased BChE activity observed in individuals with AD has been suggested to deplete acetylcholine levels. To reduce acetylcholine degradation and help restore the pool of the neurotransmitter, specific and potent BChE inhibitors are sought. Our previous findings have identified 9-fluorenylmethoxycarbonyl (Fmoc) amino acid-based inhibitors as effective BChE inhibitors. The amino acid-based compounds offered the opportunity to survey a range of structural features to enhance interactions with the enzyme active site. As enzymes interact with features of their substrates, incorporation of substrate-like features was predicted to lead to better inhibitors. Specifically, incorporation of a trimethylammonium moiety to mimic the cationic group of acetylcholine may lead to increased potency and selectivity. To test this model, a series of inhibitors bearing a cationic trimethylammonium group were synthesized, purified, and characterized. While the Fmoc-ester derivatives inhibited the enzyme, additional experiments showed the compounds acted as substrates and were enzymatically hydrolyzed. Inhibition studies with the Fmoc-amide derivatives showed that the compounds do not act as substrates and selectively inhibit BChE with IC50 values in the 0.06-10.0 µM range. Computational docking studies suggest that the inhibitors can interact with cholinyl binding site and peripheral site. Overall, the results suggest that introducing substrate-like characteristics within the Fmoc-amino acid-based background increases their potency. The versatile and ready access to amino acid-based compounds offers an attractive system to further our understanding of the relative importance of protein-small molecule interactions while guiding the development of better inhibitors.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Structural principles of B cell antigen receptor assembly

The B cell antigen receptor (BCR) is composed of a membrane-bound class M, D, G, E or A immunoglobulin for antigen recognition1-3 and a disulfide-linked Igα (also known as CD79A) and Igβ (also known as CD79B) heterodimer (Igα/β) that functions as the signalling entity through intracellular immunoreceptor tyrosine-based activation motifs (ITAMs)4,5. The organizing principle of the BCR remains unknown. Here we report cryo-electron microscopy structures of mouse full-length IgM BCR and its Fab-deleted form. At the ectodomain (ECD), the Igα/β heterodimer mainly uses Igα to associate with Cµ3 and Cµ4 domains of one heavy chain (µHC) while leaving the other heavy chain (µHC') unbound. The transmembrane domain (TMD) helices of µHC and µHC' interact with those of the Igα/β heterodimer to form a tight four-helix bundle. The asymmetry at the TMD prevents the recruitment of two Igα/β heterodimers. Notably, the connecting peptide between the ECD and TMD of µHC intervenes in between those of Igα and Igβ to guide TMD assembly through charge complementarity. Weaker but distinct density for the Igβ ITAM nestles next to the TMD, suggesting potential autoinhibition of ITAM phosphorylation. Interfacial analyses suggest that all BCR classes utilize a general organizational architecture. Our studies provide a structural platform for understanding B cell signalling and designing rational therapies against BCR-mediated diseases.

Anion-π Catalysis Enabled by the Mechanical Bond

We report a series of rotaxane-based anion-π catalysts in which the mechanical bond between a bipyridine macrocycle and an axle containing an NDI unit is intrinsic to the activity observed, including a [3]rotaxane that catalyses an otherwise disfavoured Michael addition in >60 fold selectivity over a competing decarboxylation pathway that dominates under Brønsted base conditions. The results are rationalized by detailed experimental investigations, electrochemical and computational analysis.

Mining anion-aromatic interactions in the Protein Data Bank

Mutual positioning and non-covalent interactions in anion–aromatic motifs are crucial for functional performance of biological systems. In this context, regular, comprehensive Protein Data Bank (PDB) screening that involves various scientific points of view and individual critical analysis is of utmost importance. Analysis of anions in spheres with radii of 5 Å around all 5- and 6-membered aromatic rings allowed us to distinguish 555 259 unique anion–aromatic motifs, including 92 660 structures out of the 171 588 structural files in the PDB. The use of a scarcely exploited (x, h) coordinate system led to (i) identification of three separate areas of motif accumulation: A – over the ring, B – over the ring-substituent bonds, and C – roughly in the plane of the aromatic ring, and (ii) unprecedented simultaneous comparative description of various anion–aromatic motifs located in these areas. Of the various residues considered, i.e. aminoacids, nucleotides, and ligands, the latter two exhibited a considerable tendency to locate in region Avia archetypal anion–π contacts. The applied model not only enabled statistical quantitative analysis of space around the ring, but also enabled discussion of local intermolecular arrangements, as well as detailed sequence and secondary structure analysis, e.g. anion–π interactions in the GNRA tetraloop in RNA and protein helical structures. As a purely practical issue of this work, the new code source for the PDB research was produced, tested and made freely available at https://github.com/chemiczny/PDB_supramolecular_search.

The comprehensive analysis of non-redundant PDB macromolecular structures investigating anion distributions around all aromatic molecules in available biosystems is presented.

E1784K, the most common Brugada syndrome and long-QT syndrome type 3 mutant, disrupts sodium channel inactivation through two separate mechanisms

Inheritable and de novo variants in the cardiac voltage-gated sodium channel, Nav1.5, are responsible for both long-QT syndrome type 3 (LQT3) and Brugada syndrome type 1 (BrS1). Interestingly, a subset of Nav1.5 variants can cause both LQT3 and BrS1. Many of these variants are found in channel structures that form the channel fast inactivation machinery, altering the rate, voltage dependence, and completeness of the fast inactivation process. We used a series of mutants at position 1784 to show that the most common inheritable Nav1.5 variant, E1784K, alters fast inactivation through two separable mechanisms: (1) a charge-dependent interaction that increases the noninactivating current characteristic of E1784K; and (2) a hyperpolarized voltage dependence and accelerated rate of fast inactivation that decreases the peak sodium current. Using a homology model built on the NavPaS structure, we find that the charge-dependent interaction is between E1784 and K1493 in the DIII-DIV linker of the channel, five residues downstream of the putative inactivation gate. This interaction can be disrupted by a positive charge at position 1784 and rescued with the K1493E/E1784K double mutant that abolishes the noninactivating current. However, the double mutant does not restore either the voltage dependence or rates of fast inactivation. Conversely, a mutant at the bottom of DIVS4, K1641D, causes a hyperpolarizing shift in the voltage dependence of fast inactivation and accelerates the rate of fast inactivation without causing an increase in noninactivating current. These findings provide novel mechanistic insights into how the most common inheritable arrhythmogenic mixed syndrome variant, E1784K, simultaneously decreases transient sodium currents and increases noninactivating currents, leading to both BrS1 and LQT3.

Impact of Ligands on Structural and Optical Properties of Ag29 Nanoclusters

A ligand exchange strategy has been employed to understand the role of ligands on the structural and optical properties of atomically precise 29 atom silver nanoclusters (NCs). By ligand optimization, ∼44-fold quantum yield (QY) enhancement of Ag₂₉(BDT)_12-x(DHLA)_x NCs (x = 1-6) was achieved, where BDT and DHLA refer to 1,3-benzene-dithiol and dihydrolipoic acid, respectively. High-resolution mass spectrometry was used to monitor ligand exchange, and structures of the different NCs were obtained through density functional theory (DFT). The DFT results from Ag₂₉(BDT)₁₁(DHLA) NCs were further experimentally verified through collisional cross-section (CCS) analysis using ion mobility mass spectrometry (IM MS). An excellent match in predicted CCS values and optical properties with the respective experimental data led to a likely structure of Ag₂₉(DHLA)₁₂ NCs consisting of an icosahedral core with an Ag₁₆S₂₄ shell. Combining the experimental observation with DFT structural analysis of a series of atomically precise NCs, Ag_29-yAu_y(BDT)_12-x(DHLA)_x (where y, x = 0,0; 0,1; 0,12 and 1,12; respectively), it was found that while the metal core is responsible for the origin of photoluminescence (PL), ligands play vital roles in determining their resultant PLQY.

About the relevance of anion-π interactions in water

Anion-π interactions are emerging as exotic features with potential applications in chemistry. In the last years, their relevance in living systems has been outlined, and so far there is no concluding significant evidence recognized about the participation of anion-π interactions in water because anion-π sensors contain large aromatic hydrophobic surfaces with limited solubility. By transforming a neutral heterocycle (for example quinoline) into its corresponding salt (quinolinium), we have been able to overcome these solubility issues, and new cationic water-soluble fluorophores have been prepared. Herein, we used N-alkylated heterocycles as π-acidic surfaces to shed light on the nature of anion-π in water by the direct measurement of the fluorescence and UV/Vis spectra in combination with DFT and X-ray analyses.

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray-derived conformational ensembles

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme's reaction cycle. To examine positioning and motions and test catalytic proposals, we collected "room temperature" X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Ensemble analyses indicated limited change through KSI's reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI's general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure-function relationships, providing the basis for a new era of "ensemble-function" interrogation of enzymes.

Kemp Elimination Reaction Catalyzed by Electric Fields

The Kemp elimination reaction is the most widely used in the de novo design of new enzymes. The effect of two different kinds of electric fields in the reactions of acetate as a base with benzisoxazole and 5-nitrobenzisoxazole as substrates have been theoretically studied. The effect of the solvent reaction field has been calculated using the SMD continuum model for several solvents; we have shown that solvents inhibit both reactions, the decrease of the reaction rate being larger as far as the dielectric constant is increased. The diminution of the reaction rate is especially remarkable between aprotic organic solvents and protic solvents as water, the electrostatic term of the hydrogen bonds being the main factor for the large inhibitory effect of water. The presence of an external electric field oriented in the direction of the charge transfer (z axis) increases it and, so, the reaction rate. In the reaction of the nitro compound, if the electric field is oriented in an orthogonal direction (x axis) the charge transfer to the NO₂ group is favored and there is a subsequent increase of the reaction rate. However, this increase is smaller than the one produced by the field in the z axis. It is worthwhile mentioning that one of the main effects of external electric fields of intermediate intensity is the reorientation of the reactants. Finally, the implications of our results in the de novo design of enzymes are discussed.

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.

Discovery of novel anaplastic lymphoma kinase inhibitors: Structure and energy-based pharmacophore strategy

The clinical outcomes in patients with non-small cell lung cancer have improved, as a result of anaplastic lymphoma kinase (ALK) inhibition. Therefore in the current study, substantial effort has been made to identify ALK inhibitors through systematic virtual screening experiment consisting of e-pharmacophore and pharmacophore perception techniques. Initially, a pharmacophore model (AAAHP.193) and an e-pharmacophore model (DDRRR) encompassing the whole dataset of 12 known ALK inhibitors were developed. The hypotheses could retrieve effective compounds from DrugBank database (8621 molecules), which were then subjected to molecular docking and ADME prediction. These approaches resulted in the identification of five hits, namely, nebivolol, HDY, D42, 796, and LZE having higher Glide docking scores and promising ADME properties with augmented CNS involvement. Moreover, molecular dynamics simulations were performed to validate the inhibitory activity of the hit compounds, and density functional theory calculations were carried out to scrutinize the chemical reactivity of the hits. Subsequent interaction and scaffold analysis identified prominent interactions of the hits with ALK kinase domain and scaffolds with anti-tumor activity against lung cancer cell lines. We strongly believe that the study provides an outlook for the sighting of novel and potent ALK inhibitors in the near future.

The Emergence of Anion-π Catalysis

The objective of this Account is to summarize the first five years of anion-π catalysis. The general idea of anion-π catalysis is to stabilize anionic transition states on aromatic surfaces. This is complementary to the stabilization of cationic transition states on aromatic surfaces, a mode of action that occurs in nature and is increasingly used in chemistry. Anion-π catalysis, however, rarely occurs in nature and has been unexplored in chemistry. Probably because the attraction of anions to π surfaces as such is counterintuitive, anion-π interactions in general are much younger than cation-π interactions and remain under-recognized until today. Anion-π catalysis has emerged from early findings that anion-π interactions can mediate the transport of anions across lipid bilayer membranes. With this evidence for stabilization in the ground state secured, there was no reason to believe that anion-π interactions could not also stabilize anionic transition states. As an attractive reaction to develop anion-π catalysis, the addition of malonic acid half thioesters to enolate acceptors was selected. This choice was also made because without enzymes decarboxylation is preferred and anion-π interactions promised to catalyze selectively the disfavored but relevant enolate addition. Concerning anion-π catalysts, we started with naphthalene diimides (NDIs) because their intrinsic quadrupole moment is highly positive. The NDI scaffold was used to address questions such as the positioning of substrates on the catalytic π surface or the dependence of activity on the π acidity of this π surface. With the basics in place, the next milestone was the creation of anion-π enzymes, that is, enzymes that operate with an interaction rarely used in biology, at least on intrinsically π-acidic or highly polarizable aromatic amino-acid side chains. Electric-field-assisted anion-π catalysis addresses topics such as heterogeneous catalysis on electrodes and remote control of activity by voltage. On π-stacked foldamers, anion-(π) _n-π catalysis was discovered. Fullerenes emerged as the scaffold of choice to explore contributions from polarizability. On fullerenes, anionic transition states are stabilized by large macrodipoles that appear only in response to their presence. With this growing collection of anion-π catalysts, several reactions beyond enolate addition have been explored so far. Initial efforts focused on asymmetric anion-π catalysis. Increasing enantioselectivity with increasing π acidity of the active π surface has been exemplified for enamine and iminium chemistry and for anion-π transaminase mimics. However, the delocalized nature of anion-π interactions calls for the stabilization of charge displacements over longer distances. The first step in this direction was the formation of cyclohexane rings with five stereogenic centers from achiral acyclic substrates on π-acidic surfaces. Moreover, the intrinsically disfavored exo transition state of anionic Diels-Alder reactions is stabilized selectively on π-acidic surfaces; endo products and otherwise preferred Michael addition products are completely suppressed. Taken together, we hope that these results on catalyst design and reaction scope will establish anion-π catalysis as a general principle in catalysis in the broadest sense.

Rapid assessment of conformational preferences in biaryl and aryl carbonyl fragments

The ability to rapidly assess the preferred conformation of key fragments in a structure “by visual inspection” is a very useful starting point in the process of drug design. With the ability to do so, one could address questions like: “How could we avoid planarity in a molecule?”, “Will a molecule change its conformational preference if we make it more or less basic?” or “How does this electronic repulsion affect the conformational preference in the system?” in timely fashion. In this paper, we describe how the conformational energy profile (CEP, plot of energy as a function of dihedral bond angle) of a fragment can be interpreted through the understanding the interplay between resonance stabilization, steric effects and electrostatic interactions. Fifty-nine biaryl and aryl carbonyl fragments present in oral drugs or which are close derivatives thereof were selected. Calculation of their CEPs using ab initio methodology allowed us to conclude the relative importance of these factors in the conformational preference of these fragments as follows: “steric repulsion > lone pair—lone pair repulsion > lone pair—fluorine repulsion > resonance stabilization” and to formulate “rules of thumb” that the practicing medicinal/organic chemist can apply when analysing molecules that contain these fragments.

Impact of Homologous Resistance Mutations from Pathogenic Yeast on Saccharomyces cerevisiae Lanosterol 14α-Demethylase

Fungal infections frequently affect immunodeficient individuals and are estimated to kill 1.35 million people per annum. Azole antifungals target the membrane-bound cytochrome P450 monooxygenase lanosterol 14α-demethylase (CYP51; Erg11p). Mutations in CYP51 can render the widely used triazole drugs less effective. The Candida albicans CYP51 mutation G464S and the double mutation Y132F G464S (Y140F and G464S by Saccharomyces cerevisiae numbering) as well as the CYP51A G54E/R/W mutations of Aspergillus fumigatus (G73E/R/W by S. cerevisiae numbering) have been reproduced in a recombinant C-terminal hexahistidine-tagged version of S. cerevisiae CYP51 (ScErg11p6×His). Phenotypes and X-ray crystal structures were determined for the mutant enzymes. Liquid microdilution assays showed that the G464S mutation in ScErg11p6×His conferred no difference in the susceptibility of yeast to triazole drugs but in combination with the Y140F mutation gave a 4-fold reduction in susceptibility to the short-tailed triazole fluconazole. The ScErg11p6×His Y140F G464S mutant was unstable during purification and was not crystallized. The ScErg11p6×His G73E/R/W mutations conferred increased susceptibly to all triazoles tested in liquid microdilution assays. High-resolution X-ray crystal structures reveal two different conformations of the ligand itraconazole, including a previously unseen conformation, as well as interactions between the tail of this triazole and the E/W73 residue. This study shows that S. cerevisiae CYP51 adequately represents some but not all mutations in CYP51s of pathogenic fungi. Insight into the molecular mechanisms of resistance mutations in CYP51 will assist the development of inhibitors that will overcome antifungal resistance.

The structures and properties of anionic tryptophan complexes

IRMPD spectroscopy and electronic structure calculations are employed to identify π–π interactions in ionic tryptophan clusters.

Characterization of Ionizable Groups' Environments in Proteins and Protein-Ligand Complexes through a Statistical Analysis of the Protein Data Bank

We conduct a statistical analysis of the molecular environment of common ionizable functional groups in both protein-ligand complexes and inside proteins from the Protein Data Bank (PDB). In particular, we characterize the frequency, type, and density of the interacting atoms as well as the presence of a potential counterion. We found that for ligands, most guanidinium groups, half of primary and secondary amines, and one-fourth of imidazole neighbor a carboxylate group. Tertiary amines bind more rarely near carboxylate groups, which may be explained by a crowded neighborhood and hydrophobic character. In comparison to the environment seen by the ligands, inside proteins, an environment enriched in main-chain atoms is found, and the prevalence of direct charge neutralization by carboxylate groups is different. When the ionizable character of water molecules and phenolic or hydroxyl groups is accounted, considering a high-resolution dataset (less than 1.5 Å), charge neutralization could occur for well above 80% of the ligand functional groups considered, but for tertiary amines.

Measuring inter-protein pairwise interaction energies from a single native mass spectrum by double-mutant cycle analysis

The strength and specificity of protein complex formation is crucial for most life processes and is determined by interactions between residues in the binding partners. Double-mutant cycle analysis provides a strategy for studying the energetic coupling between amino acids at the interfaces of such complexes. Here we show that these pairwise interaction energies can be determined from a single high-resolution native mass spectrum by measuring the intensities of the complexes formed by the two wild-type proteins, the complex of each wild-type protein with a mutant protein, and the complex of the two mutant proteins. This native mass spectrometry approach, which obviates the need for error-prone measurements of binding constants, can provide information regarding multiple interactions in a single spectrum much like nuclear Overhauser effects (NOEs) in nuclear magnetic resonance. Importantly, our results show that specific inter-protein contacts in solution are maintained in the gas phase.

Double mutant cycle (DMC) analyses can provide the interaction energies between amino acids at the interface of protein complexes. Here, the authors determine pairwise interaction energies using high-resolution native mass spectroscopy, offering a straightforward route for the DMC methodology.

An Activator-Blocker Pair Provides a Controllable On-Off Switch for a Ketosteroid Isomerase Active Site Mutant

Control of enzyme activity is fundamental to biology and represents a long-term goal in bioengineering and precision therapeutics. While several powerful molecular strategies have been developed, limitations remain in their generalizability and dynamic range. We demonstrate a control mechanism via separate small molecules that turn on the enzyme (activator) and turn off the activation (blocker). We show that a pocket created near the active site base of the enzyme ketosteriod isomerase (KSI) allows efficient and saturable base rescue when the enzyme's natural general base is removed. Binding a small molecule with similar properties but lacking general-base capability in this pocket shuts off rescue. The ability of small molecules to directly participate in and directly block catalysis may afford a broad controllable dynamic range. This approach may be amenable to numerous enzymes and to engineering and screening approaches to identify activators and blockers with strong, specific binding for engineering and therapeutic applications.

Systematic Functional Analysis of Active-Site Residues in l-Threonine Dehydrogenase from Thermoplasma volcanium

Enzymes have been through millions of years of evolution during which their active-site microenvironments are fine-tuned. Active-site residues are commonly conserved within protein families, indicating their importance for substrate recognition and catalysis. In this work, we systematically mutated active-site residues of l-threonine dehydrogenase from Thermoplasma volcanium and characterized the mutants against a panel of substrate analogs. Our results demonstrate that only a subset of these residues plays an essential role in substrate recognition and catalysis and that the native enzyme activity can be further enhanced roughly 4.6-fold by a single point mutation. Kinetic characterization of mutants on substrate analogs shows that l-threonine dehydrogenase possesses promiscuous activities toward other chemically similar compounds not previously observed. Quantum chemical calculations on the hydride-donating ability of these substrates also reveal that this enzyme did not evolve to harness the intrinsic substrate reactivity for enzyme catalysis. Our analysis provides insights into connections between the details of enzyme active-site structure and specific function. These results are directly applicable to rational enzyme design and engineering.

Highly Dynamic Anion-Quadrupole Networks in Proteins

The dynamics of anion-quadrupole (or anion-π) interactions formed between negatively charged (Asp/Glu) and aromatic (Phe) side chains are for the first time computationally characterized in RmlC (Protein Data Bank entry 1EP0 ), a homodimeric epimerase. Empirical force field-based molecular dynamics simulations predict anion-quadrupole pairs and triplets (anion-anion-π and anion-π-π) are formed by the protein during the simulated trajectory, which suggests that the anion-quadrupole interactions may provide a significant contribution to the overall stability of the protein, with an average of -1.6 kcal/mol per pair. Some anion-π interactions are predicted to form during the trajectory, extending the number of anion-quadrupole interactions beyond those predicted from crystal structure analysis. At the same time, some anion-π pairs observed in the crystal structure exhibit marginal stability. Overall, most anion-π interactions alternate between an "on" state, with significantly stabilizing energies, and an "off" state, with marginal or null stabilizing energies. The way proteins possibly compensate for transient loss of anion-quadrupole interactions is characterized in the RmlC aspartate 84-phenylalanine 112 anion-quadrupole pair observed in the crystal structure. A double-mutant cycle analysis of the thermal stability suggests a possible loss of anion-π interactions compensated by variations of hydration of the residues and formation of compensating electrostatic interactions. These results suggest that near-planar anion-quadrupole pairs can exist, sometimes transiently, which may play a role in maintaining the structural stability and function of the protein, in an otherwise very dynamic interplay of a nonbonded interaction network as well as solvent effects.

Evaluation of the Catalytic Contribution from a Positioned General Base in Ketosteroid Isomerase

Proton transfer reactions are ubiquitous in enzymes and utilize active site residues as general acids and bases. Crystal structures and site-directed mutagenesis are routinely used to identify these residues, but assessment of their catalytic contribution remains a major challenge. In principle, effective molarity measurements, in which exogenous acids/bases rescue the reaction in mutants lacking these residues, can estimate these catalytic contributions. However, these exogenous moieties can be restricted in reactivity by steric hindrance or enhanced by binding interactions with nearby residues, thereby resulting in over- or underestimation of the catalytic contribution, respectively. With these challenges in mind, we investigated the catalytic contribution of an aspartate general base in ketosteroid isomerase (KSI) by exogenous rescue. In addition to removing the general base, we systematically mutated nearby residues and probed each mutant with a series of carboxylate bases of similar pKa but varying size. Our results underscore the need for extensive and multifaceted variation to assess and minimize steric and positioning effects and determine effective molarities that estimate catalytic contributions. We obtained consensus effective molarities of ∼5 × 10(4) M for KSI from Comamonas testosteroni (tKSI) and ∼10(3) M for KSI from Pseudomonas putida (pKSI). An X-ray crystal structure of a tKSI general base mutant showed no additional structural rearrangements, and double mutant cycles revealed similar contributions from an oxyanion hole mutation in the wild-type and base-rescued reactions, providing no indication of mutational effects extending beyond the general base site. Thus, the high effective molarities suggest a large catalytic contribution associated with the general base. A significant portion of this effect presumably arises from positioning of the base, but its large magnitude suggests the involvement of additional catalytic mechanisms as well.

Mechanistic signs of double-barreled structure in a fluoride ion channel

The Fluc family of F⁻ ion channels protects prokaryotes and lower eukaryotes from the toxicity of environmental F⁻. In bacteria, these channels are built as dual-topology dimers whereby the two subunits assemble in antiparallel transmembrane orientation. Recent crystal structures suggested that Fluc channels contain two separate ion-conduction pathways, each with two F⁻ binding sites, but no functional correlates of this unusual architecture have been reported. Experiments here fill this gap by examining the consequences of mutating two conserved F⁻-coordinating phenylalanine residues. Substitution of each phenylalanine specifically extinguishes its associated F⁻ binding site in crystal structures and concomitantly inhibits F⁻ permeation. Functional analysis of concatemeric channels, which permit mutagenic manipulation of individual pores, show that each pore can be separately inactivated without blocking F⁻ conduction through its symmetry-related twin. The results strongly support dual-pathway architecture of Fluc channels.

DOI:http://dx.doi.org/10.7554/eLife.18767.001

Role of conserved Met112 residue in the catalytic activity and stability of ketosteroid isomerase

Ketosteroid isomerase (3-oxosteroid Δ(5)-Δ(4)-isomerase, KSI) from Pseudomonas putida catalyzes allylic rearrangement of the 5,6-double bond of Δ(5)-3-ketosteroid to 4,5-position by stereospecific intramolecular transfer of a proton. The active site of KSI is formed by several hydrophobic residues and three catalytic residues (Tyr14, Asp38, and Asp99). In this study, we investigated the role of a hydrophobic Met112 residue near the active site in the catalysis, steroid binding, and stability of KSI. Replacing Met112 with alanine (yields M112A) or leucine (M112L) decreased the kcat by 20- and 4-fold, respectively. Compared with the wild type (WT), M112A and M112L KSIs showed increased KD values for equilenin, an intermediate analogue; these changes suggest that loss of packing at position 112 might lead to unfavorable steroid binding, thereby resulting in decreased catalytic activity. Furthermore, M112A and M112L mutations reduced melting temperature (Tm) by 6.4°C and 2.5°C, respectively. These changes suggest that favorable packing in the core is important for the maintenance of stability in KSI. The M112K mutation decreased kcat by 2000-fold, compared with the WT. In M112K KSI structure, a new salt bridge was formed between Asp38 and Lys112. This bridge could change the electrostatic potential of Asp38, and thereby contribute to the decreased catalytic activity. The M112K mutation also decreased the stability by reducing Tm by 4.1°C. Our data suggest that the Met112 residue may contribute to the catalytic activity and stability of KSI by providing favorable hydrophobic environments and compact packing in the catalytic core.

X-ray Structure Analysis of Indazolium trans-[Tetrachlorobis(1H-indazole)ruthenate(III)] (KP1019) Bound to Human Serum Albumin Reveals Two Ruthenium Binding Sites and Provides Insights into the Drug Binding Mechanism

Ruthenium(III) complexes are promising candidates for anticancer drugs, especially the clinically studied indazolium trans-[tetrachlorobis(1H-indazole)ruthenate(III)] (KP1019) and its analogue sodium trans-[tetrachlorobis(1H-indazole)ruthenate(III)] (NKP-1339). Several studies have emphasized the likely role of human serum proteins in the transportation and accumulation of ruthenium(III) complexes in tumors. Therefore, the interaction between KP1019 and human serum albumin was investigated by means of X-ray crystallography and inductively coupled plasma mass spectrometry (ICP-MS). The structural data unambiguously reveal the binding of two ruthenium atoms to histidine residues 146 and 242, which are both located within well-known hydrophobic binding pockets of albumin. The ruthenium centers are octahedrally coordinated by solvent molecules revealing the dissociation of both indazole ligands from the ruthenium-based drug. However, a binding mechanism is proposed indicating the importance of the indazole ligands for binding site recognition and thus their indispensable role for the binding of KP1019.

PICI: A web server with a multi-parametric algorithm for identifying interaction sites within protein complexes

Identifying interactions between proteins in a complex is essential to analyze their role in various cellular activities and structural stability. Therefore, effective algorithms and computer programs are required for the better understanding of various interactions exist at the protein-protein interface. The protein inter-chain interaction (PICI) server was developed to identify these interaction sites and various interactions that contribute to the specificity and strength of the protein complex by using Protein Interaction Identification Algorithm (PIIA). The multi-parametric approach of this algorithm involves the interface area between the subunits, the atomic coordinate distance, the linearity of bonds, the orientation of aromatic side-chains, and physicochemical properties of the residues. Particular advantages of PICI include its ability to identify various strong and weak interactions, including those which have not been considered before by any other server. Representation of interface data in text tables, 3D interaction visualizer, colored sequence viewer, and a dynamic colored graphical matrix display makes it distinct from similar web servers. Users can analyze interactions within multi-chain proteins by their PDBids or by uploading a structure atomic coordinate file and can download the result in plain text or XML format.

Molecular Mechanism Underlying ATP-Induced Conformational Changes in the Nucleoprotein Filament of Mycobacterium smegmatis RecA

RecA plays a central role in bacterial DNA repair, homologous recombination, and restoration of stalled replication forks by virtue of its active extended nucleoprotein filament. Binding of ATP and its subsequent recognition by the carboxamide group of a highly conserved glutamine (Gln196 in MsRecA) have been implicated in the formation of active RecA nucleoprotein filaments. Although the mechanism of ATP-dependent structural transitions in RecA has been proposed on the basis of low-resolution electron microscopic reconstructions, the precise sequence of events that constitute these transitions is poorly understood. On the basis of biochemical and crystallographic analyses of MsRecA variants carrying mutations in highly conserved Gln196 and Arg198 residues, we propose that the disposition of the interprotomer interface is the structural basis of allosteric activation of RecA. Furthermore, this study accounts for the contributions of several conserved amino acids to ATP hydrolysis and to the transition from collapsed to extended filament forms in Mycobacterium smegmatis RecA (MsRecA). In addition to their role in the inactive compressed state, the study reveals a role for Gln196 and Arg198 along with Phe219 in ATP hydrolysis in the active extended nucleoprotein filament. Finally, our data suggest that the primary, but not secondary, nucleotide binding site in MsRecA isomerizes into the ATP binding site present in the extended nucleoprotein filament.

σ-Hole Bond vs π-Hole Bond: A Comparison Based on Halogen Bond

The σ-hole and π-hole are the regions with positive surface electrostatic potential on the molecule entity; the former specifically refers to the positive region of a molecular entity along extension of the Y-Ge/P/Se/X covalent σ-bond (Y = electron-rich group; Ge/P/Se/X = Groups IV-VII), while the latter refers to the positive region in the direction perpendicular to the σ-framework of the molecular entity. The directional noncovalent interactions between the σ-hole or π-hole and the negative or electron-rich sites are named σ-hole bond or π-hole bond, respectively. The contributions from electrostatic, charge transfer, and other terms or Coulombic interaction to the σ-hole bond and π-hole bond were reviewed first followed by a brief discussion on the interplay between the σ-hole bond and the π-hole bond as well as application of the two types of noncovalent interactions in the field of anion recognition. It is expected that this review could stimulate further development of the σ-hole bond and π-hole bond in theoretical exploration and practical application in the future.

BIOPHYSICS. Response to Comments on "Extreme electric fields power catalysis in the active site of ketosteroid isomerase"

Natarajan et al. and Chen and Savidge comment that comparing the electric field in ketosteroid isomerase's (KSI's) active site to zero overestimates the catalytic effect of KSI's electric field because the reference reaction occurs in water, which itself exerts a sizable electrostatic field. To compensate, Natarajan et al. argue that additional catalytic weight arises from positioning of the general base, whereas Chen and Savidge propose a separate contribution from desolvation of the general base. We note that the former claim is not well supported by published results, and the latter claim is intriguing but lacks experimental basis. We also take the opportunity to clarify some of the more conceptually subtle aspects of electrostatic catalysis.

Understanding the self-assembly of Fmoc-phenylalanine to hydrogel formation

Hydrogels of low molecular weight molecules are important in biomedical applications. Multiple factors are responsible for hydrogel formation, but their role in governing self-assembly to hydrogel formation is poorly understood. Herein, we report the hydrogel formation of fluorenylmethyloxycarbonyl phenylalanine (FmocF) molecule. We used physical and thermal stimuli for solubilizing FmocF above the critical concentration to induce gel formation. The key role of Fmoc, Fmoc and phenylalanine covalent linkage, flexibility of phe side chain, pH, and buffer ions in self-assembly of FmocF to gel formation is described. We found that the collective action of different non-covalent interactions play a role in making FmocF hydrogel. Using powder diffraction and infrared spectroscopy, we also report a new polymorphic form of FmocF after transitioning to hydrogel. In addition, we are proposing a model for drug release from FmocF hydrogel.

Landscape of π-π and sugar-π contacts in DNA-protein interactions

There were 1765 contacts identified between DNA nucleobases or deoxyribose and cyclic (W, H, F, Y) or acyclic (R, E, D) amino acids in 672 X-ray structures of DNA-protein complexes. In this first study to compare π-interactions between the cyclic and acyclic amino acids, visual inspection was used to categorize amino acid interactions as nucleobase π-π (according to biological edge) or deoxyribose sugar-π (according to sugar edge). Overall, 54% of contacts are nucleobase π-π interactions, which involve all amino acids, but are more common for Y, F, and R, and involve all DNA nucleobases with similar frequencies. Among binding arrangements, cyclic amino acids prefer more planar (stacked) π-systems than the acyclic counterparts. Although sugar-π interactions were only previously identified with the cyclic amino acids and were found to be less common (38%) than nucleobase-cyclic amino acid contacts, sugar-π interactions are more common than nucleobase π-π contacts for the acyclic series (61% of contacts). Similar to DNA-protein π-π interactions, sugar-π contacts most frequently involve Y and R, although all amino acids adopt many binding orientations relative to deoxyribose. These DNA-protein π-interactions stabilize biological systems, by up to approximately -40 kJ mol(-1) for neutral nucleobase or sugar-amino acid interactions, but up to approximately -95 kJ mol(-1) for positively or negatively charged contacts. The high frequency and strength, despite variation in structure and composition, of these π-interactions point to an important function in biological systems.

Experimental and computational mutagenesis to investigate the positioning of a general base within an enzyme active site

The positioning of catalytic groups within proteins plays an important role in enzyme catalysis, and here we investigate the positioning of the general base in the enzyme ketosteroid isomerase (KSI). The oxygen atoms of Asp38, the general base in KSI, were previously shown to be involved in anion–aromatic interactions with two neighboring Phe residues. Here we ask whether those interactions are sufficient, within the overall protein architecture, to position Asp38 for catalysis or whether the side chains that pack against Asp38 and/or the residues of the structured loop that is capped by Asp38 are necessary to achieve optimal positioning for catalysis. To test positioning, we mutated each of the aforementioned residues, alone and in combinations, in a background with the native Asp general base and in a D38E mutant background, as Glu at position 38 was previously shown to be mispositioned for general base catalysis. These double-mutant cycles reveal positioning effects as large as 10³-fold, indicating that structural features in addition to the overall protein architecture and the Phe residues neighboring the carboxylate oxygen atoms play roles in positioning. X-ray crystallography and molecular dynamics simulations suggest that the functional effects arise from both restricting dynamic fluctuations and disfavoring potential mispositioned states. Whereas it may have been anticipated that multiple interactions would be necessary for optimal general base positioning, the energetic contributions from positioning and the nonadditive nature of these interactions are not revealed by structural inspection and require functional dissection. Recognizing the extent, type, and energetic interconnectivity of interactions that contribute to positioning catalytic groups has implications for enzyme evolution and may help reveal the nature and extent of interactions required to design enzymes that rival those found in biology.

Enzyme architecture: on the importance of being in a protein cage

Substrate binding occludes water from the active sites of many enzymes. There is a correlation between the burden to enzymatic catalysis of deprotonation of carbon acids and the substrate immobilization at solvent-occluded active sites for ketosteroid isomerase (KSI--small burden, substrate pKa=13), triosephosphate isomerase (TIM, substrate pKa≈18) and diaminopimelate epimerase (DAP epimerase, large burden, substrate pKa≈29) catalyzed reaction. KSI binds substrates at a surface cleft, TIM binds substrate at an exposed 'cage' formed by closure of flexible loops; and, DAP epimerase binds substrate in a tight cage formed by an 'oyster-like' clamping motion of protein domains. Directed evolution of a solvent-occluded active site at a designed protein catalyst of the Kemp elimination reaction is discussed.

Mechanistic insight into the reaction catalysed by bacterial type II dehydroquinases1

The present study identified the residue that deprotonates the essential tyrosine that triggers the catalytic process and provides details of the required motions for the catalytic turnover. A previously unknown key role for the essential arginine and two conserved arginines is also reported.

Uncovering the determinants of a highly perturbed tyrosine pKa in the active site of ketosteroid isomerase

Within the idiosyncratic enzyme active-site environment, side chain and ligand pKa values can be profoundly perturbed relative to their values in aqueous solution. Whereas structural inspection of systems has often attributed perturbed pKa values to dominant contributions from placement near charged groups or within hydrophobic pockets, Tyr57 of a Pseudomonas putida ketosteroid isomerase (KSI) mutant, suggested to have a pKa perturbed by nearly 4 units to 6.3, is situated within a solvent-exposed active site devoid of cationic side chains, metal ions, or cofactors. Extensive comparisons among 45 variants with mutations in and around the KSI active site, along with protein semisynthesis, (13)C NMR spectroscopy, absorbance spectroscopy, and X-ray crystallography, was used to unravel the basis for this perturbed Tyr pKa. The results suggest that the origin of large energetic perturbations are more complex than suggested by visual inspection. For example, the introduction of positively charged residues near Tyr57 raises its pKa rather than lowers it; this effect, and part of the increase in the Tyr pKa from the introduction of nearby anionic groups, arises from accompanying active-site structural rearrangements. Other mutations with large effects also cause structural perturbations or appear to displace a structured water molecule that is part of a stabilizing hydrogen-bond network. Our results lead to a model in which three hydrogen bonds are donated to the stabilized ionized Tyr, with these hydrogen-bond donors, two Tyr side chains, and a water molecule positioned by other side chains and by a water-mediated hydrogen-bond network. These results support the notion that large energetic effects are often the consequence of multiple stabilizing interactions rather than a single dominant interaction. Most generally, this work provides a case study for how extensive and comprehensive comparisons via site-directed mutagenesis in a tight feedback loop with structural analysis can greatly facilitate our understanding of enzyme active-site energetics. The extensive data set provided may also be a valuable resource for those wishing to extensively test computational approaches for determining enzymatic pKa values and energetic effects.