A Search for Energy Minimized Sequences of Proteins

In this paper, we present numerical evidence that supports the notion of minimization in the sequence space of proteins for a target conformation. We use the conformations of the real proteins in the Protein Data Bank (PDB) and present computationally efficient methods to identify the sequences with minimum energy. We use edge-weighted connectivity graph for ranking the residue sites with reduced amino acid alphabet and then use continuous optimization to obtain the energy-minimizing sequences. Our methods enable the computation of a lower bound as well as a tight upper bound for the energy of a given conformation. We validate our results by using three different inter-residue energy matrices for five proteins from protein data bank (PDB), and by comparing our energy-minimizing sequences with 80 million diverse sequences that are generated based on different considerations in each case. When we submitted some of our chosen energy-minimizing sequences to Basic Local Alignment Search Tool (BLAST), we obtained some sequences from non-redundant protein sequence database that are similar to ours with an E-value of the order of 10(-7). In summary, we conclude that proteins show a trend towards minimizing energy in the sequence space but do not seem to adopt the global energy-minimizing sequence. The reason for this could be either that the existing energy matrices are not able to accurately represent the inter-residue interactions in the context of the protein environment or that Nature does not push the optimization in the sequence space, once it is able to perform the function.

Molecular characterisation of ABC transporter type FtsE and FtsX proteins of Mycobacterium tuberculosis.

Elicitation of drug resistance and various survival strategies inside host macrophages have been the hallmarks of Mycobacterium tuberculosis as a successful pathogen. ATP Binding Cassette (ABC) transporter type proteins are known to be involved in the efflux of drugs in bacterial and mammalian systems. FtsE, an ABC transporter type protein, in association with the integral membrane protein FtsX, is involved in the assembly of potassium ion transport proteins and probably of cell division proteins as well, both of which being relevant to tubercle bacillus. In this study, we cloned ftsE gene of M. tuberculosis, overexpressed and purified. The recombinant MtFtsE-6xHis protein and the native MtFtsE protein were found localized on the membrane of E. coli and M. tuberculosis cells, respectively. MtFtsE-6xHis protein showed ATP binding in vitro, for which the K42 residue in the Walker A motif was found essential. While MtFtsE-6xHis protein could partially complement growth defect of E. coli ftsE temperature-sensitive strain MFT1181, co-expression of MtFtsE and MtFtsX efficiently complemented the growth defect, indicating that the MtFtsE and MtFtsX proteins might be performing an associated function. MtFtsE and MtFtsX-6xHis proteins were found to exist as a complex on the membrane of E. coli cells co-expressing the two proteins.

The crystal and molecular structure of benzyloxycarbonyl-a-amino-isobutyryl-L-prolyl-methylamide. The observation of an X2-Pro3 Type III b-turn

The crystal and molecular structure of N-benzyloxycarbonyl-a-aminoisobutyryl-L-prolyl methylamide, the amino terminal dipeptide fragment of alamethicin, has been determined using direct methods. The compound crystallizes in the orthorhombic system with the space group P212-21. Cell dimensions are a = 7.705 A, b = 11.365 A, and c = 21.904 A. The structure has been refined using conventional procedures to a final R factor of 0.054. The molecular structure possesses a 4 - 1 intramolecular N-H--0 hydrogen bond formed between the CO group of the urethane moiety and the NH group of the methylamide function. The peptide backbone adopts the type 111 P-turn conformation, with 42 = -51.0°, +* = -39.7",&j = -65.0', $3 = -25.4'. An unusual feature is the occurrence of the proline residue at position 3 of the P-turn. The observed structure supports the view that Aib residues initiate the formation of type 111 @-turn conformations. The pyrrolidine ring is puckered in Cy-exo fashion.

An incipient 310 helix in Piv-Pro-Pro-Ala-NHMe as a model for peptide folding

The molecular mechanism of helix nucleation in peptides and proteins is not yet understood and the question of whether sharp turns in the polypeptide backbone serve as nuclei for protein folding has evoked controversy1,2. A recent study of the conformation of a tetrapeptide containing the stereochemically constrained residue alpha-aminoisobutyric acid, both in solution and the solid state, yielded a structure consisting of two consecutive beta-turns, leading to an incipient 310 helical conformation3,4. This led us to speculate that specific tri- and tetra-peptide sequences may indeed provide a helical twist to the amino-terminal segment of helical regions in proteins and provide a nucleation site for further propagation. The transformation from a 310 helical structure to an alpha-helix should be facile and requires only small changes in the phi and psi conformational angles and a rearrangement of the hydrogen bonding pattern5. If such a mechanism is involved then it should be possible to isolate an incipient 310 helical conformation in a tripeptide amide or tetrapeptide sequence, based purely on the driving force derived from short-range interactions. We have synthesised and studied the model peptide pivaloyl-Pro-Pro-Ala-NHMe (compound I) and provide here spectroscopic evidence for a 310 helical conformation in compound I.

An Insight to the Dynamics of Conserved Water Molecular Triad in IMPDH II (Human): Recognition of Cofactor and Substrate to Catalytic Arg 322

Inosine 5' monophosphate dehydrogenase (IMPDH II) is a key enzyme involved in the de novo biosynthesis pathway of purine nucleotides and is also considered to be an excellent target for cancer inhibitor design. The conserve R 322 residue (in human) is thought to play some role in the recognition of inhibitor and cofactor through the catalytic D 364 and N 303. The 15 ns simulation and the water dynamics of the three different PDB structures (1B3O, 1NF7, and 1NFB) of human IMPDH by CHARMM force field have clearly indicated the involvement of three conserved water molecules (W-L, W-M, and W-C) in the recognition of catalytic residues (R 322, D 364, and N 303) to inhibitor and cofactor. Both the guanidine nitrogen atoms (NH1 and NH 2) of the R 322 have anchored the di- and mono-nucleotide (cofactor and inhibitor) binding domains via the conserved W-C and W-L water molecules. Another conserved water molecule W-M seems to bridge the two domains including the R 322 and also the W-C and W-L through seven centers H-bonding coordination. The conserved water molecular triad (W-C - W-M - W-L) in the protein complex may thought to play some important role in the recognition of inhibitor and cofactor to the protein through R 322 residue.

Expanding the Peptide beta-Turn in alpha gamma Hybrid Sequences: 12 Atom Hydrogen Bonded Helical and Hairpin Turns

Hybrid peptide segments containing contiguous alpha and gamma amino acid residues can form C-12 hydrogen bonded turns which may be considered as backbone expanded analogues of C-10 beta-turns) found in alpha alpha segments. Exploration of the regular hydrogen bonded conformations accessible for hybrid alpha gamma sequences is facilitated by the use of a stereochemically constrained gamma amino acid residue gabapentin (1-aminomethylcyclohexaneacetic acid, Gpn), in which the two torsion angles about C-gamma-C-beta (theta(1)) and C-beta-C-alpha (theta(2)) are predominantly restricted to gauche conformations. The crystal structures of the octapeptides Boc-Gpn-Aib-Gpn-Aib-Gpn-Aib-Gpn-Aib-OMe (1) and Boc-Leu-Phe-Val-Aib-Gpn-Leu-Phe-Val-OMe (2) reveal two distinct conformations for the Aib-Gpn segment. Peptide 1 forms a continuous helix over the Aib(2)-Aib(6) segment, while the peptide 2 forms beta-hairpin structure stabilized by four cross-strand hydrogen bonds with the Aib-Gpn segment forming a nonhelical C-12 turn. The robustness of the helix in peptide 1 in solution is demonstrated by NMR methods. Peptide 2 is conformationally fragile in solution with evidence of beta-hairpin conformations being obtained in methanol. Theoretical calculations permit delineation of the various C-12 hydrogen bonded structures which are energetically feasible in alpha gamma and gamma alpha sequences.

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I-turn forming tetrapeptide and an -helical tetradecapeptide

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt -helical conformations stabilized by 11 successive 5 1 hydrogen bonds. In addition, a single 4 1 hydrogen bond is also observed at the N-terminus. All five Dpg residues adopt backbone torsion angles (, ) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle NCC() and the observed backbone , values. For > 106° , helices are observed, while fully extended structures are characterized by < 106° . The mean values for extended and folded conformations for the Dpg residue are 103.6° ± 1.7° and 109.9° ± 2.6° , respectively.

Cloning, overexpression, folding and purification of a biosynthetically derived three disulfide scorpion toxin (BTK-2) from Mesobuthus tamulus

BTK-2, a 32 residue scorpion toxin initially identified in the venom of red Indian scorpion Mesobuthus tamulus was cloned, overexpressed and purified using Cytochrome 155 fusion protein system developed in our laboratory. The synthetic gene coding for the peptide was designed taking into account optimal codon usage by Escherichia coli. High expression levels of the fusion protein enabled facile purification of this peptide. The presence of disulfide bonded isomers, occurring as distinctly populated states even in the fusion protein, were separated by gel filtration chromatography. The target peptide was liberated from the host protein by Tev protease cleavage and subsequent purification was achieved using RP-HPLC methods. Reverse phase HPLC clearly showed the presence of at least two isomeric forms of the peptide that were significantly populated. The oxidative folding of BTK-2 was achieved under ambient conditions during the course of purification. Structural characterization of the two forms, by solution homonuclear and heteronuclear NMR methods, has shown that these two forms exhibit significantly different structural properties, and represent the natively folded and a "misfolded" form of the peptide. The formation of properly folded BTK-2 as a major fraction without the use of in vitro oxidative refolding methods clearly indicate the versatility of the Cytochrome b(5) fusion protein system for the efficient production of peptides for high resolution NMR studies.

Quantitation and characterization of glutathionyl haemoglobin as an oxidative stress marker in chronic renal failure by mass spectrometry

Objectives: Glutathionyl haemoglobin (GS-Hb) belonging to the class of glutathionylated proteins has been investigated as a possible marker of oxidative stress in different chronic diseases. The purpose of this study was to examine whether glutathionyl haemoglobin can serve as an oxidative stress marker in non-diabetic chronic renal failure patients on different renal replacement therapies (RRT) through its quantitation, and characterization of the specific binding site of glutathione in haemoglobin molecule by mass spectrometric analysis. Design and methods: The study group consisted of non-diabetic chronic renal failure patients on renal replacement therapy (RRT): hemodialysis (HD), continuous ambulatory peritoneal dialysis (CAPD) and renal allograft transplant (Txp) patients. Haemoglobin samples of these subjects were analyzed by liquid chromatography electrospray ionization mass spectrometry for GS-Hb quantitation. Characterization of GS-Hb was done by tandem mass spectrometry. Levels of erythrocyte glutathione (GSH) and lipid peroxidation (as thiobarbituric acid reacting substances) were measured spectrophotometrically, while glycated baernoglobin (HbA1c) was measured by HPLC. Results: GS-Hb levels were markedly elevated in the dialysis group and marginally in the transplant group as compared to the controls. GS-Hb levels correlated positively with lipid peroxidation and negatively with the erythrocyte glutathione levels in RRT groups indicating enhanced oxidative stress. De novo sequencing of the chymotryptic fragment of GS-Hb established that glutathione is attached to Cys-93 of the beta globin chain. Mass spectrometric quantitation of total glycated haemoglobin showed good agreement with HbA1c estimation by conventional HPLC method. Conclusions: Glutathionyl haemoglobin can serve as a clinical marker of oxidative stress in chronic debilitating therapies like RRT. Mass spectrometry provides a reliable analytical tool for quantitation and residue level characterization of different post-translational modifications of haemoglobin. (c) 2007 The Canadian Society of Clinical Chemists. Published by Elsevier Inc. All rights reserved.

Role of histidine residues in EcoP15I DNA methyltransferase activity as probed by chemical modification and site-directed mutagenesis

Towards understanding the catalytic mechanism of M.EcoP15I [EcoP15I MTase (DNA methyltransferase); an adenine methyltransferase], we investigated the role of histidine residues in catalysis. M.EcoP15I, when incubated with DEPC (diethyl pyrocarbonate), a histidine-specific reagent, shows a time- and concentration-dependent inactivation of methylation of DNA containing its recognition sequence of 5'-CAGCAG-3'. The loss of enzyme activity was accompanied by an increase in absorbance at 240 nm. A difference spectrum of modified versus native enzyme shows the formation of N-carbethoxyhistidine that is diminished by hydroxylamine. This, along with other experiments, strongly suggests that the inactivation of the enzyme by DEPC was specific for histidine residues. Substrate protection experiments show that pre-incubating the methylase with DNA was able to protect the enzyme from DEPC inactivation. Site-directed mutagenesis experiments in which the 15 histidine residues in the enzyme were replaced individually with alanine corroborated the chemical modification studies and established the importance of His-335 in the methylase activity. No gross structural differences were detected between the native and H335A mutant MTases, as evident from CD spectra, native PAGE pattern or on gel filtration chromatography. Replacement of histidine with alanine residue at position 335 results in a mutant enzyme that is catalytically inactive and binds to DNA more tightly than the wild-type enzyme. Thus we have shown in the present study, through a combination of chemical modification and site-directed mutagenesis experiments, that His-335 plays an essential role in DNA methylation catalysed by M.EcoP15I.

Structures and molecular-dynamics studies of three active-site mutants of bovine pancreatic phospholipase A(2)

Phospholipase A(2) hydrolyzes phospholipids at the sn-2 position to cleave the fatty-acid ester bond of L-glycerophospholipids. The catalytic dyad (Asp99 and His48) along with a nucleophilic water molecule is responsible for enzyme hydrolysis. Furthermore, the residue Asp49 in the calcium-binding loop is essential for controlling the binding of the calcium ion and the catalytic action of phospholipase A2. To elucidate the structural role of His48 and Asp49, the crystal structures of three active-site single mutants H48N, D49N and D49K have been determined at 1.9 angstrom resolution. Although the catalytically important calcium ion is present in the H48N mutant, the crystal structure shows that proton transfer is not possible from the catalytic water to the mutated residue. In the case of the Asp49 mutants, no calcium ion was found in the active site. However, the tertiary structures of the three active-site mutants are similar to that of the trigonal recombinant enzyme. Molecular-dynamics simulation studies provide a good explanation for the crystallographic results.

Systematic study on crystal-contact engineering of diphthine synthase: influence of mutations at crystal-packing regions on X-ray diffraction quality

It is well known that protein crystallizability can be influenced by site-directed mutagenesis of residues on the molecular surface of proteins, indicating that the intermolecular interactions in crystal-packing regions may play a crucial role in the structural regularity at atomic resolution of protein crystals. Here, a systematic examination was made of the improvement in the diffraction resolution of protein crystals on introducing a single mutation of a crystal-packing residue in order to provide more favourable packing interactions, using diphthine synthase from Pyrococcus horikoshii OT3 as a model system. All of a total of 21 designed mutants at 13 different crystal-packing residues yielded almost isomorphous crystals from the same crystallization conditions as those used for the wild-type crystals, which diffracted X-rays to 2.1 angstrom resolution. Of the 21 mutants, eight provided crystals with an improved resolution of 1.8 angstrom or better. Thus, it has been clarified that crystal quality can be improved by introducing a suitable single mutation of a crystal-packing residue. In the improved crystals, more intimate crystal-packing interactions than those in the wild-type crystal are observed. Notably, the mutants K49R and T146R yielded crystals with outstandingly improved resolutions of 1.5 and 1.6 angstrom, respectively, in which a large-scale rearrangement of packing interactions was unexpectedly observed despite the retention of the same isomorphous crystal form. In contrast, the mutants that provided results that were in good agreement with the designed putative structures tended to achieve only moderate improvements in resolution of up to 1.75 angstrom. These results suggest a difficulty in the rational prediction of highly effective mutations in crystal engineering.

Self-acylation properties of type II fatty acid biosynthesis acyl carrier protein

Acyl carrier protein (ACIP) plays a central role in many metabolic processes inside the cell, and almost 4% of the total enzymes inside the cell require it as a cofactor. Here, we report self-acylation properties in ACPs from Plasmodium falciparum and Brassica napus that are essential components of type II fatty acid biosynthesis (FAS II), disproving the existing notion that this phenomenon is restricted only to ACPs involved in polyketide biosynthesis. We also provide strong evidence to suggest that catalytic self-acylation is intrinsic to the individual ACP. Mutational analysis of these ACPs revealed the key residue(s) involved in this phenomenon. We also demonstrate that these FAS 11 ACPs exhibit a high degree of selectivity for self-acylation employing only dicarboxylic acids as substrates. A plausible mechanism for the self-acylation reaction is also proposed.

Disulfide conformation and design at helix N-termini

To understand structural and thermodynamic features of disulfides within an alpha-helix, a non-redundant dataset comprising of 5025 polypeptide chains containing 2311 disulfides was examined. Thirty-five examples were found of intrahelical disulfides involving a CXXC motif between the N-Cap and third helical positions. GLY and PRO were the most common amino acids at positions 1 and 2, respectively. The N-Cap residue for disulfide bonded CXXC motifs had average values of (-112 +/- 25.2 degrees, 106 +/- 25.4 degrees). To further explore conformational requirements for intrahelical disulfides, CYS pairs were introduced at positions N-Cap-3; 1,4; 7,10 in two helices of an Escherichia coli thioredoxin mutant lacking its active site disulfide (nSS Trx). In both helices, disulfides formed spontaneously during purification only at positions N-Cap-3. Mutant stabilities were characterized by chemical denaturation studies (in both oxidized and reduced states) and differential scanning calorimetry (oxidized state only). All oxidized as well as reduced mutants were destabilized relative to nSS Trx. All mutants were redox active, but showed decreased activity relative to wild-type thioredoxin. Such engineered disulfides can be used to probe helix start sites in proteins of unknown structure and to introduce redox activity into proteins. Conversely, a protein with CYS residues at positions N-Cap and 3 of an alpha-helix is likely to have redox activity.

Gasifier stoves - science, technology and field outreach

Development of a new class of single pan high efficiency, low emission stoves, named gasifier stoves, that promise constant power that can be controlled using any solid biomass fuel in the form of pellets is reported here. These stoves use battery-run fan-based air supply for gasification (primary air) and for combustion (secondary air).Design with the correct secondary air flow ensures near-stoichiometric combustion that allows attainment of peak combustion temperatures with accompanying high water boiling efficiencies (up to 50% for vessels of practical relevance) and very low emissions (of carbon monoxide, particulate matter and oxides of nitrogen). The use of high density agro-residue based pellets or coconut shell pieces ensures operational duration of about an hour or more at power levels of 3 kWth (similar to 12 g/min). The principles involved and the optimization aspects of the design are outlined. The dependence of efficiency and emissions on the design parameters are described. The field imperatives that drive the choice of the rechargeable battery source and the fan are brought out. The implications of developments of Oorja-Plus and OorjaSuper stoves to the domestic cooking scenario of India are briefly discussed. The process development, testing and internal qualification tasks were undertaken by Indian Institute of Science. Product development and the fuel pellet production were dealt with by First Energy Private Ltd.Close interaction at several times during this period has helped progress the project from the laboratory to large scale commercial operation. At this time, over four hundred thousand stoves and 30 kilotonnes fuel have been sold in four states in India.