NMR of hydrogen bonding in cold-shock protein A and an analysis of the influence of crystallographic resolution on comparisons of hydrogen bond lengths

Hydrogen bonding in cold-shock protein A of Escherichia coli has been investigated using long-range HNCO spectroscopy. Nearly half of the amide protons involved in hydrogen bonds in solution show no measurable protection from exchange in water, cautioning against a direct correspondence between hydrogen bonding and hydrogen exchange protection. The N to O atom distance across a hydrogen bond, R(NO), is related to the size of the (3h)J(NC') trans hydrogen bond coupling constant and the amide proton chemical shift. Both NMR parameters show poorer agreement with the 2.0-A resolution X-ray structure of the cold-shock protein studied by NMR than with a 1.2-A resolution X-ray structure of a homologous cold-shock protein from the thermophile B. caldolyticus. The influence of crystallographic resolution on comparisons of hydrogen bond lengths was further investigated using a database of 33 X-ray structures of ribonuclease A. For highly similar structures, both hydrogen bond R(NO) distance and Calpha coordinate root mean square deviations (RMSD) show systematic increases as the resolution of the X-ray structure used for comparison decreases. As structures diverge, the effects of coordinate errors on R(NO) distance and Calpha coordinate root mean square deviations become progressively smaller. The results of this study are discussed with regard to the influence of data precision on establishing structure similarity relationships between proteins.     

Homology modeling and molecular dynamics simulation of odonthubuthus doriae (Od1) scorpion toxin in comparison to the BmK M1

All of the α-subgroups share similarity in their sequence and structure but different in the toxicity to various voltage-gated sodium channels (VGSCs). We modeled the first 3D structural model of the Od1 based on BmK M1 using homology modeling. The reliability of model for more investigation and compare to BmK M1 has been examined and confirmed. Then the model structure is further refined by energy minimization and molecular dynamics methods. The purpose of this modeling and simulation is comparison toxicity of two mentioned toxins by investigation structural feature of functional regions including core domain, 5-turn and C-terminal which make NC domain. In the one hand, it is intriguing that Od1 in comparison to BmK M1 shows same solvent accessible surface area (SASA) in 5-turn region but a little more exposed and feasibility (more SASA) in C-terminal region and key functional residues of C-terminal such as positive residues Arg58, lys62 and Arg (His)64. These data suggested that Od1 has similarity with BmK M1 but has more toxicity to sodium channel. In the other hand 5-turn proximity of C-terminal to 5-turn in BmK M1with cis peptide bond is less than Od1 without cis peptide bond which is a confirmation with experimental data about BmK M1.A better understanding of the 3-D structure of Od1and comparison to BmK M1 will be helpful for more investigation of functional characters action of natural toxins with a specialized role for VGSCs.     

Antifreeze protein from shorthorn sculpin: Identification of the ice-binding surface

Shorthorn sculpins, Myoxocephalus scorpius, are protected from freezing in icy seawater by alanine-rich, alpha-helical antifreeze proteins (AFPs). The major serum isoform (SS-8) has been reisolated and analyzed to establish its correct sequence. Over most of its length, this 42 amino acid protein is predicted to be an amphipathic alpha-helix with one face entirely composed of Ala residues. The other side of the helix, which is more heterogeneous and hydrophilic, contains several Lys. Computer simulations had suggested previously that these Lys residues were involved in binding of the peptide to the [11-20] plane of ice in the <-1102> direction. To test this hypothesis, a series of SS-8 variants were generated with single Ala to Lys substitutions at various points around the helix. All of the peptides retained significant alpha-helicity and remained as monomers in solution. Substitutions on the hydrophilic helix face at position 16, 19, or 22 had no obvious effect, but those on the adjacent Ala-rich surface at positions 17, 21, and 25 abolished antifreeze activity. These results, with support from our own modeling and docking studies, show that the helix interacts with the ice surface via the conserved alanine face, and lend support to the emerging idea that the interaction of fish AFPs with ice involves appreciable hydrophobic interactions. Furthermore, our modeling suggests a new N terminus cap structure, which helps to stabilize the helix, whereas the role of the lysines on the hydrophilic face may be to enhance solubility of the protein.     

ExSer: A standalone tool to mine protein data bank (PDB) for secondary structural elements

Detailed structural analysis of protein necessitates investigation at primary, secondary and tertiary levels, respectively. Insight into protein secondary structures pave way for understanding the type of secondary structural elements involved (α-helices, β-strands etc.), the amino acid sequence that encode the secondary structural elements, number of residues, length and, percentage composition of the respective elements in the protein. Here we present a standalone tool entitled "ExSer" which facilitate an automated extraction of the amino acid sequence that encode for the secondary structural regions of a protein from the protein data bank (PDB) file. AVAILABILITY: ExSer is freely downloadable from http://code.google.com/p/tool-exser/     

ICP34.5 protein of herpes simplex virus facilitates the initiation of protein translation by bridging eukaryotic initiation factor 2alpha (eIF2alpha) and protein phosphatase 1

The ICP34.5 protein of herpes simplex virus type 1 is a neurovirulence factor that plays critical roles in viral replication and anti-host responses. One of its functions is to recruit protein phosphatase 1 (PP1) that leads to the dephosphorylation of the α subunit of translation initiation factor eIF2 (eIF2α), which is inactivated by infection-induced phosphorylation. As PP1 is a protein phosphatase with a wide range of substrates, the question remains to be answered how ICP34.5 directs PP1 to specifically dephosphorylate eIF2α. Here we report that ICP34.5 not only binds PP1 but also associates with eIF2α by in vitro and in vivo assays. The binding site of eIF2α is identified at amino acids 233-248 of ICP34.5, which falls in the highly homologous region with human gene growth arrest and DNA damage 34. The interaction between ICP34.5 and eIF2α is independent of the phosphorylation status of eIF2α at serine 51. Deletion mutation of this region results in the failure of dephosphorylation of eIF2α by PP1 and, consequently, interrupts viral protein synthesis and replication. Our data illustrated that the binding between viral protein ICP34.5 and the host eIF2α is crucial for the specific dephosphorylation of eIF2α by PP1. We propose that herpes simplex virus protein ICP34.5 bridges PP1 and eIF2α via their binding motifs and thereby facilitates the protein synthesis and viral replication.     

Protein kinase Ctheta is a specific target for inhibition of the HIV type 1 replication in CD4+ T lymphocytes

Integration of HIV-1 genome in CD4(+) T cells produces latent reservoirs with long half-life that impedes the eradication of the infection. Control of viral replication is essential to reduce the size of latent reservoirs, mainly during primary infection when HIV-1 infects CD4(+) T cells massively. The addition of immunosuppressive agents to highly active antiretroviral therapy during primary infection would suppress HIV-1 replication by limiting T cell activation, but these agents show potential risk for causing lymphoproliferative disorders. Selective inhibition of PKC, crucial for T cell function, would limit T cell activation and HIV-1 replication without causing general immunosuppression due to PKC being mostly expressed in T cells. Accordingly, the effect of rottlerin, a dose-dependent PKC inhibitor, on HIV-1 replication was analyzed in T cells. Rottlerin was able to reduce HIV-1 replication more than 20-fold in MT-2 (IC(50) = 5.2 μM) and Jurkat (IC(50) = 2.2 μM) cells and more than 4-fold in peripheral blood lymphocytes (IC(50) = 4.4 μM). Selective inhibition of PKC, but not PKCδ or -ζ, was observed at <6.0 μM, decreasing the phosphorylation at residue Thr(538) on the kinase catalytic domain activation loop and avoiding PKC translocation to the lipid rafts. Consequently, the main effector at the end of PKC pathway, NF-κB, was repressed. Rottlerin also caused a significant inhibition of HIV-1 integration. Recently, several specific PKC inhibitors have been designed for the treatment of autoimmune diseases. Using these inhibitors in combination with highly active antiretroviral therapy during primary infection could be helpful to avoid massive viral infection and replication from infected CD4(+) T cells, reducing the reservoir size at early stages of the infection.     

Derivation of rules for comparative protein modeling from a database of protein structure alignments

We describe a database of protein structure alignments as well as methods and tools that use this database to improve comparative protein modeling. The current version of the database contains 105 alignments of similar proteins or protein segments. The database comprises 416 entries, 78,495 residues, 1,233 equivalent entry pairs, and 230,396 pairs of equivalent alignment positions. At present, the main application of the database is to improve comparative modeling by satisfaction of spatial restraints implemented in the program MODELLER (?ali A, Blundell TL, 1993, J Mol Biol 234:779–815). To illustrate the usefulness of the database, the restraints on the conformation of a disulfide bridge provided by an equivalent disulfide bridge in a related structure are derived from the alignments; the prediction success of the disulfide dihedral angle classes is increased to approximately 80%, compared to approximately 55% for modeling that relies on the stereochemistry of disulfide bridges alone. The second example of the use of the database is the derivation of the probability density function for comparative modeling of the cis/trans isomerism of the proline residues; the prediction success is increased from 0% to 82.9% for cis-proline and from 93.3% to 96.2% for trans-proline. The database is available via electronic mail.     

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on Cα atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 Å, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low. © 1995 Wiley-Liss, Inc.     

Surfactant protein D (Sp-D) binds to membrane-proximal domain (D3) of signal regulatory protein α (SIRPα), a site distant from binding domain of CD47, while also binding to analogous region on signal regulatory protein β (SIRPβ)

Signal regulatory protein α (SIRPα), a highly glycosylated type-1 transmembrane protein, is composed of three immunoglobulin-like extracellular loops as well as a cytoplasmic tail containing three classical tyrosine-based inhibitory motifs. Previous reports indicate that SIRPα binds to humoral pattern recognition molecules in the collectin family, namely surfactant proteins D and A (Sp-D and Sp-A, respectively), which are heavily expressed in the lung and constitute one of the first lines of innate immune defense against pathogens. However, little is known about molecular details of the structural interaction of Sp-D with SIRPs. In the present work, we examined the molecular basis of Sp-D binding to SIRPα using domain-deleted mutant proteins. We report that Sp-D binds to the membrane-proximal Ig domain (D3) of SIRPα in a calcium- and carbohydrate-dependent manner. Mutation of predicted N-glycosylation sites on SIRPα indicates that Sp-D binding is dependent on interactions with specific N-glycosylated residues on the membrane-proximal D3 domain of SIRPα. Given the remarkable sequence similarity of SIRPα to SIRPβ and the lack of known ligands for the latter, we examined Sp-D binding to SIRPβ. Here, we report specific binding of Sp-D to the membrane-proximal D3 domain of SIRPβ. Further studies confirmed that Sp-D binds to SIRPα expressed on human neutrophils and differentiated neutrophil-like cells. Because the other known ligand of SIRPα, CD47, binds to the membrane-distal domain D1, these findings indicate that multiple, distinct, functional ligand binding sites are present on SIRPα that may afford differential regulation of receptor function.     

Factors contributing to decreased protein stability when aspartic acid residues are in beta-sheet regions

Asp residues are significantly under represented in beta-sheet regions of proteins, especially in the middle of beta-strands, as found by a number of studies using statistical, modeling, or experimental methods. To further understand the reasons for this under representation of Asp, we prepared and analyzed mutants of a beta-domain. Two Gln residues of the immunoglobulin light-chain variable domain (V(L)) of protein Len were replaced with Asp, and then the effects of these changes on protein stability and protein structure were studied. The replacement of Q38D, located at the end of a beta-strand, and that of Q89D, located in the middle of a beta-strand, reduced the stability of the parent immunoglobulin V(L) domain by 2.0 kcal/mol and 5.3 kcal/mol, respectively. Because the Q89D mutant of the wild-type V(L)-Len domain was too unstable to be expressed as a soluble protein, we prepared the Q89D mutant in a triple mutant background, V(L)-Len M4L/Y27dD/T94H, which was 4.2 kcal/mol more stable than the wild-type V(L)-Len domain. The structures of mutants V(L)-Len Q38D and V(L)-Len Q89D/M4L/Y27dD/T94H were determined by X-ray diffraction at 1.6 A resolution. We found no major perturbances in the structures of these Q-->D mutant proteins relative to structures of the parent proteins. The observed stability changes have to be accounted for by cumulative effects of the following several factors: (1) by changes in main-chain dihedral angles and in side-chain rotomers, (2) by close contacts between some atoms, and, most significantly, (3) by the unfavorable electrostatic interactions between the Asp side chain and the carbonyls of the main chain. We show that the Asn side chain, which is of similar size but neutral, is less destabilizing. The detrimental effect of Asp within a beta-sheet of an immunoglobulin-type domain can have very serious consequences. A somatic mutation of a beta-strand residue to Asp could prevent the expression of the domain both in vitro and in vivo, or it could contribute to the pathogenic potential of the protein in vivo.     

Why proteins without an α-crystallin domain should not be included in the human small heat shock protein family HSPB

The presence of an α-crystallin domain documents the evolutionary relatedness of the ubiquitous family of small heat shock proteins. Sequence and three-dimensional structure provide no evidence for the presence of such a domain in HSPC034, recently proposed as the 11th member of the human HSPB family. Also, phylogenetic analyses detect no relationship between HSPC034 and the human HSPB1–10 sequences. Arguments are provided as to why inclusion in the HSPB family of proteins like HSPC034, which resemble small heat shock proteins in being heat-inducible and having chaperone-like properties and a low monomeric mass, but are evolutionarily unrelated, is misleading and confusing.     

Toward the solution of the protein structure prediction problem

Since Anfinsen demonstrated that the information encoded in a protein’s amino acid sequence determines its structure in 1973, solving the protein structure prediction problem has been the Holy Grail of structural biology. The goal of protein structure prediction approaches is to utilize computational modeling to determine the spatial location of every atom in a protein molecule starting from only its amino acid sequence. Depending on whether homologous structures can be found in the Protein Data Bank (PDB), structure prediction methods have been historically categorized as template-based modeling (TBM) or template-free modeling (FM) approaches. Until recently, TBM has been the most reliable approach to predicting protein structures, and in the absence of reliable templates, the modeling accuracy sharply declines. Nevertheless, the results of the most recent community-wide assessment of protein structure prediction experiment (CASP14) have demonstrated that the protein structure prediction problem can be largely solved through the use of end-to-end deep machine learning techniques, where correct folds could be built for nearly all single-domain proteins without using the PDB templates. Critically, the model quality exhibited little correlation with the quality of available template structures, as well as the number of sequence homologs detected for a given target protein. Thus, the implementation of deep-learning techniques has essentially broken through the 50-year-old modeling border between TBM and FM approaches and has made the success of high-resolution structure prediction significantly less dependent on template availability in the PDB library.     

Nature of the amyloid-beta monomer and the monomer-oligomer equilibrium

The monomer to oligomer transition initiates the aggregation and pathogenic transformation of Alzheimer amyloid-β (Aβ) peptide. However, the monomeric state of this aggregation-prone peptide has remained beyond the reach of most experimental techniques, and a quantitative understanding of this transition is yet to emerge. Here, we employ single-molecule level fluorescence tools to characterize the monomeric state and the monomer-oligomer transition at physiological concentrations in buffers mimicking the cerebrospinal fluid (CSF). Our measurements show that the monomer has a hydrodynamic radius of 0.9 ± 0.1 nm, which confirms the prediction made by some of the in silico studies. Surprisingly, at equilibrium, both Aβ(40) and Aβ(42) remain predominantly monomeric up to 3 μm, above which it forms large aggregates. This concentration is much higher than the estimated concentrations in the CSF of either normal or diseased brains. If Aβ oligomers are present in the CSF and are the key agents in Alzheimer pathology, as is generally believed, then these must be released in the CSF as preformed entities. Although the oligomers are thermodynamically unstable, we find that a large kinetic barrier, which is mostly entropic in origin, strongly impedes their dissociation. Thermodynamic principles therefore allow the development of a pharmacological agent that can catalytically convert metastable oligomers into nontoxic monomers.     

Computational protein structure modeling and analysis of UV-B stress protein in Synechocystis PCC 6803

This study focuses on Ultra Violet stress (UVS) gene product which is a UV stress induced protein from cyanobacteria, Synechocystis PCC 6803. Three dimensional structural modeling of target UVS protein was carried out by homology modeling method. 3F2I pdb from Nostoc sp. PCC 7120 was selected as a suitable template protein structure. Ultimately, the detection of active binding regions was carried out for characterization of functional sites in modeled UV-B stress protein. The top five probable ligand binding sites were predicted and the common binding residues between target and template protein was analyzed. It has been validated for the first time that modeled UVS protein structure from Synechocystis PCC 6803 was structurally and functionally similar to well characterized UVS protein of another cyanobacterial species, Nostoc sp PCC 7120 because of having same structural motif and fold with similar protein topology and function. Investigations revealed that UVS protein from Synechocystis sp. might play significant role during ultraviolet resistance. Thus, it could be a potential biological source for remediation for UV induced stress.     

New tricks of an old pattern: structural versatility of scorpion toxins with common cysteine spacing

Scorpion venoms are a rich source of K(+) channel-blocking peptides. For the most part, they are structurally related small disulfide-rich proteins containing a conserved pattern of six cysteines that is assumed to dictate their common three-dimensional folding. In the conventional pattern, two disulfide bridges connect an α-helical segment to the C-terminal strand of a double- or triple-stranded β-sheet, conforming a cystine-stabilized α/β scaffold (CSα/β). Here we show that two K(+) channel-blocking peptides from Tityus scorpions conserve the cysteine spacing of common scorpion venom peptides but display an unconventional disulfide pattern, accompanied by a complete rearrangement of the secondary structure topology into a CS helix-loop-helix fold. Sequence and structural comparisons of the peptides adopting this novel fold suggest that it would be a new elaboration of the widespread CSα/β scaffold, thus revealing an unexpected structural versatility of these small disulfide-rich proteins. Acknowledgment of such versatility is important to understand how venom structural complexity emerged on a limited number of molecular scaffolds.     

Refinement of protein structures by iterative comparative modeling and CryoEM density fitting

We developed a method for structure characterization of assembly components by iterative comparative protein structure modeling and fitting into cryo-electron microscopy (cryoEM) density maps. Specifically, we calculate a comparative model of a given component by considering many alternative alignments between the target sequence and a related template structure while optimizing the fit of a model into the corresponding density map. The method relies on the previously developed Moulder protocol that iterates over alignment, model building, and model assessment. The protocol was benchmarked using 20 varied target-template pairs of known structures with less than 30% sequence identity and corresponding simulated density maps at resolutions from 5A to 25A. Relative to the models based on the best existing sequence profile alignment methods, the percentage of C(alpha) atoms that are within 5A of the corresponding C(alpha) atoms in the superposed native structure increases on average from 52% to 66%, which is half-way between the starting models and the models from the best possible alignments (82%). The test also reveals that despite the improvements in the accuracy of the fitness function, this function is still the bottleneck in reducing the remaining errors. To demonstrate the usefulness of the protocol, we applied it to the upper domain of the P8 capsid protein of rice dwarf virus that has been studied by cryoEM at 6.8A. The C(alpha) root-mean-square deviation of the model based on the remotely related template, bluetongue virus VP7, improved from 8.7A to 6.0A, while the best possible model has a C(alpha) RMSD value of 5.3A. Moreover, the resulting model fits better into the cryoEM density map than the initial template structure. The method is being implemented in our program MODELLER for protein structure modeling by satisfaction of spatial restraints and will be applicable to the rapidly increasing number of cryoEM density maps of macromolecular assemblies.     

Multiple mapping method: a novel approach to the sequence-to-structure alignment problem in comparative protein structure modeling

A major bottleneck in comparative protein structure modeling is the quality of input alignment between the target sequence and the template structure. A number of alignment methods are available, but none of these techniques produce consistently good solutions for all cases. Alignments produced by alternative methods may be superior in certain segments but inferior in others when compared to each other; therefore, an accurate solution often requires an optimal combination of them. To address this problem, we have developed a new approach, Multiple Mapping Method (MMM). The algorithm first identifies the alternatively aligned regions from a set of input alignments. These alternatively aligned segments are scored using a composite scoring function, which determines their fitness within the structural environment of the template. The best scoring regions from a set of alternative segments are combined with the core part of the alignments to produce the final MMM alignment. The algorithm was tested on a dataset of 1400 protein pairs using 11 combinations of two to four alignment methods. In all cases MMM showed statistically significant improvement by reducing alignment errors in the range of 3 to 17%. MMM also compared favorably over two alignment meta-servers. The algorithm is computationally efficient; therefore, it is a suitable tool for genome scale modeling studies.     

Molecular modeling of protein function regions

Experimental protein structures often provide extensive insight into the mode and specificity of small molecule binding, and this information is useful for understanding protein function and for the design of drugs. We have performed an analysis of the reliability with which ligand-binding information can be deduced from computer model structures, as opposed to experimentally derived ones. Models produced as part of the CASP experiments are used. The accuracy of contacts between protein model atoms and experimentally determined ligand atom positions is the main criterion. Only comparative models are included (i.e., models based on a sequence relationship between the protein of interest and a known structure). We find that, as expected, contact errors increase with decreasing sequence identity used as a basis for modeling. Analysis of the causes of errors shows that sequence alignment errors between model and experimental template have the most deleterious effect. In general, good, but not perfect, insight into ligand binding can be obtained from models based on a sequence relationship, providing there are no alignment errors in the model. The results support a structural genomics strategy based on experimental sampling of structure space so that all protein domains can be modeled on the basis of 30% or higher sequence identity.     

N-terminal half of the cGMP phosphodiesterase gamma-subunit contributes to stabilization of the GTPase-accelerating protein complex

In the visual signal terminating transition state, the cyclic GMP phosphodiesterase (PDE6) inhibitory γ-subunit (PDEγ) stimulates GTPase activity of the α-subunit of transducin (αt) by enhancing the interaction between αt and its regulator of G protein signaling (RGS9), which is constitutively bound to the type 5 G protein β-subunit (β5). Although it is known from a crystal structure of partial molecules that the PDEγ C terminus contacts with both αt and RGS9, contributions from the intrinsically disordered PDEγ N-terminal half remain unclear. In this study, we were able to investigate this issue using a photolabel transfer strategy that allows for mapping the interface of full-length proteins. We observed label transfer from PDEγ N-terminal positions 50, 30, and 16 to RGS9·β5 in the GTPase-accelerating protein (GAP) complex composed of PDEγ·αt·RGS9·β5. In support of a direct PDEγ N-terminal interaction with RGS9·β5, the PDEγ N-terminal peptide PDEγ(1-61) abolished label transfer to RGS9·β5, and another N-terminal peptide, PDEγ(10-30), disassembled the GAP complex in label transfer and pulldown experiments. Furthermore, we determined that the PDEγ C-terminal interaction with αt was enhanced whereas the N-terminal interaction was weakened upon changing the αt conformation from the signaling state to the transition state. This "rearrangement" of PDEγ domain interactions with αt appears to facilitate the interaction of the PDEγ N-terminal half with RGS9·β5 and hence its contribution to optimal stabilization of the GAP complex.