Analysis of multiple crystal forms of Bacillus subtilis BacB suggests a role for a metal ion as a nucleant for crystallization

Bacillus subtilis BacB is an oxidase that is involved in the production of the antibiotic bacilysin. This protein contains two double-stranded beta-helix (cupin) domains fused in a compact arrangement. BacB crystallizes in three crystal forms under similar crystallization conditions. An interesting observation was that a slight perturbation of the crystallization droplet resulted in the nucleation of a different crystal form. An X-ray absorption scan of BacB suggested the presence of cobalt and iron in the crystal. Here, a comparative analysis of the different crystal forms of BacB is presented in an effort to identify the basis for the different lattices. It is noted that metal ions mediating interactions across the asymmetric unit dominate the different packing arrangements. Furthermore, a normalized B-factor analysis of all the crystal structures suggests that the solvent-exposed metal ions decrease the flexibility of a loop segment, perhaps influencing the choice of crystal form. The residues coordinating the surface metal ion are similar in the triclinic and monoclinic crystal forms. The coordinating ligands for the corresponding metal ion in the tetragonal crystal form are different, leading to a tighter packing arrangement. Although BacB is a monomer in solution, a dimer of BacB serves as a template on which higher order symmetrical arrangements are formed. The different crystal forms of BacB thus provide experimental evidence for metal-ion-mediated lattice formation and crystal packing.

Structures of apo and GTP-bound molybdenum cofactor biosynthesis protein MoaC from Thermus thermophilus HB8

The first step in the molybdenum cofactor (Moco) biosynthesis pathway involves the conversion of guanosine triphosphate (GTP) to precursor Z by two proteins (MoaA and MoaC). MoaA belongs to the S-adenosylmethioninedependent radical enzyme superfamily and is believed to generate protein and/or substrate radicals by reductive cleavage of S-adenosylmethionine using an Fe-S cluster. MoaC has been suggested to catalyze the release of pyrophosphate and the formation of the cyclic phosphate of precursor Z. However, structural evidence showing the binding of a substrate-like molecule to MoaC is not available. Here, apo and GTP-bound crystal structures of MoaC from Thermus thermophilus HB8 are reported. Furthermore, isothermal titration calorimetry experiments have been carried out in order to obtain thermodynamic parameters for the protein-ligand interactions. In addition, molecular-dynamics (MD) simulations have been carried out on the protein-ligand complex of known structure and on models of relevant complexes for which X-ray structures are not available. The biophysical, structural and MD results reveal the residues that are involved in substrate binding and help in speculating upon a possible mechanism.

Drosophila indirect flight muscle specific Act88F actin mutants as a model system for studying congenital myopathies of the human ACTA1 skeletal muscle actin gene

Most human ACTA1 skeletal actin gene mutations cause dominant, congenital myopathies often with severely reduced muscle function and neonatal mortality. High sequence conservation of actin means many mutated ACTA1 residues are identical to those in the Drosophila Act88F, an indirect flight muscle specific sarcomeric actin. Four known Act88F mutations occur at the same actin residues mutated in ten ACTA1 nemaline mutations, A138D/P, R256H/L, G268C/D/R/S and R372C/S. These Act88F mutants were examined for similar muscle phenotypes. Mutant homozygotes show phenotypes ranging from a lack of myofibrils to almost normal sarcomeres at eclosion. Aberrant Z-disc-like structures and serial Z-disc arrays, ‘zebra bodies’, are observed in homozygotes and heterozygotes of all four Act88F mutants. These electron-dense structures show homologies to human nemaline bodies/rods, but are much smaller than those typically found in the human myopathy. We conclude that the Drosophila indirect flight muscles provide a good model system for studying ACTA1 mutations.

HARMONY: a server for the assessment of protein structures

Protein structure validation is an important step in computational modeling and structure determination. Stereochemical assessment of protein structures examine internal parameters such as bond lengths and Ramachandran (phi, psi) angles. Gross structure prediction methods such as inverse folding procedure and structure determination especially at low resolution can sometimes give rise to models that are incorrect due to assignment of misfolds or mistracing of electron density maps. Such errors are not reflected as strain in internal parameters. HARMONY is a procedure that examines the compatibility between the sequence and the structure of a protein by assigning scores to individual residues and their amino acid exchange patterns after considering their local environments. Local environments are described by the backbone conformation, solvent accessibility and hydrogen bonding patterns. We are now providing HARMONY through a web server such that users can submit their protein structure files and, if required, the alignment of homologous sequences. Scores are mapped on the structure for subsequent examination that is useful to also recognize regions of possible local errors in protein structures. HARMONY server is located at http://caps.ncbs.res.in/harmony/

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Oligosaccharides modulate the apoptotic activity of glycodelin

GlycodelinA (GdA), a multifunctional glycoprotein secreted at high concentrations by the uterine endometrium during the early phases of pregnancy, carries glycan chains on asparagines at positions N28 and N63. GdA purified from amniotic fluid is known to be a suppressor of T-cell proliferation, an inducer of T-cell apoptosis, and an inhibitorof sperm-zona binding in contrast to its glycoform, glycodelinS (GdS), which is secreted by the seminal vesicles into the seminal plasma. The oligosaccharide chains of GdA terminate in sialic acid residues, whereas those of GdS are not sialylated but are heavily fucosylated. Our previous work has shown that the apoptogenic activity of GdA resides in the protein backbone, and we have also demonstrated the importance of sialylation for the manifestation of GdA-induced apoptosis. Recombinant glycodelin (Gd) expressed in the Sf21 insec cell line yielded an apoptotically active Gd; however, the same geneexpressed in the insect cell line Tni produced apoptotically inactive Gd, as observed with the gene expressed in the Chinese hamster ovary(CHO) cell line and earlier in Pichia pastoris. Glycan analysis of the Tni and Sf21 cell line-expressed Gd proteins reveals differences in their glycan structures, which modulate the manifestation of apoptogenic activity of Gd. Through apoptotic assays carried out with the wild-type (WT) and glycosylation mutants of Gd expressed in Sf21 and Tni cells before and after mannosidase digestion, we conclude that the accessibility to the apoptogenic region of Gd is influenced by the size of the glycans.

Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. There quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Angstrom from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Nai counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. in the absence of any coordinated ion. due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.

Solution structure of O-glycosylated C-terminal leucine zipper domain of human salivary mucin (MUC7)

Solution structures of a 23 residue glycopeptide II (KIS* RFLLYMKNLLNRIIDDMVEQ, where * denotes the glycan Gal-beta-(1-3)-alpha-GalNAc) and its deglycosylated counterpart I derived from the C-terminal leucine zipper domain of low molecular weight human salivary mucin (MUC7) were studied using CD, NMR spectroscopy and molecular modeling. The peptide I was synthesized using the Fmoc chemistry following the conventional procedure and the glycopeptide II was synthesized incorporating the O-glycosylated building block (N alpha-Fmoc-Ser-[Ac-4,-beta-D-Gal-(1,3)-Ac(2)alpha-D-GalN(3)]-OPfp) at the appropriate position in stepwise assembly of peptide chain. Solution structures of these glycosylated and nonglycosylated peptides were studied in water and in the presence of 50% of an organic cosolvent, trifluoroethanol (TFE) using circular dichroism (CD), and in 50% TFE using two-dimensional proton nuclear magnetic resonance (2D H-1 NMR) spectroscopy. CD spectra in aqueous medium indicate that the apopeptide I adapts, mostly, a beta-sheet conformation whereas the glycopeptide II assumes helical structure. This transition in the secondary structure, upon glycosylation, demonstrates that the carbohydrate moiety exerts significant effect on the peptide backbone conformation. However, in 50% TFE both the peptides show pronounced helical structure. Sequential and medium range NOEs, C alpha H chemical shift perturbations, (3)J(NH:C alpha H) couplings and deuterium exchange rates of the amide proton resonances in water containing 50% TFE indicate that the peptide I adapts alpha-helical structure from Ile2-Val21 and the glycopeptide II adapts alpha-helical structure from Ser3-Glu22. The observation of continuous stretch of helix in both the peptides as observed by both NMR and CD spectroscopy strongly suggests that the C-terminal domain of MUC7 with heptad repeats of leucines or methionine residues may be stabilized by dimeric leucine zipper motif. The results reported herein may be invaluable in understanding the aggregation (or dimerization) of MUC7 glycoprotein which would eventually have implications in determining its structure-function relationship.

Computer modelling approach to study the modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein

The modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein (ABP) have been studied by energy minimization using the low resolution (2.4 A) X-ray data of the protein. These studies suggest that these sugars preferentially bind in the alpha-form to ABP, unlike L-arabinose where both alpha- and beta-anomers bind almost equally. The best modes of binding of alpha- and beta-anomers of D-galactose and D-fucose differ slightly in the nature of the possible hydrogen bonds with the protein. The residues Arg 151 and Asn 232 of ABP from bidentate hydrogen bonds with both L-arabinose and D-galactose, but not with D-fucose or D-glucose. However in the case of L-arabinose, Arg 151 forms hydrogen bonds with the hydroxyl group at the C-4 atom and the ring oxygen, whereas in case of D-galactose it forms bonds with the hydroxyl groups at the C-4 and C-6 atoms of the pyranose ring. The calculated conformational energies also predict that D-galactose is a better inhibitor than D-fucose and D-glucose, in agreement with kinetic studies. The weak inhibitor D-glucose binds preferentially to one domain of ABP leading to the formation of a weaker complex. Thus these studies provide information about the most probable binding modes of these sugars and also provide a theoretical explanation for the observed differences in their binding affinities.

Crystal structure of L-lysyl-L-glutamic acid dihydrate

L-Lysyl-L-glutamic acid dihydrate, C11N3O5H21·2H2O, crystallizes in the monoclinic space group P21 with a = 12.474(2), b = 5.020(1), c = 13.157(2) Å, β= 114.69(1)° and Z = 2. The crystal structure was solved by direct methods and refined to an R value of 0.037 using full matrix least-squares method. The molecule exists as a double zwitterion with both the amino and carboxyl groups ionised. The peptide has a folded conformation with its Lys residue trans and Glu residue gauche−gauche+. The side chains of the Lys and Glu residues correspond to all trans and folded (g−g−g−) conformations respectively. The terminal carboxyl group forms hydrogen bonds with the ξ-amino group of the lysine side chain. The head-to-tail interaction often seen in peptide crystals is absent in the present structure. In the extended crystal structure water molecules form channels along the b direction and are enclosed within helically arranged hydrogen bonds formed by the lysine side chain and the peptide backbone.

Purification and properties of β-glucosidase from a thermophilic fungus Humicola lanuginosa (Griffon and Maublanc) Bunce

An extracellular β-glucosidase (EC 3.2.1.21) has been purified to homogeneity from the culture filtrate of a thermophilic fungus, Humicola lanuginosa (Griffon and Maublanc) Bunce, using duplicating paper as the carbon source. The enzyme was purified 82-fold with a 43% yield by ion-exchange chromatography and gel filtration. The molecular weight of the protein was estimated to be 135,000 by gel filtration and 110,000 by electrophoresis. The sedimentation coefficient was 10.5 S. It was an acidic protein containing high amounts of acidic amino acid residues. It was poor in sulphur-containing amino acids. It also contained 9% carbohydrate. The enzyme activity was optimum at pH 4.5 and at 60°C. The enzyme was stable in the pH range 6–9 for 24 h at 25°C. The enzyme had similar affinities towards cellobiose and p-nitrophenyl-β-d-glucoside with Km values of 0.44 mM and 0.50 mM, respectively. The enzyme was capable of hydrolysing larchwood xylan, xylobiose and p-nitrophenyl-β-d-xyloside, though to a lesser extent. The enzyme was specific for the β-configuration and glucose moiety in the substrate.

An unusual C-H center dot center dot center dot O hydrogen bond mediated reversal of polypeptide chain direction in a synthetic peptide helix

An unusual C-terminal conformation has been detected in a synthetic decapeptide designed to analyze the stereochemistry of helix termination in polypeptides. The crystal structure of the decapeptide Boc-Leu-Aib-Val-Ala-Leu-Aib-Val-(D)Ala-(D)Leu-Aib-OMe reveals a helical segment spanning residues 1-7 and helix termination by formation of a Schellman motif, generated by (D)Ala(8) adopting the left-handed helical (alpha(L)) conformation. The extended conformation at (D)Leu(9) results in a compact folded structure, stabilized by a potentially strong C-H ... O hydrogen bond between Ala(4) (CH)-H-alpha and (D)Leu(9)CO. The parameters for C-H ... O interaction are Ala(4) (CH)-H-alpha .. O=C (D)Leu(9) distance 3.27 Angstrom C-alpha-H .. O angle 176 degrees, and O .. H-alpha distance 2.29 Angstrom. This structure suggests that insertion of contiguous D-residues may provide a handle for the generation of designed structures containing more than one helical segment folded in a compact manner. (C) 2000 Academic Press.

Primary structure of belladonna mottle virus coat protein.

The coat protein of belladonna mottle virus (a tymovirus) was cleaved by trypsin and chymotrypsin, and the peptides were separated by high performance liquid chromatography using a combination of gel permeation, reverse phase, and ion pair chromatography. The peptides were sequenced manually using the 4-N, N-dimethylaminoazobenzene-4'-isothiocyanate/phenyl isothiocyanate double-coupling method. The chymotryptic peptides were aligned by overlapping sequences of tryptic peptides and by homology with another tymovirus, eggplant mosaic virus. The belladonna mottle virus is more closely related to eggplant mosaic virus than to turnip yellow mosaic virus, the type member of this group, as evident from the sequence homologies of 57 and 32%, respectively. The accumulation of basic residues at the amino terminus implicated in RNA-protein interactions in many spherical plant viruses was absent in all the three sequences. Interestingly, the amino-terminal region is the least conserved among the tymoviruses. The longest stretch of conserved sequence between belladonna mottle virus and eggplant mosaic virus was residues 34-44, whereas it was residues 96-102 in the case of belladonna mottle virus and turnip yellow mosaic virus. A tetrapeptide in the region (residues 154-157) was found to be common for all the three sequences. It is possible that these conserved regions (residues 34-44, 96-102, 154-157) are involved in either intersubunit or RNA-protein interactions.

Interaction of 3-hydroxybenzoate-6-hydroxylase with cibacron blue

Cibacron blue is a potent inhibitor of 3-HBA-6-hydroxylase at a concentration < 1 mu M. Kinetic analyses revealed that at a concentration below 0.5 mu M the dye behaves as an uncompetitive inhibitor with respect to 3-HBA and competes with NADH for the same site on the enzyme. The alteration of the near-UV CD spectrum and quenching of the emission fluorescence of the enzyme by cibacron blue indicates a significant alteration in the environment of aromatic amino acid residues due to a stacking interaction and subtle conformatiodnal changes in the enzyme. The concentration-dependent quenching of the intrinsic fluorescence of the enzyme by cibacron blue was employed to determine the binding parameters such as association constant (K-a) and stoichiometry (r) for the enzyme-dye complex.

Structure of uracil-DNA glycosylase from Mycobacterium tuberculosis: insights into interactions with ligands

Uracil N-glycosylase (Ung) is the most thoroughly studied of the group of uracil DNA-glycosylase (UDG) enzymes that catalyse the first step in the uracil excision-repair pathway. The overall structure of the enzyme from Mycobacterium tuberculosis is essentially the same as that of the enzyme from other sources. However, differences exist in the N- and C-terminal stretches and some catalytic loops. Comparison with appropriate structures indicate that the two-domain enzyme closes slightly when binding to DNA, while it opens slightly when binding to the proteinaceous inhibitor Ugi. The structural changes in the catalytic loops on complexation reflect the special features of their structure in the mycobacterial protein. A comparative analysis of available sequences of the enzyme from different sources indicates high conservation of amino-acid residues in the catalytic loops. The uracil-binding pocket in the structure is occupied by a citrate ion. The interactions of the citrate ion with the protein mimic those of uracil, in addition to providing insights into other possible interactions that inhibitors could be involved in.