Gene recognition via spliced sequence alignment.

Gene recognition is one of the most important problems in computational molecular biology. Previous attempts to solve this problem were based on statistics, and applications of combinatorial methods for gene recognition were almost unexplored. Recent advances in large-scale cDNA sequencing open a way toward a new approach to gene recognition that uses previously sequenced genes as a clue for recognition of newly sequenced genes. This paper describes a spliced alignment algorithm and software tool that explores all possible exon assemblies in polynomial time and finds the multiexon structure with the best fit to a related protein. Unlike other existing methods, the algorithm successfully recognizes genes even in the case of short exons or exons with unusual codon usage; we also report correct assemblies for genes with more than 10 exons. On a test sample of human genes with known mammalian relatives, the average correlation between the predicted and actual proteins was 99%. The algorithm correctly reconstructed 87% of genes and the rare discrepancies between the predicted and real exon-intron structures were caused either by short (less than 5 amino acids) initial/terminal exons or by alternative splicing. Moreover, the algorithm predicts human genes reasonably well when the homologous protein is nonvertebrate or even prokaryotic. The surprisingly good performance of the method was confirmed by extensive simulations: in particular, with target proteins at 160 accepted point mutations (PAM) (25% similarity), the correlation between the predicted and actual genes was still as high as 95%.

Involvement of Asn-293 in stereospecific agonist recognition and in activation of the beta 2-adrenergic receptor.

To investigate the molecular mechanism for stereospecific binding of agonists to beta 2-adrenergic receptors we used receptor models to identify potential binding sites for the beta-OH-group of the ligand, which defines the chiral center. Ser-165, located in transmembrane helix IV, and Asn-293, situated in the upper half of transmembrane helix VI, were identified as potential binding sites. Mutation of Ser-165 to Ala did not change the binding of either isoproterenol isomer as revealed after transient expression in human embryonic kidney (HEK)-293 cells. In contrast, a receptor mutant in which Asn-293 was replaced by Leu showed substantial loss of stereospecific isoproterenol binding. Adenylyl cyclase stimulation by this mutant after stable expression in CHO cells confirmed the substantial loss of stereospecificity for isoproterenol. In a series of agonists the loss of affinity in the Leu-293 mutant receptor was strongly correlated with the intrinsic activity of the compounds. Full agonists showed a 10-30-fold affinity loss, whereas partial agonists had almost the same affinity for both receptors. Stereospecific recognition of antagonists was unaltered in the Leu-293 mutant receptor. These data indicate a relationship between stereospecificity and intrinsic activity of agonists and suggest that Asn-293 is important for both properties of the agonist-receptor interaction.

Molecular cloning of the first metazoan beta-1,3 glucanase from eggs of the sea urchin Strongylocentrotus purpuratus.

We report the molecular cloning of the first beta-1,3 glucanase from animal tissue. Three peptide sequences were obtained from beta-1,3 glucanase that had been purified from eggs of the sea urchin Strongylocentrotus purpuratus and the gene was cloned by PCR using oligonucleotides deduced from the peptide sequences. The full-length cDNA shows a predicted enzyme structure of 499 aa with a hydrophobic signal sequence. A 3.2-kb message is present in eggs, during early embryogenesis, and in adult gut tissue. A polyclonal antibody to the native 68-kDa enzyme recognizes a single band during early embryogenesis that reappears in the adult gut, and recognizes a 57-kDa fusion protein made from a full-length cDNA clone for beta-1,3 glucanase. The identity of this molecule as beta-1,3 glucanase is confirmed by sequence homology, by the presence of all three peptide sequences in the deduced amino acid sequence, and by the recognition of the bacterial fusion protein by the antibody directed against the native enzyme. Data base searches show significant homology at the amino acid level to beta-1,3 glucanases from two species of bacteria and a clotting factor from the horseshoe crab. The homology with the bacteria is centered in a 304-aa region in which there are seven scattered regions of high homology between the four divergent species. These four species were also found to have two homologous regions in common with more distantly related plant, fungal, and bacterial proteins. A global phylogeny based on these regions strongly suggests that the glucanases are a very ancient family of genes. In particular, there is an especially deep split within genes taken from the bacterial genus Bacillus.

Calicheamicin-DNA complexes: warhead alignment and saccharide recognition of the minor groove.

The solution structures of calicheamicin gamma 1I, its cycloaromatized analog (calicheamicin epsilon), and its aryl tetrasaccharide complexed to a common DNA hairpin duplex have been determined by NMR and distance-refined molecular dynamics computations. Sequence specificity is associated with carbohydrate-DNA recognition that places the aryl tetrasaccharide component of all three ligands in similar orientations in the minor groove at the d(T-C-C-T).d(A-G-G-A) segment. The complementary fit of the ligands and the DNA minor groove binding site creates numerous van der Waals contacts as well as hydrogen bonding interactions. Notable are the iodine and sulfur atoms of calicheamicin that hydrogen bond with the exposed amino proton of the 5'- and 3'-guanines, respectively, of the d(A-G-G-A) segment. The sequence-specific carbohydrate binding orients the enediyne aglycone of calicheamicin gamma 1I such that its C3 and C6 proradical centers are adjacent to the cleavage sites. While the enediyne aglycone of calicheamicin gamma 1I is tilted relative to the helix axis and spans the minor groove, the cycloaromatized aglycone is aligned approximately parallel to the helix axis in the respective complexes. Specific localized conformational perturbations in the DNA have been identified from imino proton complexation shifts and changes in specific sugar pucker patterns on complex formation. The helical parameters for the carbohydrate binding site are comparable with corresponding values in B-DNA fibers while a widening of the groove is observed at the adjacent aglycone binding site.

Cyclic polyamides for recognition in the minor groove of DNA.

Small molecules that specifically bind with high affinity to any designated DNA sequence in the human genome would be useful tools in molecular biology and potentially in human medicine. Simple rules have been developed to rationally alter the sequence specificity of minor groove-binding polyamides containing N-methylimidazole and N-methylpyrrole amino acids. Crescent-shaped polyamides bind as antiparallel dimers with each polyamide making specific contacts with each strand on the floor of the minor groove. Cyclic polyamides have now been synthesized that bind designated DNA sequences at subnanomolar concentrations.

Mimics of the binding sites of opioid receptors obtained by molecular imprinting of enkephalin and morphine.

Molecular imprinting of morphine and the endogenous neuropeptide [Leu5]enkephalin (Leu-enkephalin) in methacrylic acid-ethylene glycol dimethacrylate copolymers is described. Such molecular imprints possess the capacity to mimic the binding activity of opioid receptors. The recognition properties of the resultant imprints were analyzed by radioactive ligand binding analysis. We demonstrate that imprinted polymers also show high binding affinity and selectivity in aqueous buffers. This is a major breakthrough for molecular imprinting technology, since the binding reaction occurs under conditions relevant to biological systems. The antimorphine imprints showed high binding affinity for morphine, with Kd values as low as 10(-7) M, and levels of selectivity similar to those of antibodies. Preparation of imprints against Leu-enkephalin was greatly facilitated by the use of the anilide derivative rather than the free peptide as the print molecule, due to improved solubility in the polymerization mixture. Free Leu-enkephalin was efficiently recognized by this polymer (Kd values as low as 10(-7) M were observed). Four tetra- and pentapeptides, with unrelated amino acid sequences, were not bound. The imprints showed only weak affinity for two D-amino acid-containing analogues of Leu-enkephalin. Enantioselective recognition of the L-enantiomer of phenylalanylglycine anilide, a truncated analogue of the N-terminal end of enkephalin, was observed.

A chloroplast homologue of the signal recognition particle subunit SRP54 is involved in the posttranslational integration of a protein into thylakoid membranes.

The mechanisms involved in the integration of proteins into the thylakoid membrane are largely unknown. However, many of the steps of this process for the light-harvesting chlorophyll a/b protein (LHCP) have been described and reconstituted in vitro. LHCP is synthesized as a precursor in the cytosol and posttranslationally imported into chloroplasts. Upon translocation across the envelope membranes, the N-terminal transit peptide is cleaved, and the apoprotein is assembled into a soluble "transit complex" and then integrated into the thylakoid membrane via three transmembrane helices. Here we show that 54CP, a chloroplast homologue of the 54-kDa subunit of the mammalian signal recognition particle (SRP54), is essential for transit complex formation, is present in the complex, and is required for LHCP integration into the thylakoid membrane. Our data indicate that 54CP functions posttranslationally as a molecular chaperone and potentially pilots LHCP to the thylakoids. These results demonstrate that one of several pathways for protein routing to the thylakoids is homologous to the SRP pathway and point to a common evolutionary origin for the protein transport systems of the endoplasmic reticulum and the thylakoid membrane.

Gene cloning, molecular modeling, and phylogenetics of serine protease P32 and serine carboxypeptidase SCP1 from nematophagous fungi Pochonia rubescens and Pochonia chlamydosporia

The fungi Pochonia chlamydosporia and Pochonia rubescens are parasites of nematode eggs and thus are biocontrol agents of nematodes. Proteolytic enzymes such as the S8 proteases VCP1 and P32, secreted during the pathogenesis of nematode eggs, are major virulence factors in these fungi. Recently, expression of these enzymes and of SCP1, a new putative S10 carboxypeptidase, was detected during endophytic colonization of barley roots by these fungi. In our study, we cloned the genomic and mRNA sequences encoding P32 from P. rubescens and SCP1 from P. chlamydosporia. P32 showed a high homology with the serine proteases Pr1A from the entomopathogenic fungus Metarhizium anisopliae and VCP1 from P. chlamydosporia (86% and 76% identity, respectively). However, the catalytic pocket of P32 showed differences in the amino acids of the substrate-recognition sites compared with the catalytic pockets of Pr1A and VCP1 proteases. Phylogenetic analysis of P32 suggests a common ancestor with protease Pr1A. SCP1 displays the characteristic features of a member of the S10 family of serine proteases. Phylogenetic comparisons show that SCP1 and other carboxypeptidases from filamentous fungi have an origin different from that of yeast vacuolar serine carboxypeptidases. Understanding protease genes from nematophagous fungi is crucial for enhancing the biocontrol potential of these organisms.

The Structure of theCyprinid herpesvirus 3ORF112-Zα·Z-DNA Complex Reveals a Mechanism of Nucleic Acids Recognition Conserved with E3L, a Poxvirus Inhibitor of Interferon Response

In vertebrate species, the innate immune system down-regulates protein translation in response to viral infection through the action of the double-stranded RNA (dsRNA)-activated protein kinase (PKR). In some teleost species another protein kinase, Z-DNA-dependent protein kinase (PKZ), plays a similar role but instead of dsRNA binding domains, PKZ has Zα domains. These domains recognize the left-handed conformer of dsDNA and dsRNA known as Z-DNA/Z-RNA. Cyprinid herpesvirus 3 infects common and koi carp, which have PKZ, and encodes the ORF112 protein that itself bears a Zα domain, a putative competitive inhibitor of PKZ. Here we present the crystal structure of ORF112-Zα in complex with an 18-bp CpG DNA repeat, at 1.5 Å. We demonstrate that the bound DNA is in the left-handed conformation and identify key interactions for the specificity of ORF112. Localization of ORF112 protein in stress granules induced in Cyprinid herpesvirus 3-infected fish cells suggests a functional behavior similar to that of Zα domains of the interferon-regulated, nucleic acid surveillance proteins ADAR1 and DAI.

hnRNP A2, a potential ssDNA/RNA molecular adapter at the telomere

The heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multi-tasking protein that acts in the cytoplasm and nucleus. We have explored the possibility that this protein is associated with telomeres and participates in their maintenance. Rat brain hnRNP A2 was shown to have two nucleic acid binding sites. In the presence of heparin one site binds single-stranded oligodeoxyribonucleotides irrespective of sequence but not the corresponding oligoribonucleotides. Both the hnRNP A2-binding cis-acting element for the cytoplasmic RNA trafficking element, A2RE, and the ssDNA telomere repeat match a consensus sequence for binding to a second sequence-specific site identified by mutational analysis. hnRNP A2 protected the telomeric repeat sequence, but not the complementary sequence, against DNase digestion: the glycine-rich domain was found to be necessary, but not sufficient, for protection. The N-terminal RRM (RNA recognition motif) and tandem RRMs of hnRNP A2 also bind the single-stranded, template-containing segment of telomerase RNA. hnRNP A2 colocalizes with telomeric chromatin in the subset of PML bodies that are a hallmark of ALT cells, reinforcing the evidence for hnRNPs having a role in telomere maintenance. Our results support a model in which hnRNP A2 acts as a molecular adapter between single-stranded telomeric repeats, or telomerase RNA, and another segment of ssDNA.

Modelling the structure of latexin-carboxypeptidase A complex based on chemical cross-linking and molecular docking

We have determined the three-dimensional structure of the protein complex between latexin and carboxypeptidase A using a combination of chemical cross-linking, mass spectrometry and molecular docking. The locations of three intermolecular cross-links were identified using mass spectrometry and these constraints were used in combination with a speed-optimised docking algorithm allowing us to evaluate more than 3 x 10(11) possible conformations. While cross-links represent only limited structural constraints, the combination of only three experimental cross-links with very basic molecular docking was sufficient to determine the complex structure. The crystal structure of the complex between latexin and carboxypeptidase A4 determined recently allowed us to assess the success of this structure determination approach. Our structure was shown to be within 4 angstrom r.m.s. deviation of C alpha atoms of the crystal structure. The study demonstrates that cross-linking in combination with mass spectrometry can lead to efficient and accurate structural modelling of protein complexes.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

Peptide binding prediction for the human class II MHC allele HLA-DP2:a molecular docking approach

MHC class II proteins bind oligopeptide fragments derived from proteolysis of pathogen antigens, presenting them at the cell surface for recognition by CD4+ T cells. Human MHC class II alleles are grouped into three loci: HLA-DP, HLA-DQ and HLA-DR. In contrast to HLA-DR and HLA-DQ, HLA-DP proteins have not been studied extensively, as they have been viewed as less important in immune responses than DRs and DQs. However, it is now known that HLA-DP alleles are associated with many autoimmune diseases. Quite recently, the X-ray structure of the HLA-DP2 molecule (DPA*0103, DPB1*0201) in complex with a self-peptide derived from the HLA-DR a-chain has been determined. In the present study, we applied a validated molecular docking protocol to a library of 247 modelled peptide-DP2 complexes, seeking to assess the contribution made by each of the 20 naturally occurred amino acids at each of the nine binding core peptide positions and the four flanking residues (two on both sides).