Lessons on protein structure from interleukin-4: All disulfides are not created equal.

Interleukin-4 (IL-4) is a hematopoietic cytokine composed by a four-helix bundle stabilized by an antiparallel beta-sheet and three disulfide bonds: Cys3-Cys127, Cys24-Cys65, and Cys46-Cys99. IL-4 is involved in several immune responses associated to infection, allergy, autoimmunity, and cancer. Besides its physiological relevance, IL-4 is often used as a "model" for protein design and engineering. Hence, to understand the role of each disulfide in the structure and dynamics of IL-4, we carried out several spectroscopic analyses (circular dichroism [CD], fluorescence, nuclear magnetic resonance [NMR]), and molecular dynamics (MD) simulations on wild-type IL-4 and four IL-4 disulfide mutants. All disulfide mutants showed loss of structure, altered interhelical angles, and looser core packings, showing that all disulfides are relevant for maintaining the overall fold and stability of the four-helix bundle motif, even at very low pH. In the absence of the disulfide connecting both protein termini Cys3-Cys127, C3T-IL4 showed a less packed protein core, loss of secondary structure (~9%) and fast motions on the sub-nanosecond time scale (lower S2 order parameters and larger τc correlation time), especially at the two protein termini, loops, beginning of helix A and end of helix D. In the absence of Cys24-Cys65, C24T-IL4 presented shorter alpha-helices (14% loss in helical content), altered interhelical angles, less propensity to form the small anti-parallel beta-sheet and increased dynamics. Simultaneously deprived of two disulfides (Cys3-Cys127 and Cys24-Cys65), IL-4 formed a partially folded "molten globule" with high 8-anilino-1-naphtalenesulphonic acid-binding affinity and considerable loss of secondary structure (~50%decrease), as shown by the far UV-CD, NMR, and MD data.

The structural basis for high affinity binding of α1-acid glycoprotein to the potent antitumor compound UCN-01.

The α1-acid glycoprotein (AGP) is an abundant blood plasma protein with important immunomodulatory functions coupled to endogenous and exogenous ligand-binding properties. Its affinity for many drug-like structures, however, means AGP can have a significant effect on the pharmokinetics and pharmacodynamics of numerous small molecule therapeutics. Staurosporine, and its hydroxylated forms UCN-01 and UCN-02, are kinase inhibitors that have been investigated at length as antitumour compounds. Despite their potency, these compounds display poor pharmokinetics due to binding to both AGP variants, AGP1 and AGP2. The recent renewed interest in UCN-01 as a cytostatic protective agent prompted us to solve the structure of the AGP2-UCN-01 complex by X-ray crystallography, revealing for the first time the precise binding mode of UCN-01. The solution NMR suggests AGP2 undergoes a significant conformational change upon ligand binding, but also that it uses a common set of sidechains with which it captures key groups of UCN-01 and other small molecule ligands. We anticipate that this structure and the supporting NMR data will facilitate rational redesign of small molecules that could evade AGP and therefore improve tissue distribution.

Quantitative High-Field NMR- and Mass Spectrometry-Based Fatty Acid Sequencing Reveals Internal Structure in Ru-Catalyzed Deuteration of Docosahexaenoic Acid.

Ru-based catalysis results in highly unsaturated fatty acid (HUFA) ethyl esters (EE) deuterated to various extents. The products carry 2H (D) mainly at their bis-allylic positions, where they are resistant to autoxidation compared to natural HUFA and are promising as neurological and retinal drugs. We characterized the extent of deuteration at each allylic position of docosa-4,7,10,13,16,19-hexaenoic acid deuterated to completion at bis-allylic and allylic positions (D-DHA) by two-dimensional (2D) and high-field (600 and 950 MHz) NMR. In separate experiments, the kinetics of docosahexaenoic acid (DHA) EE deuteration was evaluated using Paternò-Büchi (PB) reaction tandem mass spectrometry (MS/MS) analysis, enabling deuteration to be quantitatively characterized for isotopologues (D0-D14 DHA) at each internal allylic position. NMR analysis shows that the net deuteration of the isotopologue mixture is about 94% at the bis-allylic positions, and less than 1% remained as the protiated -CH2-. MS analysis shows that deuteration kinetics follow an increasing curve at bis-allylic positions with higher rate for internal bis-allylic positions. Percent D of bis-allylic positions increases linearly from D1 to D9 in which all internal bis-allylic positions (C9, C12, C15) deuterate uniformly and more rapidly than external bis-allylic positions (C6, C18). The mono-allylic positions near the methyl end (C21) show a steep increase of D only after the D10 isotopologue has been deuterated to >90%, while the mono-allylic position near the carboxyl position, C3, deuterates last and least. These data establish detailed methods for the characterization of Ru-catalyzed deuteration of HUFA as well as the phenomenological reaction kinetics as net product is formed.

An NMR and MD study of complexes of bacteriophage lambda lysozyme with tetra- and hexa-N-acetylchitohexaose.

The X-ray structure of lysozyme from bacteriophage lambda (λ lysozyme) in complex with the inhibitor hexa-N-acetylchitohexaose (NAG6) (PDB: 3D3D) has been reported previously showing sugar units from two molecules of NAG6 bound in the active site. One NAG6 is bound with four sugar units in the ABCD sites and the other with two sugar units in the E'F' sites potentially representing the cleavage reaction products; each NAG6 cross links two neighboring λ lysozyme molecules. Here we use NMR and MD simulations to study the interaction of λ lysozyme with the inhibitors NAG4 and NAG6 in solution. This allows us to study the interactions within the complex prior to cleavage of the polysaccharide. 1 HN and 15 N chemical shifts of λ lysozyme resonances were followed during NAG4/NAG6 titrations. The chemical shift changes were similar in the two titrations, consistent with sugars binding to the cleft between the upper and lower domains; the NMR data show no evidence for simultaneous binding of a NAG6 to two λ lysozyme molecules. Six 150 ns MD simulations of λ lysozyme in complex with NAG4 or NAG6 were performed starting from different conformations. The simulations with both NAG4 and NAG6 show stable binding of sugars across the D/E active site providing low energy models for the enzyme-inhibitor complexes. The MD simulations identify different binding subsites for the 5th and 6th sugars consistent with the NMR data. The structural information gained from the NMR experiments and MD simulations have been used to model the enzyme-peptidoglycan complex.

Exploitation of an iron transporter for bacterial protein antibiotic import.

Unlike their descendants, mitochondria and plastids, bacteria do not have dedicated protein import systems. However, paradoxically, import of protein bacteriocins, the mechanisms of which are poorly understood, underpins competition among pathogenic and commensal bacteria alike. Here, using X-ray crystallography, isothermal titration calorimetry, confocal fluorescence microscopy, and in vivo photoactivatable cross-linking of stalled translocation intermediates, we demonstrate how the iron transporter FpvAI in the opportunistic pathogen Pseudomonas aeruginosa is hijacked to translocate the bacteriocin pyocin S2 (pyoS2) across the outer membrane (OM). FpvAI is a TonB-dependent transporter (TBDT) that actively imports the small siderophore ferripyoverdine (Fe-Pvd) by coupling to the proton motive force (PMF) via the inner membrane (IM) protein TonB1. The crystal structure of the N-terminal domain of pyoS2 (pyoS2NTD) bound to FpvAI (Kd = 240 pM) reveals that the pyocin mimics Fe-Pvd, inducing the same conformational changes in the receptor. Mimicry leads to fluorescently labeled pyoS2NTD being imported into FpvAI-expressing P. aeruginosa cells by a process analogous to that used by bona fide TBDT ligands. PyoS2NTD induces unfolding by TonB1 of a force-labile portion of the plug domain that normally occludes the central channel of FpvAI. The pyocin is then dragged through this narrow channel following delivery of its own TonB1-binding epitope to the periplasm. Hence, energized nutrient transporters in bacteria also serve as rudimentary protein import systems, which, in the case of FpvAI, results in a protein antibiotic 60-fold bigger than the transporter's natural substrate being translocated across the OM.

The CcmC-CcmE interaction during cytochrome c maturation by System I is driven by protein-protein and not protein-heme contacts.

Cytochromes c are ubiquitous proteins, essential for life in most organisms. Their distinctive characteristic is the covalent attachment of heme to their polypeptide chain. This post-translational modification is performed by a dedicated protein system, which in many Gram-negative bacteria and plant mitochondria is a nine-protein apparatus (CcmA-I) called System I. Despite decades of study, mechanistic understanding of the protein-protein interactions in this highly complex maturation machinery is still lacking. Here, we focused on the interaction of CcmC, the protein that sources the heme cofactor, with CcmE, the pivotal component of System I responsible for the transfer of the heme to the apocytochrome. Using in silico analyses, we identified a putative interaction site between these two proteins (residues Asp47, Gln50, and Arg55 on CcmC; Arg73, Asp101, and Glu105 on CcmE), and we validated our findings by in vivo experiments in Escherichia coli Moreover, employing NMR spectroscopy, we examined whether a heme-binding site on CcmE contributes to this interaction and found that CcmC and CcmE associate via protein-protein rather than protein-heme contacts. The combination of in vivo site-directed mutagenesis studies and high-resolution structural techniques enabled us to determine at the residue level the mechanism for the formation of one of the key protein complexes for cytochrome c maturation by System I.

Structural Insights into Notch Receptor-Ligand Interactions.

Pioneering cell aggregation experiments from the Artavanis-Tsakonas group in the late 1980's localized the core ligand recognition sequence in the Drosophila Notch receptor to epidermal growth factor-like (EGF) domains 11 and 12. Since then, advances in protein expression, structure determination methods and functional assays have enabled us to define the molecular basis of the core receptor/ligand interaction and given new insights into the architecture of the Notch complex at the cell surface. We now know that Notch EGF11 and 12 interact with the Delta/Serrate/LAG-2 (DSL) and C2 domains of ligand and that membrane-binding, together with additional protein-protein interactions outside the core recognition domains, are likely to fine-tune generation of the Notch signal. Furthermore, structure determination of O-glycosylated variants of Notch alone or in complex with receptor fragments, has shown that these sugars contribute directly to the binding interface, as well as to stabilizing intra-molecular domain structure, providing some mechanistic insights into the observed modulatory effects of O-glycosylation on Notch activity.Future challenges lie in determining the complete extracellular architecture of ligand and receptor in order to understand (i) how Notch/ligand complexes may form at the cell surface in response to physiological cues, (ii) the role of lipid binding in stabilizing the Notch/ligand complex, (iii) the impact of O-glycosylation on binding and signalling and (iv) to dissect the different pathologies that arise as a consequence of mutations that affect proteins involved in the Notch pathway.

The Role of Active Site Flexible Loops in Catalysis and of Zinc in Conformational Stability of Bacillus cereus 569/H/9 ß-Lactamase.

Metallo-ß-lactamases catalyze the hydrolysis of most ß-lactam antibiotics and hence represent a major clinical concern. The development of inhibitors for these enzymes is complicated by the diversity and flexibility of their substrate-binding sites, motivating research into their structure and function. In this study, we examined the conformational properties of the Bacillus cereus ß-lactamase II in the presence of chemical denaturants using a variety of biochemical and biophysical techniques. The apoenzyme was found to unfold cooperatively, with a Gibbs free energy of stabilization (ΔG(0)) of 32 ± 2 kJ·mol(-1) For holoBcII, a first non-cooperative transition leads to multiple interconverting native-like states, in which both zinc atoms remain bound in an apparently unaltered active site, and the protein displays a well organized compact hydrophobic core with structural changes confined to the enzyme surface, but with no catalytic activity. Two-dimensional NMR data revealed that the loss of activity occurs concomitantly with perturbations in two loops that border the enzyme active site. A second cooperative transition, corresponding to global unfolding, is observed at higher denaturant concentrations, with ΔG(0) value of 65 ± 1.4 kJ·mol(-1) These combined data highlight the importance of the two zinc ions in maintaining structure as well as a relatively well defined conformation for both active site loops to maintain enzymatic activity.

An extended active-site motif controls the reactivity of the thioredoxin fold.

Proteins belonging to the thioredoxin (Trx) superfamily are abundant in all organisms. They share the same structural features, arranged in a seemingly simple fold, but they perform a multitude of functions in oxidative protein folding and electron transfer pathways. We use the C-terminal domain of the unique transmembrane reductant conductor DsbD as a model for an in-depth analysis of the factors controlling the reactivity of the Trx fold. We employ NMR spectroscopy, x-ray crystallography, mutagenesis, in vivo functional experiments applied to DsbD, and a comparative sequence analysis of Trx-fold proteins to determine the effect of residues in the vicinity of the active site on the ionization of the key nucleophilic cysteine of the -CXXC- motif. We show that the function and reactivity of Trx-fold proteins depend critically on the electrostatic features imposed by an extended active-site motif.

Structure of the terminal PCP domain of the non-ribosomal peptide synthetase in teicoplanin biosynthesis.

The biosynthesis of the glycopeptide antibiotics, of which teicoplanin and vancomycin are representative members, relies on the combination of non-ribosomal peptide synthesis and modification of the peptide by cytochrome P450 (Oxy) enzymes while the peptide remains bound to the peptide synthesis machinery. We have structurally characterized the final peptidyl carrier protein domain of the teicoplanin non-ribosomal peptide synthetase machinery: this domain is believed to mediate the interactions with tailoring Oxy enzymes in addition to its function as a shuttle for intermediates between multiple non-ribosomal peptide synthetase domains. Using solution state NMR, we have determined structures of this PCP domain in two states, the apo and the post-translationally modified holo state, both of which conform to a four-helix bundle assembly. The structures exhibit the same general fold as the majority of known carrier protein structures, in spite of the complex biosynthetic role that PCP domains from the final non-ribosomal peptide synthetase module must play in glycopeptide antibiotic biosynthesis. These structures thus support the hypothesis that it is subtle rearrangements, rather than dramatic conformational changes, which govern carrier protein interactions and selectivity during non-ribosomal peptide synthesis.

Consequences of inducing intrinsic disorder in a high-affinity protein-protein interaction.

The kinetic and thermodynamic consequences of intrinsic disorder in protein-protein recognition are controversial. We address this by inducing one partner of the high-affinity colicin E3 rRNase domain-Im3 complex (K(d) ≈ 10(-12) M) to become an intrinsically disordered protein (IDP). Through a variety of biophysical measurements, we show that a single alanine mutation at Tyr507 within the hydrophobic core of the isolated colicin E3 rRNase domain causes the enzyme to become an IDP (E3 rRNase(IDP)). E3 rRNase(IDP) binds stoichiometrically to Im3 and forms a structure that is essentially identical to the wild-type complex. However, binding of E3 rRNase(IDP) to Im3 is 4 orders of magnitude weaker than that of the folded rRNase, with thermodynamic parameters reflecting the disorder-to-order transition on forming the complex. Critically, pre-steady-state kinetic analysis of the E3 rRNase(IDP)-Im3 complex demonstrates that the decrease in affinity is mostly accounted for by a drop in the electrostatically steered association rate. Our study shows that, notwithstanding the advantages intrinsic disorder brings to biological systems, this can come at severe kinetic and thermodynamic cost.

Comparison of the backbone dynamics of wild-type Hydrogenobacter thermophilus cytochrome c(552) and its b-type variant.

Cytochrome c552 from the thermophilic bacterium Hydrogenobacter thermophilus is a typical c-type cytochrome which binds heme covalently via two thioether bonds between the two heme vinyl groups and two cysteine thiol groups in a CXXCH sequence motif. This protein was converted to a b-type cytochrome by substitution of the two cysteine residues by alanines (Tomlinson and Ferguson in Proc Natl Acad Sci USA 97:5156-5160, 2000a). To probe the significance of the covalent attachment of the heme in the c-type protein, (15)N relaxation and hydrogen exchange studies have been performed for the wild-type and b-type proteins. The two variants share very similar backbone dynamic properties, both proteins showing high (15)N order parameters in the four main helices, with reduced values in an exposed loop region (residues 18-21), and at the C-terminal residue Lys80. Some subtle changes in chemical shift and hydrogen exchange protection are seen between the wild-type and b-type variant proteins, not only for residues at and neighbouring the mutation sites, but also for some residues in the heme binding pocket. Overall, the results suggest that the main role of the covalent linkages between the heme group and the protein chain must be to increase the stability of the protein.

Molecular basis for Jagged-1/Serrate ligand recognition by the Notch receptor.

We have mapped a Jagged/Serrate-binding site to specific residues within the 12th EGF domain of human and Drosophila Notch. Two critical residues, involved in a hydrophobic interaction, provide a ligand-binding platform and are adjacent to a Fringe-sensitive residue that modulates Notch activity. Our data suggest that small variations within the binding site fine-tune ligand specificity, which may explain the observed sequence heterogeneity in mammalian Notch paralogues, and should allow the development of paralogue-specific ligand-blocking antibodies. As a proof of principle, we have generated a Notch-1-specific monoclonal antibody that blocks binding, thus paving the way for antibody tools for research and therapeutic applications.

NMR study of the structure and dynamics of the BRCT domain from the kinetochore protein KKT4.

KKT4 is a multi-domain kinetochore protein specific to kinetoplastids, such as Trypanosoma brucei. It lacks significant sequence similarity to known kinetochore proteins in other eukaryotes. Our recent X-ray structure of the C-terminal region of KKT4 shows that it has a tandem BRCT (BRCA1 C Terminus) domain fold with a sulfate ion bound in a typical binding site for a phosphorylated serine or threonine. Here we present the 1H, 13C and 15N resonance assignments for the BRCT domain of KKT4 (KKT4463-645) from T. brucei. We show that the BRCT domain can bind phosphate ions in solution using residues involved in sulfate ion binding in the X-ray structure. We have used these assignments to characterise the secondary structure and backbone dynamics of the BRCT domain in solution. Mutating the residues involved in phosphate ion binding in T. brucei KKT4 BRCT results in growth defects confirming the importance of the BRCT phosphopeptide-binding activity in vivo. These results may facilitate rational drug design efforts in the future to combat diseases caused by kinetoplastid parasites.

The dynamics of lysozyme from bacteriophage lambda in solution probed by NMR and MD simulations.

(15) N NMR relaxation studies, analyses of NMR data to include chemical shifts, residual dipolar couplings (RDC), NOEs and H(N) -H(α) coupling constants, and molecular dynamics (MD) simulations have been used to characterise the behaviour of lysozyme from bacteriophage lambda (λ lysozyme) in solution. The lower and upper lip regions in λ lysozyme (residues 51-60 and 128-141, respectively) show reduced (1) H-(15) N order parameters indicating mobility on a picosecond timescale. In addition, residues in the lower and upper lips also show exchange contributions to T2 indicative of slower timescale motions. The chemical shift, RDC, coupling constant and NOE data for λ lysozyme indicate that two fluctuating ß-strands (ß3 and ß4) are populated in the lower lip region while the N terminus of helix α6 (residues 136-139) forms dynamic helical turns in the upper lip region. This behaviour is confirmed by MD simulations that show hydrogen bonds, indicative of the ß-sheet and helical secondary structure in the lip regions, with populations of 40-60 %. Thus in solution λ lysozyme adopts a conformational ensemble that will contain both the open and closed forms observed in the crystal structures of the protein.

Water mediation is essential to nucleation of ß-turn formation in peptide folding motifs.

Water-mediated bond formation: The structure of the peptide GPG-NH2 has been investigated in aqueous solution to understand the role of water in the formation of a ß-turn. Using a combination of neutron diffraction enhanced by isotopic substitution, NMR spectroscopy, and computer simulations, it was found that water is an essential component to initiate folding in solution.

Ion binding with charge inversion combined with screening modulates DEAD box helicase phase transitions.

Membraneless organelles, or biomolecular condensates, enable cells to compartmentalize material and processes into unique biochemical environments. While specific, attractive molecular interactions are known to stabilize biomolecular condensates, repulsive interactions, and the balance between these opposing forces, are largely unexplored. Here, we demonstrate that repulsive and attractive electrostatic interactions regulate condensate stability, internal mobility, interfaces, and selective partitioning of molecules both in vitro and in cells. We find that signaling ions, such as calcium, alter repulsions between model Ddx3 and Ddx4 condensate proteins by directly binding to negatively charged amino acid sidechains and effectively inverting their charge, in a manner fundamentally dissimilar to electrostatic screening. Using a polymerization model combined with generalized stickers and spacers, we accurately quantify and predict condensate stability over a wide range of pH, salt concentrations, and amino acid sequences. Our model provides a general quantitative treatment for understanding how charge and ions reversibly control condensate stability.

Oxidation state-dependent protein-protein interactions in disulfide cascades.

Bacterial growth and pathogenicity depend on the correct formation of disulfide bonds, a process controlled by the Dsb system in the periplasm of Gram-negative bacteria. Proteins with a thioredoxin fold play a central role in this process. A general feature of thiol-disulfide exchange reactions is the need to avoid a long lived product complex between protein partners. We use a multidisciplinary approach, involving NMR, x-ray crystallography, surface plasmon resonance, mutagenesis, and in vivo experiments, to investigate the interaction between the two soluble domains of the transmembrane reductant conductor DsbD. Our results show oxidation state-dependent affinities between these two domains. These observations have implications for the interactions of the ubiquitous thioredoxin-like proteins with their substrates, provide insight into the key role played by a unique redox partner with an immunoglobulin fold, and are of general importance for oxidative protein-folding pathways in all organisms.

Residual dipolar couplings: are multiple independent alignments always possible?

RDCs for the 14 kDa protein hen egg-white lysozyme (HEWL) have been measured in eight different alignment media. The elongated shape and strongly positively charged surface of HEWL appear to limit the protein to four main alignment orientations. Furthermore, low levels of alignment and the protein's interaction with some alignment media increases the experimental error. Together with heterogeneity across the alignment media arising from constraints on temperature, pH and ionic strength for some alignment media, these data are suitable for structure refinement, but not the extraction of dynamic parameters. For an analysis of protein dynamics the data must be obtained with very low errors in at least three or five independent alignment media (depending on the method used) and so far, such data have only been reported for three small 6-8 kDa proteins with identical folds: ubiquitin, GB1 and GB3. Our results suggest that HEWL is likely to be representative of many other medium to large sized proteins commonly studied by solution NMR. Comparisons with over 60 high-resolution crystal structures of HEWL reveal that the highest resolution structures are not necessarily always the best models for the protein structure in solution.

Structure and interdomain interactions of a hybrid domain: a disulphide-rich module of the fibrillin/LTBP superfamily of matrix proteins.

The fibrillins and latent transforming growth factor-beta binding proteins (LTBPs) form a superfamily of structurally-related proteins consisting of calcium-binding epidermal growth factor-like (cbEGF) domains interspersed with 8-cysteine-containing transforming growth factor beta-binding protein-like (TB) and hybrid (hyb) domains. Fibrillins are the major components of the extracellular 10-12 nm diameter microfibrils, which mediate a variety of cell-matrix interactions. Here we present the crystal structure of a fibrillin-1 cbEGF9-hyb2-cbEGF10 fragment, solved to 1.8 A resolution. The hybrid domain fold is similar, but not identical, to the TB domain fold seen in previous fibrillin-1 and LTBP-1 fragments. Pairwise interactions with neighboring cbEGF domains demonstrate extensive interfaces, with the hyb2-cbEGF10 interface dependent on Ca(2+) binding. These observations provide accurate constraints for models of fibrillin organization within the 10-12 nm microfibrils and provide further molecular insights into how Ca(2+) binding influences the intermolecular interactions and biomechanical properties of fibrillin-1.