Structure of bovine prothrombin fragment 1 refined at 2.25 Å resolution

The structure of bovine prothrombin fragment 1 has been refined at 2.25 Å resolution using high resolution measurements made with the synchrotron beam at CHESS. The synchrotron data were collected photographically by oscillation methods (R-merge = 0.08). These were combined with lower order diffractometer data for refinement purposes. The structure was refined using restrained least-squares methods with the program PROLSQ to a crystallographic R-value of 0.175. The structure includes 105 water molecules with occupancies of >0·6. The first 35 residues (Ala1-Leu35) of the N-terminal ?-carboxy glutamic acid-domain (Ala1-Cys48) of fragment 1 are disordered as are two carbohydrate chains of Mr ? 5000; the latter two combine to render 40% of the structure disordered. The folding of the kringle of fragment 1 is related to the close intramolecular contact between the inner loop disulfide groups. Half of the conserved sequence of the kringle forms an inner core surrounding these disulfide groups. The remainder of the sequence conservation is associated with the many turns of the main chain. The Pro95 residue of the kringle has a cis conformation and Tyr74 is ordered in fragment 1, although nuclear magnetic resonance studies indicate that the comparable residue of plasminogen kringle 4 has two positions. Surface accessibility calculations indicate that none of the disulfide groups of fragment 1 is accessible to solvent.

A Deterministic Optimization Approach to Protein Sequence Design Using Continuous Models

Determining the sequence of amino acid residues in a heteropolymer chain of a protein with a given conformation is a discrete combinatorial problem that is not generally amenable for gradient-based continuous optimization algorithms. In this paper we present a new approach to this problem using continuous models. In this modeling, continuous "state functions" are proposed to designate the type of each residue in the chain. Such a continuous model helps define a continuous sequence space in which a chosen criterion is optimized to find the most appropriate sequence. Searching a continuous sequence space using a deterministic optimization algorithm makes it possible to find the optimal sequences with much less computation than many other approaches. The computational efficiency of this method is further improved by combining it with a graph spectral method, which explicitly takes into account the topology of the desired conformation and also helps make the combined method more robust. The continuous modeling used here appears to have additional advantages in mimicking the folding pathways and in creating the energy landscapes that help find sequences with high stability and kinetic accessibility. To illustrate the new approach, a widely used simplifying assumption is made by considering only two types of residues: hydrophobic (H) and polar (P). Self-avoiding compact lattice models are used to validate the method with known results in the literature and data that can be practically obtained by exhaustive enumeration on a desktop computer. We also present examples of sequence design for the HP models of some real proteins, which are solved in less than five minutes on a single-processor desktop computer Some open issues and future extensions are noted.

Functional Mimics of Glutathione Peroxidase: Bioinspired Synthetic Antioxidants

Oxidative stress is caused by an imbalance between the production of reactive oxygen species (ROS) and the biological system's ability to detoxify these reactive intermediates. Mammalian cells have elaborate antioxidant defense mechanisms to control the damaging effects of ROS. Glutathione peroxidase (GPx), a selenoenzyme, plays a key role in protecting the organism from oxidative damage by catalyzing the reduction of harmful hydroperoxides with thiol a ``catalytic triad'' with tryptophan and glutamine, which cofactors. The selenocysteine residue at the active site forms activates the selenium moiety for an efficient reduction of peroxides. After the discovery that ebselen, a synthetic organoselenium compound, mimics the catalytic activity of GPx both in vitro and in vivo, several research groups developed a number of small-molecule selenium compounds as functional mimics of GPx, either by modifying the basic structure of ebselen or by incorporating some structural features of the native enzyme. The synthetic mimics reported in the literature can be classified in three major categories: (i) cyclic selenenyl amides having a Se-N bond, (ii) diaryl diselenides, and (iii) aromatic or aliphatic monoselenides. Recent studies show that ebselen exhibits very poor GPx activity when aryl or benzylic thiols such as PhSH or BnSH are used as cosubstrates. Because the catalytic activity of each GPx mimic largely depends on the thiol cosubstrates used, the difference in the thiols causes the discrepancies observed in different studies. In this Account, we demonstrate the effect of amide and amine substituents on the GPx activity of various organoselenium compounds. The existence of strong Se ... O/N interactions in the selenenyl sulfide intermediates significantly reduces the GPx activity. These interactions facilitate an attack of thiol at selenium rather than at sulfur, leading to thiol exchange reactions that hamper the formation of catalytically active selenol. Therefore, any substituent capable of enhancing the nucleophilic attack of thiol at sulfur in the selenenyl sulfide state would enhance the antioxidant potency of organoselenium compounds. Interestingly, replacement of the sec-amide substituent by a tert-amide group leads to a weakening of Se ... 0 interactions in the selenenyl sulfide intermediates. This modification results in 10- to 20-fold enhancements in the catalytic activities. Another strategy involving the replacement of tert-amide moieties by tert-amino substituents further increases the activity by 3- to 4-fold. The most effective modification so far in benzylamine-based GPx mimics appears to be either the replacement of a tert-amino substituent by a sec-amino group or the introduction of an additional 6-methoxy group in the phenyl ring. These strategies can contribute to a remarkable enhancement in the GPx activity. In addition to enhancing catalytic activity, a change in the substituents near the selenium moiety alters the catalytic mechanisms. The mechanistic investigations of functional mimics are useful not only for understanding the complex chemistry at the active site of GPx but also for designing and synthesizing novel antioxidants and anti-inflammatory agents.

HELANAL: A program to characterise helix geometry in proteins,

A detailed analysis of structural and position dependent characteristic features of helices will give a better understanding of the secondary structure formation in globular proteins. Here we describe an algorithm that quantifies the geometry of helices in proteins on the basis of their C-alpha atoms alone. The Fortran program HELANAL can extract the helices from the PDB files and then characterises the overall geometry of each helix as being linear, curved or kinked, in terms of its local structural features, viz. local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes. Even helices with large radius of curvature are unambiguously identified as being linear or curved. The program can also be used to differentiate a kinked helix and other motifs, such as helix-loop-helix or a helix-turn-helix (with a single residue linker) with the help of local bending angles. In addition to these, the program can also be used to characterise the helix start and end as well as other types of secondary structures.

Error-free matrix symmetrizers and equivalent symmetric matrices

A symmetrizer of the matrix A is a symmetric solution X that satisfies the matrix equation XA=AprimeX. An exact matrix symmetrizer is computed by obtaining a general algorithm and superimposing a modified multiple modulus residue arithmetic on this algorithm. A procedure based on computing a symmetrizer to obtain a symmetric matrix, called here an equivalent symmetric matrix, whose eigenvalues are the same as those of a given real nonsymmetric matrix is presented.

Conformational studies on beta-bend containing a cis peptide unit

Conformational studies have been carried out on the X-cis-Pro tripeptide system (a system of three linked peptide units, in the trans-cis-trans configuration) using energy minimization techniques. For X, residues Gly, L-Ala, D-Ala and L-Pro have been used. The energy minima have been classified into different groups based upon the conformational similarity. There are 15, 20, 18 and 6 minima that are possible for the four cases respectively and these fall into 11 different groups. A study of these minima shows that, (i) some minima contain hydrogen bonds - either 4-->1 or 1-->2 type, (ii) the low energy minima qualify themselves as bend conformations, (iii) cis' and trans' conformations are possible for the prolyl residue as also the C(gamma)-endo and C(gamma)-exo puckerings, and (iv) for Pro-cis-Pro, cis' at the first prolyl residue is ruled out, due to the high energy. The available crystal structure data on proteins and peptides, containing cis-Pro segment have been examined with a view to find the minima that occur in solid state. The data from protein show that they fall under two groups. The conformation at X in X-cis-Pro is near extended when it is a non-glycyl residue. In both peptides and proteins there exists a preference for trans' conformation at prolyl residue over cis' when X is a non-glycyl residue. The minima obtained can be useful in modelling studies.

(2E)-3-(4-Ethoxyphenyl)-1-(2-methyl-4-phenylquinolin-3-yl)prop-2-en-1-on e monohydrate

The title hydrate, C27H23NO2 center dot H2O, features an almost planar quinoline residue (r.m.s. deviation = 0.015 angstrom) with the benzene dihedral angle = 63.80 (7)degrees] and chalcone C-C-C-O torsion angle = -103.38 (18)degrees] substituents twisted significantly out of its plane. The configuration about the C=C bond 1.340 (2) angstrom] is E. In the crystal, molecules related by the 21 symmetry operation are linked along the b axis via water molecules that form O-H center dot center dot center dot O-c and O-H center dot center dot center dot N-q hydrogen bonds (c = carbonyl and q = quinoline). A C-H center dot center dot center dot O interaction also occurs.

Effect of constraints by threonine on proline containing αa-helix—A molecular dynamics approach

Proline plays an important role in the secondary structure of proteins. In the pursuit of understanding its structural role, Proline containing helices with constraints have been studied by employing molecular dynamics (MD) technique. In the present study, the constraint introduced is a threonine residue, whose sidechain has intramolecular hydrogen bond interaction with the backbone oxygen atom. The three systems that have been chosen for characterization are: (1) Ace-(Ala)12−Thr-Pro-(Ala)10−NHMe, (2) Ace-(Ala)13-Pro-Ala-Thr- (Ala)8-NHMe and (3) Ace-(Ala)13-Pro-(Ala)3-Thr-(Ala)6-NHMe. The equilibrium structures and structural transitions have been identified by monitoring the backbone dihedral angles, bend related parameters and the hydrogen bond interactions. The MD averages and root mean square (r.m.s.) fluctuations are compared and discussed. Energy minimization has been carried out on selected MD simulated points in order to analyze the characteristics of different conformations.

Identification of Critical Residues of the Mycobacterial Dephosphocoenzyme A Kinase by Site-Directed Mutagenesis

Dephosphocoenzyme A kinase performs the transfer of the c-phosphate of ATP to dephosphocoenzyme A, catalyzing the last step of coenzyme A biosynthesis. This enzyme belongs to the P-loop-containing NTP hydrolase superfamily, all members of which posses a three domain topology consisting of a CoA domain that binds the acceptor substrate, the nucleotide binding domain and the lid domain. Differences in the enzymatic organization and regulation between the human and mycobacterial counterparts, have pointed out the tubercular CoaE as a high confidence drug target (HAMAP database). Unfortunately the absence of a three-dimensional crystal structure of the enzyme, either alone or complexed with either of its substrates/regulators, leaves both the reaction mechanism unidentified and the chief players involved in substrate binding, stabilization and catalysis unknown. Based on homology modeling and sequence analysis, we chose residues in the three functional domains of the enzyme to assess their contributions to ligand binding and catalysis using site-directed mutagenesis. Systematically mutating the residues from the P-loop and the nucleotide-binding site identified Lys14 and Arg140 in ATP binding and the stabilization of the phosphoryl intermediate during the phosphotransfer reaction. Mutagenesis of Asp32 and Arg140 showed catalytic efficiencies less than 5-10% of the wild type, indicating the pivotal roles played by these residues in catalysis. Non-conservative substitution of the Leu114 residue identifies this leucine as the critical residue from the hydrophobic cleft involved in leading substrate, DCoA binding. We show that the mycobacterial enzyme requires the Mg2+ for its catalytic activity. The binding energetics of the interactions of the mutant enzymes with the substrates were characterized in terms of their enthalpic and entropic contributions by ITC, providing a complete picture of the effects of the mutations on activity. The properties of mutants defective in substrate recognition were consistent with the ordered sequential mechanism of substrate addition for CoaE.

Amino acid selective unlabeling for sequence specific resonance assignments in proteins

Sequence specific resonance assignment constitutes an important step towards high-resolution structure determination of proteins by NMR and is aided by selective identification and assignment of amino acid types. The traditional approach to selective labeling yields only the chemical shifts of the particular amino acid being selected and does not help in establishing a link between adjacent residues along the polypeptide chain, which is important for sequential assignments. An alternative approach is the method of amino acid selective `unlabeling' or reverse labeling, which involves selective unlabeling of specific amino acid types against a uniformly C-13/N-15 labeled background. Based on this method, we present a novel approach for sequential assignments in proteins. The method involves a new NMR experiment named, {(CO)-C-12 (i) -N-15 (i+1)}-filtered HSQC, which aids in linking the H-1(N)/N-15 resonances of the selectively unlabeled residue, i, and its C-terminal neighbor, i + 1, in HN-detected double and triple resonance spectra. This leads to the assignment of a tri-peptide segment from the knowledge of the amino acid types of residues: i - 1, i and i + 1, thereby speeding up the sequential assignment process. The method has the advantage of being relatively inexpensive, applicable to H-2 labeled protein and can be coupled with cell-free synthesis and/or automated assignment approaches. A detailed survey involving unlabeling of different amino acid types individually or in pairs reveals that the proposed approach is also robust to misincorporation of N-14 at undesired sites. Taken together, this study represents the first application of selective unlabeling for sequence specific resonance assignments and opens up new avenues to using this methodology in protein structural studies.

Analysis of the effects of amino-acid-sequence on the structure of proteins

The conformation of amino acid side chains as observed in well-determined structures of globular proteins has earlier been extensively investigated. In contrast, the structural features of the polypeptide backbone that result from the occurrence of specific amino acids along the polypeptide have not been analysed. In this article, we present the statistically significant features in the backbone geometry that appear to be a consequence of the occurrence of rotamers of different amino acid side chains by analysing 102 well-refined structures that form a random collection of proteins. It is found that the persistence of helical segments around each residue is influenced by the residue type. Several residues exert asymmetrical influence between the carboxyl and amino terminal polypeptide segments. The degree to which secondary structures depart from an average geometry also appears to depend on residue type. These departures are correlated to the corresponding Chou and Fasman parameters of amino acid residues. The frequency distribution of the side chain rotamers is influenced by polypeptide secondary structure. In turn, the rotamer conformation of side chain affects the extension of the secondary structure of the backbone. The strongest correlation is found between the occurrence of g+ conformation and helix propagation on the carboxyl side of many residues.

Diversity in Functional Organization of Class I and Class II Biotin Protein Ligase

The cell envelope of Mycobacterium tuberculosis (M. tuberculosis) is composed of a variety of lipids including mycolic acids, sulpholipids, lipoarabinomannans, etc., which impart rigidity crucial for its survival and pathogenesis. Acyl CoA carboxylase (ACC) provides malonyl-CoA and methylmalonyl-CoA, committed precursors for fatty acid and essential for mycolic acid synthesis respectively. Biotin Protein Ligase (BPL/BirA) activates apo-biotin carboxyl carrier protein (BCCP) by biotinylating it to an active holo-BCCP. A minimal peptide (Schatz), an efficient substrate for Escherichia coli BirA, failed to serve as substrate for M. tuberculosis Biotin Protein Ligase (MtBPL). MtBPL specifically biotinylates homologous BCCP domain, MtBCCP87, but not EcBCCP87. This is a unique feature of MtBPL as EcBirA lacks such a stringent substrate specificity. This feature is also reflected in the lack of self/promiscuous biotinylation by MtBPL. The N-terminus/HTH domain of EcBirA has the selfbiotinable lysine residue that is inhibited in the presence of Schatz peptide, a peptide designed to act as a universal acceptor for EcBirA. This suggests that when biotin is limiting, EcBirA preferentially catalyzes, biotinylation of BCCP over selfbiotinylation. R118G mutant of EcBirA showed enhanced self and promiscuous biotinylation but its homologue, R69A MtBPL did not exhibit these properties. The catalytic domain of MtBPL was characterized further by limited proteolysis. Holo-MtBPL is protected from proteolysis by biotinyl-59 AMP, an intermediate of MtBPL catalyzed reaction. In contrast, apo-MtBPL is completely digested by trypsin within 20 min of co-incubation. Substrate selectivity and inability to promote self biotinylation are exquisite features of MtBPL and are a consequence of the unique molecular mechanism of an enzyme adapted for the high turnover of fatty acid biosynthesis.

Effect of Substrate Binding Loop Mutations on the Structure, Kinetics, and Inhibition of Enoyl Acyl Carrier Protein Reductase from Plasmodium falciparum

Enoyl acyl carrier protein reductase (ENR), which catalyzes the final and rate limiting step of fatty acid elongation, has been validated as a potential drug target. Triclosan is known to be an effective inhibitor for this enzyme. We mutated the substrate binding site residue Ala372 of the ENR of Plasmodium falciparum (PfENR) to Methionine and Valine which increased the affinity of the enzyme towards triclosan to almost double, close to that of Escherichia coli ENR (EcENR) which has a Methionine at the structurally similar position of Ala372 of PfENR. Kinetic studies of the mutants of PfENR and the crystal structure analysis of the A372M mutant revealed that a more hydrophobic environment enhances the affinity of the enzyme for the inhibitor. A triclosan derivative showed a threefold increase in the affinity towards the mutants compared to the wild type, due to additional interactions with the A372M mutant as revealed by the crystal structure. The enzyme has a conserved salt bridge which stabilizes the substrate binding loop and appears to be important for the active conformation of the enzyme. We generated a second set of mutants to check this hypothesis. These mutants showed loss of function, except in one case, where the crystal structure showed that the substrate binding loop is stabilized by a water bridge network. (C) 2011 IUBMB mum Life, 63(1): 30-41,2011

Modification of the sugar specificity of a plant lectin: structural studies on a point mutant of Erythrina corallodendron lectin

A mutant of Erythrina corallodendron lectin was generated with the aim of enhancing its affinity for N-acetylgalactosamine. A tyrosine residue close to the binding site of the lectin was mutated to a glycine in order to facilitate stronger interactions between the acetamido group of the sugar and the lectin which were prevented by the side chain of the tyrosine in the wild-type lectin. The crystal structures of this Y106G mutant lectin in complex with galactose and N-acetylgalactosamine have been determined. A structural rationale has been provided for the differences in the relative binding affinities of the wild-type and mutant lectins towards the two sugars based on the structures. A hydrogen bond between the O6 atom of the sugars and the variable loop of the carbohydrate-binding site of the lectin is lost in the mutant complexes owing to a conformational change in the loop. This loss is compensated by an additional hydrogen bond that is formed between the acetamido group of the sugar and the mutant lectin in the complex with N-acetylgalactosamine, resulting in a higher affinity of the mutant lectin for N-acetylgalactosamine compared with that for galactose, in contrast to the almost equal affinity of the wild-type lectin for the two sugars. The structure of a complex of the mutant with a citrate ion bound at the carbohydrate-binding site that was obtained while attempting to crystallize the complexes with sugars is also presented.

Effect of the valine-threonine constraint on the dynamics of the proline helix — A molecular dynamics study

Proline residues in helices play an important role in the structure of proteins. The proline residue introduces a kink in the helix which varies from about 5-degrees to 50-degrees. The presence of other residues such as threonine or valine near the proline region can influence the flexibility exhibited by the kinked helix, which can have an important biological role. In the present paper, the constraint introduced by threonine and valine on a proline helix is investigated by molecular dynamics studies. The systems considered am (1) a poly-alanine helix with threonine-proline residues (TP) and (2) a poly-alanine helix with valine-threonine-proline residues (VTP), in the middle. Molecular dynamics simulations are carried out on these two systems for 500 ps. The results are analyzed in terms of structural transitions, bend-related parameters and sidechain orientations.