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.

Assessing the quality of solution nuclear magnetic resonance structures by complete cross-validation.

Structure determination of macromolecules in solution by nuclear magnetic resonance (NMR) spectroscopy involves the fitting of atomic models to the observed nuclear Overhauser effect (NOE) data. Complete cross-validation has been used to define reliable and unbiased criteria for the quality of solution NMR structures. The method is based on the partitioning of NOE data into test sets and the cross-validation of statistical quantities for each of the test sets. A high correlation between cross-validated measures of fit, such as distance bound violations and NMR R values, and the quality of solution NMR structures was observed. Less complete data resulted in poorer satisfaction of the cross-validated measures of fit. Optimization of cross-validated measures of fit will likely produce solution NMR structures with maximal information content.

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 folding catalyst protein disulfide isomerase is constructed of active and inactive thioredoxin modules.

BACKGROUND: Protein disulfide isomerase (PDI), a multifunctional protein of the endoplasmic reticulum, catalyzes the formation, breakage and rearrangement of disulfide bonds during protein folding. Dissection of this protein into its individual domains has confirmed the presence of the a and a' domains, which are homologous to thioredoxin, having related structures and activities. The a and a' domains both contain a -Cys-Gly-His-Cys- active-site sequence motif. The remainder of the molecule consists primarily of two further domains, designated b and b' which are thought to be sequence repeats on the basis of a limited sequence similarity. The functions of the b and b' domains are unknown and, until now, the structure of neither domain was known. RESULTS: Heteronuclear nuclear magnetic resonance (NMR) methods have been used to determine the global fold of the PDI b domain. The protein has an alpha/beta fold with the order of the elements of secondary structure being beta1-alpha1-beta2-alpha2-beta3-alpha3-beta4-beta5+ ++-alpha4. The strands are all in a parallel arrangement with respect to each other, except for beta4 which is antiparallel. The arrangement of the secondary structure elements of the b domain is identical to that found in the a domain of PDI and in the ubiquitous redox protein thioredoxin; the three-dimensional folding topology of the b domain is also very similar to that of these proteins. CONCLUSIONS: Our determination of the global fold of the b domain of PDI by NMR reveals that, like the a domain, the b domain contains the thioredoxin motif, even though the b domain has no significant amino-acid sequence similarities to any members of the thioredoxin family. This observation, together with indications that the b' domain adopts a similar fold, suggests that PDI consists of active and inactive thioredoxin modules. These modules may have been adapted during evolution to provide PDI with its complete spectrum of enzymatic activities.

Structure calculation from NMR data.

NMR calculation methods have kept pace with the rapid extension of NMR experiments to larger molecules. By including additional data and effects of local dynamics in the refinement, we can obtain a more complete picture of the molecule in solution. The structure determination process is being aided by new methods to solve some aspects of spectral assignment during the structure calculation.

The three-dimensional structure of a type I module from titin: a prototype of intracellular fibronectin type III domains.

BACKGROUND: Titin is a huge protein ( approximately 3 MDa) that is present in the contractile unit (sarcomere) of striated muscle and has a key role in muscle assembly and elasticity. Titin is mainly composed of two types of module (type I and II). Type I modules are found exclusively in the region of titin localised in the A band, where they are arranged in a super-repeat pattern that correlates with the ultrastructure of the thick filament. No structure of a titin type I module has been reported so far. RESULTS: We have determined the structure of a representative type I module, A71, using nuclear magnetic resonance (NMR) spectroscopy. The structure has the predicted fibronectin type III fold. Titin-specific conserved residues are either located at the putative module-module interfaces or along one side of the protein surface. Several proline residues that contribute to two stretches in a polyproline II helix conformation are solvent-exposed and line up as a continuous ribbon extending over more than two-thirds of the module surface. Homology models of the type I module N-terminal to A71 (A70) and the double module A70-A71 were used to discuss possible intermodule interactions and their role in module-module orientation. CONCLUSIONS: As residues at the module-module interfaces are highly conserved, we speculate that similar interactions govern all of the interfaces between type I modules in titin. This conservation would lead to a regular multiple array of similar surface structures. Such an arrangement would allow arrays of contiguous type I modules to expose multiple proline stretches in a highly regular way and these may act as binding sites for other thick filament proteins.

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.

The three-dimensional structure of the HRDC domain and implications for the Werner and Bloom syndrome proteins.

BACKGROUND: The HRDC (helicase and RNaseD C-terminal) domain is found at the C terminus of many RecQ helicases, including the human Werner and Bloom syndrome proteins. RecQ helicases have been shown to unwind DNA in an ATP-dependent manner. However, the specific functional roles of these proteins in DNA recombination and replication are not known. An HRDC domain exists in both of the human RecQ homologues that are implicated in human disease and may have an important role in their function. RESULTS: We have determined the three-dimensional structure of the HRDC domain in the Saccharomyces cerevisiae RecQ helicase Sgs1p by nuclear magnetic resonance (NMR) spectroscopy. The structure resembles auxiliary domains in bacterial DNA helicases and other proteins that interact with nucleic acids. We show that a positively charged region on the surface of the Sgs1p HRDC domain can interact with DNA. Structural similarities to bacterial DNA helicases suggest that the HRDC domain functions as an auxiliary domain in RecQ helicases. Homology models of the Werner and Bloom HRDC domains show different surface properties when compared with Sgs1p. CONCLUSIONS: The HRDC domain represents a structural scaffold that resembles auxiliary domains in proteins that are involved in nucleic acid metabolism. In Sgs1p, the HRDC domain could modulate the helicase function via auxiliary contacts to DNA. However, in the Werner and Bloom syndrome helicases the HRDC domain may have a role in their functional differences by mediating diverse molecular interactions.

Calculation of protein structures with ambiguous distance restraints. Automated assignment of ambiguous NOE crosspeaks and disulphide connectivities.

The distances derived from nuclear Overhauser effect (NOE) spectra are usually converted into three-dimensional structures by computer algorithms loosely termed distance geometry. To a varying degree, these methods require that the distance data is unambiguously assigned to pairs of atoms. Typically, however, there are many NOE crosspeaks that cannot be assigned without some knowledge of the structure. These crosspeaks have to be assigned in an iterative manner, using preliminary structures calculated from the unambiguous crosspeaks. In this paper, I present an alternative to this iterative approach. The ambiguity of an NOE crosspeak is correctly described in terms of the distances between all pairs of protons that may be involved. A simple restraining term is defined in terms of "ambiguous" distance restraints that can allow all possible assignments. A new minimization procedure based on simulated annealing is described that is capable of using highly ambiguous data for ab initio structure calculations. In particular, it is feasible to specify the restraint list directly in terms of the proton chemical shift assignment and the NOE peak table, without having assigned NOE crosspeaks to proton pairs. While the primary aim of this paper is determining the global fold of proteins from NMR data, similar strategies can be used for other types of ambiguous distance data. The application to one example, disulphide bridges with unknown connectivity, is described. Model NOE data were generated from the X-ray crystal structure of a small protein with known chemical shift assignments. Varying degrees of ambiguity in the data were assumed. The method obtained the correct polypeptide fold even when all distance restraints were ambiguous. Thus, the new approach may facilitate structure calculations with data derived from very overlapped spectra. It is also a step towards automating the calculation of structures from NMR data. This could prove especially valuable for data derived from three- and four-dimensional experiments. The approach may also prove useful for model building studies and tertiary structure prediction.

Are there non-trivial dynamic cross-correlations in proteins?

The analysis of internal motion in ensembles of flexible molecules in coordinate space requires the removal of overall motion by a least-squares fitting procedure of the Cartesian coordinates. It has been demonstrated that the choice of the atom set used for fitting influences the picture of the internal motion of BPTI. We have performed essential dynamics analyses of a 1 ns molecular dynamics trajectory of the single-stranded DNA-binding protein from the Pf3 phage using either all alpha-carbon atoms or the least mobile ones for fitting the trajectory prior to the analysis. We found that covariances of atoms separated by long distances were significantly reduced in the latter case; the overall overlap of essential spaces was still high. In the second part we present a method that does not introduce and bias caused by overall motion: principal component analysis in distance space. Non-trivial dynamic cross-correlations were preserved in distance space, which answers the question posed in the title in the affirmative. However, cross-correlations were throughout smaller than those detected by standard essential dynamics analyses.

Influence of internal dynamics on accuracy of protein NMR structures: derivation of realistic model distance data from a long molecular dynamics trajectory.

In order to study the effect of internal dynamics on the accuracy of NMR structures in detail, we generated NOE distance data from a long molecular dynamics trajectory of BPTI. Cross-relaxation rates were calculated from the trajectory by analysis of the appropriate proton-proton vector autocorrelation functions. A criterion for the convergence of correlation functions was developed, and the analysis was restricted to those correlation functions that had converged within the simulation time. Effective distances were determined from the calculated cross-relaxation rates. Internal dynamics affected the derived distances in a realistic way, since they were subject both to radial averaging (which increases the cross-relaxation rate) and angular averaging (which decreases the cross-relaxation rate). The comparison of the effective distances with average distance between the protons during the trajectory showed that for most the effects of angular and distance averaging essentially cancel out. For these distances, the effective distance derived from an NOE is therefore a very good estimate of the average distance, or the distance in the average structure. However, for about 10% of the distances, the effective distance was more than 10% larger than the average distance, while for about 5%, it was more than 10% smaller, in some cases by more than 2 A. Little correlation is observed between the effects on cross-relaxation rates to different protons of the same residue. The results of this analysis have implications for the way structures are calculated from NOE distance data. For many distances, the assumption of a rigid structure is valid, and large error bounds would result in the loss of too much information content. On the other hand, the error bounds very often employed are not wide enough for some of the effects seen in our study.

Relaxation matrix refinement of the solution structure of squash trypsin inhibitor.

The structure of the small squash trypsin inhibitor CMTI-I is refined by directly minimizing the difference between the observed two-dimensional nuclear Overhauser enhancement (NOE) intensities and those calculated by the full relaxation matrix approach. To achieve this, a term proportional to this difference was added to the potential energy function of the molecular dynamics program X-PLOR. Derivatives with respect to atomic co-ordinates are calculated analytically. Spin diffusion effects are thus accounted for fully during the refinement. Initial structures for the refinement were those determined recently by solution nuclear magnetic resonance using the isolated two-spin approximation to derive distance range estimates. The fits to the nuclear magnetic resonance data improve significantly with only small shifts in the refined structures during a few cycles of conjugate gradient minimization. However, larger changes (approximately 1 A) in the conformation occur during simulated annealing, which is accompanied by a further reduction of the difference between experimental and calculated two-dimensional NOE intensities. The refined structures are closer to the X-ray structure of the inhibitor complexed with trypsin than the initial structures. The root-mean-square difference for backbone atoms between the initial structures and the X-ray structure is 0.96 A, and that between the refined structures and the X-ray structure 0.61 A.

Automated NOESY interpretation with ambiguous distance restraints: the refined NMR solution structure of the pleckstrin homology domain from beta-spectrin.

We have used a novel, largely automated, calculation method to refine the NMR solution structure of the pleckstrin homology domain of beta-spectrin. The method is called ARIA for Ambiguous Restraints for Iterative Assignment. The starting point for ARIA is an almost complete assignment of the proton chemical shifts, and a list of partially assigned NOEs, mostly sequential and secondary structure NOEs. The restraint list is then augmented by automatically interpreting peak lists generated by automated peak-picking. The central task of ARIA is the assignment of ambiguous NOEs during the structure calculation using a combination of ambiguous distance restraints and an iterative assignment strategy. In addition, ARIA calibrates ambiguous NOEs to derive distance restraints, merges overlapping data sets to remove duplicate information, and uses empirical rules to identify erroneous peaks. While the distance restraints for the structure calculations were exclusively extracted from homonuclear 2D experiments, ARIA is especially suited for the analysis of multidimensional spectra. Applied to the pleckstrin homology domain, ARIA generated structures of good quality, and of sufficiently high accuracy to solve the X-ray crystal structure of the same domain by molecular replacement. The comparison of the free NMR solution structure to the X-ray structure, which is complexed to D-myo-inositol-1,4,5-triphosphate, shows that the ligand primarily induces a disorder-order transition in the binding loops, which are disordered in the NMR ensemble but well ordered in the crystal. The structural core of the protein is unaffected, as evidenced by a backbone root-mean-square difference between the average NMR coordinates and the X-ray crystal structure for the secondary structure elements of less than 0.6 A.

Solution structure of the spectrin repeat: a left-handed antiparallel triple-helical coiled-coil.

Cytoskeletal proteins belonging to the spectrin family have an elongated structure composed of repetitive units. The three-dimensional solution structure of the 16th repeat from chicken brain alpha-spectrin (R16) has been determined by NMR spectroscopy and distance geometry-simulated annealing calculations. We used a total of 1035 distance restraints, which included 719 NOE-based values obtained by applying the ambiguous restraints for iterative assignment (ARIA) method. In addition, we performed a direct refinement against 1H-chemical shifts. The final ensemble of 20 structures shows an average RMSD of 1.52 A from the mean for the backbone atoms, excluding loops and N and C termini. R16 is made up of three antiparallel alpha-helices separated by two loops, and folds into a left-handed coiled-coil. The basic unit of spectrin is an antiparallel heterodimer composed of two homologous chains, beta and alpha. These assemble a tetramer via a mechanism that relies on the completion of a single repeat by association of the partial repeats located at the C terminus of the beta-chain (two helices) and at the N terminus of the alpha-chain (one helix). This tetramer is the assemblage able to cross-link actin filaments. Model building by homology of the "tetramerization" repeat from human erythrocyte spectrin illuminates the possible role of point mutations which cause hemolytic anemias.

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.

Heteronuclear relaxation study of the PH domain of beta-spectrin: restriction of loop motions upon binding inositol trisphosphate.

The structural dynamics of protein ligand binding sites is one factor determining the specificity towards related ligands. In this context, the spectrin PH domain, which binds to a number of phosphatidylinositol lipid head groups, was investigated with respect to the dynamics of the binding loops. The latter were found to be of intermediate flexibility on a picosecond to nanosecond time-scale in the free protein and become more rigid upon ligand binding. Significant 15N and proton chemical shift changes occur in the binding loops. The internal correlation time, determined from 15N heteronuclear relaxation data using the standard model-free approach, decreases upon ligand binding. For several residues a concomitant rise in the generalized order parameter is observed. This is interpreted as a dampening effect of the ligand on a slow loop motion, while a fast component is not affected. Molecular dynamics simulations were performed to further investigate this situation. In fact, two time-scales of loop motions in the free state are observed in a 9 ns molecular dynamics trajectory. Agreement with generalized order parameters obtained from the experiment improves when a subtrajectory is analyzed that excludes rare dihedral transitions.

A model of the complex between single-stranded DNA and the single-stranded DNA binding protein encoded by gene V of filamentous bacteriophage M13.

A contact analysis and a series of restrained molecular dynamics simulations were employed to derive a model of the complex between single-stranded DNA and the single-stranded DNA-binding protein encoded by gene V of the filamentous phage M13. The study is based on the recently elucidated solution structure of the Tyr41-->His mutant of the protein. Electron microscopy studies, indicating that the complex forms a flexible, left-handed helical coil with a diameter of 8 to 9 nm and an average pitch of 9 nm, were taken into consideration. The contact analysis served to determine the helix parameters that permit the energetically most favourable packing of protein molecules. Then a protein super-helix was built, into which two extended strands of DNA were modelled using restrained molecular dynamics. Specific constraints were included to ensure that the DNA would position itself into the binding groove of the protein. These constraints are based on recent NMR spin label experiments which offered a direct identification of the amino acids of the protein present in the DNA-binding domain. We present a model for the complex which is in full agreement with the existing reliable biophysical and biochemical data. A description of the protein-protein interface is given and the protein-DNA interaction is discussed in view of the derived model. In addition, we demonstrate that, on the basis of the available experimental data, and not imposing the left-handedness of the nucleoprotein complex, it is feasible to build also a plausible model for the complex which exhibits the opposite, i.e. right-handed, helical sense. This nucleoprotein structure features characteristics highly similar to those of the left-handed helix.

The solution structure of the Tyr41-->His mutant of the single-stranded DNA binding protein encoded by gene V of the filamentous bacteriophage M13.

The solution structure of mutant Tyr41-->His of the single-stranded DNA binding protein encoded by gene V of the filamentous bacteriophage M13 has been investigated by nuclear magnetic resonance spectroscopy. Two- and three-dimensional NMR experiments have been employed with a variety of NMR samples of gene V protein, some of which were uniformly enriched with either 15N or 13C. The structure of mutant Tyr41-->His of the M13 gene V protein which occurs in solution as a symmetric dimer was calculated using a two-stage procedure. The first step of the procedure involved the calculation of a set of individual monomer structures using the distance geometry program DIANA. This was then followed by the calculation of dimer structures employing "simulated annealing" protocols with the program X-PLOR. Hereby, the problem of assignment of intra- and inter-subunit NOEs of the symmetric dimer was circumvented through use of a target function that correctly deals with the intra- and inter-subunit contributions to the NOE peaks. Furthermore, a pseudo energy term was employed to restrain the symmetry of the dimer. In addition to this novel calculation strategy, we have incorporated distance information for a set of NOEs which were unambiguously identified as inter-subunit NOEs using an NMR strategy based on asymmetric labelling. A total of 20 structures were calculated for the M13 gene V protein mutant Tyr41-->His based on approximately 1000 experimental restraints derived from the NMR data. The structure of residues 1 to 15 and 29 to 87 of both monomers is reasonably well determined with an average atomic r.m.s. difference between the individual structures and the respective mean structure of approximately 0.9 A for the backbone atoms and approximately 1.4 A for all atoms. The orientation of the exposed anti-parallel beta-loop (residues 16 to 28) with respect to the core could not be determined. The molecular architecture of each of the monomers includes a five-stranded beta-barrel enclosing a hydrophobic core and two-antiparallel beta-loops. The dimer structure is stabilized predominantly by hydrophobic residues primarily involving the symmetry-related dyad domains (residues 64 to 82) of the monomers. Residues which are close to bound single-stranded DNA were identified previously from binding experiments with spin-labelled oligonucleotides. The solution structure of mutant Tyr41-->His of the M13 gene V protein is consistent with these binding data and provides a clear view of the protein's single-stranded DNA binding path.

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.