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

Structures of synthetic nanobody-SARS-CoV-2 receptor-binding domain complexes reveal distinct sites of interaction

Combating the worldwide spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the emergence of new variants demands understanding of the structural basis of the interaction of antibodies with the SARS-CoV-2 receptor-binding domain (RBD). Here, we report five X-ray crystal structures of sybodies (synthetic nanobodies) including those of binary and ternary complexes of Sb16-RBD, Sb45-RBD, Sb14-RBD-Sb68, and Sb45-RBD-Sb68, as well as unliganded Sb16. These structures reveal that Sb14, Sb16, and Sb45 bind the RBD at the angiotensin-converting enzyme 2 interface and that the Sb16 interaction is accompanied by a large conformational adjustment of complementarity-determining region 2. In contrast, Sb68 interacts at the periphery of the SARS-CoV-2 RBD-angiotensin-converting enzyme 2 interface. We also determined cryo-EM structures of Sb45 bound to the SARS-CoV-2 spike protein. Superposition of the X-ray structures of sybodies onto the trimeric spike protein cryo-EM map indicates that some sybodies may bind in both "up" and "down" configurations, but others may not. Differences in sybody recognition of several recently identified RBD variants are explained by these structures.

Synthetic nanobody-SARS-CoV-2 receptor-binding domain structures identify distinct epitopes

The worldwide spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) demands unprecedented attention. We report four X-ray crystal structures of three synthetic nanobodies (sybodies) (Sb16, Sb45 and Sb68) bind to the receptor-binding domain (RBD) of SARS-CoV-2: binary complexes of Sb16-RBD and Sb45-RBD; a ternary complex of Sb45-RBD-Sb68; and Sb16 unliganded. Sb16 and Sb45 bind the RBD at the ACE2 interface, positioning their CDR2 and CDR3 loops diametrically. Sb16 reveals a large CDR2 shift when binding the RBD. Sb68 interacts peripherally at the ACE2 interface; steric clashes with glycans explain its mechanism of viral neutralization. Superposing these structures onto trimeric spike (S) protein models indicates these sybodies bind conformations of the mature S protein differently, which may aid therapeutic design.

ONE SENTENCE SUMMARY

X-ray structures of synthetic nanobodies complexed with the receptor-binding domain of the spike protein of SARS-CoV-2 reveal details of CDR loop interactions in recognition of distinct epitopic sites.

Prediction of Protein Tertiary Structure via Regularized Template Classification Techniques

We discuss the use of the regularized linear discriminant analysis (LDA) as a model reduction technique combined with particle swarm optimization (PSO) in protein tertiary structure prediction, followed by structure refinement based on singular value decomposition (SVD) and PSO. The algorithm presented in this paper corresponds to the category of template-based modeling. The algorithm performs a preselection of protein templates before constructing a lower dimensional subspace via a regularized LDA. The protein coordinates in the reduced spaced are sampled using a highly explorative optimization algorithm, regressive–regressive PSO (RR-PSO). The obtained structure is then projected onto a reduced space via singular value decomposition and further optimized via RR-PSO to carry out a structure refinement. The final structures are similar to those predicted by best structure prediction tools, such as Rossetta and Zhang servers. The main advantage of our methodology is that alleviates the ill-posed character of protein structure prediction problems related to high dimensional optimization. It is also capable of sampling a wide range of conformational space due to the application of a regularized linear discriminant analysis, which allows us to expand the differences over a reduced basis set.

Computational Redesign of PD-1 Interface for PD-L1 Ligand Selectivity

Chronic or persistent stimulation of the programmed cell death-1 (PD-1) pathway prevents T cells from mounting anti-tumor and anti-viral immune responses. Blockade of this inhibitory checkpoint pathway has shown therapeutic importance by rescuing T cells from their exhausted state. Cognate ligands of the PD-1 receptor include the tissue-specific PD-L1 and PD-L2 proteins. Engineering a human PD-1 interface specific for PD-L1 or PD-L2 can provide a specific reagent and therapeutic advantage for tissue-specific disruption of the PD-1 pathway. We utilized ProtLID, a computational framework, which constitutes a residue-based pharmacophore approach, to custom-design a human PD-1 interface specific to human PD-L1 without any significant affinity to PD-L2. In subsequent cell assay experiments, half of all single-point mutant designs proved to introduce a statistically significant selectivity, with nine of these maintaining a close to wild-type affinity to PD-L1. This proof-of-concept study suggests a general approach to re-engineer protein interfaces for specificity.

Substitutions that lock and unlock the proton-coupled folate transporter (PCFT-SLC46A1) in an inward-open conformation

The proton-coupled folate transporter (PCFT) mediates intestinal absorption of folates and their transport from blood to cerebrospinal fluid across the choroid plexus. Substitutions at Asp-109 in the first intracellular loop between the first and second transmembrane domains (TMDs) abolish PCFT function, but protein expression and trafficking to the cell membrane are retained. Here, we used site-directed mutagenesis, the substituted-cysteine accessibility method, functional analyses, and homology modeling to determine whether the D109A substitution locks PCFT in one of its conformational states. Cys-substituted residues lining the PCFT aqueous translocation pathway and accessible in WT PCFT to the membrane-impermeable cysteine-biotinylation reagent, MTSEA-biotin, lost accessibility when introduced into the D109A scaffold. Substitutions at Gly-305 located exofacially within the eighth TMD, particularly with bulky residues, when introduced into the D109A scaffold largely restored function and MTSEA-biotin accessibility to Cys-substituted residues within the pathway. Likewise, Ser-196 substitution in the fifth TMD, predicted by homology modeling to be in proximity to Gly-305, also partially restored function found in solute transporters, is critical to oscillation of the carrier among its conformational states. Substitutions at Asp-109 and Gly-112 lock PCFT in an inward-open conformation, resulting in the loss of function. However, the integrity of the locked protein is preserved, indicated by the restoration of function after insertion of a second "unlocking" mutation. and accessibility. Similarly, the inactivating G112K substitution within the first intracellular loop was partially reactivated by introducing the G305L substitution. These data indicate that the first intracellular loop, with a sequence identical to "motif A" (GXXXDXXGR(R/K)).

Chemical disruption of the pyroptotic pore-forming protein gasdermin D inhibits inflammatory cell death and sepsis

Dysregulation of inflammatory cell death is a key driver of many inflammatory diseases. Pyroptosis, a highly inflammatory form of cell death, uses intracellularly generated pores to disrupt electrolyte homeostasis and execute cell death. Gasdermin D, the pore-forming effector protein of pyroptosis, coordinates membrane lysis and the release of highly inflammatory molecules, such as interleukin-1β, which potentiate the overactivation of the innate immune response. However, to date, there is no pharmacologic mechanism to disrupt pyroptosis. Here, we identify necrosulfonamide as a direct chemical inhibitor of gasdermin D, the pyroptotic pore-forming protein, which binds directly to gasdermin D to inhibit pyroptosis. Pharmacologic inhibition of pyroptotic cell death by necrosulfonamide is efficacious in sepsis models and suggests that gasdermin D inhibitors may be efficacious clinically in inflammatory diseases.

Hereditary folate malabsorption due to a mutation in the external gate of the proton-coupled folate transporter SLC46A1

Hereditary folate malabsorption (HFM) is an autosomal recessive disorder characterized by impaired intestinal folate absorption and impaired folate transport across the choroid plexus due to loss of function of the proton-coupled folate transporter (PCFT-SLC46A1). We report a novel mutation, causing HFM, affecting a residue located in the 11th transmembrane helix within the external gate. The mutant N411K-PCFT was stable, trafficked to the cell membrane, and had sufficient residual activity to characterize the transport defect and the structural requirements at this site for gate function. The influx V_max of the N411K mutant was markedly decreased, as was the affinity for most, but not all, folate/antifolate substrates. The greatest loss of activity was for 5-methyltetrahydrofolate. Substitutions with positive charged residues resulted in a loss of activity (arginine > lysine > histidine). Function was retained for the negative charged aspartate, but not the larger glutamate substitutions, whereas the bulky hydrophobic (leucine), or polar (glutamine) substitutions, were tolerated. Homology models of PCFT, in the inward and outward open conformations, based upon the mammalian Glut5 fructose transporter structures, localize Asn411 protruding into the aqueous pathway. This is most prominent when the carrier is in the inward open conformation when the external gate is closed. Mutations at this site likely result in highly specific steric and electrostatic interactions between the Asn411-substituted, and other, residues in the gate region that impede carrier function. The substrate specificity of the N411K mutant may be due to alterations of substrate flows through the external gate, downstream allosteric alterations in the folate-binding pocket, or both.

Substituted-cysteine accessibility and cross-linking identify an exofacial cleft in the 7th and 8th helices of the proton-coupled folate transporter (SLC46A1)

The proton-coupled folate transporter (PCFT-SLC46A1) is required for folate transport across the apical membrane of the small intestine and across the choroid plexus. This study focuses on the structure/function of the 7th transmembrane domain (TMD), and its relationship to the 8th TMD as assessed by the substituted cysteine accessibility method (SCAM) and dicysteine cross-linking. Nine exofacial residues (I278C; H281C-L288C) of 23 residues in the 7th TMD were accessible to 2-((biotinoyl)amino)ethyl methanethiosulfonate (MTSEA-biotin). Pemetrexed, a high-affinity substrate for PCFT, decreased or abolished biotinylation of seven of these residues consistent with their location in or near the folate binding pocket. Homology models of PCFT based on Glut5 fructose transporter structures in both inward- and outward- open conformations were constructed and predicted that two pairs of residues (T289-I304C and Q285-Q311C) from the 7th and 8th TMDs should be in sufficiently close proximity to form a disulfide bond when substituted with cysteines. The single Cys-substituted mutants were accessible to MTSEA-biotin and functional with and without pretreatment with dithiotreitol. However, the double mutants were either not accessible at all, or accessibility was markedly reduced and function markedly impaired. This occurred spontaneously without inclusion of an oxidizing agent. Dithiotreitol restored accessibility and function consistent with disulfide bond disruption. The data establish the proximity of exofacial regions of the 7th and 8th TMDs and their role in defining the aqueous translocation pathway and suggest that these helices may be a component of an exofacial cleft through which substrates enter the protein binding pocket in its outward-open conformation.

Residues in the eighth transmembrane domain of the proton-coupled folate transporter (SLC46A1) play an important role in defining the aqueous translocation pathway and in folate substrate binding

The proton-coupled folate transporter (PCFT-SLC46A1) is required for intestinal folate absorption and folate transport across the choroid plexus. This report addresses the structure/function of the 8th transmembrane helix. Based upon biotinylation of cysteine-substituted residues by MTSEA-biotin, 14 contiguous exofacial residues to Leu316 were accessible to the extracellular compartment of the 23 residues in this helix (Leu303-Leu325). Pemetrexed blocked biotinylation of six Cys-substituted residues deep within the helix implicating an important role for this region in folate binding. Accessibility decreased at 4°C vs RT. The influx K_t, K_i and V_max were markedly increased for the P314C mutant, similar to what was observed for Y315A and Y315P mutants. However, the K_t, alone, was increased for the P314Y mutant. To correlate these observations with PCFT structural changes during the transport cycle, homology models were built for PCFT based upon the recently reported structures of bovine and rodent GLUT5 fructose transporters in the inward-open and outward- open conformations, respectively. The models predict substantial structural alterations in the exofacial region of the eighth transmembrane helix as it cycles between its conformational states that can account for the extended and contiguous aqueous accessibility of this region of the helix. Further, a helix break in one of the two conformations can account for the critical roles Pro314 and Tyr315, located in this region, play in PCFT function. The data indicates that the 8th transmembrane helix of PCFT plays an important role in defining the aqueous channel and the folate binding pocket.

Live-cell visualization of gasdermin D-driven pyroptotic cell death

Pyroptosis is a form of cell death important in defenses against pathogens that can also result in a potent and sometimes pathological inflammatory response. During pyroptosis, GSDMD (gasdermin D), the pore-forming effector protein, is cleaved, forms oligomers, and inserts into the membranes of the cell, resulting in rapid cell death. However, the potent cell death induction caused by GSDMD has complicated our ability to understand the biology of this protein. Studies aimed at visualizing GSDMD have relied on expression of GSDMD fragments in epithelial cell lines that naturally lack GSDMD expression and also lack the proteases necessary to cleave GSDMD. In this work, we performed mutagenesis and molecular modeling to strategically place tags and fluorescent proteins within GSDMD that support native pyroptosis and facilitate live-cell imaging of pyroptotic cell death. Here, we demonstrate that these fusion proteins are cleaved by caspases-1 and -11 at Asp-276. Mutations that disrupted the predicted p30-p20 autoinhibitory interface resulted in GSDMD aggregation, supporting the oligomerizing activity of these mutations. Furthermore, we show that these novel GSDMD fusions execute inflammasome-dependent pyroptotic cell death in response to multiple stimuli and allow for visualization of the morphological changes associated with pyroptotic cell death in real time. This work therefore provides new tools that not only expand the molecular understanding of pyroptosis but also enable its direct visualization.

cDNA cloning, molecular modeling and docking calculations of L-type lectins from Swartzia simplex var. grandiflora (Leguminosae, Papilionoideae), a member of the tribe Swartzieae

The genus Swartzia is a member of the tribe Swartzieae, whose genera constitute the living descendants of one of the early branches of the papilionoid legumes. Legume lectins comprise one of the main families of structurally and evolutionarily related carbohydrate-binding proteins of plant origin. However, these proteins have been poorly investigated in Swartzia and to date, only the lectin from S. laevicarpa seeds (SLL) has been purified. Moreover, no sequence information is known from lectins of any member of the tribe Swartzieae. In the present study, partial cDNA sequences encoding L-type lectins were obtained from developing seeds of S. simplex var. grandiflora. The amino acid sequences of the S. simplex grandiflora lectins (SSGLs) were only averagely related to the known primary structures of legume lectins, with sequence identities not greater than 50-52%. The SSGL sequences were more related to amino acid sequences of papilionoid lectins from members of the tribes Sophoreae and Dalbergieae and from the Cladratis and Vataireoid clades, which constitute with other taxa, the first branching lineages of the subfamily Papilionoideae. The three-dimensional structures of 2 representative SSGLs (SSGL-A and SSGL-E) were predicted by homology modeling using templates that exhibit the characteristic β-sandwich fold of the L-type lectins. Molecular docking calculations predicted that SSGL-A is able to interact with D-galactose, N-acetyl-D-galactosamine and α-lactose, whereas SSGL-E is probably a non-functional lectin due to 2 mutations in the carbohydrate-binding site. Using molecular dynamics simulations followed by density functional theory calculations, the binding free energies of the interaction of SSGL-A with GalNAc and α-lactose were estimated as -31.7 and -47.5 kcal/mol, respectively. These findings gave insights about the carbohydrate-binding specificity of SLL, which binds to immobilized lactose but is not retained in a matrix containing D-GalNAc as ligand.

Distinct domains of the β1-subunit cytosolic N terminus control surface expression and functional properties of large-conductance calcium-activated potassium (BK) channels

The properties and function of large-conductance calcium- and voltage-activated potassium (BK) channels are modified by the tissue-specific expression of regulatory β1-subunits. Although the short cytosolic N-terminal domain of the β1-subunit is important for controlling both BK channel trafficking and function, whether the same, or different, regions of the N terminus control these distinct processes remains unknown. Here we demonstrate that the first six N-terminal residues including Lys-3, Lys-4, and Leu-5 are critical for controlling functional regulation, but not trafficking, of BK channels. This membrane-distal region has features of an amphipathic helix that is predicted to control the orientation of the first transmembrane-spanning domain (TM1) of the β1-subunit. In contrast, a membrane-proximal leucine residue (Leu-17) controls trafficking without affecting functional coupling, an effect that is in part dependent on controlling efficient endoplasmic reticulum exit of the pore-forming α-subunit. Thus cell surface trafficking and functional coupling with BK channels are controlled by distinct domains of the β1-subunit N terminus.

Polyglutamine Tract Expansion Increases S-Nitrosylation of Huntingtin and Ataxin-1

Expansion of the polyglutamine (polyQ) tract in the huntingtin (Htt) protein causes Huntington’s disease (HD), a fatal inherited movement disorder linked to neurodegeneration in the striatum and cortex. S-nitrosylation and S-acylation of cysteine residues regulate many functions of cytosolic proteins. We therefore used a resin-assisted capture approach to identify these modifications in Htt. In contrast to many proteins that have only a single S-nitrosylation or S-acylation site, we identified sites along much of the length of Htt. Moreover, analysis of cells expressing full-length Htt or a large N-terminal fragment of Htt shows that polyQ expansion strongly increases Htt S-nitrosylation. This effect appears to be general since it is also observed in Ataxin-1, which causes spinocerebellar ataxia type 1 (SCA1) when its polyQ tract is expanded. Overexpression of nitric oxide synthase increases the S-nitrosylation of normal Htt and the frequency of conspicuous juxtanuclear inclusions of Htt N-terminal fragments in transfected cells. Taken together with the evidence that S-nitrosylation of Htt is widespread and parallels polyQ expansion, these subcellular changes show that S-nitrosylation affects the biology of this protein in vivo.

The PP-motif in luminal loop 2 of ZnT transporters plays a pivotal role in TNAP activation

Secretory and membrane-bound zinc-requiring enzymes are thought to be activated by binding zinc in the early secretory pathway. One such enzyme, tissue-non-specific alkaline phosphatase (TNAP), is activated through a two-step mechanism, via protein stabilization and subsequent enzyme activation through metalation, by ZnT5-ZnT6 heterodimers or ZnT7 homodimers. However, little is known about the molecular basis underlying the activation process. In the present study, we found that the di-proline motif (PP-motif) in luminal loop 2 of ZnT5 and ZnT7 is important for TNAP activation. TNAP activity was significantly reduced in cells lacking ZnT5-ZnT6 heterodimers and ZnT7 homodimers [triple knockout (TKO) cells]. The decreased TNAP activity was restored by expressing hZnT5 with hZnT6 or hZnT7, but significantly less so (almost 90% less) by expressing mutants thereof in which the PP-motif was mutated to alanine (PP-AA). In TKO cells, overexpressed hTNAP was not completely activated, and it was converted less efficiently into the holo form by expressing a PP-AA mutant of hZnT5 with hZnT6, whose defects were not restored by zinc supplementation. The zinc transport activity of hZnT7 was not significantly impaired by the PP-AA mutation, indicating that the PP-motif is involved in the TNAP maturation process, although it does not control zinc transport activity. The PP-motif is highly conserved in ZnT5 and ZnT7 orthologues, and its importance for TNAP activation is conserved in the Caenorhabditis elegans hZnT5 orthologue CDF5. These results provide novel molecular insights into the TNAP activation process in the early secretory pathway.

An Entamoeba histolytica ADP-ribosyl transferase from the diphtheria toxin family modifies the bacterial elongation factor Tu

ADP-ribosyl transferases are enzymes involved in the post-translational modification of proteins; they participate in multiple physiological processes, pathogenesis and host-pathogen interactions. Several reports have characterized the functions of these enzymes in viruses, prokaryotes and higher eukaryotes, but few studies have reported ADP-ribosyl transferases in lower eukaryotes, such as parasites. The locus EHI_155600 from Entamoeba histolytica encodes a hypothetical protein that possesses a domain from the ADP-ribosylation superfamily; this protein belongs to the diphtheria toxin family according to a homology model using poly-ADP-ribosyl polymerase 12 (PARP12 or ARTD12) as a template. The recombinant protein expressed in Escherichia coli exhibited in vitro ADP-ribosylation activity that was dependent on the time and temperature. Unlabeled βNAD(+), but not ADP-ribose, competed in the enzymatic reaction using biotin-βNAD(+) as the ADP-ribose donor. The recombinant enzyme, denominated EhToxin-like, auto-ADP-ribosylated and modified an acceptor from E. coli that was identified by MS/MS as the elongation factor Tu (EF-Tu). To the best of our knowledge, this is the first report to identify an ADP-ribosyl transferase from the diphtheria toxin family in a protozoan parasite. The known toxins from this family (i.e., the diphtheria toxin, the Pseudomonas aeruginosa toxin Exo-A, and Cholix from Vibrio cholerae) modify eukaryotic elongation factor two (eEF-2), whereas the amoeba EhToxin-like modified EF-Tu, which is another elongation factor involved in protein synthesis in bacteria and mitochondria.

Identification of an Extracellular Gate for the Proton-coupled Folate Transporter (PCFT-SLC46A1) by Cysteine Cross-linking

The proton-coupled folate transporter (PCFT, SLC46A1) is required for intestinal folate absorption and folate homeostasis in humans. A homology model of PCFT, based upon theEscherichia coliglycerol 3-phosphate transporter structure, predicted that PCFT transmembrane domains (TMDs) 1, 2, 7, and 11 form an extracellular gate in the inward-open conformation. To assess this model, five residues (Gln(45)-TMD1, Asn(90)-TMD2, Leu(290)-TMD7, Ser(407)-TMD11 and Asn(411)-TMD11) in the predicted gate were substituted with Cys to generate single and nine double mutants. Transport function of the mutants was assayed in transient transfectants by measurement of [(3)H]substrate influx as was accessibility of the Cys residues to biotinylation. Pairs of Cys residues were assessed for spontaneous formation of a disulfide bond, induction of a disulfide bond by oxidization with dichloro(1,10-phenanthroline)copper (II) (CuPh), or the formation of a Cd(2+)complex. The data were consistent with the formation of a spontaneous disulfide bond between the N90C/S407C pair and a CuPh- and Cd(2+)-induced disulfide bond and complex, respectively, for the Q45C/L290C and L290C/N411C pairs. The decrease in activity induced by cross-linkage of the Cys residue pairs was due to a decrease in the influxVmaxconsistent with restriction in the mobility of the transporter. The presence of folate substrate decreased the CuPh-induced inhibition of transport. Hence, the data support the glycerol 3-phosphate transporter-based homology model of PCFT and the presence of an extracellular gate formed by TMDs 1, 2, 7, and 11.

Interaction of TAPBPR, a tapasin homolog, with MHC-I molecules promotes peptide editing

Peptide loading of major histocompatibility complex class I (MHC-I) molecules is central to antigen presentation, self-tolerance, and CD8(+) T-cell activation. TAP binding protein, related (TAPBPR), a widely expressed tapasin homolog, is not part of the classical MHC-I peptide-loading complex (PLC). Using recombinant MHC-I molecules, we show that TAPBPR binds HLA-A*02:01 and several other MHC-I molecules that are either peptide-free or loaded with low-affinity peptides. Fluorescence polarization experiments establish that TAPBPR augments peptide binding by MHC-I. The TAPBPR/MHC-I interaction is reversed by specific peptides, related to their affinity. Mutational and small-angle X-ray scattering (SAXS) studies confirm the structural similarities of TAPBPR with tapasin. These results support a role of TAPBPR in stabilizing peptide-receptive conformation(s) of MHC-I, permitting peptide editing.

Structural and sequencing analysis of local target DNA recognition by MLV integrase

Target-site selection by retroviral integrase (IN) proteins profoundly affects viral pathogenesis. We describe the solution nuclear magnetic resonance structure of the Moloney murine leukemia virus IN (M-MLV) C-terminal domain (CTD) and a structural homology model of the catalytic core domain (CCD). In solution, the isolated MLV IN CTD adopts an SH3 domain fold flanked by a C-terminal unstructured tail. We generated a concordant MLV IN CCD structural model using SWISS-MODEL, MMM-tree and I-TASSER. Using the X-ray crystal structure of the prototype foamy virus IN target capture complex together with our MLV domain structures, residues within the CCD α2 helical region and the CTD β1-β2 loop were predicted to bind target DNA. The role of these residues was analyzed in vivo through point mutants and motif interchanges. Viable viruses with substitutions at the IN CCD α2 helical region and the CTD β1-β2 loop were tested for effects on integration target site selection. Next-generation sequencing and analysis of integration target sequences indicate that the CCD α2 helical region, in particular P187, interacts with the sequences distal to the scissile bonds whereas the CTD β1-β2 loop binds to residues proximal to it. These findings validate our structural model and disclose IN-DNA interactions relevant to target site selection.

Identification of Tyr residues that enhance folate substrate binding and constrain oscillation of the proton-coupled folate transporter (PCFT-SLC46A1)

The proton-coupled folate transporter (PCFT) mediates intestinal folate absorption and transport of folates across the choroid plexus. This study focuses on the role of Tyr residues in PCFT function. The substituted Cys-accessibility method identified four Tyr residues (Y291, Y362, Y315, and Y414) that are accessible to the extracellular compartment; three of these (Y291, Y362, and Y315) are located within or near the folate binding pocket. When the Tyr residues were replaced with Cys or Ala, these mutants showed similar (up to 6-fold) increases in influx Vmax and Kt/Ki for [(3)H]methotrexate and [(3)H]pemetrexed. When the Tyr residues were replaced with Phe, these changes were moderated or absent. When Y315A PCFT was used as representative of the mutants and [(3)H]pemetrexed as the transport substrate, this substitution did not increase the efflux rate constant. Furthermore, neither influx nor efflux mediated by Y315A PCFT was transstimulated by the presence of substrate in the opposite compartment; however, substantial bidirectional transstimulation of transport was mediated by wild-type PCFT. This resulted in a threefold greater efflux rate constant for cells that express wild-type PCFT than for cells that express Y315 PCFT under exchange conditions. These data suggest that these Tyr residues, possibly through their rigid side chains, secure the carrier in a high-affinity state for its folate substrates. However, this may be achieved at the expense of constraining the carrier's mobility, thereby decreasing the rate at which the protein oscillates between its conformational states. The Vmax generated by these Tyr mutants may be so rapid that further augmentation during transstimulation may not be possible.

Computational methods in drug discovery

Computer-aided drug discovery/design methods have played a major role in the development of therapeutically important small molecules for over three decades. These methods are broadly classified as either structure-based or ligand-based methods. Structure-based methods are in principle analogous to high-throughput screening in that both target and ligand structure information is imperative. Structure-based approaches include ligand docking, pharmacophore, and ligand design methods. The article discusses theory behind the most important methods and recent successful applications. Ligand-based methods use only ligand information for predicting activity depending on its similarity/dissimilarity to previously known active ligands. We review widely used ligand-based methods such as ligand-based pharmacophores, molecular descriptors, and quantitative structure-activity relationships. In addition, important tools such as target/ligand data bases, homology modeling, ligand fingerprint methods, etc., necessary for successful implementation of various computer-aided drug discovery/design methods in a drug discovery campaign are discussed. Finally, computational methods for toxicity prediction and optimization for favorable physiologic properties are discussed with successful examples from literature.

Structural insights into a novel interkingdom signaling circuit by cartography of the ligand-binding sites of the homologous quorum sensing LuxR-family

Recent studies have identified a novel interkingdom signaling circuit, via plant signaling molecules, and a bacterial sub-family of LuxR proteins, bridging eukaryotes and prokaryotes. Indeed pivotal plant-bacteria interactions are regulated by the so called Plant Associated Bacteria (PAB) LuxR solo regulators that, although closely related to the quorum sensing (QS) LuxR family, do not bind or respond to canonical quorum sensing N-acyl homoserine lactones (AHLs), but only to specific host plant signal molecules. The large body of structural data available for several members of the QS LuxR family complexed with different classes of ligands (AHLs and other compounds), has been exploited to dissect the cartography of their regulatory domains through structure-based multiple sequence alignments, structural superimposition and a comparative analysis of the contact residues involved in ligand binding. In the absence of experimentally determined structures of members of the PAB LuxR solos subfamily, an homology model of its prototype OryR is presented, aiming to elucidate the architecture of its ligand-binding site. The obtained model, in combination with the cartography of the regulatory domains of the homologous QS LuxRs, provides novel insights into the 3D structure of its ligand-binding site and unveils the probable molecular determinants responsible for differences in selectivity towards specific host plant signal molecules, rather than to canonical QS compounds.

Role of the fourth transmembrane domain in proton-coupled folate transporter function as assessed by the substituted cysteine accessibility method

The proton-coupled folate transporter (PCFT, SLC46A1) mediates folate transport across the apical brush-border membrane of the proximal small intestine and the basolateral membrane of choroid plexus ependymal cells. Two loss-of-function mutations in PCFT, which are the basis for hereditary folate malabsorption, have been identified within the fourth transmembrane domain (TMD4) in subjects with this disorder. We have employed the substituted Cys accessibility method (SCAM) to study the accessibilities of all residues in TMD4 and their roles in folate substrate binding to the carrier. When residues 146-167 were replaced by Cys, all except R148C were expressed at the cell surface. Modification of five of these substituted Cys residues (positions 147, 152, 157, 158, and 161) by methanethiosulfonate (MTS) reagents led to reduction of PCFT function. All five residues could be labeled with N-biotinylaminoethyl-MTS, and this could be blocked by the high-affinity PCFT substrate pemetrexed. Pemetrexed also protected PCFT mutant function from inhibitory modification of the substituted Cys at positions 157, 158, and 161 by a MTS. The findings indicate that these five residues in TMD4 are accessible to the aqueous translocation pathway, play a role in folate substrate binding, and are likely located within or near the folate binding pocket. A homology model of PCFT places three of these residues, Phe¹⁵⁷, Gly¹⁵⁸, and Leu¹⁶¹, within a breakpoint in the midportion of TMD4, a region that likely participates in alterations in the PCFT conformational state during carrier cycling.

Structural model for tubulin recognition and deformation by kinesin-13 microtubule depolymerases

To elucidate the structural basis of the mechanism of microtubule depolymerization by kinesin-13s, we analyzed complexes of tubulin and the Drosophila melanogaster kinesin-13 KLP10A by electron microscopy (EM) and fluorescence polarization microscopy. We report a nanometer-resolution (1.1 nm) cryo-EM three-dimensional structure of the KLP10A head domain (KLP10AHD) bound to curved tubulin. We found that binding of KLP10AHD induces a distinct tubulin configuration with displacement (shear) between tubulin subunits in addition to curvature. In this configuration, the kinesin-binding site differs from that in straight tubulin, providing an explanation for the distinct interaction modes of kinesin-13s with the microtubule lattice or its ends. The KLP10AHD-tubulin interface comprises three areas of interaction, suggesting a crossbow-type tubulin-bending mechanism. These areas include the kinesin-13 family conserved KVD residues, and as predicted from the crossbow model, mutating these residues changes the orientation and mobility of KLP10AHDs interacting with the microtubule.

eThread: a highly optimized machine learning-based approach to meta-threading and the modeling of protein tertiary structures

Template-based modeling that employs various meta-threading techniques is currently the most accurate, and consequently the most commonly used, approach for protein structure prediction. Despite the evident progress in this field, accurate structure models cannot be constructed for a significant fraction of gene products, thus the development of new algorithms is required. Here, we describe the development, optimization and large-scale benchmarking of eThread, a highly accurate meta-threading procedure for the identification of structural templates and the construction of corresponding target-to-template alignments. eThread integrates ten state-of-the-art threading/fold recognition algorithms in a local environment and extensively uses various machine learning techniques to carry out fully automated template-based protein structure modeling. Tertiary structure prediction employs two protocols based on widely used modeling algorithms: Modeller and TASSER-Lite. As a part of eThread, we also developed eContact, which is a Bayesian classifier for the prediction of inter-residue contacts and eRank, which effectively ranks generated multiple protein models and provides reliable confidence estimates as structure quality assessment. Excluding closely related templates from the modeling process, eThread generates models, which are correct at the fold level, for >80% of the targets; 40–50% of the constructed models are of a very high quality, which would be considered accurate at the family level. Furthermore, in large-scale benchmarking, we compare the performance of eThread to several alternative methods commonly used in protein structure prediction. Finally, we estimate the upper bound for this type of approach and discuss the directions towards further improvements.

Functional roles of the A335 and G338 residues of the proton-coupled folate transporter (PCFT-SLC46A1) mutated in hereditary folate malabsorption

The proton-coupled folate transporter (PCFT-SLC46A1) mediates intestinal folate absorption and folate transport across the choroid plexus, processes defective in hereditary folate malabsorption (HFM). This paper characterizes the functional defect, and the roles of two mutated PCFT residues, associated with HFM (G338R and A335D). The A335D-PCFT and other mutations at this residue result in an unstable protein; when expression of a mutant protein was preserved, function was always retained. The G338R and other charged mutants resulted in an unstable protein; substitutions with small neutral and polar amino acids preserved protein but with impaired function. Pemetrexed and methotrexate (MTX) influx kinetics mediated by the G338C mutant PCFT revealed marked (15- to 20-fold) decreases in K(t) and V(max) compared with wild-type PCFT. In contrast, there was only a small (∼2-fold) decrease in the MTX influx K(i) and an increase (∼3-fold) in the pemetrexed influx K(i) for the G338C-PCFT mutant. Neither a decrease in pH to 4.5, nor an increase to 7.4, restored function of the G338C mutant relative to wild-type PCFT excluding a role for this residue in proton binding or proton coupling. Homology modeling localized the A335 and G338 residues embedded in the 9th transmembrane, consistent with the inaccessibility of the A335C and G338C proteins to MTS reagents. Hence, the loss of intrinsic G338C-PCFT function was due solely to impaired oscillation of the carrier between its conformational states. The data illustrate how alterations in carrier cycling can impact influx K(t) without comparable alterations in substrate binding to the carrier.

Identification of a functionally critical GXXG motif and its relationship to the folate binding site of the proton-coupled folate transporter (PCFT-SLC46A1)

The proton-coupled folate transporter (PCFT) mediates intestinal folate absorption, and loss-of-function mutations in this gene result in the autosomal recessive disorder hereditary folate malabsorption. The current study, focused on a structure-functional analysis of this transporter, identified Gly-189 and Gly-192 (a GxxG motif) located in the fifth transmembrane domain as residues that could not be replaced with alanine without a loss of function. In contrast, function was preserved when Gly-56 and Gly-59 (the other conservative GXXG motif in human PCFT) were replaced with alanine. Similarly, Gly-93 and Gly-97, which constitute the only conserved GXXXG dimerization motif in human PCFT, tolerated alanine substitution. To explore the role of this region in folate binding, the residues around Gly-189 and Gly-192 were analyzed by the substituted cysteine accessibility method. Both I188C and M193C mutants were functional and were inhibited by membrane-impermeable sulfhydryl-reactive reagents; this could be prevented with PCFT substrate, but the protection was sustained at 0°C only for the I188C mutant, consistent with localization of Ile-188 in the PCFT folate binding pocket. The functional role of residues around Gly-189 and Gly-192 is consistent with a molecular structural model in which these two residues along with Ieu-188 are accessible to the PCFT aqueous translocation pathway.

A P425R mutation of the proton-coupled folate transporter causing hereditary folate malabsorption produces a highly selective alteration in folate binding

Proton-coupled folate transporter (PCFT) mediates folate intestinal absorption and transport across the choroid plexus, processes defective in subjects with hereditary folate malabsorption (HFM). PCFT is also widely expressed in human solid tumors where it contributes to the transport of pemetrexed and other antifolates. This study defines the basis for the functional changes due to a P425R mutation detected in a subject with HFM. Among various substitutions, only positively charged mutants (P425R and P425K) lost function but in a highly selective manner. Transport of reduced folates mediated by P425R-PCFT was virtually abolished; the methotrexate influx K(t) was increased fivefold (from 2 to 10 μM). In contrast, the pemetrexed influx K(t) mediated by P425R-PCFT was decreased 30% compared with wild-type (WT)-PCFT. Methotrexate inhibition of pemetrexed influx was competitive with a K(i) for WT-PCFT comparable to its influx K(t). However, the methotrexate influx K(i) for P425R-PCFT was ∼15-fold higher than the WT-PCFT influx K(t) and threefold higher than the methotrexate influx K(t) for the P425R-PCFT mutant. The confirmed secondary structure and homology modeling place the P425 residue at the junction of the 6th external loop and 12th transmembrane domain, remote from the aqueous translocation pathway, a prediction confirmed by the failure to label P425C-PCFT with N-biotinylaminoethyl methanethiosulfonate-biotin and the absence of inhibition of P425C-PCFT function by water-soluble sulfhydryl reagents. Hence, despite its location, the P425R-PCFT mutation produces a conformational change that fully preserves pemetrexed binding but markedly impairs binding of methotrexate and other folates to the carrier.

Protein structure prediction using a docking-based hierarchical folding scheme

The pathways by which proteins fold into their specific native structure are still an unsolved mystery. Currently, many methods for protein structure prediction are available, and most of them tackle the problem by relying on the vast amounts of data collected from known protein structures. These methods are often not concerned with the route the protein follows to reach its final fold. This work is based on the premise that proteins fold in a hierarchical manner. We present FOBIA, an automated method for predicting a protein structure. FOBIA consists of two main stages: the first finds matches between parts of the target sequence and independently folding structural units using profile-profile comparison. The second assembles these units into a 3D structure by searching and ranking their possible orientations toward each other using a docking-based approach. We have previously reported an application of an initial version of this strategy to homology based targets. Since then we have considerably enhanced our method's abilities to allow it to address the more difficult template-based target category. This allows us to now apply FOBIA to the template-based targets of CASP8 and to show that it is both very efficient and promising. Our method can provide an alternative for template-based structure prediction, and in particular, the docking-basedranking technique presented here can be incorporated into any profile-profile comparison based method.

VISTA, a novel mouse Ig superfamily ligand that negatively regulates T cell responses

The immunoglobulin (Ig) superfamily consists of many critical immune regulators, including the B7 family ligands and receptors. In this study, we identify a novel and structurally distinct Ig superfamily inhibitory ligand, whose extracellular domain bears homology to the B7 family ligand PD-L1. This molecule is designated V-domain Ig suppressor of T cell activation (VISTA). VISTA is primarily expressed on hematopoietic cells, and VISTA expression is highly regulated on myeloid antigen-presenting cells (APCs) and T cells. A soluble VISTA-Ig fusion protein or VISTA expression on APCs inhibits T cell proliferation and cytokine production in vitro. A VISTA-specific monoclonal antibody interferes with VISTA-induced suppression of T cell responses by VISTA-expressing APCs in vitro. Furthermore, anti-VISTA treatment exacerbates the development of the T cell-mediated autoimmune disease experimental autoimmune encephalomyelitis in mice. Finally, VISTA overexpression on tumor cells interferes with protective antitumor immunity in vivo in mice. These findings show that VISTA, a novel immunoregulatory molecule, has functional activities that are nonredundant with other Ig superfamily members and may play a role in the development of autoimmunity and immune surveillance in cancer.

3D models of lamprey progesterone receptor complexed with progesterone, 7α-hydroxy-progesterone and 15α-hydroxy-progesterone

Sea lamprey, a basal vertebrate, contains a progesterone receptor [PR]. An unusual property of lamprey is that gonadotropin-releasing hormone induces synthesis of 15α-hydroxy-progesterone [15α-OH-P] instead of progesterone. There also is indirect evidence for 7α-OH-P in lamprey serum. To determine if there is a structural basis for the binding of 7α-OH-P and 15α-OH-P to lamprey PR, we constructed 3D models of the lamprey PR complexed with progesterone, 7α-OH-P and 15α-OH-P. These 3D models reveal that Met-277 in lamprey PR has a specific interaction with the 15α-hydroxyl on 15α-OH-P and with Met-192, which also contacts the 15α-hydroxyl group. We also find that 7α-OH-P has favorable contacts with side-chains in lamprey PR. BLAST searches reveal that Met-277 on lamprey PR is unique among vertebrate PRs. This unique site on lamprey PR could be a target for compounds to control reproduction in sea lamprey, an environmental pest in Lake Michigan.

A biochemical mechanism for the oncogenic potential of the p110beta catalytic subunit of phosphoinositide 3-kinase

Class I PI3-kinases signal downstream of receptor tyrosine kinases and G protein-coupled receptors and have been implicated in tumorigenesis. Although the oncogenic potential of the PI3-kinase subunit p110α requires its mutational activation, other p110 isoforms can induce transformation when overexpressed in the wild-type state. In wild-type p110α, N345 in the C2 domain forms hydrogen bonds with D560 and N564 in the inter-SH2 (iSH2) domain of p85, and mutations of p110α or p85 that disrupt this interface lead to increased basal activity and transformation. Sequence analysis reveals that N345 in p110α aligns with K342 in p110β. This difference makes wild-type p110β analogous to a previously described oncogenic mutant, p110α-N345K. We now show that p110β is inhibited by p85 to a lesser extent than p110α and is not differentially inhibited by wild-type p85 versus p85 mutants that disrupt the C2-iSH2 domain interface. Similar results were seen in soft agar and focus-formation assays, where p110β was similar to p110α-N345K in transforming potential. Inhibition of p110β by p85 was enhanced by a K342N mutation in p110β, which led to decreased activity in vitro, decreased basal Akt and ribosomal protein S6 kinase (S6K1) activation, and decreased transformation in NIH 3T3 cells. Moreover, unlike wild-type p110β, p110β-K342N was differentially regulated by wild-type and mutant p85, suggesting that the inhibitory C2-iSH2 interface is functional in this mutant. This study shows that the enhanced transforming potential of p110β is the result of its decreased inhibition by p85, due to the disruption of an inhibitory C2-iSH2 domain interface.

Structure-function activity of the human sodium-dependent multivitamin transporter: role of His¹¹⁵ and His²⁵⁴

Intestinal absorption of biotin occurs via a Na(+)-dependent carrier-mediated process that involves the sodium-dependent multivitamin transporter (SMVT; product of the Slc5a6 gene). The SMVT system is exclusively expressed at the apical membrane domain of the polarized intestinal epithelial cells. Whereas previous studies from our laboratory and others have characterized different physiological and biological aspects of SMVT, little is currently known about its structure-function activity relationship. Using site-directed mutagenesis approach, we examined the role of the positively charged histidine (His) residues of the human SMVT (hSMVT) in transporting the negatively charged biotin. Of the seven conserved (across species) His residues in the hSMVT polypeptide, only His¹¹⁵ and His²⁵⁴ were found to be important for the function of hSMVT as their mutation led to a significant reduction in carrier-mediated biotin uptake. This inhibition was mediated via a significant reduction in the maximal velocity (V(max)), but not the apparent Michaelis constant (K(m)), of the biotin uptake process and was not related to the charge of the His residue. The inhibition was also not due to changes in transcriptional or translational efficiency of the mutated hSMVT compared with wild-type carrier. However, surface biotinylation assay showed a significant reduction in the level of expression of the mutated hSMVT at the cell surface, a finding that was further confirmed by confocal imaging. Our results show important role for His¹¹⁵ and His²⁵⁴ residues in hSMVT function, which is most probably mediated via an effect on level of hSMVT expression at the cell membrane.

Insights into proteolytic processing of the major peanut allergen Ara h 2 by endogenous peanut proteases

BACKGROUND

The major peanut allergens are Ara h 1, Ara h 2 and Ara h 6. Proteolytic processing has been shown to be required for the maturation process of Ara h 6. The aim of this study was to examine whether Ara h 2 undergoes proteolytic processing and, if so, whether proteolytic processing influences its ability to bind human immunoglobulin E (IgE).

RESULTS

Ara h 2 isolated from peanut extract under conditions of protease inhibition revealed a single additional peak for its two known isoforms (Ara h 2.01 and Ara h 2.02), corresponding to a C-terminally truncated form lacking a dipeptide (RY). Ara h 2 isolated in the absence of protease inhibition, however, yielded two additional peaks, identified as C-terminally truncated forms lacking either a dipeptide (RY) or a single tyrosine residue. The IgE-binding capacity of the Ara h 2 truncated forms was not altered.

CONCLUSION

Ara h 2 undergoes proteolytic processing by peanut proteases that involves C-terminal removal of a dipeptide. Hence Ara h 2 isolated from peanut extract is a complex mixture of two isoforms expressed by different genes, Ara h 2.01 and Ara h 2.02, as well as truncated forms generated by the proteolytic processing of these isoforms.

Template-based protein structure modeling

Functional characterization of a protein is often facilitated by its 3D structure. However, the fraction of experimentally known 3D models is currently less than 1% due to the inherently time-consuming and complicated nature of structure determination techniques. Computational approaches are employed to bridge the gap between the number of known sequences and that of 3D models. Template-based protein structure modeling techniques rely on the study of principles that dictate the 3D structure of natural proteins from the theory of evolution viewpoint. Strategies for template-based structure modeling will be discussed with a focus on comparative modeling, by reviewing techniques available for all the major steps involved in the comparative modeling pipeline.

Identification of linear epitopes in Bacillus anthracis protective antigen bound by neutralizing antibodies

Protective antigen (PA), the binding subunit of anthrax toxin, is the major component in the current anthrax vaccine, but the fine antigenic structure of PA is not well defined. To identify linear neutralizing epitopes of PA, 145 overlapping peptides covering the entire sequence of the protein were synthesized. Six monoclonal antibodies (mAbs) and antisera from mice specific for PA were tested for their reactivity to the peptides by enzyme-linked immunosorbent assays. Three major linear immunodominant B-cell epitopes were mapped to residues Leu(156) to Ser(170), Val(196) to Ile(210), and Ser(312) to Asn(326) of the PA protein. Two mAbs with toxin-neutralizing activity recognized two different epitopes in close proximity to the furin cleavage site in domain 1. The three-dimensional complex structure of PA and its neutralizing mAbs 7.5G and 19D9 were modeled using the molecular docking method providing models for the interacting epitope and paratope residues. For both mAbs, LeTx neutralization was associated with interference with furin cleavage, but they differed in effectiveness depending on whether they bound on the N- or C-terminal aspect of the cleaved products. The two peptides containing these epitopes that include amino acids Leu(156)-Ser(170) and Val(196)-Ile(210) were immunogenic and elicited neutralizing antibody responses to PA. These results identify the first linear neutralizing epitopes of PA and show that peptides containing epitope sequences can elicit neutralizing antibody responses, a finding that could be exploited for vaccine design.

The functional roles of the His247 and His281 residues in folate and proton translocation mediated by the human proton-coupled folate transporter SLC46A1

This report addresses the functional role of His residues in the proton-coupled folate transporter (PCFT; SLC46A1), which mediates intestinal folate absorption. Of ten His residues, only H247A and H281A mutations altered function. The folic acid influx Kt at pH 5.5 for H247A was downward arrow 8.4-fold. Although wild type (WT)-PCFT Ki values varied among the folates, Ki values were much lower and comparable for H247-A, -R, -Q, or -E mutants. Homology modeling localized His247 to the large loop separating transmembrane domains 6 and 7 at the cytoplasmic entrance of the translocation pathway in hydrogen-bond distance to Ser172. The folic acid influx Kt for S172A-PCFT was decreased similar to H247A. His281 faces the extracellular region in the seventh transmembrane domain. H281A-PCFT results in loss-of-function due to approximately 12-fold upward arrow in the folic acid influx Kt. When the pH was decreased from 5.5 to 4.5, the WT-PCFT folic acid influx Kt was unchanged, but the Kt decreased 4-fold for H281A. In electrophysiological studies in Xenopus oocytes, both WT-PCFT- and H281A-PCFT-mediated folic acid uptake produced current and acidification, and both exhibited a low level of folate-independent proton transport (slippage). Slippage was markedly increased for the H247A-PCFT mutant. The data suggest that disruption of the His247 to Ser172 interaction results in a PCFT conformational alteration causing a loss of selectivity, increased substrate access to a high affinity binding pocket, and proton transport in the absence of a folate gradient. The His281 residue is not essential for proton coupling but plays an important role in PCFT protonation, which, in turn, augments folate binding to the carrier.

The sodium-dependent D-glucose transport protein of Helicobacter pylori

Helicobacter pylori is a gram-negative pathogenic microaerophile with a particular tropism for the mucosal surface of the gastric epithelium. Despite its obligatory microaerophilic character, it can metabolize D-glucose and/or D-galactose in both oxidative and fermentative pathways via a Na(+)-dependent secondary active transport, a glucokinase and enzymes of the pentose phosphate pathway. We have assigned the Na(+)-dependent transport of glucose to the protein product of the H. pylori 1174 gene. The gene was heterologously expressed in a glucose transport-deficient Escherichia coli strain, where transport activities of radiolabelled D-glucose, D-galactose and 2-deoxy-D-glucose were restored, consistent with the expected specificity of the hexose uptake system in H. pylori. D-mannose was also identified as a substrate. The HP1174 transport protein was purified and reconstituted into proteoliposomes, where sodium dependence of sugar transport activity was demonstrated. Additionally the tryptophan/tyrosine fluorescence of the purified protein showed quenching by 2-deoxy-D-glucose, D-mannose, D-glucose or D-galactose in the presence of sodium ions. This is the first reported purification and characterization of an active glucose transport protein member of the TC 2.1.7 subgroup of the Major Facilitator Superfamily, constituting the route for entry of sugar nutrients into H. pylori. A model is derived of its three-dimensional structure as a paradigm of the family.

Refining homology models by combining replica-exchange molecular dynamics and statistical potentials

A protocol is presented for the global refinement of homology models of proteins. It combines the advantages of temperature-based replica-exchange molecular dynamics (REMD) for conformational sampling and the use of statistical potentials for model selection. The protocol was tested using 21 models. Of these 14 were models of 10 small proteins for which high-resolution crystal structures were available, the remainder were targets of the recent CASPR exercise. It was found that REMD in combination with currently available force fields could sample near-native conformational states starting from high-quality homology models. Conformations in which the backbone RMSD of secondary structure elements (SSE-RMSD) was lower than the starting value by 0.5-1.0 A were found for 15 out of the 21 cases (average 0.82 A). Furthermore, when a simple scoring function consisting of two statistical potentials was used to rank the structures, one or more structures with SSE-RMSD of at least 0.2 A lower than the starting value was found among the five best ranked structures in 11 out of the 21 cases. The average improvement in SSE-RMSD for the best models was 0.42 A. However, none of the scoring functions tested identified the structures with the lowest SSE-RMSD as the best models although all identified the native conformation as the one with lowest energy. This suggests that while the proposed protocol proved effective for the refinement of high-quality models of small proteins scoring functions remain one of the major limiting factors in structure refinement. This and other aspects by which the methodology could be further improved are discussed.

Boosting Protein Threading Accuracy

Protein threading is one of the most successful protein structure prediction methods. Most protein threading methods use a scoring function linearly combining sequence and structure features to measure the quality of a sequence-template alignment so that a dynamic programming algorithm can be used to optimize the scoring function. However, a linear scoring function cannot fully exploit interdependency among features and thus, limits alignment accuracy.This paper presents a nonlinear scoring function for protein threading, which not only can model interactions among different protein features, but also can be efficiently optimized using a dynamic programming algorithm. We achieve this by modeling the threading problem using a probabilistic graphical model Conditional Random Fields (CRF) and training the model using the gradient tree boosting algorithm. The resultant model is a nonlinear scoring function consisting of a collection of regression trees. Each regression tree models a type of nonlinear relationship among sequence and structure features. Experimental results indicate that this new threading model can effectively leverage weak biological signals and improve both alignment accuracy and fold recognition rate greatly.

Mutations in the Arabidopsis peroxisomal ABC transporter COMATOSE allow differentiation between multiple functions in planta: insights from an allelic series

COMATOSE (CTS), the Arabidopsis homologue of human Adrenoleukodystrophy protein (ALDP), is required for import of substrates for peroxisomal beta-oxidation. A new allelic series and a homology model based on the bacterial ABC transporter, Sav1866, provide novel insights into structure-function relations of ABC subfamily D proteins. In contrast to ALDP, where the majority of mutations result in protein absence from the peroxisomal membrane, all CTS mutants produced stable protein. Mutation of conserved residues in the Walker A and B motifs in CTS nucleotide-binding domain (NBD) 1 resulted in a null phenotype but had little effect in NBD2, indicating that the NBDs are functionally distinct in vivo. Two alleles containing mutations in NBD1 outside the Walker motifs (E617K and C631Y) exhibited resistance to auxin precursors 2,4-dichlorophenoxybutyric acid (2,4-DB) and indole butyric acid (IBA) but were wild type in all other tests. The homology model predicted that the transmission interfaces are domain-swapped in CTS, and the differential effects of mutations in the conserved "EAA motif" of coupling helix 2 supported this prediction, consistent with distinct roles for each NBD. Our findings demonstrate that CTS functions can be separated by mutagenesis and the structural model provides a framework for interpretation of phenotypic data.

Improved scoring function for comparative modeling using the M4T method

Improvements in comparative protein structure modeling for the remote target-template sequence similarity cases are possible through the optimal combination of multiple template structures and by improving the quality of target-template alignment. Recently developed MMM and M4T methods were designed to address these problems. Here we describe new developments in both the alignment generation and the template selection parts of the modeling algorithms. We set up a new scoring function in MMM to deliver more accurate target-template alignments. This was achieved by developing and incorporating into the composite scoring function a novel statistical pairwise potential that combines local and non-local terms. The non-local term of the statistical potential utilizes a shuffled reference state definition that helped to eliminate most of the false positive signal from the background distribution of pairwise contacts. The accuracy of the scoring function was further increased by using BLOSUM mutation table scores.

How well can the accuracy of comparative protein structure models be predicted?

Comparative structure models are available for two orders of magnitude more protein sequences than are experimentally determined structures. These models, however, suffer from two limitations that experimentally determined structures do not: They frequently contain significant errors, and their accuracy cannot be readily assessed. We have addressed the latter limitation by developing a protocol optimized specifically for predicting the Calpha root-mean-squared deviation (RMSD) and native overlap (NO3.5A) errors of a model in the absence of its native structure. In contrast to most traditional assessment scores that merely predict one model is more accurate than others, this approach quantifies the error in an absolute sense, thus helping to determine whether or not the model is suitable for intended applications. The assessment relies on a model-specific scoring function constructed by a support vector machine. This regression optimizes the weights of up to nine features, including various sequence similarity measures and statistical potentials, extracted from a tailored training set of models unique to the model being assessed: If possible, we use similarly sized models with the same fold; otherwise, we use similarly sized models with the same secondary structure composition. This protocol predicts the RMSD and NO3.5A errors for a diverse set of 580,317 comparative models of 6174 sequences with correlation coefficients (r) of 0.84 and 0.86, respectively, to the actual errors. This scoring function achieves the best correlation compared to 13 other tested assessment criteria that achieved correlations ranging from 0.35 to 0.71.

Benchmarking of TASSER_2.0: an improved protein structure prediction algorithm with more accurate predicted contact restraints

To improve tertiary structure predictions of more difficult targets, the next generation of TASSER, TASSER_2.0, has been developed. TASSER_2.0 incorporates more accurate side-chain contact restraint predictions from a new approach, the composite-sequence method, based on consensus restraints generated by an improved threading algorithm, PROSPECTOR_3.5, which uses computationally evolved and wild-type template sequences as input. TASSER_2.0 was tested on a large-scale, benchmark set of 2591 nonhomologous, single domain proteins < or =200 residues that cover the Protein Data Bank at 35% pairwise sequence identity. Compared with the average fraction of accurately predicted side-chain contacts of 0.37 using PROSPECTOR_3.5 with wild-type template sequences, the average accuracy of the composite-sequence method increases to 0.60. The resulting TASSER_2.0 models are closer to their native structures, with an average root mean-square deviation of 4.99 A compared to the 5.31 A result of TASSER. Defining a successful prediction as a model with a root mean-square deviation to native <6.5 A, the success rate of TASSER_2.0 (TASSER) for Medium targets (targets with good templates/poor alignments) is 74.3% (64.7%) and 40.8% (35.5%) for the Hard targets (incorrect templates/alignments). For Easy targets (good templates/alignments), the success rate slightly increases from 86.3% to 88.4%.

Development and benchmarking of TASSER(iter) for the iterative improvement of protein structure predictions

To improve the accuracy of TASSER models especially in the limit where threading provided template alignments are of poor quality, we have developed the TASSER(iter) algorithm which uses the templates and contact restraints from TASSER generated models for iterative structure refinement. We apply TASSER(iter) to a large benchmark set of 2,773 nonhomologous single domain proteins that are < or = 200 in length and that cover the PDB at the level of 35% pairwise sequence identity. Overall, TASSER(iter) models have a smaller global average RMSD of 5.48 A compared to 5.81 A RMSD of the original TASSER models. Classifying the targets by the level of prediction difficulty (where Easy targets have a good template with a corresponding good threading alignment, Medium targets have a good template but a poor alignment, and Hard targets have an incorrectly identified template), TASSER(iter) (TASSER) models have an average RMSD of 4.15 A (4.35 A) for the Easy set and 9.05 A (9.52 A) for the Hard set. The largest reduction of average RMSD is for the Medium set where the TASSER(iter) models have an average global RMSD of 5.67 A compared to 6.72 A of the TASSER models. Seventy percent of the Medium set TASSER(iter) models have a smaller RMSD than the TASSER models, while 63% of the Easy and 60% of the Hard TASSER models are improved by TASSER(iter). For the foldable cases, where the targets have a RMSD to the native <6.5 A, TASSER(iter) shows obvious improvement over TASSER models: For the Medium set, it improves the success rate from 57.0 to 67.2%, followed by the Hard targets where the success rate improves from 32.0 to 34.8%, with the smallest improvement in the Easy targets from 82.6 to 84.0%. These results suggest that TASSER(iter) can provide more reliable predictions for targets of Medium difficulty, a range that had resisted improvement in the quality of protein structure predictions.

Exchanging murine and human immunoglobulin constant chains affects the kinetics and thermodynamics of antigen binding and chimeric antibody autoreactivity

Mouse-human chimeric antibodies composed of murine variable (V) and human (C) chains are useful therapeutic reagents. Consequently, we investigated whether heterologous C-regions from mice and humans affected specificity and affinity, and determined the contribution of C(H) glycosylation to antigen binding. The interaction of a 12-mer peptide mimetic with monoclonal antibody (mAb) 18B7 to Cryptococcus neoformans glucuronoxylomannan, and its chimeric (ch) and deglycosylated forms were studied by surface plasmon resonance. The equilibrium and rate association constants for the chAb were higher than for mAb 18B7. V region affinity was not affected by C(H) region glycosylation whereas heterologous C region of the same isotype altered the Ab binding affinity and the specificity for self-antigens. Structural models displayed local differences that implied changes on the connectivity of residues. These findings suggest that V region conformational changes can be dictated by the C(H) domains through an allosteric effect involving networks of highly connected amino acids.

Predicting and improving the protein sequence alignment quality by support vector regression

BACKGROUND

For successful protein structure prediction by comparative modeling, in addition to identifying a good template protein with known structure, obtaining an accurate sequence alignment between a query protein and a template protein is critical. It has been known that the alignment accuracy can vary significantly depending on our choice of various alignment parameters such as gap opening penalty and gap extension penalty. Because the accuracy of sequence alignment is typically measured by comparing it with its corresponding structure alignment, there is no good way of evaluating alignment accuracy without knowing the structure of a query protein, which is obviously not available at the time of structure prediction. Moreover, there is no universal alignment parameter option that would always yield the optimal alignment.

RESULTS

In this work, we develop a method to predict the quality of the alignment between a query and a template. We train the support vector regression (SVR) models to predict the MaxSub scores as a measure of alignment quality. The alignment between a query protein and a template of length n is transformed into a (n + 1)-dimensional feature vector, then it is used as an input to predict the alignment quality by the trained SVR model. Performance of our work is evaluated by various measures including Pearson correlation coefficient between the observed and predicted MaxSub scores. Result shows high correlation coefficient of 0.945. For a pair of query and template, 48 alignments are generated by changing alignment options. Trained SVR models are then applied to predict the MaxSub scores of those and to select the best alignment option which is chosen specifically to the query-template pair. This adaptive selection procedure results in 7.4% improvement of MaxSub scores, compared to those when the single best parameter option is used for all query-template pairs.

CONCLUSION

The present work demonstrates that the alignment quality can be predicted with reasonable accuracy. Our method is useful not only for selecting the optimal alignment parameters for a chosen template based on predicted alignment quality, but also for filtering out problematic templates that are not suitable for structure prediction due to poor alignment accuracy. This is implemented as a part in FORECAST, the server for fold-recognition and is freely available on the web at http://pbil.kaist.ac.kr/forecast.