Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Interaction Energy Based Protein Structure Networks

The three dimensional structure of a protein is formed and maintained by the noncovalent interactions among the amino acid residues of the polypeptide chain These interactions can be represented collectively in the form of a network So far such networks have been investigated by considering the connections based on distances between the amino acid residues Here we present a method of constructing the structure network based on interaction energies among the amino acid residues in the protein We have investigated the properties of such protein energy based networks (PENs) and have shown correlations to protein structural features such as the clusters of residues involved in stability formation of secondary and super secondary structural units Further we demonstrate that the analysis of PENs in terms of parameters such as hubs and shortest paths can provide a variety of biologically important information such as the residues crucial for stabilizing the folded units and the paths of communication between distal residues in the protein Finally the energy regimes for different levels of stabilization in the protein structure have clearly emerged from the PEN analysis

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Allostery in tRNA Synthetases Elucidated from MD Simulations and Protein Structure Networks

tRNA synthetases (aaRS) are enzymes crucial in the translation of genetic code. The enzyme accylates the acceptor stem of tRNA by the congnate amino acid bound at the active site, when the anti-codon is recognized by the anti-codon site of aaRS. In a typical aaRS, the distance between the anti-codon region and the amino accylation site is approximately 70 Å. We have investigated this allosteric phenomenon at molecular level by MD simulations followed by the analysis of protein structure networks (PSN) of non-covalent interactions. Specifically, we have generated conformational ensembles by performing MD simulations on different liganded states of methionyl tRNA synthetase (MetRS) from Escherichia coli and tryptophenyl tRNA synthetase (TrpRS) from Human. The correlated residues during the MD simulations are identified by cross correlation maps. We have identified the amino acids connecting the correlated residues by the shortest path between the two selected members of the PSN. The frequencies of paths have been evaluated from the MD snapshots[1]. The conformational populations in different liganded states of the protein have been beautifully captured in terms of network parameters such as hubs, cliques and communities[2]. These parameters have been associated with the rigidity and plasticity of the protein conformations and can be associated with free energy landscape. A comparison of allosteric communication in MetRS and TrpRS [3] elucidated in this study highlights diverse means adopted by different enzymes to perform a similar function. The computational method described for these two enzymes can be applied to the investigation of allostery in other systems.

Photo-induced double-strand DNA and site-specific protein cleavage activity of L-histidine (mu-oxo)diiron(III) complexes of heterocyclic bases

Three oxo-bridged diiron(III) complexes of L-histidine and heterocyclic bases [Fe-2(mu-O)(L-his)(2)(B)(2)](ClO4)(2) (1-3), where B is 2,2'-bipyridine (bpy),1,10-phenanthroline (phen), dipyrido[3,2-d:2',3'-f]quinoxaline (dpq), were prepared and characterized. The bpy complex 1 was structurally characterized by X-ray crystallography. The molecular structure showed a {Fe-2(mu-O)} core in which iron(III) in a FeN4O2 coordination is bound to tridentate monoanionic L-histidine and bidentate bpy ligands. The Fe center dot center dot center dot Fe distance is similar to 3.5 angstrom. The Fe-O-Fe unit is essentially linear, giving a bond angle of similar to 172 degrees. The complexes showed irreversible cyclic voltammetric cathodic response near -0.1 V vs. SCE in H2O-0.1 M KCl. The binuclear units displayed antiferromagnetic interaction between two high-spin (S = 5/2) iron(III) centers giving a -J value of -110 cm(-1). The complexes showed good DNA binding propensity giving a binding constant value of similar to 10(5) M-1. Isothermal titration calorimetric data indicated single binding mode to the DNA. The binding was found to be driven by negative free energy change and enthalpy. The dpq complex 3 showed oxidative double-strand DNA cleavage on exposure to UV-A and visible light. The phen complex 2 displayed single-strand photocleavage of DNA. The DNA double-strand breaks were rationalized from theoretical molecular docking calculations. Mechanistic investigations showed formation of hydroxyl radicals as the reactive species through photodecarboxylation of the L-histidine ligand. The complexes exhibited good binding propensity to bovine serum albumin (BSA) protein in Tris-HCl/NaCl buffer medium. The dpq complex 3 showed UV-A light-induced site-specific oxidative BSA cleavage forming fragments of similar to 45 kDa and similar to 20 kDa molecular weights via SOH pathway.

Stereochemical studies on cyclic peptides. XIII. Energy minimization studies on cyclic hexapeptides having hydrogen bonds

The basic cyclic hexapeptide conformations which accommodate hydrogen bonded β and γ turns in the backbone have been worked out using stereochemical criteria and energy minimization procedures. It was found that cyclic hexapeptides can be made up of all possible combinations of 4 ± 1 hydrogen bonded types I, I', II and II' β turns, giving rise to symmetric conformations having twofold and inversion symmetries as well as nonsymmetric structures. Conformations having exclusive features of 3 ± 1 hydrogen bonded γ turns were found to be possible in threefold and S6 symmetric cyclic hexapeptides. The results show that the cyclic hexapeptides formed by the linking of two β turn tripeptide fragments differ mainly in (a) the hydrogen bonding scheme present in the β turn tripeptides and (b) the conformation at the α-carbon atoms where the two tripeptide fragments link. The different hydrogen bonding schemes found in the component β turns are: 1) a β turn with only a 4 ± 1 hydrogen bond, 2) a type I or I' β turn with 4 ± 1 and 3 ± 1 hydrogen bonds occurring in a bifurcated form and 3) a type II or II' β turn having both the 4 ± 1 and the 3 ± 1 hydrogen bonds with the same acceptor oxygen atom. The conformation at the linking α-carbon atoms was found to lie either in the extended region or in the 3 ± 1 hydrogen bonded γ turn or inverse γ turn regions. Further, the threefold and the S6 symmetric conformations have three γ turns interleaved by three extended regions or three inverse γ turns, respectively. The feasibility of accommodating alanyl residues of both isomeric forms in the CHP minima has been explored. Finally, the available experimental data are reviewed in the light of the present results.

Effect of Crowding Agents, Signal Peptide, and Chaperone SecB on the Folding and Aggregation of E. coli Maltose Binding Protein

SecB, a soluble cytosolic chaperone component of the Secexport pathway, binds to newly synthesized precursor proteins and prevents their premature aggregation and folding and subsequently targets them to the translocation machinery on the membrane. PreMBP, the precursor form of maltose binding protein, has a 26-residue signal sequence attached to the N-terminus of MBP and is a physiological substrate of SecB. We examine the effect of macromolecular crowding and SecB on the stability and refolding of denatured preMBP and MBP. PreMBP was less stable than MBP (ΔTm =7( 0.5 K) in both crowded and uncrowded solutions. Crowding did not cause any substantial changes in the thermal stability ofMBP(ΔTm=1(0.4 K) or preMBP (ΔTm=0(0.6 K), as observed in spectroscopically monitored thermal unfolding experiments. However, both MBP and preMBP were prone to aggregation while refolding under crowded conditions. In contrast to MBP aggregates, which were amorphous, preMBP aggregates form amyloid fibrils.Under uncrowded conditions, a molar excess of SecB was able to completely prevent aggregation and promote disaggregation of preformed aggregates of MBP. When a complex of the denatured protein and SecB was preformed, SecB could completely prevent aggregation and promote folding of MBP and preMBP even in crowded solution. Thus, in addition to maintaining substrates in an unfolded, export-competent conformation, SecB also suppresses the aggregation of its substrates in the crowded intracellular environment. SecB is also able to promote passive disaggregation of macroscopic aggregates of MBP in the absence of an energy source such as ATP or additional cofactors. These experiments also demonstrate that signal peptide can reatly influence protein stability and aggregation propensity.

Mode of inhibition of mitochondrial energy transduction by chlorophenoxyisobutyrate

1. a-p-Chlorophenoxyisobutyric acid, the ethyl ester of which is widely used as an antihypercholesterolaemic drug, is an inhibitor of energy-transfer reactions in isolated rat liver mitochondria. 2. The compound at lower concentrations (<4.0mmol/mg of mitochondrial protein) inhibits state 3 oxidation, stimulates state 4 oxidation, abolishes respiratory control and stimulates the latent adenosine triphosphatase activity of mitochondria. The inhibition imposed on state 3 oxidation is relieved by dinitrophenol. 3. At higher concentrations it inhibits coupled phosphorylation as well as dinitrophenol-stimulated adenosine triphosphatase activity. The inhibition of state 3 oxidation under these conditions is not reversed by uncouplers. 4. The three coupling sites of phosphorylation exhibit differential susceptibility to inactivation by this compound. Coupled phosphorylation at the first site is abolished at a drug concentration of 3.0mmol/mg of protein. The third site is inactivated when the concentration of the drug reaches 5.0mmol/mg of protein. The second site is the most refractory and drug concentrations of the order of 10.0mmol/mg of protein are required effectively to inhibit phosphorylation at this site. 5. The compound also inhibits ATP-dependent reversal of electron transport as well as the adenosine triphosphatase activity in submitochondrial particles. 6. The oxidation of NADH and succinate in these particles is not inhibited. 7. These properties indicate that the compound acts as an `inhibitory uncoupler' of energy-transfer reactions in isolated mitochondria.

Succinate and malate improve development of hamster eight-cell embryos in vitro: Confirmation of viability by embryo transfer

The in vitro development of hamster preimplantation embryos is supported by non-glucose energy substrates. To investigate the importance of embryonic metabolism, influence of succinate and malate on the development of hamster 8-cell embryos to blastocysts was examined using a chemically defined protein-free modified hamster embryo culture medium-2 (HECM-2m). There was a dose-dependent influence of succinate on blastocyst development; 0.5 mM succinate was optimal (85.1% ± 3.9 vs. 54.5% ± 3.5). In succinate-supplemented HECM-2m, blastocyst development was reduced by omission of lactate (68.5% ± 7.2), but not pyruvate (85.8% ± 6.2) or glutamine (84.1% ± 2.1). Succinate along with either glutamine or lactate or pyruvate poorly supported blastocyst development (28%-58%). Malate also stimulated blastocyst development; 0.01 mM malate was optimal (86.3% ± 2.8). Supplementation of both succinate and malate to HECM-2m supported maximal (100%) blastocyst development, which was inhibited 4-fold by the addition of glucose/phosphate. The mean cell numbers (MCN) of blastocysts cultured in succinate-supplemented HECM-2m was higher (28.3 ± 1.1) than it was for those cultured in the absence of glutamine or pyruvate (range 20-24). The MCN was the highest (33.4 ± 1.6) for blastocysts cultured in succinate-malate-supplemented HECM-2m followed by those in succinate (28.3 ± 1.1) or malate (24.7 ± 0.5) supplemented HECM-2m. Embryo transfer experiments showed that 29.8% (±4.5) of transferred blastocysts cultured in succinate-malate-supplemented HECM-2m produced live births, similar (P > 0.1) to the control transfers of freshly recovered 8-cells (33.5% ± 2.0) or blastocysts (28.9% ± 3.0). These data show that supplementation of succinate and malate to HECM-2m supports 100% development of hamster 8-cell embryos to high quality viable blastocysts and that non-glucose oxidizable energy substrates are the most preferred components in hamster embryo culture medium. Mol. Reprod. Dev. 47:440-447, 1997. © 1997 Wiley-Liss, Inc.

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.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

Antithyroid Drugs and their Analogues Protect Against Peroxynitrite-Mediated Protein Tyrosine Nitration-A Mechanistic Study

In this paper, the effect of some commonly used antithyroid drugs and their analogues on peroxynitrite-mediated nitration of proteins is described. The nitration of tyrosine residues in bovine serum albumin (BSA) and cytochromec was studied by Western blot analysis. These studies reveal that the antithyroid drugs methimazole (MMI), 6-n-propyl-2-thiouracil (PTU), and 6-methyl-2-thiouracil (MTU), which contain thione moieties, significantly reduce the tyrosine nitration of both BSA and cytochrome c. While MMI exhibits good peroxynitrite (PN) scavenging activity, the thiouracil compounds PTU and MTU are slightly less effective than MMI. The S- and Se-methylated compounds show a weak inhibitory effect in the nitration of tyrosine, indicating that the presence of a thione or selone moiety is important for an efficient inhibition. Similarly, the replacement of N-H moiety in MMI by N-methyl or N-m-methoxybenzyl substituents dramatically reduces the antioxidant activity of the parent compound. Theoretical studies indicate that the substitution of N-H moiety by N-Me significantly increases the energy required for the oxidation of sulfur center by PN. However, such substitution in the selenium analogue of MMI increases the activity of parent compound. This is due to the facile oxidation of the selone moiety to the corresponding selenenic and seleninic acids. Unlike N,N'-disubstituted thiones, the corresponding selones efficiently scavenge PN, as they predominantly exist in their zwitterionic forms in which the selenium atom carries a large negative charge.

Studies on the conformation of amino acids XII. Energy calculations on prolyl residue

The favoured conformations of the prolyl residue have been obtained by calculating their potential energies arising from bond-angle strain, torsion-angle strain, non-bonded and electrostatic interatomic energies. In addition to the five membered ring, the peptide unit at the amino end (with ω = 180°) and the C′ atom at the carboxyl end have been taken into account. It is found that there are two local minima in the configurational space of the parameters defining the conformation, as is actually observed-one (denoted by B) with Cγ displaced on the same side as C′, which is lower in energy than the other (denoted by A) with Cγ displaced on the opposite side of C′. The other four atoms Cδ, N, Cα, Cβ are nearly in a plane. The conformations of minimum energy (for both A and B) have bond angles very close to the mean observed values while the torsion angles are well within the range observed in various structures for each type. Taking into account the fact that the influence of neighbouring molecules in a crystal structure may make the conformation of a molecule different from the minimal one, the ranges of the conformational parameters for which the energy is within 0.6 kcal/mole above the minimum value (called the "most probable range") and within 1.2 kcal/mole (called the "probable range") have been determined. The ranges thus obtained, agree well with observation, and most of the observed data lie within the most probable ranges, although differing appreciably from the conformation of minimum energy. The study has been extended, in a limited way, to the conformation of the ring in the amino acid proline. Since the nitrogen is tetrahedral in this (as contrasted with being planar in the prolyl residue), it is found that any one of the five atoms can be out of plane (either way), with the other four lying nearly in a plane. These correspond to low energy conformations (up to 1.2 kcal/mole above the minimum). One such example, in which the Cα atom is out of plane is known for dl-proline · HCl. It is also shown that in these calculations energies due to bond length distortions can be neglected to a good degree of approximation, provided the 'best' values of the bond lengths for the particular compound are used in the theoretical calculations.

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.