Cation–π, amino–π, π–π, and H‐bond interactions stabilize antigen–antibody interfaces

The identification of immunogenic regions on the surface of antigens, which are able to stimulate an immune response, is a major challenge for the design of new vaccines. Computational immunology aims at predicting such regions—in particular B‐cell epitopes—but is far from being reliably applicable on a large scale. To gain understanding into the factors that contribute to the antigen–antibody affinity and specificity, we perform a detailed analysis of the amino acid composition and secondary structure of antigen and antibody surfaces, and of the interactions that stabilize the complexes, in comparison with the composition and interactions observed in other heterodimeric protein interfaces. We make a distinction between linear and conformational B‐cell epitopes, according to whether they consist of successive residues along the polypeptide chain or not. The antigen–antibody interfaces were shown to differ from other protein–protein interfaces by their smaller size, their secondary structure with less helices and more loops, and the interactions that stabilize them: more H‐bond, cation–π, amino–π, and π–π interactions, and less hydrophobic packing; linear and conformational epitopes can clearly be distinguished. Often, chains of successive interactions, called cation/amino–π and π–π chains, are formed. The amino acid composition differs significantly between the interfaces: antigen–antibody interfaces are less aliphatic and more charged, polar and aromatic than other heterodimeric protein interfaces. Moreover, paratopes and epitopes—albeit to a lesser extent—have amino acid compositions that are distinct from general protein surfaces. This specificity holds promise for improving B‐cell epitope prediction. Proteins 2014; 82:1734–1746. © 2014 Wiley Periodicals, Inc.     

Stabilization of protein structure through π–π interaction in the second coordination sphere of pseudoazurin

Noncovalent, weak interactions in the second coordination sphere of the copper active site of Pseudoazurin (PAz) from Achromobacter cycloclastes were examined using a series of Met16X variants. In this study, the differences in protein stability due to the changes in the nature of the 16th amino acid (Met, Phe, Val, Ile) were investigated by electrospray ionization mass spectrometry (ESI‐MS) and far‐UV circular dichroism (CD) as a result of acid denaturation. The percentage of native states (folded holo forms) of Met16Phe variants was estimated to be 75% at pH 2.9 although the wild‐type (WT), Met16Val and Met16Ile PAz, became completely unfolded. The high stability under acidic conditions is correlated with the result of the active site being stabilized by the aromatic substitution of the Met16 residue. The π–π interaction in the second coordination sphere makes a significant contribution to the stability of active site and the protein matrix.     

Evidence of π‐stacking interactions in the self‐assembly of hIAPP22‐29

The role aromatic amino acids play in the formation of amyloid is a subject of controversy. In an effort to clarify the contribution of aromaticity to the self‐assembly of human islet amyloid polypeptide (hIAPP)_22‐29, peptide analogs containing electron donating groups (EDGs) or electron withdrawing groups (EWGs) as substituents on the aromatic ring of Phe‐23 at the para position have been synthesized and characterized using turbidity measurements in conjunction with Raman and fluorescence spectroscopy. Results indicate the incorporation of EDGs on the aromatic ring of Phe‐23 virtually abolish the ability of hIAPP_22‐29 to form amyloid. Peptides containing EWGs were still capable of forming aggregates. These aggregates were found to be rich in β‐sheet secondary structure. Transmission electron microscopy images of the aggregates confirm the presence of amyloid fibrils. The observed difference in amyloidogenic propensity between peptides containing EDGs and those with EWGs appears not to be based on differences in peptide hydrophobicity. Fluorescence and Raman spectroscopic investigations reveal that the environment surrounding the aromatic ring becomes more hydrophobic and ordered upon aggregation. Furthermore, Raman measurements of peptide analogs containing EWGs, conclusively demonstrate a distinct downshift in the ? C?C? ring mode (ca. 1600 cm^?1) upon aggregation that has previously been shown to be indicative of π‐stacking. While previous work has demonstrated that π‐stacking is not an absolute requirement for fibrillization, our findings indicate that Phe‐23 also contributes to fibril formation through π‐stacking interactions and that it is not only the hydrophobic nature of this residue that is relevant in the self‐assembly of hIAPP_22‐29. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

The β‐sheet breakers and π‐stacking

Fibrillation of β‐amyloid is recognized as a key process leading to the development of Alzheimer's disease. Small peptides called β‐sheet breakers were found to inhibit the process of β‐amyloid fibrillation and to dissolve amyloid fibrils in vitro, in vivo, and in cell culture studies [1,2]. The mechanism by which peptide inhibition takes place remains elusive and a detailed model needs to be established. Here, we present new insights into the possible role of consecutive Phe residues, present in the structure of β‐sheet breakers, supported by the results obtained by means of MD simulations. We performed a 30‐ns MD of two β‐sheet breakers: iAβ5 (LPFFD) and iAβ6 (LPFFFD) which have two and three consecutive Phe residues, respectively. We have found that Phe rings in these peptides tend to form stacked conformations. For one of the peptides – iAβ6 – the calculated electrostatic contribution to free energy of one of the conformers with three rings stacked (c2) is significantly lower than that corresponding to the unstacked one (c1), two rings stacked (c0) and second conformer with three rings stacked (c3). This may favor the interaction of the c2 conformer with the target on amyloid fibril. We hypothesize that the mechanism of inhibition of amyloidogenesis by β‐sheet breaker involves competition among π‐stacked Phe residues of the inhibitor and π‐stacking within the β‐amyloid fibril. iAβ6 may be a promising candidate for a lead compound of amyloidogenesis inhibitors. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Essential requirements for substrate binding affinity and selectivity toward human CYP2 family enzymes

A detailed analysis of substrate selectivity within the cytochrome P450 2 (CYP2) family is reported. From a consideration of specific interactions between drug substrates for human CYP2 family enzymes and the putative active sites of CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, and CYP2E1, it is likely that the number and disposition of hydrogen bond donor/acceptors and aromatic rings within the various P450 substrate molecules determines their enzyme selectivity and binding affinity, together with directing their preferred routes of metabolism by the CYP2 enzymes concerned. Although many aliphatic residues are present in most P450 active sites, it would appear that their main contribution centers around hydrophobic interactions and desolvation processes accompanying substrate binding. Molecular modeling studies based on the recent CYP2C5 crystal structure appear to show close agreement with site-directed mutagenesis experiments and with information on substrate metabolism and selectivity within the CYP2 family.     

Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine

The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure–function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.     

Structures of human cytochrome P-450 2E1. Insights into the binding of inhibitors and both small molecular weight and fatty acid substrates

Human microsomal cytochrome P-450 2E1 (CYP2E1) monooxygenates > 70 low molecular weight xenobiotic compounds, as well as much larger endogenous fatty acid signaling molecules such as arachidonic acid. In the process, CYP2E1 can generate toxic or carcinogenic compounds, as occurs with acetaminophen overdose, nitrosamines in cigarette smoke, and reactive oxygen species from uncoupled catalysis. Thus, the diverse roles that CYP2E1 has in normal physiology, toxicity, and drug metabolism are related to its ability to metabolize diverse classes of ligands, but the structural basis for this was previously unknown. Structures of human CYP2E1 have been solved to 2.2 angstroms for an indazole complex and 2.6 angstroms for a 4-methylpyrazole complex. Both inhibitors bind to the heme iron and hydrogen bond to Thr303 within the active site. Complementing its small molecular weight substrates, the hydrophobic CYP2E1 active site is the smallest yet observed for a human cytochrome P-450. The CYP2E1 active site also has two adjacent voids: one enclosed above the I helix and the other forming a channel to the protein surface. Minor repositioning of the Phe478 aromatic ring that separates the active site and access channel would allow the carboxylate of fatty acid substrates to interact with conserved 216QXXNN220 residues in the access channel while positioning the hydrocarbon terminus in the active site, consistent with experimentally observed omega-1 hydroxylation of saturated fatty acids. Thus, these structures provide insights into the ability of CYP2E1 to effectively bind and metabolize both small molecule substrates and fatty acids.     

Side‐chain hydrophobicity and the stability of Aβ16–22 aggregates

Recent mutagenesis studies using the hydrophobic segment of Aβ suggest that aromatic π‐stacking interactions may not be critical for fibril formation. We have tested this conjecture by probing the effect of Leu, Ile, and Ala mutation of the aromatic Phe residues at positions 19 and 20, on the double‐layer hexametric chains of Aβ fragment Aβ_16–22 using explicit solvent all‐atom molecular dynamics. As these simulations rely on the accuracy of the utilized force fields, we first evaluated the dynamic and stability dependence on various force fields of small amyloid aggregates. These initial investigations led us to choose AMBER99SB‐ILDN as force field in multiple long molecular dynamics simulations of 100 ns that probe the stability of the wild‐type and mutants oligomers. Single‐point and double‐point mutants confirm that size and hydrophobicity are key for the aggregation and stability of the hydrophobic core region (Aβ_16–22). This suggests as a venue for designing Aβ aggregation inhibitors the substitution of residues (especially, Phe 19 and 20) in the hydrophobic region (Aβ_16–22) with natural and non‐natural amino acids of similar size and hydrophobicity.     

Effects of active site cleft residues on oligosaccharide binding,hydrolysis, and glycosynthase activities of rice BGlu1 and its mutants

Rice BGlu1 (Os3BGlu7) is a glycoside hydrolase family 1 β‐glucosidase that hydrolyzes cellooligosaccharides with increasing efficiency as the degree of polymerization (DP) increases from 2 to 6, indicating six subsites for glucosyl residue binding. Five subsites have been identified in X‐ray crystal structures of cellooligosaccharide complexes with its E176Q acid‐base and E386G nucleophile mutants. X‐ray crystal structures indicate that cellotetraose binds in a similar mode in BGlu1 E176Q and E386G, but in a different mode in the BGlu1 E386G/Y341A variant, in which glucosyl residue 4 (Glc4) interacts with Q187 instead of the eliminated phenolic group of Y341. Here, we found that the Q187A mutation has little effect on BGlu1 cellooligosaccharide hydrolysis activity or oligosaccharide binding in BGlu1 E176Q, and only slight effects on BGlu1 E386G glycosynthase activity. X‐ray crystal structures showed that cellotetraose binds in a different position in BGlu1 E176Q/Y341A, in which it interacts directly with R178 and W337, and the Q187A mutation had little effect on cellotetraose binding. Mutations of R178 and W337 to A had significant and nonadditive effects on oligosaccharide hydrolysis by BGlu1, pNPGlc cleavage and cellooligosaccharide inhibition of BGlu1 E176Q and BGlu1 E386G glycosynthase activity. Hydrolysis activity was partially rescued by Y341 for longer substrates, suggesting stacking of Glc4 on Y341 stabilizes binding of cellooligosaccharides in the optimal position for hydrolysis. This analysis indicates that complex interactions between active site cleft residues modulate substrate binding and hydrolysis.     

Modeling interactions between a β‐O‐4 type lignin model compound and 1‐allyl‐3‐methylimidazolium chloride ionic liquid

Aiming at understanding the molecular mechanism of the lignin dissolution in imidazolium‐based ionic liquids (ILs), this work presents a combined quantum chemistry (QC) calculation and molecular dynamics (MD) simulation study on the interaction of the lignin model compound, veratrylglycerol‐β‐guaiacyl ether (VG) with 1‐allyl‐3‐methylimidazolium chloride ([Amim]Cl). The monomer of VG is shown to feature a strong intramolecular hydrogen bond, and its dimer is indicated to present important π‐π stacking and intermolecular hydrogen bonding interactions. The interactions of both the cation and anion of [Amim]Cl with VG are shown to be stronger than that between the two monomers, indicating that [Amim]Cl is capable of dissolving lignin. While Cl^– anion forms a hydrogen‐bonded complex with VG, the imidazolium cation interacts with VG via both the π‐π stacking and intermolecular hydrogen bonding. The calculated interaction energies between VG and the IL or its components (the cation, anion, and ion pair) indicate the anion plays a more important role than the cation for the dissolution of lignin in the IL. Theoretical results provide help for understanding the molecular mechanism of lignin dissolution in imidazolium‐based IL. The theoretical calculations on the interaction between the lignin model compound and [Amim]Cl ionic liquid indicate that the anion of [Amim]Cl plays a more important role for lignin dissolution although the cation also makes a substantial contribution.     

Substrate determinants in the C99 juxtamembrane domains differentially affect γ–secretase cleavage specificity and modulator pharmacology

The molecular mechanisms governing γ‐secretase cleavage specificity are not fully understood. Herein, we demonstrate that extending the transmembrane domain of the amyloid precursor protein‐derived C99 substrate in proximity to the cytosolic face strongly influences γ–secretase cleavage specificity. Sequential insertion of leucines or replacement of membrane‐anchoring lysines by leucines elevated the production of Aβ42, whilst lowering production of Aβ40. A single insertion or replacement was sufficient to produce this phenotype, suggesting that the helical length distal to the ε–site is a critical determinant of γ‐secretase cleavage specificity. Replacing the lysine at the luminal membrane border (K28) with glutamic acid (K28E) increased Aβ37 and reduced Aβ42 production. Maintaining a positive charge with an arginine replacement, however, did not alter cleavage specificity. Using two potent and structurally distinct γ–secretase modulators (GSMs), we elucidated the contribution of K28 to the modulatory mechanism. Surprisingly, whilst lowering the potency of the non‐steroidal anti‐inflammatory drug‐type GSM, the K28E mutation converted a heteroaryl‐type GSM to an inverse GSM. This result implies the proximal lysine is critical for the GSM mechanism and pharmacology. This region is likely a major determinant for substrate binding and we speculate that modulation of substrate binding is the fundamental mechanism by which GSMs exert their action.     

Molecular modelling of human CYP1B1 substrate interactions and investigation of allelic variant effects on metabolism

Molecular modelling of human CYP1B1 based on homology with the mammalian P450, CYP2C5, of known three-dimensional structure is reported. The enzyme model has been used to investigate the likely mode of binding for selected CYP1B1 substrates, particularly with regard to the possible effects of allelic variants of CYP1B1 on metabolism. In general, it appears that the CYP1B1 model is consistent with known substrate selectivity for the enzyme, and the sites of metabolism can be rationalized in terms of specific contacts with key amino acid residues within the CYP1B1 heme locus. Furthermore, a mode of binding interaction for the inhibitor, alpha-naphthoflavone, is presented which accords with currently available information. The current paper shows that a combination of molecular modelling and experimental determinations on the substrate metabolism for CYP1B1 allelic variants can aid in the understanding of structure-function relationships within P450 enzymes.     

Substrate recruitment of γ‐secretase and mechanism of clinical presenilin mutations revealed by photoaffinity mapping

Intramembrane proteases execute fundamental biological processes ranging from crucial signaling events to general membrane proteostasis. Despite the availability of structural information on these proteases, it remains unclear how these enzymes bind and recruit substrates, particularly for the Alzheimer's disease‐associated γ‐secretase. Systematically scanning amyloid precursor protein substrates containing a genetically inserted photocrosslinkable amino acid for binding to γ‐secretase allowed us to identify residues contacting the protease. These were primarily found in the transmembrane cleavage domain of the substrate and were also present in the extramembranous domains. The N‐terminal fragment of the catalytic subunit presenilin was determined as principal substrate‐binding site. Clinical presenilin mutations altered substrate binding in the active site region, implying a pathogenic mechanism for familial Alzheimer's disease. Remarkably, PEN‐2 was identified besides nicastrin as additional substrate‐binding subunit. Probing proteolysis of crosslinked substrates revealed a mechanistic model of how these subunits interact to mediate a stepwise transfer of bound substrate to the catalytic site. We propose that sequential binding steps might be common for intramembrane proteases to sample and select cognate substrates for catalysis.     

Biochemical characterization and structural insight into aliphatic β‐amino acid adenylation enzymes IdnL1 and CmiS6

Macrolactam antibiotics such as incednine and cremimycin possess an aliphatic β‐amino acid as a starter unit of their polyketide chain. In the biosynthesis of incednine and cremimycin, unique stand‐alone adenylation enzymes IdnL1 and CmiS6 select and activate the proper aliphatic β‐amino acid as a starter unit. In this study, we describe the enzymatic characterization and the structural basis of substrate specificity of IdnL1 and CmiS6. Functional analysis revealed that IdnL1 and CmiS6 recognize 3‐aminobutanoic acid and 3‐aminononanoic acid, respectively. We solved the X‐ray crystal structures of IdnL1 and CmiS6 to understand the recognition mechanism of these aliphatic β‐amino acids. These structures revealed that IdnL1 and CmiS6 share a common recognition motif that interacts with the β‐amino group of the substrates. However, the hydrophobic side‐chains of the substrates are accommodated differently in the two enzymes. IdnL1 has a bulky Leu220 located close to the terminal methyl group of 3‐aminobutanoate of the trapped acyl‐adenylate intermediate to construct a shallow substrate‐binding pocket. In contrast, CmiS6 possesses Gly220 at the corresponding position to accommodate 3‐aminononanoic acid. This structural observation was supported by a mutational study. Thus, the size of amino acid residue at the 220 position is critical for the selection of an aliphatic β‐amino acid substrate in these adenylation enzymes. Proteins 2017; 85:1238–1247. © 2017 Wiley Periodicals, Inc.     

Homology modeling and substrate binding study of human CYP2C9 enzyme

The CYP2C subfamily of human liver P450 isozymes is of major importance in drug metabolism. The most abundant 2C isozyme, CYP2C9, regioselectively hydroxylates a wide variety of substrates. A major obstacle to understanding this specificity in human CYP2C9 is the absence of a 3D structure. A 3D model of CYP2C9 was built, assessed, and used to characterize explicit enzyme-substrate complexes using methods previously developed in our laboratory. The 3D model was assessed by determining its stability to unconstrained molecular dynamics and by comparison of specific properties with those of known protein structures. The CYP2C9 model was then used to characterize explicit enzyme complexes with three structurally and chemically diverse substrates: (S)-naproxen, phenytoin, and progesterone. Each substrate was found to bind to the enzyme with a favorable interaction energy and to remain in the binding site during unconstrained molecular dynamics. Moreover, the mode of binding of each substrate led to calculated preferred hydroxylation sites consistent with experiment. Binding-site residues identified for the models included Arg 105 and Arg97 as key cationic residues, as well as Asn 202, Asp 293, Pro 101, Leu 102, Gly 296, and Phe 476. Site-specific mutations are proposed for further integrated computational and experimental study.     

Crystal structure of Saccharomyces cerevisiae Aro8, a putative α‐aminoadipate aminotransferase

α‐Aminoadipate aminotransferase (AAA‐AT) catalyzes the amination of 2‐oxoadipate to α‐aminoadipate in the fourth step of the α‐aminoadipate pathway of lysine biosynthesis in fungi. The aromatic aminotransferase Aro8 has recently been identified as an AAA‐AT in Saccharomyces cerevisiae. This enzyme displays broad substrate selectivity, utilizing several amino acids and 2‐oxo acids as substrates. Here we report the 1.91Å resolution crystal structure of Aro8 and compare it to AAA‐AT LysN from Thermus thermophilus and human kynurenine aminotransferase II. Inspection of the active site of Aro8 reveals asymmetric cofactor binding with lysine‐pyridoxal‐5‐phosphate bound within the active site of one subunit in the Aro8 homodimer and pyridoxamine phosphate and a HEPES molecule bound to the other subunit. The HEPES buffer molecule binds within the substrate‐binding site of Aro8, yielding insights into the mechanism by which it recognizes multiple substrates and how this recognition differs from other AAA‐AT/kynurenine aminotransferases.     

Design of anti‐thyroid drugs: Binding studies and structure determination of the complex of lactoperoxidase with 2‐mercaptoimidazole at 2.30 Å resolution

Lactoperoxidase (LPO) belongs to mammalian heme peroxidase superfamily, which also includes myeloperoxidase (MPO), eosinophil peroxidase (EPO), and thyroid peroxidase (TPO). LPO catalyzes the oxidation of a number of substrates including thiocyanate while TPO catalyzes the biosynthesis of thyroid hormones. LPO is also been shown to catalyze the biosynthesis of thyroid hormones indicating similar functional and structural properties. The binding studies showed that 2‐mercaptoimidazole (MZY) bound to LPO with a dissociation constant of 0.63 µM. The inhibition studies showed that the value of IC₅₀ was 17 µM. The crystal structure of the complex of LPO with MZY showed that MZY bound to LPO in the substrate‐binding site on the distal heme side. MZY was oriented in the substrate‐binding site in such a way that the sulfur atom is at a distance of 2.58 Å from the heme iron. Previously, a similar compound, 3‐amino‐1,2,4‐triazole (amitrole) was also shown to bind to LPO in the substrate‐binding site on the distal heme side. The amino nitrogen atom of amitrole occupied the same position as that of sulfur atom in the present structure indicating a similar mode of binding. Recently, the structure of the complex of LPO with a potent antithyroid drug, 1‐methylimidazole‐2‐thiol (methimazole, MMZ) was also determined. It showed that MMZ bound to LPO in the substrate‐binding site on the distal heme side with 2 orientations. The position of methyl group was same in the 2 orientations while the positions of sulfur atom differed indicating a higher preference for a methyl group.     

Hsp74/14‐3‐3σ Complex Mediates Centrosome Amplification by High Glucose,Insulin, and Palmitic Acid

It has been reported recently that type 2 diabetes promotes centrosome amplification via 14‐3‐3σ/ROCK1 complex. In the present study, 14‐3‐3σ interacting proteins are characterized and their roles in the centrosome amplification by high glucose, insulin, and palmitic acid are investigated. Co‐immunoprecipitation in combination with MS analysis identified 134 proteins that interact with 14‐3‐3σ, which include heat shock 70 kDa protein 4 (Hsp74). Gene ontology analyses reveal that many of them are enriched in binding activity. Kyoto Encyclopedia of Genes and Genomes analysis shows that the top three enriched pathways are ribosome, carbon metabolism, and biosynthesis of amino acids. Molecular and functional investigations show that the high glucose, insulin, and palmitic acid increase the expression and binding of 14‐3‐3σ and Hsp74 as well as centrosome amplification, all of which are inhibited by knockdown of 14‐3‐3σ or Hsp74. Moreover, molecular docking analysis shows that the interaction between the 14‐3‐3σ and the Hsp74 is mainly through hydrophobic contacts and a lesser degree ionic interactions and hydrogen bond by different amino acids residues. In conclusion, the results suggest that the experimental treatment triggers centrosome amplification via upregulations of expression and binding of 14‐3‐3σ and Hsp74.