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.

Sequence and Structure Based Protein Folding Studies With Implications

As the expression of the genetic blueprint, proteins are at the heart of all biological systems. The ever increasing set of available protein structures has taught us that diversity is the hallmark of their architecture, a fundamental characteristic that enables them to perform the vast array of functionality upon which all of life depends. This diversity, however, is central to one of the most challenging problems in molecular biology: how does a folding polypeptide chain navigate its way through all of the myriad of possible conformations to find its own particular biologically active form? With few overarching structural principles to draw upon that can be applied to all protein architecture, the search for a solution to the protein folding problem has yet to produce an algorithm that can explain and duplicate this fundamental biological process. In this thesis, we take a two-pronged approach for investigating the protein folding process. Our initial statistical studies of the distributions of hydrophobic and hydrophilic residues within α-helices and β-sheets suggest (i) that hydrophobicity plays a critical role in helix and sheet formation; and (ii) that the nucleation of these motifs may result in largely unidirectional growth. Most tellingly, from an examination of the amino acids found in the smallest β-sheets, we do not find any evidence of a β-nucleating code in the primary protein sequence. Complementing these statistical analyses, we have analyzed the structural environments of several ever-widening aspects of protein topology. Our examination of the gaps between strands in the smallest β-sheets reveals a common organizational principle underlying β-formation involving strands separated by large sequential gaps: with very few exceptions, these large gaps fold into single, compact structural modules, bringing the β-strands that are otherwise far apart in the sequence close together in space. We conclude, therefore, that β-nucleation in the smallest sheets results from the co-location of two strands that are either local in sequence, or local in space following prior folding events. A second study of larger β-sheets both corroborates and extends these findings: virtually all large sequential gaps between pairs of β-strands organize themselves into an hierarchical arrangement, creating a bread-crumb model of go-and-come-back structural organization that ultimately juxtaposes two strands of a parental β-structure that are far apart in the sequence in close spatial proximity. In a final study, we have formalized this go-and-come-back notion into the concept of anti-parallel double-strandedness (DS), and measure this property across protein architecture in general. With over 90% of all residues in a large, non-redundant set of protein structures classified as DS, we conclude that DS is a unifying structural principle that underpins all globular proteins. We postulate, moreover, that this one simple principle, anti-parallel double-strandedness, unites protein structure, protein folding and protein evolution.

Unraveling the Structure and Function of G Protein-Coupled Receptors Through NMR Spectroscopy

G protein-coupled receptors (GPCRs) are a large superfamily of signaling proteins expressed on the plasma membrane. They are involved in a wide range of physiological processes and, therefore, are exploited as drug targets in a multitude of therapeutic areas. In this extent, knowledge of structural and functional properties of GPCRs may greatly facilitate rational design of modulator compounds. Solution and solid-state nuclear magnetic resonance (NMR) spectroscopy represents a powerful method to gather atomistic insights into protein structure and dynamics. In spite of the difficulties inherent the solution of the structure of membrane proteins through NMR, these methods have been successfully applied, sometimes in combination with molecular modeling, to the determination of the structure of GPCR fragments, the mapping of receptor-ligand interactions, and the study of the conformational changes associated with the activation of the receptors. In this review, we provide a summary of the NMR contributions to the study of the structure and function of GPCRs, also in light of the published crystal structures.

The structure of an unconventional HD-GYP protein from Bdellovibrio reveals the roles of conserved residues in this class of cyclic-di-GMP phosphodiesterases

UNLABELLED: Cyclic-di-GMP is a near-ubiquitous bacterial second messenger that is important in localized signal transmission during the control of various processes, including virulence and switching between planktonic and biofilm-based lifestyles. Cyclic-di-GMP is synthesized by GGDEF diguanylate cyclases and hydrolyzed by EAL or HD-GYP phosphodiesterases, with each functional domain often appended to distinct sensory modules. HD-GYP domain proteins have resisted structural analysis, but here we present the first structural representative of this family (1.28 Å), obtained using the unusual Bd1817 HD-GYP protein from the predatory bacterium Bdellovibrio bacteriovorus. Bd1817 lacks the active-site tyrosine present in most HD-GYP family members yet remains an excellent model of their features, sharing 48% sequence similarity with the archetype RpfG. The protein structure is highly modular and thus provides a basis for delineating domain boundaries in other stimulus-dependent homologues. Conserved residues in the HD-GYP family cluster around a binuclear metal center, which is observed complexed to a molecule of phosphate, providing information on the mode of hydroxide ion attack on substrate. The fold and active site of the HD-GYP domain are different from those of EAL proteins, and restricted access to the active-site cleft is indicative of a different mode of activity regulation. The region encompassing the GYP motif has a novel conformation and is surface exposed and available for complexation with binding partners, including GGDEF proteins.

IMPORTANCE: It is becoming apparent that many bacteria use the signaling molecule cyclic-di-GMP to regulate a variety of processes, most notably, transitions between motility and sessility. Importantly, this regulation is central to several traits implicated in chronic disease (adhesion, biofilm formation, and virulence gene expression). The mechanisms of cyclic-di-GMP synthesis via GGDEF enzymes and hydrolysis via EAL enzymes have been suggested by the analysis of several crystal structures, but no information has been available to date for the unrelated HD-GYP class of hydrolases. Here we present the multidomain structure of an unusual member of the HD-GYP family from the predatory bacterium Bdellovibrio bacteriovorus and detail the features that distinguish it from the wider structural family of general HD fold hydrolases. The structure reveals how a binuclear iron center is formed from several conserved residues and provides a basis for understanding HD-GYP family sequence requirements for c-di-GMP hydrolysis.

Magnitude differences in bioactive compounds, chemical functional groups, fatty acid profiles, nutrient degradation and digestion, molecular structure and metabolic characteristics of protein in newly developed yellow-seeded and black-seeded canola lines.

Recently, new lines of yellow-seeded (CS-Y) and black-seeded canola (CS-B) have been developed with chemical and structural alteration through modern breeding technology. However, no systematic study was found on the bioactive compounds, chemical functional groups, fatty acid profiles, inherent structure, nutrient degradation and absorption, or metabolic characteristics between the newly developed yellow- and black-seeded canola lines. This study aimed to systematically characterize chemical, structural, and nutritional features in these canola lines. The parameters accessed include bioactive compounds and antinutrition factors, chemical functional groups, detailed chemical and nutrient profiles, energy value, nutrient fractions, protein structure, degradation kinetics, intestinal digestion, true intestinal protein supply, and feed milk value. The results showed that the CS-Y line was lower (P ≤ 0.05) in neutral detergent fiber (122 vs 154 g/kg DM), acid detergent fiber (61 vs 99 g/kg DM), lignin (58 vs 77 g/kg DM), nonprotein nitrogen (56 vs 68 g/kg DM), and acid detergent insoluble protein (11 vs 35 g/kg DM) than the CS-B line. There was no difference in fatty acid profiles except C20:1 eicosenoic acid content (omega-9) which was in lower in the CS-Y line (P < 0.05) compared to the CS-B line. The glucosinolate compounds differed (P < 0.05) in terms of 4-pentenyl, phenylethyl, 3-CH3-indolyl, and 3-butenyl glucosinolates (2.9 vs 1.0 μmol/g) between the CS-Y and CS-B lines. For bioactive compounds, total polyphenols tended to be different (6.3 vs 7.2 g/kg DM), but there were no differences in erucic acid and condensed tannins with averages of 0.3 and 3.1 g/kg DM, respectively. When protein was portioned into five subfractions, significant differences were found in PA, PB1 (65 vs 79 g/kg CP), PB2, and PC fractions (10 vs 33 g/kg CP), indicating protein degradation and supply to small intestine differed between two new lines. In terms of protein structure spectral profile, there were no significant differences in functional groups of amides I and II, α helix, and β-sheet structure as well as their ratio between the two new lines, indicating no difference in protein structure makeup and conformation between the two lines. In terms of energy values, there were significant differences in total digestible nutrient (TDN; 149 vs 133 g/kg DM), metabolizable energy (ME; 58 vs 52 MJ/kg DM), and net energy for lactation (NEL; 42 vs 37 MJ/kg DM) between CS-Y and CS-B lines. For in situ rumen degradation kinetics, the two lines differed in soluble fraction (S; 284 vs 341 g/kg CP), potential degradation fraction (D; 672 vs 590 g/kg CP), and effective degraded organic matter (EDOM; 710 vs 684 g/kg OM), but no difference in degradation rate. CS-Y had higher digestibility of rumen bypass protein in the intestine than CS-B (566 vs 446 g/kg of RUP, P < 0.05). Modeling nutrient supply results showed that microbial protein synthesis (MCP; 148 vs 171 g/kg DM) and rumen protein degraded balance (DPB; 108 vs 127 g/kg DM) were lower in the CS-Y line, but there were no differences in total truly digested protein in small intestine (DVE) and feed milk value (FMV) between the two lines. In conclusion, the new yellow line had different nutritional, chemical, and structural features compared to the black line. CS-Y provided better nutrient utilization and availability.

Nuclear magnetic resonance structure of the nucleic acid-binding domain of severe acute respiratory syndrome coronavirus nonstructural protein 3

The nuclear magnetic resonance (NMR) structure of a globular domain of residues 1071 to 1178 within the previously annotated nucleic acid-binding region (NAB) of severe acute respiratory syndrome coronavirus nonstructural protein 3 (nsp3) has been determined, and N- and C-terminally adjoining polypeptide segments of 37 and 25 residues, respectively, have been shown to form flexibly extended linkers to the preceding globular domain and to the following, as yet uncharacterized domain. This extension of the structural coverage of nsp3 was obtained from NMR studies with an nsp3 construct comprising residues 1066 to 1181 [ nsp3(1066-1181)] and the constructs nsp3(1066-1203) and nsp3(1035-1181). A search of the protein structure database indicates that the globular domain of the NAB represents a new fold, with a parallel four-strand beta-sheet holding two alpha-helices of three and four turns that are oriented antiparallel to the beta-strands. Two antiparallel two-strand beta-sheets and two 3(10)-helices are anchored against the surface of this barrel-like molecular core. Chemical shift changes upon the addition of single-stranded RNAs (ssRNAs) identified a group of residues that form a positively charged patch on the protein surface as the binding site responsible for the previously reported affinity for nucleic acids. This binding site is similar to the ssRNA-binding site of the sterile alpha motif domain of the Saccharomyces cerevisiae Vts1p protein, although the two proteins do not share a common globular fold.

Structure, Dynamics, and RNA Interaction Analysis of the Human SBDS Protein

Shwachman-Bodian-Diamond syndrome is an autosomal recessive genetic syndrome with pleiotropic phenotypes, including pancreatic deficiencies, bone marrow dysfunctions with increased risk of myelodysplasia or leukemia, and skeletal abnormalities. This syndrome has been associated with mutations in the SBDS gene, which encodes a conserved protein showing orthologs in Archaea and eukaryotes. The Shwachman-Bodian-Diamond syndrome pleiotropic phenotypes may be an indication of different cell type requirements for a fully functional SBDS protein. RNA-binding activity has been predicted for archaeal and yeast SBDS orthologs, with the latter also being implicated in ribosome biogenesis. However, full-length SBDS orthologs function in a species-specific manner, indicating that the knowledge obtained from model systems may be of limited use in understanding major unresolved issues regarding SBDS function, namely, the effect of mutations in human SBDS on its biochemical function and the specificity of RNA interaction. We determined the solution structure and backbone dynamics of the human SBDS protein and describe its RNA binding site using NMR spectroscopy. Similarly to the crystal structures of Archaea, the overall structure of human SBDS comprises three well-folded domains. However, significant conformational exchange was observed in NMR dynamics experiments for the flexible linker between the N-terminal domain and the central domain, and these experiments also reflect the relative motions of the domains. RNA titrations monitored by heteronuclear correlation experiments and chemical shift mapping analysis identified a classic RNA binding site at the N-terminal FYSH (fungal, Yhr087wp, Shwachman) domain that concentrates most of the mutations described for the human SBDS. (C) 2010 Elsevier Ltd. All rights reserved.

X-ray structure of the human hyperplastic discs protein: An ortholog of the C-terminal domain of poly(A)-binding protein

The poly(A)-binding protein (PABP) recognizes the 3′ mRNA poly(A) tail and plays an essential role in eukaryotic translation initiation and mRNA stabilization/degradation. PABP is a modular protein, with four N-terminal RNA-binding domains and an extensive C terminus. The C-terminal region of PABP is essential for normal growth in yeast and has been implicated in mediating PABP homo-oligomerization and protein–protein interactions. A small, proteolytically stable, highly conserved domain has been identified within this C-terminal segment. Remarkably, this domain is also present in the hyperplastic discs protein (HYD) family of ubiquitin ligases. To better understand the function of this conserved region, an x-ray structure of the PABP-like segment of the human HYD protein has been determined at 1.04-Å resolution. The conserved domain adopts a novel fold resembling a right-handed supercoil of four α-helices. Sequence profile searches and comparative protein structure modeling identified a small ORF from the Arabidopsis thaliana genome that encodes a structurally similar but distantly related PABP/HYD domain. Phylogenetic analysis of the experimentally determined (HYD) and homology modeled (PABP) protein surfaces revealed a conserved feature that may be responsible for binding to a PABP interacting protein, Paip1, and other shared interaction partners.

Fourier Transform Infrared Microspectroscopy Detects Changes in Protein Secondary Structure Associated with Desiccation Tolerance in Developing Maize Embryos1

Isolated immature maize (Zea mays L.) embryos have been shown to acquire tolerance to rapid drying between 22 and 25 d after pollination (DAP) and to slow drying from 18 DAP onward. To investigate adaptations in protein profile in association with the acquisition of desiccation tolerance in isolated, immature maize embryos, we applied in situ Fourier transform infrared microspectroscopy. In fresh, viable, 20- and 25-DAP embryo axes, the shapes of the different amide-I bands were identical, and this was maintained after flash drying. On rapid drying, the 20-DAP axes had a reduced relative proportion of α-helical protein structure and lost viability. Rapidly dried 25-DAP embryos germinated (74%) and had a protein profile similar to the fresh control axes. On slow drying, the α-helical contribution in both the 20- and 25-DAP embryo axes increased compared with that in the fresh control axes, and survival of desiccation was high. The protein profile in dry, mature axes resembled that after slow drying of the immature axes. Rapid drying resulted in an almost complete loss of membrane integrity in the 20-DAP embryo axes and much less so in the 25-DAP axes. After slow drying, low plasma membrane permeability ensued in both the 20- and 25-DAP axes. We conclude that slow drying of excised, immature embryos leads to an increased proportion of α-helical protein structures in their axes, which coincides with additional tolerance of desiccation stress.

West nile virus core protein: Tetramer structure and ribbon formation

We have determined the crystal structure of the core (C) protein from the Kunjin subtype of West Nile virus (WNV), closely related to the NY99 strain of WNV, currently a major health threat in the U.S. WNV is a member of the Flaviviridae family of enveloped RNA viruses that contains many important human pathogens. The C protein is associated with the RNA genome and forms the internal core which is surrounded by the envelope in the virion. The C protein structure contains four a. helices and forms dimers that are organized into tetramers. The tetramers form extended filamentous ribbons resembling the stacked alpha helices seen in HEAT protein structures.

The 2.1 angstrom crystal structure of the far-red fluorescent protein HcRed: Inherent conformational flexibility of the chromophore

We have determined the crystal structure of HcRed, a far-red fluorescent protein isolated from Heteractis crispa, to 2.1 resolution. HcRed was observed to form a dimer, in contrast to the monomeric form of green fluorescent protein (GFP) or the tetrameric forms of the GFP-like proteins (eqFP611, Rtms5 and DsRed). Unlike the well-defined chromophore conformation observed in GFP and the GFP-like proteins, the HcRed chromophore was observed to be considerably mobile. Within the HcRed structure, the cyclic tripeptide chromophore, Glu64-Tyr65-Gly66, was observed to adopt both a cis coplanar and a tran. non-coplanar conformation. As a result of these two con formations, the hydroxyphenyl moiety of the chromophore makes distinct interactions within the interior of the b-can. These data together with a quantum chemical model of the chromophore, suggest the cis coplanar conformation to be consistent with the fluorescent properties of HcRed, and the trans non-coplanar conformation to be consistent with non-fluorescent properties of hcCP, the chromoprotein parent of HcRed. Moreover, within the GFP-like family, it appears that where conformational freedom is permissible then flexibility in the chromophore conformation is possible. 2005 Elsevier Ltd. All rights reserved.

Prediction of cis/trans isomerization in proteins using PSI-BLAST profiles and secondary structure information

Background The majority of peptide bonds in proteins are found to occur in the trans conformation. However, for proline residues, a considerable fraction of Prolyl peptide bonds adopt the cis form. Proline cis/trans isomerization is known to play a critical role in protein folding, splicing, cell signaling and transmembrane active transport. Accurate prediction of proline cis/trans isomerization in proteins would have many important applications towards the understanding of protein structure and function. Results In this paper, we propose a new approach to predict the proline cis/trans isomerization in proteins using support vector machine (SVM). The preliminary results indicated that using Radial Basis Function (RBF) kernels could lead to better prediction performance than that of polynomial and linear kernel functions. We used single sequence information of different local window sizes, amino acid compositions of different local sequences, multiple sequence alignment obtained from PSI-BLAST and the secondary structure information predicted by PSIPRED. We explored these different sequence encoding schemes in order to investigate their effects on the prediction performance. The training and testing of this approach was performed on a newly enlarged dataset of 2424 non-homologous proteins determined by X-Ray diffraction method using 5-fold cross-validation. Selecting the window size 11 provided the best performance for determining the proline cis/trans isomerization based on the single amino acid sequence. It was found that using multiple sequence alignments in the form of PSI-BLAST profiles could significantly improve the prediction performance, the prediction accuracy increased from 62.8% with single sequence to 69.8% and Matthews Correlation Coefficient (MCC) improved from 0.26 with single local sequence to 0.40. Furthermore, if coupled with the predicted secondary structure information by PSIPRED, our method yielded a prediction accuracy of 71.5% and MCC of 0.43, 9% and 0.17 higher than the accuracy achieved based on the singe sequence information, respectively. Conclusion A new method has been developed to predict the proline cis/trans isomerization in proteins based on support vector machine, which used the single amino acid sequence with different local window sizes, the amino acid compositions of local sequence flanking centered proline residues, the position-specific scoring matrices (PSSMs) extracted by PSI-BLAST and the predicted secondary structures generated by PSIPRED. The successful application of SVM approach in this study reinforced that SVM is a powerful tool in predicting proline cis/trans isomerization in proteins and biological sequence analysis.

Ozonolysis of phospholipid double bonds during electrospray ionization : A new tool for structure determination

Phospholipids are the key structural component of cell membranes, and recent advances in electrospray ionization mass spectrometry provide for the fast and efficient analysis of these compounds in biological extracts.1-3 The application of electrospray ionization tandem mass spectrometry (ESI-MS/MS) to phospholipid analysis has demonstrated several key advantages over the more traditional chromatographic methods, including speed and greater structural information.4 For example, the ESI-MS/MS spectrum of a typical phospholipidsparticularly in negative ion modesreadily identifies the carbon chain length and the degree of unsaturation of each of the fatty acids esterified to the parent molecule.5 A critical limitation of conventional ESI-MS/MS analysis, however, is the inability to uniquely identify the position of double bonds within the fatty acid chains. This is especially problematic given the importance of double bond position in determining the biological function of lipid classes.6 Previous attempts to identify double bond position in intact phospholipids using mass spectrometry employ either MS3 or offline chemical derivatization.7-11 The former method requires specialized instrumentation and is rarely applied, while the latter methods suffer from complications inherent in sample handling prior to analysis. In this communication we outline a novel on-line approach for the identification of double bond position in intact phospholipids. In our method, the double bond(s) present in unsaturated phospholipids are cleaved by ozonolysis within the ion source of a conventional ESI mass spectrometer to give two chemically induced fragment ions that may be used to unambiguously assign the position of the double bond. This is achieved by using oxygen as the electrospray nebulizing gas in combination with high electrospray voltages to initiate the formation of an ozoneproducing.

Functionalized Ag nanoparticles with tunable optical properties for selective protein analysis

We present a preparation procedure for small sized biocompatibly coated Ag nanoparticles with tunable surface plasmon resonances. The conditions were optimised with respect to the resonance Raman signal enhancement of heme proteins and to the preservation of the native protein structure....