Inhibiting Type VI Secretion System Activity with a Biomimetic Peptide Designed To Target the Baseplate Wedge Complex.

Human health is threatened by bacterial infections that are increasingly resistant to multiple drugs. A recently emerged strategy consists of disarming pathogenic bacteria by targeting and blocking their virulence factors. The type VI secretion system (T6SS) is a widespread secretion nanomachine encoded and employed by pathogenic strains to establish their virulence process during host invasion. Given the conservation of T6SS in several human bacterial pathogens, the discovery of an effective broad-spectrum T6SS virulence blocker represents an attractive target for development of antivirulence therapies. Here, we identified and validated a protein-protein interaction interface, TssK-TssG, as a key factor in the assembly of the T6SS baseplate (BP) complex in the pathogen enteroaggregative Escherichia coli (EAEC). In silico and biochemical studies revealed that the determinants of the interface are broadly conserved among pathogenic species, suggesting a role for this interface as a target for T6SS inhibition. Based on the high-resolution structure of the TssKFGE wedge complex, we rationally designed a biomimetic cyclic peptide (BCP) that blocks the assembly of the EAEC BP complex and inhibits the function of T6SS in bacterial cultures. Our BCP is the first compound completely designed from prior structural knowledge with anti-T6SS activity that can be used as a model to target human pathogens. IMPORTANCE New therapeutic options are urgently needed to fight drug-resistant and life-threatening infections. In contrast to antibiotics that inhibit the growth pathways of bacteria, the antivirulence strategy is a promising approach to disarm pathogens by interfering with bacterial virulence factors without exerting evolutionary pressure. The type VI secretion system (T6SS) is used by many pathogens, including members of the antibiotic-resistant ESKAPE bacteria (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.), to establish their virulence during the invasion of the human host. Although the T6SS is undoubtedly involved in pathogenesis, strategies targeting this virulence factor are crucially lacking. Here, we used a combination of genetics, microbiology, biochemical, biophysics, and bioinformatics approaches to rationally design a biomimetic peptide that interferes with T6SS assembly and functioning. This study represents a novel proof of concept for an antivirulence strategy which aims to interfere with the assembly of the T6SS.

A structure refinement protocol combining NMR residual dipolar couplings and small angle scattering restraints.

We present the implementation of a target function based on Small Angle Scattering data (Gabel et al. Eur Biophys J 35(4):313-327, 2006) into the Crystallography and NMR Systems (CNS) and demonstrate its utility in NMR structure calculations by simultaneous application of small angle scattering (SAS) and residual dipolar coupling (RDC) restraints. The efficiency and stability of the approach are demonstrated by reconstructing the structure of a two domain region of the 31 kDa nuclear export factor TAP (TIP-associated protein). Starting with the high resolution X-ray structures of the two individual TAP domains, the translational and orientational domain arrangement is refined simultaneously. We tested the stability of the protocol against variations of the SAS target parameters and the number of RDCs and their uncertainties. The activation of SAS restraints results in an improved translational clustering of the domain positions and lifts part of the fourfold degeneracy of their orientations (associated with a single alignment tensor). The resulting ensemble of structures reflects the conformational space that is consistent with the experimental SAS and RDC data. The SAS target function is computationally very efficient. SAS restraints can be activated at different levels of precision and only a limited SAS angular range is required. When combined with additional data from chemical shift perturbation, paramagnetic relaxation enhancement or mutational analysis the SAS refinement is an efficient approach for defining the topology of multi-domain and/or multimeric biomolecular complexes in solution based on available high resolution structures (NMR or X-ray) of the individual domains.

NMR in the SPINE Structural Proteomics project.

This paper describes the developments, role and contributions of the NMR spectroscopy groups in the Structural Proteomics In Europe (SPINE) consortium. Focusing on the development of high-throughput (HTP) pipelines for NMR structure determinations of proteins, all aspects from sample preparation, data acquisition, data processing, data analysis to structure determination have been improved with respect to sensitivity, automation, speed, robustness and validation. Specific highlights are protonless (13)C-direct detection methods and inferential structure determinations (ISD). In addition to technological improvements, these methods have been applied to deliver over 60 NMR structures of proteins, among which are five that failed to crystallize. The inclusion of NMR spectroscopy in structural proteomics pipelines improves the success rate for protein structure determinations.

Clinical and mutational investigations of tyrosinemia type II in Northern Tunisia: identification and structural characterization of two novel TAT mutations.

Tyrosinemia type II or Richner-Hanhart Syndrome (RHS) is an autosomal recessive disorder characterized by keratitis, palmoplantar keratosis, mental retardation, and elevated blood tyrosine levels. The disease is due to a deficiency of hepatic cytosolic tyrosine aminotransferase (TATc), an enzyme involved in the tyrosine catabolic pathway. Because of the high rate of consanguinity this disorder seems to be relatively common among the Arab and Mediterranean populations. RHS is characterized by inter and intrafamilial phenotypic variability. A large spectrum of mutations within TATc gene has been shown to be responsible for RHS. In the present study, we report the clinical features and the molecular investigation of RHS in three unrelated consanguineous Tunisian families including 7 patients with confirmed biochemical diagnosis of tyrosinemia type II. Mutation analyses were performed and two novel missense mutations were identified (C151Y) and (L273P) within exon 5 and exon 8, respectively. The 3D-structural characterization of these mutations provides evidence of defective folding of the mutant proteins, and likely alteration of the enzymatic activity. Phenotype variability was observed even among individuals sharing the same pathogenic mutation.

NMR studies of the sporulation protein SpoIIAA: implications for the regulation of the transcription factor sigmaF in Bacillus subtilis.

SpoIIAA participates in a four-component mechanism for phosphorylation-dependent transcription control at the outset of sporulation. We report the refinement of the solution structure of SpoIIAA by using the automated iterative NOE assignment method ARIA. To complement the structural data, the protein dynamics were determined by measuring the T1, T2 and NOE of the backbone 15N-nuclei. The refined structure permits a discussion of the structural features that are important for the function of SpoIIAA in the regulation of the sporulation sigma factor sigmaF, and for homologous regulatory pathways present in B. subtilis and in other bacilli.

Continuum solvent molecular dynamics study of flexibility in interleukin-8.

Generalized Born continuum solvent methods have been shown to provide a reasonable description of the equilibrium thermodynamics of aqueous solvation in a variety of applications to peptides, proteins, and nucleic acids. Here we study the performance of these methods in molecular dynamics simulations of interleukin-8, comparing nanosecond-length explicit solvent simulations with those using the generalized Born model. In general, the simulations show similar results, although movement away from the initial NMR-determined structure and average fluctuations about the mean are slightly higher for the continuum solvent results. In both simulations, the two helices that are packed on top of the core sheet move closer together, resulting in a structure that more closely resembles the X-ray structure. Principal-component (quasiharmonic) analysis is used to analyze the motions of these helices in both of the simulations and in the NMR ensemble of structures. Prospects for making more general use of continuum solvent models in protein dynamics simulations are discussed.

Structural and functional studies of titin's fn3 modules reveal conserved surface patterns and binding to myosin S1--a possible role in the Frank-Starling mechanism of the heart.

The A-band part of titin, a striated-muscle specific protein spanning from the Z-line to the M-line, mainly consists of a well-ordered super-repeat array of immunoglobulin-like and fibronectin-type III (fn3)-like domains. Since it has been suspected that the fn3 domains might represent titin's binding sites to myosin, we have developed structural models for all of titin's 132 fn3-like domains. A subset of eight experimentally determined fn3 structures from a range of proteins, including titin itself, was used as homology templates. After grouping the models according to their position within the super-repeat segment of the central A-band titin region, we analyzed the models with respect to side-chain conservation. This showed that conserved residues form an extensive surface pattern predominantly at one side of the domains, whereas domains outside the central C-zone super-repeat region show generally less conserved surfaces. Since the conserved surface residues may function as protein-binding sites, we experimentally studied the binding properties of expressed multi-domain fn3 fragments. This revealed that fn3 fragments specifically bind to the sub-fragment 1 of myosin. We also measured the effect of fn3 fragments on the contractile properties of single cardiac myocytes. At sub-maximal Ca(2+) concentrations, fn3 fragments significantly enhance active tension. This effect is most pronounced at short sarcomere length, and as a result the length-dependence of Ca(2+) activation is reduced. A model of how titin's fn3-like domains may influence actomyosin interaction is proposed.

Re-face stereospecificity of methylenetetrahydromethanopterin and methylenetetrahydrofolate dehydrogenases is predetermined by intrinsic properties of the substrate.

Four different dehydrogenases are known that catalyse the reversible dehydrogenation of N5,N10-methylenetetrahydromethanopterin (methylene-H4MPT) or N5,N10-methylenetetrahydrofolate (methylene-H4F) to the respective N5,N10-methenyl compounds. Sequence comparison indicates that the four enzymes are phylogenetically unrelated. They all catalyse the Re-face-stereospecific removal of the pro-R hydrogen atom of the coenzyme's methylene group. The Re-face stereospecificity is in contrast to the finding that in solution the pro-S hydrogen atom of methylene-H4MPT and of methylene-H4F is more reactive to heterolytic cleavage. For a better understanding we determined the conformations of methylene-H4MPT in solution and when enzyme-bound by using NMR spectroscopy and semiempirical quantum mechanical calculations. For the conformation free in solution we find an envelope conformation for the imidazolidine ring, with the flap at N10. The methylene pro-S C-H bond is anticlinal and the methylene pro-R C-H bond is synclinal to the lone electron pair of N10. Semiempirical quantum mechanical calculations of heats of formation of methylene-H4MPT and methylene-H4F indicate that changing this conformation into an activated one in which the pro-S C-H bond is antiperiplanar, resulting in the preformation of the leaving hydride, would require a deltadeltaH(f) of +53 kJ mol-1 for methylene-H4MPT and of +51 kJ mol-1 for methylene-H4F. This is almost twice the energy required to force the imidazolidine ring in the enzyme-bound conformation of methylene-H4MPT (+29 kJ mol-1) or of methylene-H4F (+35 kJ mol-1) into an activated conformation in which the pro-R hydrogen atom is antiperiplanar to the lone electron pair of N10. The much lower energy for pro-R hydrogen activation thus probably predetermines the Re-face stereospecificity of the four dehydrogenases. Results are also presented explaining why the chemical reduction of methenyl-H4MPT+ and methenyl-H4F+ with NaBD4 proceeds Si-face-specific, in contrast to the enzyme-catalysed reaction.

Structure of a PH domain from the C. elegans muscle protein UNC-89 suggests a novel function.

BACKGROUND: Pleckstrin homology (PH) domains constitute a structurally conserved family present in many signaling and regulatory proteins. PH domains have been shown to bind to phospholipids, and many function in membrane targeting. They generally have a strong electrostatic polarization and interact with negatively charged phospholipids via the positive pole. On the basis of electrostatic modeling, however, we have previously identified a class of PH domains with a predominantly negative charge and predicted that these domains recognize other targets. Here, we report the first experimental structure of such a PH domain. RESULTS: The structure of the PH domain from Caenorhabditis elegans muscle protein UNC-89 has been determined by heteronuclear NMR. The domain adopts the classic PH fold, but has an unusual closed conformation of the "inositol binding loops. This creates a small opening to a deep hydrophobic pocket lined with negative charges on one side, and provides a molecular explanation for the lack of association with inositol-1,4,5-triphosphate. As predicted, the PH domain of UNC-89 has a strongly negative overall electrostatic potential. Modeling the Dbl homology (DH)-linked PH domains from the C. elegans genome shows that a large proportion of these modules are negatively charged. CONCLUSIONS: We present the first structure of a PH domain with a strong negative overall electrostatic potential. The presence of a deep pocket lined with negative charges suggests that the domain binds to ligands other than acidic phospholipids. The abundance of this class of PH domain in the C. elegans genome suggests a prominent role in mediating protein-protein interactions.

A new approach for applying residual dipolar couplings as restraints in structure elucidation.

Residual dipolar couplings are useful global structural restraints. The dipolar couplings define the orientation of a vector with respect to the alignment tensor. Although the size of the alignment tensor can be derived from the distribution of the experimental dipolar couplings, its orientation with respect to the coordinate system of the molecule is unknown at the beginning of structure determination. This causes convergence problems in the simulated annealing process. We therefore propose a protocol that translates dipolar couplings into intervector projection angles, which are independent of the orientation of the alignment tensor with respect to the molecule. These restraints can be used during the whole simulated annealing protocol.

Efficient sampling in collective coordinate space.

Collective motions in biological macromolecules have been shown to be important for function. The most important collective motions occur on slow time scales, which poses a sampling problem in dynamic simulation of biomolecules. We present a novel method for efficient conformational sampling. The method combines the simulation of an ensemble of concurrent trajectories with restraints acting on the ensemble of structures as a whole. Two properties of the ensemble may be restrained: (i) the variance of the ensemble and (ii) the average position of the ensemble. Both properties are defined in a subspace of collective coordinate space spanned by an arbitrary number of modes. We show that weak restraints on the ensemble variance suffice for an increase in sampling efficiency along soft modes by two orders of magnitudes. The resulting trajectories exhibit virtually the same structural quality as trajectories generated by restraint-free-molecular dynamics simulation, as judged by standard structure validation tools. The method is used to probe the resistance of a structure against conformational changes along collective modes and clearly distinguishes soft from stiff modes. Further applications are discussed. Proteins 2000;39:82-88.

Refinement of the protein backbone angle psi in NMR structure calculations.

Cross-correlated relaxation rates involving the Calpha-Halpha dipolar interaction and the carbonyl (C') chemical shift anisotropy (CSA) have been measured using two complementary 3D experiments. We show that the protein backbone angle psi can be directly refined against such cross-correlated relaxation rates (gammaHalphaCalpha,C') and the three-bond H/D isotope effect on the Calpha chemical shifts (3 deltaCalpha(ND)). By simultaneously using both experimental parameters as restraints during NMR structure calculations, a unique value for the backbone angle psi is defined. We have applied the new refinement method to the alpha-Spectrin SH3 domain (a beta-sheet protein) and to the Sgs1p HRDC domain (an alpha-helical protein) and show that the quality of the NMR structures is substantially improved, judging from the atomic coordinate precision and the Ramachandran map. In addition, the psi-refined NMR structures of the SH3 domain deviate less from the 1.8 A crystal structure, suggesting an improved accuracy. The proposed refinement method can be used to significantly improve the quality of NMR structures and will be applicable to larger proteins.

Unraveling the symmetry ambiguity in a hexamer: calculation of the R6 human insulin structure.

Crystallographic and NMR studies of insulin have revealed a highly flexible molecule with a range of different aggregation and structural states; the importance of these states for the function of the hormone is still unclear. To address this question, we have studied the solution structure of the insulin R6 symmetric hexamer using NMR spectroscopy. Structure determination of symmetric oligomers by NMR is complicated due to 'symmetry ambiguity' between intra- and intermonomer NOEs, and between different classes of intermonomer NOEs. Hence, to date, only two symmetric tetramers and one symmetric pentamer (VTB, B subunit of verotoxin) have been solved by NMR: there has been no other symmetric hexamer or higher-order oligomer. Recently, we reported a solution structure for R6 insulin hexamer. However, in that study, a crystal structure was used as a reference to resolve ambiguities caused by the threefold symmetry; the same method was used in solving VTB. Here, we have successfully recalculated R6 insulin using the symmetry-ADR method, a computational strategy in which ambiguities are resolved using the NMR data alone. Thus the obtained structure is a refinement of the previous R6 solution structure. Correlated motions in the final structural ensemble were analysed using a recently developed principal component method; this suggests the presence of two major conformational substates. The study demonstrates that the solution structure of higher-order symmetric oligomers can be determined unambiguously from NMR data alone, using the symmetry-ADR method. This success bodes well for future NMR studies of higher-order symmetric oligomers. The correlated motions observed in the structural ensemble suggest a new insight into the mechanism of phenol exchange and the T6 <--> R6 transition of insulin in solution.

Classification of protein sequences by homology modeling and quantitative analysis of electrostatic similarity.

Protein electrostatics plays a key role in ligand binding and protein-protein interactions. Therefore, similarities or dissimilarities in electrostatic potentials can be used as indicators of similarities or dissimilarities in protein function. We here describe a method to compare the electrostatic properties within protein families objectively and quantitatively. Three-dimensional structures are built from database sequences by comparative modeling. Molecular potentials are then computed for these with a continuum solvation model by finite difference solution of the Poisson-Boltzmann equation or analytically as a multipole expansion that permits rapid comparison of very large datasets. This approach is applied to 104 members of the Pleckstrin homology (PH) domain family. The deviation of the potentials of the homology models from those of the corresponding experimental structures is comparable to the variation of the potential in an ensemble of structures from nuclear magnetic resonance data or between snapshots from a molecular dynamics simulation. For this dataset, the results for analysis of the full electrostatic potential and the analysis using only monopole and dipole terms are very similar. The electrostatic properties of the PH domains are generally conserved despite the extreme sequence divergence in this family. Notable exceptions from this conservation are seen for PH domains linked to a Db1 homology (DH) domain and in proteins with internal PH domain repeats.

The PH superfold: a structural scaffold for multiple functions.

Pleckstrin homology (PH) domains form a structurally conserved family that is associated with many regulatory pathways within the cell. Domains with a nearly identical fold are found in other families that share no sequence similarity, suggesting the existence of a stable PH superfold. The PH domains generally function as regulated membrane-binding modules that bind to inositol lipids and respond to upstream signals by targeting the host proteins to the correct cellular sites. The other domains with a similar fold, such as the phosphotyrosine binding domains, recognize protein ligands.

The structure in solution of the b domain of protein disulfide isomerase.

Protein disulfide isomerase (PDI) is a multifunctional protein of the endoplasmic reticulum, which catalyzes the formation, breakage and rearrangement of disulfide bonds during protein folding. It consists of four domains designated a, b, b and a. Both a and a domains contains an active site with the sequence motif -Cys-Gly-His-Cys-involved directly in thiol-disulfide exchange reactions. As expected these domains have structures very similar to the ubiquitous redox protein thioredoxin. A low-resolution NMR structure of the b domain revealed that this domain adopts a fold similar to the PDI a domain and thioredoxin [Kemmink, J., Darby, N.J., Dijkstra, K., Nilges, M. and Creighton, T.E. (1997) Curr. Biol. 7, 239-245]. A refined ensemble of solution structures based on the input of 1865 structural restraints shows that the structure of PDI b is well defined throughout the complete protein except for about 10 residues at the C-terminus of the sequence. 15N relaxation data show that these residues are disordered and not part of this structural domain. Therefore the domain boundaries of PDI can now be fixed with reasonable precision. Structural comparison of the PDI b domain with thioredoxin and PDI a reveals several features important for thiol-disulfide exchange activity.

Functionally important correlated motions in the single-stranded DNA-binding protein encoded by filamentous phage Pf3.

To elucidate the interplay between different parts of dimeric single-stranded DNA-binding proteins we have studied the correlated motions in the protein encoded by filamentous phage Pf3 via the combined use of 15N-NMR relaxation experiments, molecular dynamics simulations and essential dynamics calculations. These studies provide insight into the mechanism underlying the protein-DNA binding reaction. The most important motions can be described by a few essential modes. Most outstanding is the correlated symmetric motion of the DNA-binding wings, which are far apart in the structure. This motion determines the access of DNA to the DNA-binding domain. A correlation between the motion of the DNA-binding wing and the complex loop is indicated to play a role in the cooperative binding of the protein to DNA. These motions are in the nanosecond regime in correspondence with the 15N-NMR relaxation experiments.

Molecular dynamics and accuracy of NMR structures: effects of error bounds and data removal.

The effect of internal dynamics on the accuracy of nuclear magnetic resonance (NMR) structures was studied in detail using model distance restraint sets (DRS) generated from a 6.6 nanosecond molecular dynamics trajectory of bovine pancreatic trypsin inhibitor. The model data included the effects of internal dynamics in a very realistic way. Structure calculations using different error estimates were performed with iterative removal of systematically violated restraints. The accuracy of each calculated structure was measured as the atomic root mean square (RMS) difference to the optimized average structure derived from the trajectory by structure factors refinement. Many of the distance restraints were derived from NOEs that were significantly affected by internal dynamics. Depending on the error bounds used, these distance restraints seriously distorted the structure, leading to deviations from the coordinate average of the dynamics trajectory even in rigid regions. Increasing error bounds uniformly for all distance restraints relieved the strain on the structures. However, the accuracy did not improve. Significant improvement of accuracy was obtained by identifying inconsistent restraints with violation analysis, and excluding them from the calculation. The highest accuracy was obtained by setting bounds rather tightly, and removing about a third of the restraints. The limiting accuracy for all backbone atoms was between 0.6 and 0.7 A. Also, the precision of the structures increased with removal of inconsistent restraints, indicating that a high precision is not simply the consequence of tight error bounds but of the consistency of the DRS. The precision consistently overestimated the accuracy.

Exploring protein interiors: the role of a buried histidine in the KH module fold.

The K-homology (KH) module is a novel RNA-binding motif. The structures of a representative KH motif from vigilin (vig-KH6) and of the first KH domain of fmr1 have been recently solved by nuclear magnetic resonance (NMR) and automated assignment-refinement techniques (ARIA). While a hydrophobic residue is found at position 21 in most of the KH modules, a buried His is conserved in all the 15 KH repeats of vigilin. This position must therefore have a key structural role in stabilizing the hydrophobic core. In the present work, we have addressed the following questions in order to obtain a detailed description of the role of His 21: i) what is the exact role of the histidine in the hydrophobic core of vig-KH6? ii) can we define the interactions that allow a conserved buried position to be occupied by a histidine both in vig-KH6 and in the whole vigilin KH sub-family? iii) how is the structure and stability of vig-KH6 influenced by the state of protonation of this histidine? To answer these questions, we have carried out an extensive refinement of the vig-KH6 structure using both an improved ARIA protocol starting from different initial structures and successively running restrained and unrestrained trajectories in water. An analysis of the stability of secondary structural elements, solvent accessibility, and hydrogen bonding patterns allows hypothesis on the structural role of residue His 21 and on the interactions that this residue forms with the environment. The importance of the protonation state of His 21 on the stability of the KH fold was addressed and validated by experimental results.