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.

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.

Computer simulation of protein—carbohydrate complexes: application to arabinose-binding protein and pea lectin

The CCEM method (Contact Criteria and Energy Minimisation) has been developed and applied to study protein-carbohydrate interactions. The method uses available X-ray data even on the native protein at low resolution (above 2.4 Å) to generate realistic models of a variety of proteins with various ligands.The two examples discussed in this paper are arabinose-binding protein (ABP) and pea lectin. The X-ray crystal structure data reported on ABP-β-l-arabinose complex at 2.8, 2.4 and 1.7 Å resolution differ drastically in predicting the nature of the interactions between the protein and ligand. It is shown that, using the data at 2.4 Å resolution, the CCEM method generates complexes which are as good as the higher (1.7 Å) resolution data. The CCEM method predicts some of the important hydrogen bonds between the ligand and the protein which are missing in the interpretation of the X-ray data at 2.4 Å resolution. The theoretically predicted hydrogen bonds are in good agreement with those reported at 1.7 Å resolution. Pea lectin has been solved only in the native form at 3 Å resolution. Application of the CCEM method also enables us to generate complexes of pea lectin with methyl-α-d-glucopyranoside and methyl-2,3-dimethyl-α-d-glucopyranoside which explain well the available experimental data in solution.

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.

Conformational analysis of the right-hand twisted antiparallel β-Structure

The conformational analysis of a pair of two-linked peptide units in the anti-parallel arrangement is reported here with a view to study the effect of association of one chain with the other. The pair of two-linked peptide units were fixed in space through the hydrogen bonds between them, in accordance with certain hydrogen bond criteria. Model building was undertaken to ascertain whether the proximity of the side-chains could be used to eliminate any one of the right-hand twisted, left-hand twisted or regular β-structures. Stereochemically, it was found possible with all of them. The preference for a right-hand twisted β-structure, however, was indicated by the classical energy calculations. The relevance of the results thus obtained is discussed in the context of the preferential right-hand twist of the β-pleated sheets present in globular proteins. The agreement between the minimum energy conformations obtained for the pair of two-linked peptide units and the globular protein data is also indicated.

Crystal Structure of Staphylococcus aureus Metallopeptidase (Sapep) Reveals Large Domain Motions between the Manganese-bound and Apo-states

Proteases belonging to the M20 family are characterized by diverse substrate specificity and participate in several metabolic pathways. The Staphylococcus aureus metallopeptidase, Sapep, is a member of the aminoacylase-I/M20 protein family. This protein is a Mn2+-dependent dipeptidase. The crystal structure of this protein in the Mn2+-bound form and in the open, metal-free state suggests that large interdomain movements could potentially regulate the activity of this enzyme. We note that the extended inactive conformation is stabilized by a disulfide bond in the vicinity of the active site. Although these cysteines, Cys(155) and Cys(178), are not active site residues, the reduced form of this enzyme is substantially more active as a dipeptidase. These findings acquire further relevance given a recent observation that this enzyme is only active in methicillin-resistant S. aureus. The structural and biochemical features of this enzyme provide a template for the design of novel methicillin-resistant S. aureus-specific therapeutics.

An algorithm to find distant repeats in a pair of protein sequences

Distant repeats between a pair of protein sequences can be exploited to study the various aspects of proteins such as structure-function relationship, disorders due to protein malfunction, evolutionary analysis, etc. An in-depth analysis of the distant repeats would facilitate to establish a stable evolutionary relation of the repeats with respect to their three-dimensional structure. To this effect, an algorithm has been devised to identify the distant repeats in a pair of protein sequences by essentially using the scores of PAM (Percent Accepted Mutation) matrices. The proposed algorithm will be of much use to researchers involved in the comparative study of various organisms based on the amino-acid repeats in protein sequences. (C) 2010 Elsevier B.V. All rights reserved.

Inhibition of peroxidase-catalyzed protein tyrosine nitration by antithyroid drugs and their analogues

In this paper, we describe the effect of some commonly used thiourea-based antithyroid drugs and their analogues on the peroxidase-catalyzed nitration reactions. The nitration of bovine serum albumin (BSA) and cytochrome c was studied using the antibody against 3-nitro-L-tyrosine. This study reveals that the thione-based antithyroid drugs effectively inhibit lactoperoxidase (LPO)-catalyzed nitration of BSA. These compounds show very weak inhibition towards the nitration of cytochrome c. Some of these compounds also inhibit myeloperoxidase (MPO)-catalyzed nitration of L-tyrosine. A structure-activity correlation study on the peroxidase-catalyzed nitration of L-tyrosine reveals that the presence of thione/selone moiety is important for the inhibition. Although the presence of a free N-H group adjacent to C=S moiety is necessary for most of the thiones to inhibit the LPO-catalyzed nitration, the corresponding selones do not require the presence of any free N-H group for their activity. Furthermore, experiments with different concentrations of H2O2 suggest that the antithyroid drugs and related thiones inhibit the nitration reaction mainly by coordinating to the Fe(III)-center of the enzyme active site as previously proposed for the inhibition of peroxidase-catalyzed iodination. On the other hand, the selenium compounds inhibit the nitration by scavenging H2O2 without interacting with the enzyme active site. This assumption is based on the observations that catalase effectively inhibits tyrosine nitration by scavenging H2O2, which is one of the substrates for the nitration. In contrast, superoxide dismutase (SOD) does not alter the nitration reactions, indicating the absence of superoxide radical anion (O-2 center dot(-)) during the peroxidase-catalyzed nitration reactions. (C) 2010 Elsevier B.V. All rights reserved.

A helical peptide containing a majority of valyl residues. The crystal structure of t-butyloxy- carbonyl-(L-valyl-_-aminoisobutyryl)3-L-valyl methyl ester onohydrate

The monohydrate of the heptapeptide t-butyloxycarbonyl-(L-valyl-α-aminoiso-butyryl)3-L-valyl methyl ester crystallizes in the orthorhombic space group P212121 with four molecules in a unit cell with the dimensions α= 9.375, b = 19.413 and c = 25.878 ÅA. The structure has been solved by direct methods and refined to an R value of 0.059 for 3633 observed reflections. The molecule in the structure exists as a slightly distorted 310-helix stabilized by five 4 -> 1 intramolecular hydrogen bonds, indicating the overwhelming influence of α-aminoisobutyryl (Aib) residues in dictating helical fold even when a majority of residues in the peptide have a low intrinsic propensity to be in helices. Contrary to what is expected in helical structures, the valyl side chains, two of which are disordered, exhibit all three possible conformations. The molecules arrange themselves in a head-to-tail fashion along the c-axis. The columns thus generated pack nearly hexagonally in the crystal.

X-ray and molecular-dynamics studies on Mycobacterium leprae single-stranded DNA-binding protein and comparison with other eubacterial SSB structures

The crystal structures of two forms of Mycobacterium leprae single-stranded DNA-binding protein (SSB) have been determined at 2.05 and 2.8 A resolution. Comparison of these structures with the structures of other eubacterial SSBs indicates considerable variation in their quaternary association, although the DNA-binding domains in all of them exhibit the same OB-fold. This variation has no linear correlation with sequence variation, but could be related to variation in protein stability. Molecular-dynamics simulations have been carried out on tetrameric molecules derived from the two forms and the prototype Escherichia coli SSB and the individual subunits of both proteins. Together, the X-ray studies and molecular-dynamics simulations yield information on the relatively rigid and flexible regions of the molecule and on the effect of oligomerization on flexibility. The simulations provide insight into the changes in subunit structure on oligomerization. They also provide insight into the stability and time evolution of the hydrogen bonds/water bridges that connect the two pairs of monomers in the tetramer.

Crystal and molecular structure of l-valyl-l-lysine hydrochloride

l-Valyl-l-lysine hydrochloride, C11N3O3H23 HCl, rystallizes in the monoclinic space group P2, with a = 5.438(5), b = 14.188(5), c = 9.521(5) Å, β= 95.38(2)° and Z = 2. The crystal structure, solved by direct methods, refined to R = 0.036, using full matrix least-squares method. The peptide exists in a zwitterionic form, with the N atom of the lysine side-chain protonated. The two γ-carbons of the valine side-chain have positional disorder, giving rise to two conformations, χ111= -67.3 and 65.9°, one of which (65.9°) is sterically less favourable and has been found to be less popular amongst residues branching at β-C. The lysine side-chain has the geometry of g− tgt, not seen in crystal structures of the dipeptides reported so far. Interestingly, χ32 (63.6°) of lysine side-chain has a gauche+ conformation unlike in most of the other tructures, where it is trans. The neighbouring peptide molecules are hydrogen bonded in a head-to-tail fashion, a rather uncommon interaction in lysine peptide structures. The structure shows considerable similarity with that of l-Lys-l-Val HO in conformational angles and H-bond interactions [4].

Effect of amino acid substitutions at the subunit interface on the stabiity and aggregation properties of a dimeric protein: role of Arg 178 and Arg 218 at the dimer interface of thymidylate synthase

The significance of two interface arginine residues on the structural integrity of an obligatory dimeric enzyme thymidylate synthase (TS) from Lactobacillus casei was investigated by thermal and chemical denaturation. While the R178F mutant showed apparent stability to thermal denaturation by its decreased tendency to aggregate, the Tm of the R218K mutant was lowered by 5 degrees C. Equilibrium denaturation studies in guanidinium chloride (GdmCl) and urea indicate that in both the mutants, replacement of Arg residues results in more labile quaternary and tertiary interactions. Circular dichroism studies in aqueous buffer suggest that the protein interior in R218K may be less well-packed as compared to the wild type protein. The results emphasize that quaternary interactions may influence the stability of the tertiary fold of TS. The amino acid replacements also lead to notable alteration in the ability of the unfolding intermediate of TS to aggregate. The aggregated state of partially unfolded intermediate in the R178F mutant is stable over a narrower range of denaturant concentrations. In contrast, there is an exaggerated tendency on the part of R218K to aggregate in intermediate concentrations of the denaturant. The 3 A crystal structure of the R178F mutant reveals no major structural change as a consequence of amino acid substitution. The results may be rationalized in terms of mutational effects on both the folded and unfolded state of the protein. Site specific amino acid substitutions are useful in identifying specific regions of TS involved in association of non-native protein structures.

Structure of the putative mutarotase YeaD from Salmonella typhimurium in two different forms: in vivo binding of a sugar phosphate

Salmonella typhimurium YeaD (stYeaD), annotated as a putative aldose 1-epimerase, has a very low sequence identity to other well characterized mutarotases. Sequence analysis suggested that the catalytic residues and a few of the substrate-binding residues of galactose mutarotases (GalMs) are conserved in stYeaD. Determination of the crystal structure of stYeaD in an orthorhombic form at 1.9 angstrom resolution and in a monoclinic form at 2.5 angstrom resolution revealed this protein to adopt the beta-sandwich fold similar to GalMs. Structural comparison of stYeaD with GalMs has permitted the identification of residues involved in catalysis and substrate binding. In spite of the similar fold and conservation of catalytic residues, minor but significant differences were observed in the substrate- binding pocket. These analyses pointed out the possible role of Arg74 and Arg99, found only in YeaD-like proteins, in ligand anchoring and suggested that the specificity of stYeaD may be distinct from those of GalMs

Structure-function relationships of icosahedral plant viruses

X-ray diffraction studies on single crystals of a few viruses have led to the elucidation of their three dimensional structure at near atomic resolution. Both the tertiary structure of the coat protein subunit and the quaternary morganization of the icosahedral capsid in these viruses are remarkably similar. These studies have led to a critical re-examination of the structural principles in the architecture of isometric viruses and suggestions of alternative mechanisms of assembly. Apart from their role in the assembly of the virus particle, the coat proteins of certian viruses have been shown to inhibit the replication of the cognate RNA leading to cross-protection. The coat protein amino acid sequence and the genomic sequence of several spherical plant RNA viruses have been determined in the last decade. Experimental data on the mechanisms of uncoating, gene expression and replication of several classes of viruses have also become available. The function of the non-structural proteins of some viruses have been determined. This rapid progress has provided a wealth of information on several key steps in the life cycle of RNA viruses. The function of the viral coat protein, capsid architecture, assembly and disassembly and replication of isometric RNA plant viruses are discussed in the light of this accumulated knowledge.

Communication pathways from the anti-codon region to the aminoacylation site in Methionyl tRNA Synthetase: Molecular dynamics simulations and the structure network analysis

The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by ?70 in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross-correlations between the residues in MetRS have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.