Direct evidence that the carboxyl-terminal sequence of a bacterial chemoreceptor is an unstructured linker and enzyme tether

Sensory adaptation in bacterial chemotaxis involves reversible methylation of specific glutamyl residues on chemoreceptors. The reactions are catalyzed by a dedicated methyltransferase and dedicated methylesterase. In Escherichia coli and related organisms, control of these enzymes includes an evolutionarily recent addition of interaction with a pentapeptide activator located at the carboxyl terminus of the receptor polypeptide chain. Effective enzyme activation requires not only the pentapeptide but also a segment of the receptor polypeptide chain between that sequence and the coiled-coil body of the chemoreceptor. This segment has features consistent with a role as a flexible and presumably unstructured linker and enzyme tether, but there has been no direct information about its structure. We used site-directed spin labeling and electron paramagnetic resonance spectroscopy to characterize structural features of the carboxyl-terminal 40 residues of E. coli chemoreceptor Tar. Beginning ～ 35 residues from the carboxyl terminus and continuing to the end of the protein, spectra of spin-labeled Tar embedded in native membranes or in reconstituted proteoliposomes, exhibited mobilities characteristic of unstructured, disordered segments. Binding of methyltransferase substantially reduced mobility for positions in or near the pentapeptide but mobility for the linker sequence remained high, being only modestly reduced in a gradient of decreasing effects for 10-15 residues, a pattern consistent with the linker providing a flexible arm that would allow enzyme diffusion within defined limits. Thus, our data identify that the carboxyl-terminal linker between the receptor body and the pentapeptide is an unstructured, disordered segment that can serve as a flexible arm and enzyme tether.     

Modeling the transmembrane domain of bacterial chemoreceptors

Bacterial chemoreceptors signal across the membrane by conformational changes that traverse a four-helix transmembrane domain. High-resolution structures are available for the chemoreceptor periplasmic domain and part of the cytoplasmic domain but not for the transmembrane domain. Thus, we constructed molecular models of the transmembrane domains of chemoreceptors Trg and Tar, using coordinates of an unrelated four-helix coiled coil as a template and the X-ray structure of a chemoreceptor periplasmic domain to establish register and positioning. We tested the models using the extensive data for cross-linking propensities between cysteines introduced into adjacent transmembrane helices, and we found that many aspects of the models corresponded with experimental observations. The one striking disparity, the register of transmembrane helix 2 (TM2) relative to its partner transmembrane helix 1, could be corrected by sliding TM2 along its long axis toward the periplasm. The correction implied that axial sliding of TM2, the signaling movement indicated by a large body of data, was of greater magnitude than previously thought. The refined models were used to assess effects of inter-helical disulfides on the two ligand-induced conformational changes observed in alternative crystal structures of periplasmic domains: axial sliding within a subunit and subunit rotation. Analyses using a measure of disulfide potential energy provided strong support for the helical sliding model of transmembrane signaling but indicated that subunit rotation could be involved in other ligand-induced effects. Those analyses plus modeled distances between diagnostic cysteine pairs indicated a magnitude for TM2 sliding in transmembrane signaling of several angstroms.     

Molecular modeling and multi‐spectroscopic approaches to study the interaction between antibacterial drug and human immunoglobulin G

Mechanistic and conformational studies on the interaction of sulfamethoxazole (SMX) with human immunoglobulin G (HIgG) were performed by molecular modeling and multi‐spectroscopic methods. The interaction mechanism was firstly predicted through molecular modeling that confirmed the interaction between SMX and HIgG. The binding parameters and thermodynamic parameters at different temperatures had been calculated according to the Stern?Volmer, Scatchard, Sips and Van ’t Hoff equations, respectively. Experimental results showed that the fluorescence intensity of HIgG was quenched by the gradual addition of SMX. The binding constants of SMX with HIgG decreased with the increase of temperature, which meant that the quenching mechanism was a static quenching. Meanwhile, the results also confirmed that there was one independent class of binding site on HIgG for SMX during their interaction. The thermodynamic parameters of the reaction, namely standard enthalpy ΔH⁰ and entropy ΔS⁰, had been calculated to be ?14.69 kJ·mol^?1 and 22.99 J·mol^?1·K^?1, respectively, which suggested that the electrostatic and hydrophobic interactions were the predominant intermolecular forces in stabilizing the SMX?HIgG complex. Furthermore, experimental results obtained from three‐dimensional fluorescence spectroscopy, UV?vis absorption spectroscopy and circular dichroism (CD) spectroscopy confirmed that the conformational structure of HIgG was altered in the presence of SMX. Copyright © 2015 John Wiley & Sons, Ltd.     

Molecular modeling and design of regioselectively addressable functionalized templates with rigidified three‐dimensional structures

Extensive conformational analysis of a series of β‐alkyl substituted cyclopeptides—cyclo(Pro¹–Xaa²–Nle³–Ala⁴–Nle⁵–Pro⁶–Xaa⁷–Nle⁸–Ala⁹–Nle¹⁰) and cyclo[Pro¹–Xaa²–Nle³–(Cys⁴– Nle⁵–Pro⁶–Xaa⁷–Nle⁸–Cys⁹)–Nle¹⁰] as well as their corresponding unsubstituted core structures cyclo(Pro¹–Xaa²–Ala³–Ala⁴–Ala⁵–Pro⁶–Xaa⁷–Ala⁸–Ala⁹–Ala¹⁰) and cyclo(Pro¹–Xaa²–Ala³–Cys⁴– Ala⁵–Pro⁶–Xaa⁷–Ala⁸–Cys⁹–Ala¹⁰) has been performed employing both the ECEPP/2 and the MAB force fields (Xaa = Gly, L ‐Ala, D ‐Ala, Aib, and D ‐Pro). Results show that (a) possible three‐dimensional structures of the cyclo(Pro¹–Gly²–Lys³–Ala⁴–Lys⁵–Pro⁶–Gly⁷–Lys⁸–Ala⁹–Lys¹⁰) molecule are not limited to a single extended “rectangular” conformation with all Lys side chains oriented at the same side of the molecule; (b) conformational equilibrium in monocyclic analogues obtained by replacements of conformationally flexible Gly residues for L ‐Ala, D ‐Ala, Aib, or D ‐Pro is not significantly shifted towards the target “rectangular” conformational type; and (c) introduction of disulfide bridges between positions 4 and 9 is a very powerful way to stabilize the target conformations in the resulting bicyclic molecules. These findings form the basis for further design of rigidified regioselectively addressable functionalized templates with many application areas ranging from biostructural to diagnostic purposes. © 1999 John Wiley & Sons, Inc. Biopoly 50: 361–372, 1999     

Molecular modeling of pathogenesis-related proteins of family 5

The family of pathogenesis-related (PR) 5 proteins have diverse functions, and some of them are classified as thaumatins, osmotins, and inhibitors of α-amylase or trypsin. Although the specific function of many PR5 in plants is unknown, they are involved in the acquired systemic resistance and response to biotic stress, causing the inhibition of hyphal growth and reduction of spore germination, probably by a membrane permeabilization mechanism or by interaction with pathogen receptors. We have constructed three-dimensional models of four proteins belonging to the Rosaceae and Fagaceae botanical families by using the technique of comparative molecular modelling by homology. There are four main structural differences between all the PR5, corresponding to regions with replacements of amino acids. Folding and the secondary structures are very similar for all of them. However, the isoelectric point and charge distributions differ for earch protein.     

Systems‐level interactions between insulin–EGF networks amplify mitogenic signaling

Crosstalk mechanisms have not been studied as thoroughly as individual signaling pathways. We exploit experimental and computational approaches to reveal how a concordant interplay between the insulin and epidermal growth factor (EGF) signaling networks can potentiate mitogenic signaling. In HEK293 cells, insulin is a poor activator of the Ras/ERK (extracellular signal‐regulated kinase) cascade, yet it enhances ERK activation by low EGF doses. We find that major crosstalk mechanisms that amplify ERK signaling are localized upstream of Ras and at the Ras/Raf level. Computational modeling unveils how critical network nodes, the adaptor proteins GAB1 and insulin receptor substrate (IRS), Src kinase, and phosphatase SHP2, convert insulin‐induced increase in the phosphatidylinositol‐3,4,5‐triphosphate (PIP₃) concentration into enhanced Ras/ERK activity. The model predicts and experiments confirm that insulin‐induced amplification of mitogenic signaling is abolished by disrupting PIP₃‐mediated positive feedback via GAB1 and IRS. We demonstrate that GAB1 behaves as a non‐linear amplifier of mitogenic responses and insulin endows EGF signaling with robustness to GAB1 suppression. Our results show the feasibility of using computational models to identify key target combinations and predict complex cellular responses to a mixture of external cues.     

Investigations on the interactions between naphthalimide‐based anti‐tumor drugs and human serum albumin by spectroscopic and molecular modeling methods

The interactions between the three kinds of naphthalimide‐based anti‐tumor drugs (NADA, NADB, NADC) and human serum albumin (HSA) under simulated physiological conditions were investigated by fluorescence spectroscopy, circular dichroism spectroscopy and molecular modeling. The results of the fluorescence quenching spectroscopy showed that the quenching mechanisms for different drugs were static and their affinity was in a descending order of NADA > NADB > NADC. The relative thermodynamic parameters indicated that hydrophobic force was the predominant intermolecular force in the binding of NAD to HSA, while van der Waals interactions and hydrogen bonds could not be ignored. The results of site marker competitive experiment confirmed that the binding site of HSA primarily took place in site I. Furthermore, the molecular modeling study was consistent with these results. The study of circular dichroism spectra demonstrated that the presence of NADs decreased the α‐helical content of HSA and induced the change of the secondary structure of HSA. Copyright © 2015 John Wiley & Sons, Ltd.     

Molecular modeling reveals binding interface of γ‐tubulin with GCP4 and interactions with noscapinoids

The initiation of microtubule assembly within cells is guided by a cone shaped multi‐protein complex, γ‐tubulin ring complex (γTuRC) containing γ‐tubulin and atleast five other γ‐tubulin‐complex proteins (GCPs), i.e., GCP2, GCP3, GCP4, GCP5, and GCP6. The rim of γTuRC is a ring of γ‐tubulin molecules that interacts, via one of its longitudinal interfaces, with GCP2, GCP3, or GCP4 and, via other interface, with α/β?tubulin dimers recruited for the microtubule lattice formation. These interactions however, are not well understood in the absence of crystal structure of functional reconstitution of γTuRC subunits. In this study, we elucidate the atomic interactions between γ‐tubulin and GCP4 through computational techniques. We simulated two complexes of γ‐tubulin‐GCP4 complex (we called dimer1 and dimer2) for 25 ns to obtain a stable complex and calculated the ensemble average of binding free energies of ?158.82 and ?170.19 kcal/mol for dimer1 and ?79.53 and ?101.50 kcal/mol for dimer2 using MM‐PBSA and MM‐GBSA methods, respectively. These highly favourable binding free energy values points to very robust interactions between GCP4 and γ‐tubulin. From the results of the free‐energy decomposition and the computational alanine scanning calculation, we identified the amino acids crucial for the interaction of γ‐tubulin with GCP4, called hotspots. Furthermore, in the endeavour to identify chemical leads that might interact at the interface of γ‐tubulin‐GCP4 complex; we found a class of compounds based on the plant alkaloid, noscapine that binds with high affinity in a cavity close to γ‐tubulin‐GCP4 interface compared with previously reported compounds. All noscapinoids displayed stable interaction throughout the simulation, however, most robust interaction was observed for bromo‐noscapine followed by noscapine and amino‐noscapine. This offers a novel chemical scaffold for γ‐tubulin binding drugs near γ‐tubulin‐GCP4 interface. Proteins 2015; 83:827–843. © 2015 Wiley Periodicals, Inc.     

Fast and accurate modeling of protein-protein interactions by combining template-interface-based docking with flexible refinement

The similarity between folding and binding led us to posit the concept that the number of protein-protein interface motifs in nature is limited, and interacting protein pairs can use similar interface architectures repeatedly, even if their global folds completely vary. Thus, known protein-protein interface architectures can be used to model the complexes between two target proteins on the proteome scale, even if their global structures differ. This powerful concept is combined with a flexible refinement and global energy assessment tool. The accuracy of the method is highly dependent on the structural diversity of the interface architectures in the template dataset. Here, we validate this knowledge-based combinatorial method on the Docking Benchmark and show that it efficiently finds high-quality models for benchmark complexes and their binding regions even in the absence of template interfaces having sequence similarity to the targets. Compared to "classical" docking, it is computationally faster; as the number of target proteins increases, the difference becomes more dramatic. Further, it is able to distinguish binders from nonbinders. These features allow performing large-scale network modeling. The results on an independent target set (proteins in the p53 molecular interaction map) show that current method can be used to predict whether a given protein pair interacts. Overall, while constrained by the diversity of the template set, this approach efficiently produces high-quality models of protein-protein complexes. We expect that with the growing number of known interface architectures, this type of knowledge-based methods will be increasingly used by the broad proteomics community.     

Sub‐angstrom modeling of complexes between flexible peptides and globular proteins

A wide range of regulatory processes in the cell are mediated by flexible peptides that fold upon binding to globular proteins. Computational efforts to model these interactions are hindered by the large number of rotatable bonds in flexible peptides relative to typical ligand molecules, and the fact that different peptides assume different backbone conformations within the same binding site. In this study, we present Rosetta FlexPepDock, a novel tool for refining coarse peptide–protein models that allows significant changes in both peptide backbone and side chains. We obtain high resolution models, often of sub‐angstrom backbone quality, over an extensive and general benchmark that is based on a large nonredundant dataset of 89 peptide–protein interactions. Importantly, side chains of known binding motifs are modeled particularly well, typically with atomic accuracy. In addition, our protocol has improved modeling quality for the important application of cross docking to PDZ domains. We anticipate that the ability to create high resolution models for a wide range of peptide–protein complexes will have significant impact on structure‐based functional characterization, controlled manipulation of peptide interactions, and on peptide‐based drug design. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Application of biasing‐potential replica‐exchange simulations for loop modeling and refinement of proteins in explicit solvent

Comparative or homology modeling of a target protein based on sequence similarity to a protein with known structure is widely used to provide structural models of proteins. Depending on the target‐template similarity these model structures may contain regions of limited structural accuracy. In principle, molecular dynamics (MD) simulations can be used to refine protein model structures and also to model loop regions that connect structurally conserved regions but it is limited by the currently accessible simulation time scales. A recently developed biasing potential replica exchange (BP‐REMD) method was used to refine loops and complete decoy protein structures at atomic resolution including explicit solvent. In standard REMD simulations several replicas of a system are run in parallel at different temperatures allowing exchanges at preset time intervals. In a BP‐REMD simulation replicas are controlled by various levels of a biasing potential to reduce the energy barriers associated with peptide backbone dihedral transitions. The method requires much fewer replicas for efficient sampling compared with T‐REMD. Application of the approach to several protein loops indicated improved conformational sampling of backbone dihedral angle of loop residues compared to conventional MD simulations. BP‐REMD refinement simulations on several test cases starting from decoy structures deviating significantly from the native structure resulted in final structures in much closer agreement with experiment compared to conventional MD simulations. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Supramolecular interactions between β‐lapachone with cyclodextrins studied using isothermal titration calorimetry and molecular modeling

Supramolecular interactions between β‐lapachone (β‐lap) and cyclodextrins (CDs) were investigated by isothermal titration calorimetry. The most favorable host: guest interaction was characterized using X‐ray powder diffraction (XRD), differential scanning calorimetry and thermogravimetry (DSC/TG), spectroscopy (FT‐IR), spectroscopy (2D ROESY) nuclear magnetic resonance (NMR), and molecular modeling. Phase solubility diagrams showed β‐, HP‐β‐, SBE‐β‐, γ‐, and HP‐γ‐CDs at 1.5% (w/w) allowed an increase in apparent solubility of β‐lap with enhancement factors of 12.0, 10.1, 11.8, 2.4, and 2.2, respectively. β‐lap has a weak interaction with γ‐ and HP‐γ‐CDs and tends to interact more favorably with β‐CD and its derivatives, especially SBE‐β‐CD (K = 4160 M⁻¹; ΔG = −20.66 kJ·mol⁻¹). Thermodynamic analysis suggests a hydrophobic interaction associated with the displacement of water from the cavity of the CD by the β‐lap. In addition, van der Waals forces and hydrogen bonds were responsible for the formation of complexes. Taken together, the results showed intermolecular interactions between β‐lap and SBE‐β‐CD, thereby confirming the formation of the inclusion complex. Molecular docking results showed 2 main orientations in which the interaction of benzene moiety at the wider rim of the SBE‐β‐CD is the most stable (average docking energy of −7.0 kcal/mol). In conclusion, β‐lap:SBE‐β‐CD is proposed as an approach for use in drug delivery systems in cancer research.     

Comprehensive modeling and functional analysis of Toll‐like receptor ligand‐recognition domains

Toll‐like receptors (TLRs) are innate immune pattern‐recognition receptors endowed with the capacity to detect microbial pathogens based on pathogen‐associated molecular patterns. The understanding of the molecular principles of ligand recognition by TLRs has been greatly accelerated by recent structural information, in particular the crystal structures of leucine‐rich repeat‐containing ectodomains of TLR2, 3, and 4 in complex with their cognate ligands. Unfortunately, for other family members such as TLR7, 8, and 9, no experimental structural information is currently available. Methods such as X‐ray crystallography or nuclear magnetic resonance are not applicable to all proteins. Homology modeling in combination with molecular dynamics may provide a straightforward yet powerful alternative to obtain structural information in the absence of experimental (structural) data, provided that the generated three‐dimensional models adequately approximate what is found in nature. Here, we report the development of modeling procedures tailored to the structural analysis of the extracellular domains of TLRs. We comprehensively compared secondary structure, torsion angles, accessibility for glycosylation, surface charge, and solvent accessibility between published crystal structures and independently built TLR2, 3, and 4 homology models. Finding that models and crystal structures were in good agreement, we extended our modeling approach to the remaining members of the TLR family from human and mouse, including TLR7, 8, and 9.     

Design of novel ligands of CDP‐methylerythritol kinase by mimicking direct protein‐protein and solvent‐mediated interactions

The methylerythritol 4‐phosphate (MEP) pathway for the biosynthesis of the isoprenoid universal building blocks (isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP)) is present in most of human pathogens and is absent in animals, turning it into a promising therapeutic druggable pathway. Two different strategies, a pharmacophore‐directed virtual screening and a protein‐protein interaction (PPI)‐mimicking cyclic peptide were used to search for compounds that bind to the PPI surface of the 4‐(cytidine 5‐diphospho)‐2C‐methyl‐D‐erythritol kinase (CMK), which catalyzes the fourth step of the MEP pathway. A significant part of the pharmacophore hypothesis used in this study was designed by mimicking water‐mediated PPI relevant in the CMK homodimer complex stabilization. After database search and with the aid of docking and molecular dynamics (MD) simulations, a 7H‐furo[3,2‐g]chromen‐7‐one derivative and a cyclic peptide were chosen as candidates to be ligands of CMK. Their binding affinities were measured using surface plasmon resonance (SPR) technology. Copyright © 2010 John Wiley & Sons, Ltd.     

Differential backbone dynamics of companion helices in the extended helical coiled‐coil domain of a bacterial chemoreceptor

Cytoplasmic domains of transmembrane bacterial chemoreceptors are largely extended four‐helix coiled coils. Previous observations suggested the domain was structurally dynamic. We probed directly backbone dynamics of this domain of the transmembrane chemoreceptor Tar from Escherichia coli using site‐directed spin labeling and electron paramagnetic resonance (EPR) spectroscopy. Spin labels were positioned on solvent‐exposed helical faces because EPR spectra for such positions reflect primarily polypeptide backbone movements. We acquired spectra for spin‐labeled, intact receptor homodimers solubilized in detergent or inserted into native E. coli lipid bilayers in Nanodiscs, characterizing 16 positions distributed throughout the cytoplasmic domain and on both helices of its helical hairpins, one amino terminal to the membrane‐distal tight turn (N‐helix), and the other carboxyl terminal (C‐helix). Detergent solubilization increased backbone dynamics for much of the domain, suggesting that loss of receptor activities upon solubilization reflects wide‐spread destabilization. For receptors in either condition, we observed an unanticipated difference between the N‐ and C‐helices. For bilayer‐inserted receptors, EPR spectra from sites in the membrane‐distal protein‐interaction region and throughout the C‐helix were typical of well‐structured helices. In contrast, for approximately two‐thirds of the N‐helix, from its origin as the AS‐2 helix of the membrane‐proximal HAMP domain to the beginning of the membrane‐distal protein‐interaction region, spectra had a significantly mobile component, estimated by spectral deconvolution to average approximately 15%. Differential helical dynamics suggests a four‐helix bundle organization with a pair of core scaffold helices and two more dynamic partner helices. This newly observed feature of chemoreceptor structure could be involved in receptor function.     

Crystal structure of the deglycating enzyme Amadoriase I in its free form and substrate‐bound complex

Amadoriases, also known as fructosyl amine oxidases (FAOX), are enzymes that catalyze the de‐glycosylation of fructosyl amino acids. As such, they are excellent candidates for the development of enzyme‐based diagnostic and therapeutic tools against age‐ and diabetes‐induced protein glycation. However, mostly because of the lack of a complete structural characterization of the different members of the family, the molecular bases of their substrate specificity have yet to be fully understood. The high resolution crystal structures of the free and the substrate‐bound form of Amadoriase I shown herein allow for the identification of key structural features that account for the diverse substrate specificity shown by this class of enzymes. This is of particular importance in the context of the rather limited and partially incomplete structural information that has so far been available in the literature on the members of the FAOX family. Moreover, using molecular dynamics simulations, we describe the tunnel conformation and the free energy profile experienced by the ligand in going from bulk water to the catalytic cavity, showing the presence of four gating helices/loops, followed by an “L‐shaped” narrow cavity. In summary, the tridimensional architecture of Amadoriase I presented herein provides a reference structural framework for the design of novel enzymes for diabetes monitoring and protein deglycation. Proteins 2016; 84:744–758. © 2016 Wiley Periodicals, Inc.