Protein structure refinement algorithms

Protein structure prediction has remained a major challenge in structural biology for more than half a century. Accelerated and cost efficient sequencing technologies have allowed researchers to sequence new organisms and discover new protein sequences. Novel protein structure prediction technologies will allow researchers to study the structure of proteins and to determine their roles in the underlying biology processes and develop novel therapeutics.

Difficulty of the problem stems from two folds: (a) describing the energy landscape that corresponds to the protein structure, commonly referred to as force field problem; and (b) sampling of the energy landscape, trying to find the lowest energy configuration that is hypothesized to be the native state of the structure in solution. The two problems are interweaved and they have to be solved simultaneously. This thesis is composed of three major contributions. In the first chapter we describe a novel high-resolution protein structure refinement algorithm called GRID. In the second chapter we present REMCGRID, an algorithm for generation of low energy decoy sets. In the third chapter, we present a machine learning approach to ranking decoys by incorporating coarse-grain features of protein structures.

Structure Prediction of G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.

In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.

To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.

Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.

DNA-Mediated Charge Transport Devices for Protein Detection

Detection of biologically relevant targets, including small molecules, proteins, DNA, and RNA, is vital for fundamental research as well as clinical diagnostics. Sensors with biological elements provide a natural foundation for such devices because of the inherent recognition capabilities of biomolecules. Electrochemical DNA platforms are simple, sensitive, and do not require complex target labeling or expensive instrumentation. Sensitivity and specificity are added to DNA electrochemical platforms when the physical properties of DNA are harnessed. The inherent structure of DNA, with its stacked core of aromatic bases, enables DNA to act as a wire via DNA-mediated charge transport (DNA CT). DNA CT is not only robust over long molecular distances of at least 34 nm, but is also especially sensitive to anything that perturbs proper base stacking, including DNA mismatches, lesions, or DNA-binding proteins that distort the π-stack. Electrochemical sensors based on DNA CT have previously been used for single-nucleotide polymorphism detection, hybridization assays, and DNA-binding protein detection. Here, improvements to (i) the structure of DNA monolayers and (ii) the signal amplification with DNA CT platforms for improved sensitivity and detection are described.

First, improvements to the control over DNA monolayer formation are reported through the incorporation of copper-free click chemistry into DNA monolayer assembly. As opposed to conventional film formation involving the self-assembly of thiolated DNA, copper-free click chemistry enables DNA to be tethered to a pre-formed mixed alkylthiol monolayer. The total amount of DNA in the final film is directly related to the amount of azide in the underlying alkylthiol monolayer. DNA monolayers formed with this technique are significantly more homogeneous and lower density, with a larger amount of individual helices exposed to the analyte solution. With these improved monolayers, significantly more sensitive detection of the transcription factor TATA binding protein (TBP) is achieved.

Using low-density DNA monolayers, two-electrode DNA arrays were designed and fabricated to enable the placement of multiple DNA sequences onto a single underlying electrode. To pattern DNA onto the primary electrode surface of these arrays, a copper precatalyst for click chemistry was electrochemically activated at the secondary electrode. The location of the secondary electrode relative to the primary electrode enabled the patterning of up to four sequences of DNA onto a single electrode surface. As opposed to conventional electrochemical readout from the primary, DNA-modified electrode, a secondary microelectrode, coupled with electrocatalytic signal amplification, enables more sensitive detection with spatial resolution on the DNA array electrode surface. Using this two-electrode platform, arrays have been formed that facilitate differentiation between well-matched and mismatched sequences, detection of transcription factors, and sequence-selective DNA hybridization, all with the incorporation of internal controls.

For effective clinical detection, the two working electrode platform was multiplexed to contain two complementary arrays, each with fifteen electrodes. This platform, coupled with low density DNA monolayers and electrocatalysis with readout from a secondary electrode, enabled even more sensitive detection from especially small volumes (4 μL per well). This multiplexed platform has enabled the simultaneous detection of two transcription factors, TBP and CopG, with surface dissociation constants comparable to their solution dissociation constants.

With the sensitivity and selectivity obtained from the multiplexed, two working electrode array, an electrochemical signal-on assay for activity of the human methyltransferase DNMT1 was incorporated. DNMT1 is the most abundant human methyltransferase, and its aberrant methylation has been linked to the development of cancer. However, current methods to monitor methyltransferase activity are either ineffective with crude samples or are impractical to develop for clinical applications due to a reliance on radioactivity. Electrochemical detection of methyltransferase activity, in contrast, circumvents these issues. The signal-on detection assay translates methylation events into electrochemical signals via a methylation-specific restriction enzyme. Using the two working electrode platform combined with this assay, DNMT1 activity from tumor and healthy adjacent tissue lysate were evaluated. Our electrochemical measurements revealed significant differences in methyltransferase activity between tumor tissue and healthy adjacent tissue.

As differential activity was observed between colorectal tumor tissue and healthy adjacent tissue, ten tumor sets were subsequently analyzed for DNMT1 activity both electrochemically and by tritium incorporation. These results were compared to expression levels of DNMT1, measured by qPCR, and total DNMT1 protein content, measured by Western blot. The only trend detected was that hyperactivity was observed in the tumor samples as compared to the healthy adjacent tissue when measured electrochemically. These advances in DNA CT-based platforms have propelled this class of sensors from the purely academic realm into the realm of clinically relevant detection.

I. The 3.7 Å crystal structure of horse heart ferricytochrome C. II. The application of the Karle-Hauptman tangent formula to protein phasing

I. The 3.7 Å Crystal Structure of Horse Heart Ferricytochrome C.

The crystal structure of horse heart ferricytochrome c has been determined to a resolution of 3.7 Å using the multiple isomorphous replacement technique. Two isomorphous derivatives were used in the analysis, leading to a map with a mean figure of merit of 0.458. The quality of the resulting map was extremely high, even though the derivative data did not appear to be of high quality.

Although it was impossible to fit the known amino acid sequence to the calculated structure in an unambiguous way, many important features of the molecule could still be determined from the 3.7 Å electron density map. Among these was the fact that cytochrome c contains little or no α-helix. The polypeptide chain appears to be wound about the heme group in such a way as to form a loosely packed hydrophobic core in the molecule.

The heme group is located in a cleft on the molecule with one edge exposed to the solvent. The fifth coordinating ligand is His 18 and the sixth coordinating ligand is probably neither His 26 nor His 33.

The high resolution analysis of cytochrome c is now in progress and should be completed within the next year.

II. The Application of the Karle-Hauptman Tangent Formula to Protein Phasing.

The Karle-Hauptman tangent formula has been shown to be applicable to the refinement of previously determined protein phases. Tests were made with both the cytochrome c data from Part I and a theoretical structure based on the myoglobin molecule. The refinement process was found to be highly dependent upon the manner in which the tangent formula was applied. Iterative procedures did not work well, at least at low resolution.

The tangent formula worked very well in selecting the true phase from the two possible phase choices resulting from a single isomorphous replacement phase analysis. The only restriction on this application is that the heavy atoms form a non-centric cluster in the unit cell.

Pages 156 through 284 in this Thesis consist of previously published papers relating to the above two sections. References to these papers can be found on page 155.

A new protein structure of P-II class snake venom metalloproteinases: it comprises metalloproteinase and disintegrin domains

A new metalloproteinase-disintegrin, named Jerdonitin, was purified from Trimeresurus jerdonii venom with a molecular weight of 36 kDa on SDS-PAGE. It dose-dependently inhibited ADP-induced human platelet aggregation with IC50 of 120 nM. cDNA cloning and sequencing revealed that Jerdonitin belonged to the class II of snake venom metalloproteinases (SVMPs) (P-II class). Different from other P-II class SVMPs, metalloproteinase and disintegrin domains of its natural protein were not separated, confirmed by internal peptide sequencing. Compared to other P-II class SVMPs, Jerdonitin has two additional cysteines (Cys219 and Cys238) located in the spacer domain and disintegrin domain, respectively. They probably form a disulfide bond and therefore the metalloproteinase and disintegrin domains cannot be separated by posttranslationally processing. In summary, comparison of the amino acid sequences of Jerdonitin with those of other P-II class SVMPs by sequence alignment and phylogenetic analysis, in conjunction with natural protein structure data, suggested that it was a new type of P-II class SVMPs. (C) 2003 Elsevier Inc. All rights reserved.

Protein three-dimensional structural databases: Domains, structurally aligned homologues and superfamilies

This paper reports the availability of a database of protein structural domains (DDBASE), an alignment database of homologous proteins (HOMSTRAD) and a database of structurally aligned superfamilies (CAMPASS) on the World Wide Web (WWW). DDBASE contains information on the organization of structural domains and their boundaries; it includes only one representative domain from each of the homologous families. This database has been derived by identifying the presence of structural domains in proteins on the basis of inter-secondary structural distances using the program DIAL [Sowdhamini & Blundell (1995), Protein Sci. 4, 506-520]. The alignment of proteins in superfamilies has been performed on the basis of the structural features and relationships of individual residues using the program COMPARER [Sali & Blundell (1990), J. Mol. Biol. 212, 403-428]. The alignment databases contain information on the conserved structural features in homologous proteins and those belonging to superfamilies. Available data include the sequence alignments in structure-annotated formats and the provision for viewing superposed structures of proteins using a graphical interface. Such information, which is freely accessible on the WWW, should be of value to crystallographers in the comparison of newly determined protein structures with previously identified protein domains or existing families.

Identification of protein superfamily from structure-based sequence motif

The structure-based sequence motif of the distant proteins in evolution, protein tyrosine phosphatases (PTP) I and II superfamilies, as an example, has been defined by the structural comparison, structure-based sequence alignment and analyses on substitut

Generation of recombinant monoclonal antibodies to study structure-function of envelope protein VP28 of white spot syndrome virus from shrimp

White spot syndrome virus (WSSV) is a major pathogen in shrimp aquaculture. VP28 is one of the most important envelope proteins of WSSV. In this study, a recombinant antibody library, as single-chain fragment variable (scFv) format, displayed on phage was constructed using mRNA from spleen cells of mice immunized with-full-length VP28 expressed in Escherichia coli. After several rounds of panning, six scFv antibodies specifically binding to the epitopes in the N-terminal, middle, and C-terminal regions of VP28, respectively, were isolated from the library. Using these scFv antibodies as tools, the epitopes in VP28 were located on the envelope of the virion by immuno-electron Microscopy, Neutralization assay with these antibodies in vitro suggested that these epitopes may not be the attachment site of WSSV to host cell receptor. This study provides a new way to investigate the structure and function of the envelope proteins of WSSV. (c) 2008 Published by Elsevier Inc.

Three dimensional solution structure of neurotoxin CM-11 from the Ophiophagus hannah

The Ophiophagus hannah (King Cobra) neurotoxin CM-11 is a small protein with 72 amino acid residues, Based on complete assignments of H-1-NMR resonances and determination of secondary structures of CM-11, 349 distance and 27 dihedral angle constraints including 19 phi's and 8 chi's were collected from NOESY and DQF-COSY , and the chemical stereospecific assignment of beta(1)H was partially achieved, Twelve structures with lower energy was obtained via metric matrix distance geometry and refinement with simulated annealing, These structures have a low RMSD of 0.14 nm for backbone atoms and 0.20 nm for heavy atoms, with no distance constraint violation more than 0.05 nm, and no dihedral angle violation more than 3 degrees.

General expression for probabilistic estimation of multiphase structure invariants in the case of a native protein and multiple derivatives. Application to estimates of the three-phase structure invariants

Concise probabilistic formulae with definite crystallographic implications are obtained from the distribution for eight three-phase structure invariants (3PSIs) in the case of a native protein and a heavy-atom derivative [Hauptman (1982). Acta Cryst. A38, 289-294] and from the distribution for 27 3PSIs in the case of a native and two derivatives [Fortier, Weeks & Hauptman (1984). Acta Cryst. A40, 646-651]. The main results of the probabilistic formulae for the four-phase structure invariants are presented and compared with those for the 3PSIs. The analysis directly leads to a general formula of probabilistic estimation for the n-phase structure invariants in the case of a native and m derivatives. The factors affecting the estimated accuracy of the 3PSIs are examined using the diffraction data from a moderate-sized protein. A method to estimate a set of the large-modulus invariants, each corresponding to one of the eight 3PSIs, that has the largest \Delta\ values and relatively large structure-factor moduli between the native and derivative is suggested, which remarkably improves the accuracy, and thus a phasing procedure making full use of all eight 3PSIs is proposed.

The genomic structure, alternative splicing and immune response of Chlamys farreri thioester-containing protein

MEP is a member of thioester-containing protein (TEP) family found in Zhikong scallop Chlamys farreri and is involved in innate immunity against invading microbes. In the present study, the genomic DNA of CfTEP was cloned and characterized. The genomic DNA sequence of CfTEP consisted of 40 exons and 39 introns spanning 35 kb with all exon-intron junction sequences agreeing with the GT/AG consensus. The genomic organization of CfTEP was similar to human and mouse 0 rather than ciona C3-1 and Drosophila dTEP2. By RT-PCR technique, seven different cDNA variants of CfTEP (designated as CfTEP-A-CfTEP-G) were cloned from scallop gonad. CfTEP-A-CfTEP-F were produced by alternative splicing of six mutually exclusive exons (exons 19-24), respectively, which encoded the highly variable central region. While in CfTEP-G, the deletion of all the six exons introduced a new translation stop site and might trigger nonsense mediated decay (NMD). The mRNA expression and the proportion of the seven CfTEP variant transcripts were examined in the gonad of scallops after bacterial challenge. The fragments containing the highly variable central region of UTEP were amplified by RT-PCR and a 100 positive clones were sequenced randomly. The expression profiles of the seven MEP variants were different and displayed the sex and bacteria dependent manner. In the blank, sea water and Listonella anguillarum challenged subgroups of male scallops, all the transcripts detected were CfTEP-G isoform. In the Micrococcus luteus challenged subgroup, the isoforms expressed and their proportions were CfTEP-F (54%), CfTEP-B (23%), CfTEP-A (10%), CfTEP-C (7%) and CfTEP-E (6%). However, in the gonad of female scallops, only CfTEP-A were found in blank and sea water challenged subgroups. After L anguillarum or M. luteus challenge, four and five isoforms were detected, respectively, with CfTEP-F isoform being the most one in the both subgroups. These results suggested that the evolution of TEP genes was very complex, and that the diverse CfTEP transcripts generated by alternative splicing played an important role as pattern recognition receptors in the innate immune defense of scallops. (C) 2009 Elsevier Ltd. All rights reserved.

Intracellular copper/zinc superoxide dismutase from bay scallop Argopecten irradians: Its gene structure, mRNA expression and recombinant protein

Superoxide dismutases are an ubiquitous family of enzymes that function to efficiently catalyze the dismutation of superoxide anions. Two unique and highly compartmentalized bay scallop Argopecten irradians superoxide dismutases: MnSOD and ecCuZnSOD, have been molecularly characterized in our previous study. To complete characterize the SOD family in A. irradians, a novel intracellular copper/zinc SOD from the A. irradians (Ai-icCuZnSOD) was obtained and characterized. The full-length cDNA of Ai-icCuZnSOD was 1047 bp with a 459 bp open reading frame encoding 152 amino acids. The genomic length of the Ai-icCuZnSOD gene was about 4279 bp containing 4 exons and 3 introns. The promoter region containing many putative transcription factor binding sites were analyzed. Furthermore, quantitative reverse transcriptase real-time PCR (qRT-PCR) analysis indicated that the highest expression of the Ai-icCuZnSOD was detected in gill and the expression profiles in hemocytes of bay scallops challenged with bacteria Vibrio anguillarum and lipopolysaccharide (LPS) were different. The result presented an increased expression after injection with LPS whereas no significant changes were observed after V. anguillarum injection. A fusion protein containing Ai-icCuZnSOD was produced in vitro. The rAi-icCuZnSOD is a stable enzyme, retaining more than 80% of its activity between 10 and 60 degrees C and keeping above 88% of its activity at pH values between 5.8 and 9. Ai-icCuZnSOD is more stable under alkaline than acidic conditions. Crown Copyright (C) 2009 Published by Elsevier Ltd. All rights reserved.