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

Protein Hydration and Water Structure: X-ray Analysis of a Closely Packed Protein Crystal with Very Low Solvent Content

Low-humidity monoclinic lysozyme, resulting from a water-mediated transformation, has one of the lowest solvent contents (22% by volume) observed in a protein crystal. Its structure has been solved by the molecular replacement method and refined to an R value of 0.175 for 7684 observed reflections in the 10–1.75 Å resolution shell. 90% of the solvent in the well ordered crystals could be located. Favourable sites of hydration on the protein surface include side chains with multiple hydrogen-bonding centres, and regions between short hydrophilic side chains and the main-chain CO or NH groups of the same or nearby residues. Major secondary structural features are not disrupted by hydration. However, the free CO groups at the C terminii and, to a lesser extent, the NH groups at the N terminii of helices provide favourable sites for water interactions, as do reverse turns and regions which connect β-structure and helices. The hydration shell consists of discontinuous networks of water molecules, the maximum number of molecules in a network being ten. The substrate-binding cleft is heavily hydrated, as is the main loop region which is stabilized by water interactions. The protein molecules are close packed in the crystals with a molecular coordination number of 14. Arginyl residues are extensively involved in intermolecular hydrogen bonds and water bridges. The water molecules in the crystal are organized into discrete clusters. A distinctive feature of the clusters is the frequent occurrence of three-membered rings. The protein molecules undergo substantial rearrangement during the transformation from the native to the low-humidity form. The main-chain conformations in the two forms are nearly the same, but differences exist in the side-chain conformation. The differences are particularly pronounced in relation to Trp 62 and Trp 63. The shift in Trp 62 is especially interesting as it is also known to move during inhibitor binding.

Structure and peroxidase activity of ferric Streptomyces clavuligerus orf10-encoded protein P450CLA: UV-visible, CD, MCD and EPR spectroscopic characterization

The present study reports the spectroscopic characterization by UV-visible absorption spectroscopy, magnetic circular dichroism (MCD) and electron paramagnetic resonance (EPR) of the recombinant orf10-encoded P450-camphor like protein (P450CLA)of Streptomyces clavuligerus expressed in Escherichia coli Rosetta in the native form and associated to external ligands containing the β-lactam, oxazole and alkylamine-derived (alcohol) moieties of the clavulamic acid. Considering the diversity of potential applications for the enzyme, the reactivity with tert-butylhydroperoxide (tert-BuOOH) was also characterized. P450CLA presents a covalently bound heme group and exhibited the UV-visible, CD and MCD spectral features of P450CAM including the fingerprint Soret band at 450 nm generated by the ferrous CO-complex. P450CLA was converted to high valence species by tert-BuOOH and promoted homolytic scission of the O-O bond. The radical profile of the reaction was tert-butyloxyl as primary and methyl and butylperoxyl as secondary radicals. The secondary methyl and butylperoxyl radicals resulted respectively from the β-scission of the alkoxyl radical and from the reaction of methyl radical with molecular oxygen.

A link between protein structure and enzyme catalyzed hydrogen tunneling

We present evidence that the size of an active site side chain may modulate the degree of hydrogen tunneling in an enzyme-catalyzed reaction. Primary and secondary kH/kT and kD/kT kinetic isotope effects have been measured for the oxidation of benzyl alcohol catalyzed by horse liver alcohol dehydrogenase at 25°C. As reported in earlier studies, the relationship between secondary kH/kT and kD/kT isotope effects provides a sensitive probe for deviations from classical behavior. In the present work, catalytic efficiency and the extent of hydrogen tunneling have been correlated for the alcohol dehydrogenase-catalyzed hydride transfer among a group of site-directed mutants at position 203. Val-203 interacts with the opposite face of the cofactor NAD+ from the alcohol substrate. The reduction in size of this residue is correlated with diminished tunneling and a two orders of magnitude decrease in catalytic efficiency. Comparison of the x-ray crystal structures of a ternary complex of a high-tunneling (Phe-93 → Trp) and a low-tunneling (Val-203 → Ala) mutant provides a structural basis for the observed effects, demonstrating an increase in the hydrogen transfer distance for the low-tunneling mutant. The Val-203 → Ala ternary complex crystal structure also shows a hyperclosed interdomain geometry relative to the wild-type and the Phe-93 → Trp mutant ternary complex structures. This demonstrates a flexibility in interdomain movement that could potentially narrow the distance between the donor and acceptor carbons in the native enzyme and may enhance the role of tunneling in the hydride transfer reaction.

Folding of a nascent polypeptide chain in vitro: cooperative formation of structure in a protein module.

We have prepared a family of peptide fragments of the 64-residue chymotrypsin inhibitor 2, corresponding to its progressive elongation from the N terminus. The growing polypeptide chain has little tendency to form stable structure until it is largely synthesized, and what structures are formed are nonnative and lack, in particular, the native secondary structural elements of alpha-helix and beta-sheet. These elements then develop as sufficient tertiary interactions are made in the nearly full-length chain. The growth of structure in the small module is highly cooperative and does not result from the hierarchical accretion of substructures.

The accuracy of protein structure alignment servers

Background: Protein structural alignment is one of the most fundamental and crucial areas of research in the domain of computational structural biology. Comparison of a protein structure with known structures helps to classify it as a new or belonging to a known group of proteins. This, in turn, is useful to determine the function of protein, its evolutionary relationship with other protein molecules and grasping principles underlying protein architecture and folding. Results: A large number of protein structure alignment methods are available. Each protein structure alignment tool has its own strengths andweaknesses that need to be highlighted.We compared and presented results of six most popular and publically available servers for protein structure comparison. These web-based servers were compared with the respect to functionality (features provided by these servers) and accuracy (how well the structural comparison is performed). The CATH was used as a reference. The results showed that overall CE was top performer. DALI and PhyreStorm showed similar results whereas PDBeFold showed the lowest performance. In case of few secondary structural elements, CE, DALI and PhyreStorm gave 100% success rate. Conclusion: Overall none of the structural alignment servers showed 100% success rate. Studies of overall performance, effect of mainly alpha and effect of mainly beta showed consistent performance. CE, DALI, FatCat and PhyreStorm showed more than 90% success rate.

In situ preparation and protein delivery of silicate/alginate composite microspheres with core-shell structure

It is predicted that with increased life expectancy in the developed world, there will be a greater demand for synthetic materials to repair or regenerate lost, injured or diseased bone (Hench & Thompson 2010). There are still few synthetic materials having true bone inductivity, which limits their application for bone regeneration, especially in large-size bone defects. To solve this problem, growth factors, such as bone morphogenetic proteins (BMPs), have been incorporated into synthetic materials in order to stimulate de novo bone formation in the center of large-size bone defects. The greatest obstacle with this approach is that the rapid diffusion of the protein from the carrier material, leading to a precipitous loss of bioactivity; the result is often insufficient local induction or failure of bone regeneration (Wei et al. 2007). It is critical that the protein is loaded in the carrier material in conditions which maintains its bioactivity (van de Manakker et al. 2009). For this reason, the efficient loading and controlled release of a protein from a synthetic material has remained a significant challenge. The use of microspheres as protein/drug carriers has received considerable attention in recent years (Lee et al. 2010; Pareta & Edirisinghe 2006; Wu & Zreiqat 2010). Compared to macroporous block scaffolds, the chief advantage of microspheres is their superior protein-delivery properties and ability to fill bone defects with irregular and complex shapes and sizes. Upon implantation, the microspheres are easily conformed to the irregular implant site, and the interstices between the particles provide space for both tissue and vascular ingrowth, which are important for effective and functional bone regeneration (Hsu et al. 1999). Alginates are natural polysaccharides and their production does not have the implicit risk of contamination with allo or xeno-proteins or viruses (Xie et al. 2010). Because alginate is generally cytocompatible, it has been used extensively in medicine, including cell therapy and tissue engineering applications (Tampieri et al. 2005; Xie et al. 2010; Xu et al. 2007). Calcium cross-linked alginate hydrogel is considered a promising material as a delivery matrix for drugs and proteins, since its gel microspheres form readily in aqueous solutions at room temperature, eliminating the need for harsh organic solvents, thereby maintaining the bioactivity of proteins in the process of loading into the microspheres (Jay & Saltzman 2009; Kikuchi et al. 1999). In addition, calcium cross-linked alginate hydrogel is degradable under physiological conditions (Kibat PG et al. 1990; Park K et al. 1993), which makes alginate stand out as an attractive candidate material for the protein carrier and bone regeneration (Hosoya et al. 2004; Matsuno et al. 2008; Turco et al. 2009). However, the major disadvantages of alginate microspheres is their low loading efficiency and also rapid release of proteins due to the mesh-like networks of the gel (Halder et al. 2005). Previous studies have shown that a core-shell structure in drug/protein carriers can overcome the issues of limited loading efficiencies and rapid release of drug or protein (Chang et al. 2010; Molvinger et al. 2004; Soppimath et al. 2007). We therefore hypothesized that introducing a core-shell structure into the alginate microspheres could solve the shortcomings of the pure alginate. Calcium silicate (CS) has been tested as a biodegradable biomaterial for bone tissue regeneration. CS is capable of inducing bone-like apatite formation in simulated body fluid (SBF) and its apatite-formation rate in SBF is faster than that of Bioglass® and A-W glass-ceramics (De Aza et al. 2000; Siriphannon et al. 2002). Titanium alloys plasma-spray coated with CS have excellent in vivo bioactivity (Xue et al. 2005) and porous CS scaffolds have enhanced in vivo bone formation ability compared to porous β-tricalcium phosphate ceramics (Xu et al. 2008). In light of the many advantages of this material, we decided to prepare CS/alginate composite microspheres by combining a CS shell with an alginate core to improve their protein delivery and mineralization for potential protein delivery and bone repair applications

Assessment of the molecular structure of natrodufrénite – NaFe2+Fe3+5(PO4)4(OH)6•2(H2O), a secondary pegmatite phosphate mineral from Minas Gerais, Brazil

The mineral natrodufrénite a secondary pegmatite phosphate mineral from Minas Gerais, Brazil, has been studied by a combination of scanning electron microscopy and vibrational spectroscopic techniques. Electron probe analysis shows the formula of the studied mineral as (Na0.88Ca0.12)∑1.00(Mn0.11Mg0.08Ca0.04Zr0.01Cu0.01)∑0.97(Al0.02)∑4.91(PO4)3.96(OH6.15F0.07)6.22⋅2.05(H2O). Raman spectroscopy identifies an intense peak at 1003 cm−1 assigned to the ν1 symmetric stretching mode. Raman bands are observed at 1059 and 1118 cm−1 and are attributed to the ν3 antisymmetric stretching vibrations. A comparison is made with the spectral data of other hydrate hydroxy phosphate minerals including cyrilovite and wardite. Raman bands at 560, 582, 619 and 668 cm−1 are assigned to the ν4 bending modes and Raman bands at 425, 444, 477 and 507 cm−1 are due to the ν2 bending modes. Raman bands in the 2600–3800 cm−1 spectral range are attributed to water and OH stretching vibrations. Vibrational spectroscopy enables aspects of the molecular structure of natrodufrénite to be assessed.

Phenytoin metabolism by human cytochrome P450 : involvement of P450 3A and 2C forms in secondary metabolism and drug-protein adduct formation

The anticonvulsant phenytoin (5,5-diphenylhydantoin) provokes a skin rash in 5 to 10% of patients, which heralds the start of an idiosyncratic reaction that may result from covalent modification of normal self proteins by reactive drug metabolites. Phenytoin is metabolized by cytochrome P450 (P450) enzymes primarily to 5-(p-hydroxyphenyl-),5-phenylhydantoin (HPPH), which may be further metabolized to a catechol that spontaneously oxidizes to semiquinone and quinone species that covalently modify proteins. The aim of this study was to determine which P450s catalyze HPPH metabolism to the catechol, proposed to be the final enzymatic step in phenytoin bioactivation. Recombinant human P450s were coexpressed with NADPH-cytochrome P450 reductase in Escherichia coli. Novel bicistronic expression vectors were constructed for P450 2C19 and the three major variants of P450 2C9, i.e., 2C9*1, 2C9*2, and 2C9*3. HPPH metabolism and covalent adduct formation were assessed in parallel. P450 2C19 was the most effective catalyst of HPPH oxidation to the catechol metabolite and was also associated with the highest levels of covalent adduct formation. P450 3A4, 3A5, 3A7, 2C9*1, and 2C9*2 also catalyzed bioactivation of HPPH, but to a lesser extent. Fluorographic analysis showed that the major targets of adduct formation in bacterial membranes were the catalytic P450 forms, as suggested from experiments with human liver microsomes. These results suggest that P450 2C19 and other forms from the 2C and 3A subfamilies may be targets as well as catalysts of drug-protein adduct formation from phenytoin.

Spectroscopic vibrations of austinite (CaZnAsO4⋅OH) and its mineral structure : Implications for identification of secondary arsenic-containing mineral

Austinite (CaZnAsO4⋅OH) is a unique secondary mineral in arsenic-contaminated mine wastes. The infrared and Raman spectroscopies were used to characterize the austenite vibrations. The IR bands at 369, 790 and 416 cm−1 are assigned to the ν2, ν3 and ν4 vibrations of AsO43− unit, respectively. The Raman bands at 814, 779 and 403 cm−1 correspond to the ν1, ν3 and ν4 vibrations of AsO43− unit respectively. The sharp bands at 3265 cm−1 for IR and 3270 cm−1 both reveals that the structural hydroxyl units exist in the austenite structure. The IR and Raman spectra both show that some SO4 units isomorphically replace AsO4 in austinite. X-ray single crystal diffraction provides the arrangement of each atom in the mineral structure, and also confirms that the conclusions made from the vibrational spectra. Micro-powder diffraction was used to confirm our mineral identification due to the small quantity of the austenite crystals.

Similarities in the Genomic Sequence and Coat Protein-Structure of Plant Viruses

The genomic sequences of several RNA plant viruses including cucumber mosaic virus, brome mosaic virus, alfalfa mosaic virus and tobacco mosaic virus have become available recently. The former two viruses are icosahedral while the latter two are bullet and rod shaped, respectively in particle morphology. The non-structural 3a proteins of cucumber mosaic virus and brome mosaic virus have an amino acid sequence homology of 35% and hence are evolutionarily related. In contrast, the coat proteins exhibit little homology, although the circular dichroism spectrum of these viruses are similar. The non-coding regions of the genome also exhibit variable but extensive homology. Comparison of the brome mosaic virus and alfalfa mosaic virus sequences reveals that they are probably related although with a much larger evolutionary distance. The polypeptide folds of the coat protein of three biologically distinct isometric plant viruses, tomato bushy stunt virus, southern bean mosaic virus and satellite tobacco necrosis virus have been shown to display a striking resemblance. All of them consist of a topologically similar 8-standard β-barrel. The implications of these studies to the understanding of the evolution of plant viruses will be discussed.

Intra and Inter-Molecular Communications Through Protein Structure Network

Communication within and across proteins is crucial for the biological functioning of proteins. Experiments such as mutational studies on proteins provide important information on the amino acids, which are crucial for their function. However, the protein structures are complex and it is unlikely that the entire responsibility of the function rests on only a few amino acids. A large fraction of the protein is expected to participate in its function at some level or other. Thus, it is relevant to consider the protein structures as a completely connected network and then deduce the properties, which are related to the global network features. In this direction, our laboratory has been engaged in representing the protein structure as a network of non-covalent connections and we have investigated a variety of problems in structural biology, such as the identification of functional and folding clusters, determinants of quaternary association and characterization of the network properties of protein structures. We have also addressed a few important issues related to protein dynamics, such as the process of oligomerization in multimers, mechanism on protein folding, and ligand induced communications (allosteric effect). In this review we highlight some of the investigations which we have carried out in the recent past. A review on protein structure graphs was presented earlier, in which the focus was on the graphs and graph spectral properties and their implementation in the study of protein structure graphs/networks (PSN). In this article, we briefly summarize the relevant parts of the methodology and the focus is on the advancement brought out in the understanding of protein structure-function relationships through structure networks. The investigations of structural/biological problems are divided into two parts, in which the first part deals with the analysis of PSNs based on static structures obtained from x-ray crystallography. The second part highlights the changes in the network, associated with biological functions, which are deduced from the network analysis on the structures obtained from molecular dynamics simulations.

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.