Mapping the Accessible Conformational Landscape of an Insect Carboxylesterase Using Conformational Ensemble Analysis and Kinetic Crystallography

Quantifying and comparing radiation damage in the Protein Data Bank

Radiation damage remains one of the major bottlenecks to accurate structure solution in protein crystallography. It can induce structural and chemical changes in protein crystals, and is hence an important consideration when assessing the quality and biological veracity of crystal structures in repositories like the Protein Data Bank (PDB). However, detection of radiation damage artefacts has traditionally proved very challenging. To address this, here we introduce the Bnet metric. Bnet summarises in a single value the extent of damage suffered by a crystal structure by comparing the B-factor values of damage-prone and non-damage-prone atoms in a similar local environment. After validating that Bnet successfully detects damage in 23 different crystal structures previously characterised as damaged, we calculate Bnet values for 93,978 PDB crystal structures. Our metric highlights a range of damage features, many of which would remain unidentified by the other summary statistics typically calculated for PDB structures.

Exploring ligand dynamics in protein crystal structures with ensemble refinement

Understanding the dynamics of ligands bound to proteins is an important task in medicinal chemistry and drug design. However, the dominant technique for determining protein-ligand structures, X-ray crystallography, does not fully account for dynamics and cannot accurately describe the movements of ligands in protein binding sites. In this article, an alternative method, ensemble refinement, is used on six protein-ligand complexes with the aim of understanding the conformational diversity of ligands in protein crystal structures. The results show that ensemble refinement sometimes indicates that the flexibility of parts of the ligand and some protein side chains is larger than that which can be described by a single conformation and atomic displacement parameters. However, since the electron-density maps are comparable and Rfree values are slightly increased, the original crystal structure is still a better model from a statistical point of view. On the other hand, it is shown that molecular-dynamics simulations and automatic generation of alternative conformations in crystallographic refinement confirm that the flexibility of these groups is larger than is observed in standard refinement. Moreover, the flexible groups in ensemble refinement coincide with groups that give high atomic displacement parameters or non-unity occupancy if optimized in standard refinement. Therefore, the conformational diversity indicated by ensemble refinement seems to be qualitatively correct, indicating that ensemble refinement can be an important complement to standard crystallographic refinement as a tool to discover which parts of crystal structures may show extensive flexibility and therefore are poorly described by a single conformation. However, the diversity of the ensembles is often exaggerated (probably partly owing to the rather poor force field employed) and the ensembles should not be trusted in detail.

An Overview of Antennal Esterases in Lepidoptera

Lepidoptera are used as a model for the study of insect olfactory proteins. Among them, odorant degrading enzymes (ODEs), that degrade odorant molecules to maintain the sensitivity of antennae, have received less attention. In particular, antennal esterases (AEs; responsible for ester degradation) are crucial for intraspecific communication in Lepidoptera. Currently, transcriptomic and genomic studies have provided AEs in several species. However, efforts in gene annotation, classification, and functional assignment are still lacking. Therefore, we propose to combine evidence at evolutionary, structural, and functional level to update ODEs as well as key information into an easier classification, particularly of AEs. Finally, the kinetic parameters for putative inhibition of ODEs are discussed in terms of its role in future integrated pest management (IPM) strategies.

Overcoming insecticide resistance through computational inhibitor design

Insecticides allow control of agricultural pests and disease vectors and are vital for global food security and health. The evolution of resistance to insecticides, such as organophosphates (OPs), is a serious and growing concern. OP resistance often involves sequestration or hydrolysis of OPs by carboxylesterases. Inhibiting carboxylesterases could, therefore, restore the effectiveness of OPs for which resistance has evolved. Here, we use covalent virtual screening to produce nano-/picomolar boronic acid inhibitors of the carboxylesterase αE7 from the agricultural pest Lucilia cuprina as well as a common Gly137Asp αE7 mutant that confers OP resistance. These inhibitors, with high selectivity against human acetylcholinesterase and low to no toxicity in human cells and in mice, act synergistically with the OPs diazinon and malathion to reduce the amount of OP required to kill L. cuprina by up to 16-fold and abolish resistance. The compounds exhibit broad utility in significantly potentiating another OP, chlorpyrifos, against the common pest, the peach-potato aphid (Myzus persicae). These compounds represent a solution to OP resistance as well as to environmental concerns regarding overuse of OPs, allowing significant reduction of use without compromising efficacy.

The molecular basis for the neofunctionalization of the juvenile hormone esterase duplication in Drosophila

The Drosophila melanogaster enzymes juvenile hormone esterase (DmJHE) and its duplicate, DmJHEdup, present ideal examples for studying the structural changes involved in the neofunctionalization of enzyme duplicates. DmJHE is a hormone esterase with precise regulation and highly specific activity for its substrate, juvenile hormone. DmJHEdup is an odorant degrading esterase (ODE) responsible for processing various kairomones in antennae. Our phylogenetic analysis shows that the JHE lineage predates the hemi/holometabolan split and that several duplications of JHEs have been templates for the evolution of secreted β-esterases such as ODEs through the course of insect evolution. Our biochemical comparisons further show that DmJHE has sufficient substrate promiscuity and activity against odorant esters for a duplicate to evolve a general ODE function against a range of mid-long chain food esters, as is shown in DmJHEdup. This substrate range complements that of the only other general ODE known in this species, Esterase 6. Homology models of DmJHE and DmJHEdup enabled comparisons between each enzyme and the known structures of a lepidopteran JHE and Esterase 6. Both JHEs showed very similar active sites despite low sequence identity (30%). Both ODEs differed drastically from the JHEs and each other, explaining their complementary substrate ranges. A small number of amino acid changes are identified that may have been involved in the early stages of the neofunctionalization of DmJHEdup. Our results provide key insights into the process of neofunctionalization and the structural changes that can be involved.

Functional analysis of four upregulated carboxylesterase genes associated with fenpropathrin resistance in Tetranychus cinnabarinus (Boisduval)

BACKGROUND

Carboxylesterases (CarEs) are important in pesticide resistance. Four overexpressed CarE genes with inducible character were screened out in fenpropathrin-resistant Tetranychus cinnabarinus, but their functional roles remained to be further analyzed by RNAi and protein expression.

RESULTS

Feeding a single double-stranded (ds)RNA of each of four genes led to gene-specific downregulation of mRNA, decreased esterase activity and diminished resistance in T. cinnabarinus. More interestingly, feeding four dsRNAs simultaneously led to a more significant decrease in enzymatic activity and fold resistance than feeding a single dsRNA individually, suggesting that these CarE genes were involved in fenpropathrin-resistance and had cooperative roles. The gene CarE6 was regarded as the primary and representative candidate to be functionally expressed, because silencing of CarE6 led to the most significant decrease in resistance level. The activity of CarE6 protein was competitively inhibited by fenpropathrin. It could effectively decompose 41.7 ± 0.09% of fenpropathrin within 3 h, proving that CarE6 protein was capable of metabolizing fenpropathrin effectively in T. cinnabarinus.

CONCLUSION

The results confirm that four CarE genes are cooperatively involved in fenpropathrin resistance and the metabolic enzymes encoded by these overexpressed genes do indeed metabolize acaricide in resistant T. cinnabarinus in the evolution of acaricide resistance. © 2018 Society of Chemical Industry.

Integrating Experiment and Theory to Understand TCR-pMHC Dynamics

The conformational dynamism of proteins is well established. Rather than having a single structure, proteins are more accurately described as a conformational ensemble that exists across a rugged energy landscape, where different conformational sub-states interconvert. The interaction between αβ T cell receptors (TCR) and cognate peptide-MHC (pMHC) is no exception, and is a dynamic process that involves substantial conformational change. This review focuses on technological advances that have begun to establish the role of conformational dynamics and dynamic allostery in TCR recognition of the pMHC and the early stages of signaling. We discuss how the marriage of molecular dynamics (MD) simulations with experimental techniques provides us with new ways to dissect and interpret the process of TCR ligation. Notably, application of simulation techniques lags behind other fields, but is predicted to make substantial contributions. Finally, we highlight integrated approaches that are being used to shed light on some of the key outstanding questions in the early events leading to TCR signaling.

The maximum penalty criterion for ridge regression: application to the calibration of the force constant in elastic network models

Tikhonov regularization, or ridge regression, is a popular technique to deal with collinearity in multivariate regression. We unveil a formal analogy between ridge regression and statistical mechanics, where the objective function is comparable to a free energy, and the ridge parameter plays the role of temperature. This analogy suggests two novel criteria for selecting a suitable ridge parameter: specific-heat (C_v) and maximum penalty (MP). We apply these fits to evaluate the relative contributions of rigid-body and internal fluctuations, which are typically highly collinear, to crystallographic B-factors. This issue is particularly important for computational models of protein dynamics, such as the elastic network model (ENM), since the amplitude of the predicted internal motion is commonly calibrated using B-factor data. After validation on simulated datasets, our results indicate that rigid-body motions account on average for more than 80% of the amplitude of B-factors. Furthermore, we evaluate the ability of different fits to reproduce the amplitudes of internal fluctuations in X-ray ensembles from the B-factors in the corresponding single X-ray structures. The new ridge criteria are shown to be markedly superior to the commonly used two-parameter fit that neglects rigid-body rotations and to the full fits regularized under generalized cross-validation. In conclusion, the proposed fits ensure a more robust calibration of the ENM force constant and should prove valuable in other applications.

Previously Hidden Dynamics at the TCR-Peptide-MHC Interface Revealed

A structural characterization of the interaction between αβ TCRs and cognate peptide-MHC (pMHC) is central to understanding adaptive T cell-mediated immunity. X-ray crystallography, although the source of much structural data, traditionally provides only a static snapshot of the protein. Given the emerging evidence for the important role of conformational dynamics in protein function, we interrogated 309 crystallographic structures of pMHC complexes using ensemble refinement, a technique that can extract dynamic information from the x-ray data. Focusing on a subset of human pMHC class I systems, we found that in many cases, ensemble methods were able to uncover previously hidden evidence of significant conformational plasticity, thereby revealing additional information that can build upon and significantly enhance functional interpretations that are based on a single static structure. Notable examples include the interpretation of differences in the disease association of HLA subtypes, the relationship between peptide prominence and TCR recognition, the role of conformational flexibility in vaccine design, and the discrimination between induced fit and conformational selection models of TCR binding. We show that the currently widespread practice of analyzing pMHC interactions via the study of a single crystallographic structure does not make use of pertinent and easily accessible information from x-ray data concerning alternative protein conformations. This new analysis therefore not only highlights the capacity for ensemble methods to significantly enrich the interpretation of decades of structural data but also provides previously missing information concerning the dynamics of existing characterized TCR-pMHC interactions.

Functional characterization of carboxylesterase gene mutations involved in Aphis gossypii resistance to organophosphate insecticides

Carboxylesterases (CarEs) play an important role in detoxifying insecticides in insects. Over-expression and structural modification of CarEs have been implicated in the development of organophosphate (OP) insecticide resistance in insects. A previous study identified four nonsynonymous mutations (resulting in four amino acid residue substitutions) in the open reading frame of the carboxylesterase gene of resistant cotton aphids compared to the omethoate susceptible strain, which has possibly influenced the development of resistance to omethoate (a systemic OP insecticide). The current study further characterized the function of these mutations, both alone and in combination, in the hydrolysis of OP insecticides. The metabolism results suggest that the combination of four mutations, mainly existing in the laboratory-selected OP-resistant cotton aphid population, increased the OP hydrolase activity (approximately twofold) at the cost of detectable carboxylesterase activity. The functional studies of single or multiple mutations suggest the positive effect of H104R, A128V and T333P on the acquisition of OP hydrolase activity, especially the combination of H104R with A128V or T333P. K484R substitution decreased both the OP hydrolase activity and the CarE activity, indicating that this mutation primarily drives the negative effect on the acquisition of OP hydrolase activity amongst these four mutations in the resistant strain. The modelling and docking results are basically consistent with the metabolic results, which strongly suggest that the structural gene modification is the molecular basis for the OP resistance in this laboratory-selected cotton aphid strain.

Identifying and Visualizing Macromolecular Flexibility in Structural Biology

Structural biology comprises a variety of tools to obtain atomic resolution data for the investigation of macromolecules. Conventional structural methodologies including crystallography, NMR and electron microscopy often do not provide sufficient details concerning flexibility and dynamics, even though these aspects are critical for the physiological functions of the systems under investigation. However, the increasing complexity of the molecules studied by structural biology (including large macromolecular assemblies, integral membrane proteins, intrinsically disordered systems, and folding intermediates) continuously demands in-depth analyses of the roles of flexibility and conformational specificity involved in interactions with ligands and inhibitors. The intrinsic difficulties in capturing often subtle but critical molecular motions in biological systems have restrained the investigation of flexible molecules into a small niche of structural biology. Introduction of massive technological developments over the recent years, which include time-resolved studies, solution X-ray scattering, and new detectors for cryo-electron microscopy, have pushed the limits of structural investigation of flexible systems far beyond traditional approaches of NMR analysis. By integrating these modern methods with powerful biophysical and computational approaches such as generation of ensembles of molecular models and selective particle picking in electron microscopy, more feasible investigations of dynamic systems are now possible. Using some prominent examples from recent literature, we review how current structural biology methods can contribute useful data to accurately visualize flexibility in macromolecular structures and understand its important roles in regulation of biological processes.