Cooperative effects in hydrogen-bonding of protein secondary structure elements: a systematic analysis of crystal data using Secbase

A systematic analysis of the hydrogen-bonding geometry in helices and beta sheets has been performed. The distances and angles between the backbone carbonyl O and amide N atoms were correlated considering more than 1500 protein chains in crystal structures determined to a resolution better than 1.5 A. They reveal statistically significant trends in the H-bond geometry across the different secondary structural elements. The analysis has been performed using Secbase, a modular extension of Relibase (Receptor Ligand Database) which integrates information about secondary structural elements assigned to individual protein structures with the various search facilities implemented into Relibase. A comparison of the mean hydrogen-bond distances in alpha helices and 3(10) helices of increasing length shows opposing trends. Whereas in alpha helices the mean H-bond distance shrinks with increasing helix length and turn number, the corresponding mean dimension in 3(10) helices expands in a comparable series. Comparing similarly the hydrogen-bond lengths in beta sheets there is no difference to be found between the mean H-bond length in antiparallel and parallel beta sheets along the strand direction. In contrast, an interesting systematic trend appears to be given for the hydrogen bonds perpendicular to the strands bridging across an extended sheet. With increasing number of accumulated strands, which results in a growing number of back-to-back piling hydrogen bonds across the strands, a slight decrease of the mean H-bond distance is apparent in parallel beta sheets whereas such trends are obviously not given in antiparallel beta sheets. This observation suggests that cooperative effects mutually polarizing spatially well-aligned hydrogen bonds are present either in alpha helices and parallel beta sheets whereas such influences seem to be lacking in 3(10) helices and antiparallel beta sheets.     

Importance of secondary structural specificity determinants in protein folding: insertion of a native beta-sheet sequence into an alpha-helical coiled-coil

To examine how a short secondary structural element derived from a native protein folds when in a different protein environment, we inserted an 11-residue beta-sheet segment (cassette) from human immunoglobulin fold, Fab new, into an alpha-helical coiled-coil host protein (cassette holder). This de novo design protein model, the structural cassette mutagenesis (SCM) model, allows us to study protein folding principles involving both short- and long-range interactions that affect secondary structure stability and conformation. In this study, we address whether the insertion of this beta-sheet cassette into the alpha-helical coiled-coil protein would result in conformational change nucleated by the long-range tertiary stabilization of the coiled-coil, therefore overriding the local propensity of the cassette to form beta-sheet, observed in its native immunoglobulin fold. The results showed that not only did the nucleating helices of the coiled-coil on either end of the cassette fail to nucleate the beta-sheet cassette to fold with an alpha-helical conformation, but also the entire chimeric protein became a random coil. We identified two determinants in this cassette that prevented coiled-coil formation: (1) a tandem dipeptide NN motif at the N-terminal of the beta-sheet cassette, and (2) the hydrophilic Ser residue, which would be buried in the hydrophobic core if the coiled-coil structure were to fold. By amino acid substitution of these helix disruptive residues, that is, either the replacement of the NN motif with high helical propensity Ala residues or the substitution of Ser with Leu to enhance hydrophobicity, we were able to convert the random coil chimeric protein into a fully folded alpha-helical coiled-coil. We hypothesized that this NN motif is a "secondary structural specificity determinant" which is very selective for one type of secondary structure and may prevent neighboring residues from adopting an alternate protein fold. These sequences with secondary structural specificity determinants have very strong local propensity to fold into a specific secondary structure and may affect overall protein folding by acting as a folding initiation site.     

A QSAR study on some series of anti-hepatitis B virus (HBV) agents

A quantitative structure-activity relationship (QSAR) study has been made on some series of anti-hepatitis B virus (HBV) agents, namely, a series of novel bis(L-amino acid) ester prodrugs of 9-[2--(phosphonomethoxy)ethyl]adenine, a similar series of compounds comprising of 2- amino-6-arylthio-9-[2-(phosphonoethoxy)ethyl] purine bis(2,2,2- trifluoroethyl) esters, and a series of 1-isopropylsulfonyl-2-amine benzimidazoles. In each case significant correlations are found between the anti-HBV potencies and some physicochemical and steric properties of the compounds, indicating that for the first two series the activity is controlled by the hydrophobic and the bulk properties of the molecules and, for the third series, the steric and hydrogen bonding properties of compounds are crucial for their anti-HBV potency.     

Interaction of CTX-M-15 enzyme with cefotaxime: a molecular modelling and docking study

Extended-spectrum β-lactamases (ESBLs) are the bacterial enzymes that make them resistant to advanced-generation cephalosporins. CTXM enzymes (the most prevalent ESBL-type) target cefotaxime. Aims of the study were: Modelling of CTX-M enzyme from bla(CTX-M) sequences of clinical Escherichia coli isolatesDocking of cefotaxime with modelled CTX-M enzymes to identify amino acid residues crucial to their interaction To hypothesize a possible relationship between 'interaction energy of the docked enzyme-antibiotic complex' and 'minimum inhibitory concentration (MIC) of the antibiotic against the bacteria producing that enzyme'. Seven E. coli strains of clinical origin which were confirmed as PCR-positive for bla(CTX-M) were selected for the study. C600 cells harboring cloned bla(CTX-M) were tested for ESBL-production by double-disk-synergy test. BLAST analysis confirmed all the bla(CTX-M) genes as blaCTX-M-15. Four of the 7 strains were found to be clonally related. Modelling was performed using Swiss Model Server. Discovery Studio 2.0 (Accelrys) was used to prepare Ramachandran plots for the modelled structures. Ramachandran Z-scores for modelled CTX-M enzymes from E. coli strains D8, D183, D253, D281, D282, D295 and D296 were found to be -0.449, 0.096, 0.027, 0.043, 0.032, -1.249 and -1.107, respectively. Docking was performed using Hex 5.1 and the results were further confirmed by Autodock 4.0. The amino acid residues Asn 104, Asn132, Gly 227, Thr 235, Gly 236, and Ser237 were found to be responsible for positioning cefotaxime into the active site of the CTX-M-15 enzyme. It was found that cefotaxime MICs for the CTX-M-15-producers increased with the increasing negative interaction energy of the enzyme-antibiotic complex.     

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.     

Designed molecules that fold to mimic protein secondary structures

Molecules that fold to mimic protein secondary structures have emerged as important targets of bioorganic chemistry. Recently, a variety of compounds that mimic helices, turns, and sheets have been developed, with notable advances in the design of beta-peptides that mimic each of these structures. These compounds hold promise as a step toward synthetic molecules with protein-like properties and as drugs that block protein-protein interactions.     

Analysis of interactive packing of secondary structural elements in alpha/beta units in proteins

An alpha-helix and a beta-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical alpha/beta units in protein structures, defined at < or = 2.5 A resolution, shows that the interaxial distance between the alpha-helix and the beta-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The beta-sheet unit (beta u), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded beta-strands, packing with an alpha-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in alpha/beta units using the interacting residue information of equivalent alpha/beta units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.     

EstB from Burkholderia gladioli: a novel esterase with a beta-lactamase fold reveals steric factors to discriminate between esterolytic and beta-lactam cleaving activity

Esterases form a diverse class of enzymes of largely unknown physiological role. Because many drugs and pesticides carry ester functions, the hydrolysis of such compounds forms at least one potential biological function. Carboxylesterases catalyze the hydrolysis of short chain aliphatic and aromatic carboxylic ester compounds. Esterases, D-alanyl-D-alanine-peptidases (DD-peptidases) and beta-lactamases can be grouped into two distinct classes of hydrolases with different folds and topologically unrelated catalytic residues, the one class comprising of esterases, the other one of beta-lactamases and DD-peptidases. The chemical reactivities of esters and beta-lactams towards hydrolysis are quite similar, which raises the question of which factors prevent esterases from displaying beta-lactamase activity and vice versa. Here we describe the crystal structure of EstB, an esterase isolated from Burkholderia gladioli. It shows the protein to belong to a novel class of esterases with homology to Penicillin binding proteins, notably DD-peptidase and class C beta-lactamases. Site-directed mutagenesis and the crystal structure of the complex with diisopropyl-fluorophosphate suggest Ser75 within the "beta-lactamase" Ser-x-x-Lys motif to act as catalytic nucleophile. Despite its structural homology to beta-lactamases, EstB shows no beta-lactamase activity. Although the nature and arrangement of active-site residues is very similar between EstB and homologous beta-lactamases, there are considerable differences in the shape of the active site tunnel. Modeling studies suggest steric factors to account for the enzyme's selectivity for ester hydrolysis versus beta-lactam cleavage.     

Long-range side-chain-main-chain interactions play crucial roles in stabilizing the (betaalpha)8 barrel motif of the alpha subunit of tryptophan synthase

The role of hither-to-fore unrecognized long-range hydrogen bonds between main-chain amide hydrogens and polar side chains on the stability of a well-studied (betaalpha)8, TIM barrel protein, the alpha subunit of tryptophan synthase (alphaTS), was probed by mutational analysis. The F19-D46 and I97-D124 hydrogen bonds link the N terminus of a beta-strand with the C terminus of the succeeding antiparallel alpha-helix, and the A103-D130 hydrogen bond links the N terminus of an alpha-helix with the C terminus of the succeeding antiparallel beta-strand, forming clamps for the respective betaalpha or alphabeta hairpins. The individual replacement of these aspartic acid side chains with alanine leads to what appear to be closely related partially folded structures with significantly reduced far-UV CD ellipticity and thermodynamic stability. Comparisons with the effects of eliminating another main-chain-side-chain hydrogen bond, G26-S33, and two electrostatic side-chain-side-chain hydrogen bonds, D38-H92 and D112-H146, all in the same N-terminal folding unit of alphaTS, demonstrated a unique role for the clamp interactions in stabilizing the native barrel conformation. Because neither the asparagine nor glutamic acid variant at position 46 can completely reproduce the spectroscopic, thermodynamic, or kinetic folding properties of aspartic acid, both size and charge are crucial to its unique role in the clamp hydrogen bond. Kinetic studies suggest that the three clamp hydrogen bonds act in concert to stabilize the transition state leading to the fully folded TIM barrel motif.     

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.     

DichroWeb,a website for calculating protein secondary structure from circular dichroism spectroscopic data

Circular dichroism (CD) spectroscopy is a widely‐used method for characterizing the secondary structures of proteins. The well‐established and highly used analysis website, DichroWeb (located at: http://dichroweb.cryst.bbk.ac.uk/html/home.shtml) enables the facile quantitative determination of helix, sheet, and other secondary structure contents of proteins based on their CD spectra. DichroWeb includes a range of reference datasets and algorithms, plus graphical and quantitative methods for determining the quality of the analyses produced. This article describes the current website content, usage and accessibility, as well as the many upgraded features now present in this highly popular tool that was originally created nearly two decades ago.     

Crystal structure of the predicted phospholipase LYPLAL1 reveals unexpected functional plasticity despite close relationship to acyl protein thioesterases

Sequence homology indicates the existence of three human cytosolic acyl protein thioesterases, including APT1 that is known to depalmitoylate H- and N-Ras. One of them is the lysophospholipase-like 1 (LYPLAL1) protein that on the one hand is predicted to be closely related to APT1 but on the other hand might also function as a potential triacylglycerol lipase involved in obesity. However, its role remained unclear. The 1.7 Å crystal structure of LYPLAL1 reveals a fold very similar to APT1, as expected, but features a shape of the active site that precludes binding of long-chain substrates. Biochemical data demonstrate that LYPLAL1 exhibits neither phospholipase nor triacylglycerol lipase activity, but rather accepts short-chain substrates. Furthermore, extensive screening efforts using chemical array technique revealed a first small molecule inhibitor of LYPLAL1.     

PFIT and PFRIT: bioinformatic algorithms for detecting glycosidase function from structure and sequence

The identification of the enzymes involved in the metabolism of simple and complex carbohydrates presents one bioinformatic challenge in the post-genomic era. Here, we present the PFIT and PFRIT algorithms for identifying those proteins adopting the alpha/beta barrel fold that function as glycosidases. These algorithms are based on the observation that proteins adopting the alpha/beta barrel fold share positions in their tertiary structures having equivalent sets of atomic interactions. These are conserved tertiary interaction positions, which have been implicated in both structure and function. Glycosidases adopting the alpha/beta barrel fold share more conserved tertiary interactions than alpha/beta barrel proteins having other functions. The enrichment pattern of conserved tertiary interactions in the glycosidases is the information that PFIT and PFRIT use to predict whether any given alpha/beta barrel will function as a glycosidase or not. Using as a test set a database of 19 glycosidase and 45 nonglycosidase alpha/beta barrel proteins with low sequence similarity, PFIT and PFRIT can correctly predict glycosidase function for 84% of the proteins known to function as glycosidases. PFIT and PFRIT incorrectly predict glycosidase function for 25% of the nonglycosidases. The program PSI-BLAST can also correctly identify 84% of the 19 glycosidases, however, it incorrectly predicts glycosidase function for 50% of the nonglycosidases (twofold greater than PFIT and PFRIT). Overall, we demonstrate that the structure-based PFIT and PFRIT algorithms are both more selective and sensitive for predicting glycosidase function than the sequence-based PSI-BLAST algorithm.     

A high-accuracy protein structural class prediction algorithm using predicted secondary structural information

One major problem with the existing algorithm for the prediction of protein structural classes is low accuracies for proteins from α/β and α+β classes. In this study, three novel features were rationally designed to model the differences between proteins from these two classes. In combination with other rational designed features, an 11-dimensional vector prediction method was proposed. By means of this method, the overall prediction accuracy based on 25PDB dataset was 1.5% higher than the previous best-performing method, MODAS. Furthermore, the prediction accuracy for proteins from α+β class based on 25PDB dataset was 5% higher than the previous best-performing method, SCPRED. The prediction accuracies obtained with the D675 and FC699 datasets were also improved.     

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.     

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.     

Structural characterization of binding mode of smoking cessation drugs to nicotinic acetylcholine receptors through study of ligand complexes with acetylcholine-binding protein

Smoking cessation is an important aim in public health worldwide as tobacco smoking causes many preventable deaths. Addiction to tobacco smoking results from the binding of nicotine to nicotinic acetylcholine receptors (nAChRs) in the brain, in particular the α4β2 receptor. One way to aid smoking cessation is by the use of nicotine replacement therapies or partial nAChR agonists like cytisine or varenicline. Here we present the co-crystal structures of cytisine and varenicline in complex with Aplysia californica acetylcholine-binding protein and use these as models to investigate binding of these ligands binding to nAChRs. This analysis of the binding properties of these two partial agonists provides insight into differences with nicotine binding to nAChRs. A mutational analysis reveals that the residues conveying subtype selectivity in nAChRs reside on the binding site complementary face and include features extending beyond the first shell of contacting residues.