GAPDOCK: A genetic algorithm approach to protein docking in CAPRI round 1

As part of the first Critical Assessment of PRotein Interactions, round 1, we predict the structure of two protein-protein complexes, by using a genetic algorithm, GAPDOCK, in combination with surface complementarity, buried surface area, biochemical information, and human intervention. Among the five models submitted for target 1, HPr phosphocarrier protein (B. subtilis) and the hexameric HPr kinase (L. lactis), the best correctly predicts 17 of 52 interprotein contacts, whereas for target 2, bovine rotavirus VP6 protein-monoclonal antibody, the best model predicts 27 of 52 correct contacts. Given the difficult nature of the targets, these predictions are very encouraging and compare well with those obtained by other methods. Nevertheless, it is clear that there is a need for improved methods for distinguishing between correct and plausible but incorrect complexes. Proteins 2003;52:10-14

Diversity in Free Energy Landscape of Proteins with the Same Native Topology

In order to elucidate the role of the native state topology and the stability of subdomains in protein folding, we investigate free energy landscape of human lysozyme, which is composed of two subdomains, by Monte Carlo simulations. A realistic lattice model with Go-like interaction is used. We take the relative interaction strength (stability, in other word) of two subdomains as a variable parameter and study the folding process. A variety of folding process is observed and we obtained a phase diagram of folding in terms of temperature and the relative stability. Experimentally-observed diversity in folding process of c-type lysozimes is thus understood as a consequence of the difference in the relative stability of subdomains.Comment: 4 pages, 3 figure

Co-populated Conformational Ensembles of β(2)-Microglobulin Uncovered Quantitatively by Electrospray Ionization Mass Spectrometry

Ordered assembly of monomeric human β(2)-microglobulin (β(2)m) into amyloid fibrils is associated with the disorder hemodialysis-related amyloidosis. Previously, we have shown that under acidic conditions (pH <5.0 at 37 °C), wild-type β(2)m assembles spontaneously into fibrils with different morphologies. Under these conditions, β(2)m populates a number of different conformational states in vitro. However, this equilibrium mixture of conformationally different species is difficult to resolve using ensemble techniques such as nuclear magnetic resonance or circular dichroism. Here we use electrospray ionization mass spectrometry to resolve different species of β(2)m populated between pH 6.0 and 2.0. We show that by linear deconvolution of the charge state distributions, the extent to which each conformational ensemble is populated throughout the pH range can be determined and quantified. Thus, at pH 3.6, conditions under which short fibrils are produced, the conformational ensemble is dominated by a charge state distribution centered on the 9+ ions. By contrast, under more acidic conditions (pH 2.6), where long straight fibrils are formed, the charge state distribution is dominated by the 10+ and 11+ ions. The data are reinforced by investigations on two variants of β(2)m (V9A and F30A) that have reduced stability to pH denaturation and show changes in the pH dependence of the charge state distribution that correlate with the decrease in stability measured by tryptophan fluorescence. The data highlight the potential of electrospray ionization mass spectrometry to resolve and quantify complex mixtures of different conformational species, one or more of which may be important in the formation of amyloid

Structure of Tagatose-1,6-bisphosphate Aldolase. Insight into chiral discrimination, mechanism, and specificity of class II aldolases

Tagatose-1,6-bisphosphate aldolase (TBPA) is a tetrameric class II aldolase that catalyzes the reversible condensation of dihydroxyacetone phosphate with glyceraldehyde 3-phosphate to produce tagatose 1,6-bisphosphate. The high resolution (1.45 Å) crystal structure of the Escherichia coli enzyme, encoded by the agaY gene, complexed with phosphoglycolohydroxamate (PGH) has been determined. Two subunits comprise the asymmetric unit, and a crystallographic 2-fold axis generates the functional tetramer. A complex network of hydrogen bonds position side chains in the active site that is occupied by two cations. An unusual Na(+) binding site is created using a interaction with Tyr(183) in addition to five oxygen ligands. The catalytic Zn(2+) is five-coordinate using three histidine nitrogens and two PGH oxygens. Comparisons of TBPA with the related fructose-1,6-bisphosphate aldolase (FBPA) identifies common features with implications for the mechanism. Because the major product of the condensation catalyzed by the enzymes differs in the chirality at a single position, models of FBPA and TBPA with their cognate bisphosphate products provide insight into chiral discrimination by these aldolases. The TBPA active site is more open on one side than FBPA, and this contributes to a less specific enzyme. The availability of more space and a wider range of aldehyde partners used by TBPA together with the highly specific nature of FBPA suggest that TBPA might be a preferred enzyme to modify for use in biotransformation chemistry

Three-dimensional Structure of L-2-Haloacid Dehalogenase from Xanthobacter autotrophicus GJ10 Complexed with the Substrate-analogue Formate

The L-2-haloacid dehalogenase from the 1,2-dichloroethane degrading bacterium Xanthobacter autotrophicus GJ10 catalyzes the hydrolytic dehalogenation of small L-2-haloalkanoic acids to yield the corresponding D-2-hydroxyalkanoic acids. Its crystal structure was solved by the method of multiple isomorphous replacement with incorporation of anomalous scattering information and solvent flattening, and was refined at 1.95-Å resolution to an R factor of 21.3%. The three-dimensional structure is similar to that of the homologous L-2-haloacid dehalogenase from Pseudomonas sp. YL (1), but the X. autotrophicus enzyme has an extra dimerization domain, an active site cavity that is completely shielded from the solvent, and a different orientation of several catalytically important amino acid residues. Moreover, under the conditions used, a formate ion is bound in the active site. The position of this substrate-analogue provides valuable information on the reaction mechanism and explains the limited substrate specificity of the Xanthobacter L-2-haloacid dehalogenase.

A 1-acetamido derivative of 6-epi-valienamine: an inhibitor of a diverse group of β-N-acetylglucosaminidases

The synthesis of an analogue of 6-epi-valienamine bearing an acetamido group and its characterisation as an inhibitor of β-N-acetylglucosaminidases are described. The compound is a good inhibitor of both human O-GlcNAcase and human β-hexosaminidase, as well as two bacterial β-N-acetylglucosaminidases. A 3-D structure of the complex of Bacteroides thetaiotaomicron BtGH84 with the inhibitor shows the unsaturated ring is surprisingly distorted away from its favoured solution phase conformation and reveals potential for improved inhibitor potency

Secondary-Structure Design of Proteins by a Backbone Torsion Energy

We propose a new backbone-torsion-energy term in the force field for protein systems. This torsion-energy term is represented by a double Fourier series in two variables, the backbone dihedral angles phi and psi. It gives a natural representation of the torsion energy in the Ramachandran space in the sense that any two-dimensional energy surface periodic in both phi and psi can be expanded by the double Fourier series. We can then easily control secondary-structure-forming tendencies by modifying the torsion-energy surface. For instance, we can increase/decrease the alpha-helix-forming-tendencies by lowering/raising the torsion-energy surface in the alpha-helix region and likewise increase/decrease the beta-sheet-forming tendencies by lowering/raising the surface in the beta-sheet region in the Ramachandran space. We applied our approach to AMBER parm94 and AMBER parm96 force fields and demonstrated that our modifications of the torsion-energy terms resulted in the expected changes of secondary-structure-forming-tendencies by performing folding simulations of alpha-helical and beta-hairpin peptides.Comment: 13 pages, (Revtex4), 5 figure

Protein dynamics with off-lattice Monte Carlo moves

A Monte Carlo method for dynamics simulation of all-atom protein models is introduced, to reach long times not accessible to conventional molecular dynamics. The considered degrees of freedom are the dihedrals at C$_\alpha$-atoms. Two Monte Carlo moves are used: single rotations about torsion axes, and cooperative rotations in windows of amide planes, changing the conformation globally and locally, respectively. For local moves Jacobians are used to obtain an unbiased distribution of dihedrals. A molecular dynamics energy function adapted to the protein model is employed. A polypeptide is folded into native-like structures by local but not by global moves.Comment: 10 pages, 4 Postscript figures, uses epsf.sty and a4.sty; scheduled tentatively for Phys.Rev.E issue of 1 March 199

An ATP-binding cassette-type cysteine transporter in Campylobacter jejuni inferred from the structure of an extracytoplasmic solute receptor protein

Campylobacter jejuni is a Gram-negative food-borne pathogen associated with gastroenteritis in humans as well as cases of the autoimmune disease Guillain Barre syndrome. C. jejuni is asaccharolytic because it lacks an active glycolytic pathway for the use of sugars as a carbon source. This suggests an increased reliance on amino acids as nutrients and indeed the genome sequence of this organism indicates the presence of a number of amino acid uptake systems. Cj0982, also known as CjaA, is a putative extracytoplasmic solute receptor for one such uptake system as well as a major surface antigen and vaccine candidate. The crystal structure of Cj0982 reveals a two-domain protein with density in the enclosed cavity between the domains that clearly defines the presence of a bound cysteine ligand. Fluorescence titration experiments were used to demonstrate that Cj0982 binds cysteine tightly and specifically with a K-d of similar to 10(-7) M consistent with a role as a receptor for a high- affinity transporter. These data imply that Cj0982 is the binding protein component of an ABC-type cysteine transporter system and that cysteine uptake is important in the physiology of C. jejuni

The PDZ domain of the SpoIVB serine peptidase facilitates multiple functions

During spore formation in Bacillus subtilis, the SpoIVB protein is a critical component of the sigma (K) regulatory checkpoint. SpoIVB has been shown to be a serine peptidase that is synthesized in the spore chamber and which self-cleaves, releasing active forms. These forms can signal proteolytic processing of the transcription factor sigma (K) in the outer mother cell chamber of the sporulating cell. This forms the basis of the sigma (K) checkpoint and ensures accurate sigma (K)-controlled gene expression. SpoIVB has also been shown to activate a second distinct process, termed the second function, which is essential for the formation of heat-resistant spores. In addition to the serine peptidase domain, SpoIVB contains a PDZ domain. We have altered a number of conserved residues in the PDZ domain by site-directed mutagenesis and assayed the sporulation phenotype and signaling properties of mutant SpoIVB proteins. Our work has revealed that the SpoIVB PDZ domain could be used for up to four distinct processes, (i) targeting of itself for trans proteolysis, (11) binding to the protease inhibitor BofC, (iii) signaling of pro-sigma (K) processing, and (iv) signaling of the second function of SpoIVB