MMM: a sequence-to-structure alignment protocol

Fluorescent proteins generate a genetic color polymorphism and counteract oxidative stress in intertidal sea anemones

Fluorescent proteins (FPs) are ubiquitous tools in research, yet their endogenous functions in nature are poorly understood. In this work, we describe a combination of functions for FPs in a clade of intertidal sea anemones whose FPs control a genetic color polymorphism together with the ability to combat oxidative stress. Focusing on the underlying genetics of a fluorescent green "Neon" color morph, we show that allelic differences in a single FP gene generate its strong and vibrant color, by increasing both molecular brightness and FP gene expression level. Natural variation in FP sequences also produces differences in antioxidant capacity. We demonstrate that these FPs are strong antioxidants that can protect live cells against oxidative stress. Finally, based on structural modeling of the responsible amino acids, we propose a model for FP antioxidant function that is driven by molecular surface charge. Together, our findings shed light on the multifaceted functions that can co-occur within a single FP and provide a framework for studying the evolution of fluorescence as it balances spectral and physiological functions in nature.

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)).

Human DHEA sulfation requires direct interaction between PAPS synthase 2 and DHEA sulfotransferase SULT2A1

The high-energy sulfate donor 3′-phosphoadenosine-5′-phosphosulfate (PAPS), generated by human PAPS synthase isoforms PAPSS1 and PAPSS2, is required for all human sulfation pathways. Sulfotransferase SULT2A1 uses PAPS for sulfation of the androgen precursor dehydroepiandrosterone (DHEA), thereby reducing downstream activation of DHEA to active androgens. Human PAPSS2 mutations manifest with undetectable DHEA sulfate, androgen excess, and metabolic disease, suggesting that ubiquitous PAPSS1 cannot compensate for deficient PAPSS2 in supporting DHEA sulfation. In knockdown studies in human adrenocortical NCI-H295R1 cells, we found that PAPSS2, but not PAPSS1, is required for efficient DHEA sulfation. Specific APS kinase activity, the rate-limiting step in PAPS biosynthesis, did not differ between PAPSS1 and PAPSS2. Co-expression of cytoplasmic SULT2A1 with a cytoplasmic PAPSS2 variant supported DHEA sulfation more efficiently than co-expression with nuclear PAPSS2 or nuclear/cytosolic PAPSS1. Proximity ligation assays revealed protein–protein interactions between SULT2A1 and PAPSS2 and, to a lesser extent, PAPSS1. Molecular docking studies showed a putative binding site for SULT2A1 within the PAPSS2 APS kinase domain. Energy-dependent scoring of docking solutions identified the interaction as specific for the PAPSS2 and SULT2A1 isoforms. These findings elucidate the mechanistic basis for the selective requirement for PAPSS2 in human DHEA sulfation.

Structures of the Gasdermin D C-Terminal Domains Reveal Mechanisms of Autoinhibition

Pyroptosis is an inflammatory form of programmed cell death that plays important roles in immune protection against infections and in inflammatory disorders. Gasdermin D (GSDMD) is an executor of pyroptosis upon cleavage by caspases-1/4/5/11 following canonical and noncanonical inflammasome activation. GSDMD N-terminal domain assembles membrane pores to induce cytolysis, whereas its C-terminal domain inhibits cell death through intramolecular association with the N domain. The molecular mechanisms of autoinhibition for GSDMD are poorly characterized. Here we report the crystal structures of the human and murine GSDMD C-terminal domains, which differ from those of the full-length murine GSDMA3 and the human GSDMB C-terminal domain. Mutations of GSDMD C-domain residues predicted to locate at its interface with the N-domain enhanced pyroptosis. Our results suggest that GSDMDs may employ a distinct mode of intramolecular domain interaction and autoinhibition, which may be relevant to its unique role in pyroptosis downstream of inflammasome activation.

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.

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.

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.

The molecular evolution of cytochrome P450 genes within and between drosophila species

We map 114 gene gains and 74 gene losses in the P450 gene family across the phylogeny of 12 Drosophila species by examining the congruence of gene trees and species trees. Although the number of P450 genes varies from 74 to 94 in the species examined, we infer that there were at least 77 P450 genes in the ancestral Drosophila genome. One of the most striking observations in the data set is the elevated loss of P450 genes in the Drosophila sechellia lineage. The gain and loss events are not evenly distributed among the P450 genes-with 30 genes showing no gene gains or losses whereas others show as many as 20 copy number changes among the species examined. The P450 gene clades showing the fewest number of gene gain and loss events tend to be those evolving with the most purifying selection acting on the protein sequences, although there are exceptions, such as the rapid rate of amino acid replacement observed in the single copy phantom (Cyp306a1) gene. Within D. melanogaster, we observe gene copy number polymorphism in ten P450 genes including multiple cases of interparalog chimeras. Nonallelic homologous recombination (NAHR) has been associated with deleterious mutations in humans, but here we provide a second possible example of an NAHR event in insect P450s being adaptive. Specifically, we find that a polymorphic Cyp12a4/Cyp12a5 chimera correlates with resistance to an insecticide. Although we observe such interparalog exchange in our within-species data sets, we have little evidence of it between species, raising the possibility that such events may occur more frequently than appreciated but are masked by subsequent sequence change.

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.

Protein structure based prediction of catalytic residues

BACKGROUND

Worldwide structural genomics projects continue to release new protein structures at an unprecedented pace, so far nearly 6000, but only about 60% of these proteins have any sort of functional annotation.

RESULTS

We explored a range of features that can be used for the prediction of functional residues given a known three-dimensional structure. These features include various centrality measures of nodes in graphs of interacting residues: closeness, betweenness and page-rank centrality. We also analyzed the distance of functional amino acids to the general center of mass (GCM) of the structure, relative solvent accessibility (RSA), and the use of relative entropy as a measure of sequence conservation. From the selected features, neural networks were trained to identify catalytic residues. We found that using distance to the GCM together with amino acid type provide a good discriminant function, when combined independently with sequence conservation. Using an independent test set of 29 annotated protein structures, the method returned 411 of the initial 9262 residues as the most likely to be involved in function. The output 411 residues contain 70 of the annotated 111 catalytic residues. This represents an approximately 14-fold enrichment of catalytic residues on the entire input set (corresponding to a sensitivity of 63% and a precision of 17%), a performance competitive with that of other state-of-the-art methods.

CONCLUSIONS

We found that several of the graph based measures utilize the same underlying feature of protein structures, which can be simply and more effectively captured with the distance to GCM definition. This also has the added the advantage of simplicity and easy implementation. Meanwhile sequence conservation remains by far the most influential feature in identifying functional residues. We also found that due the rapid changes in size and composition of sequence databases, conservation calculations must be recalibrated for specific reference databases.

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.

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.

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.

PCRPi: Presaging Critical Residues in Protein interfaces, a new computational tool to chart hot spots in protein interfaces

Protein-protein interactions (PPIs) are ubiquitous in Biology, and thus offer an enormous potential for the discovery of novel therapeutics. Although protein interfaces are large and lack defining physiochemical traits, is well established that only a small portion of interface residues, the so-called hot spot residues, contribute the most to the binding energy of the protein complex. Moreover, recent successes in development of novel drugs aimed at disrupting PPIs rely on targeting such residues. Experimental methods for describing critical residues are lengthy and costly; therefore, there is a need for computational tools that can complement experimental efforts. Here, we describe a new computational approach to predict hot spot residues in protein interfaces. The method, called Presaging Critical Residues in Protein interfaces (PCRPi), depends on the integration of diverse metrics into a unique probabilistic measure by using Bayesian Networks. We have benchmarked our method using a large set of experimentally verified hot spot residues and on a blind prediction on the protein complex formed by HRAS protein and a single domain antibody. Under both scenarios, PCRPi delivered consistent and accurate predictions. Finally, PCRPi is able to handle cases where some of the input data is either missing or not reliable (e.g. evolutionary information).

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.

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.

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.

Subdividing repressor function: DNA binding affinity, selectivity, and allostery can be altered by amino acid substitution of nonconserved residues in a LacI/GalR homologue

Many mutations that impact protein function occur at residues that do not directly contact ligand. To understand the functional contributions from the sequence that links the DNA-binding and regulatory domains of the LacI/GalR homologues, we have created a chimeric protein (LLhP), which comprises the LacI DNA-binding domain, the LacI linker, and the PurR regulatory domain. Although DNA binding site residues are identical in LLhP and LacI, thermodynamic measurements of DNA binding affinity show that LLhP does not discriminate between alternative DNA ligands as well as LacI. In addition, small-angle scattering experiments show that LLhP is more compact than LacI. When DNA is released, LacI shows a 20 A increase in length that was previously attributed to unfolding of the linker. This change is not seen in apo-LLhP, even though the linker sequences of the two proteins are identical. Together, results indicate that long-range functional and structural changes are propagated across the interface that forms between the linker and regulatory domain. These changes could be mediated via the side chains of several linker residues that contact the regulatory domains of the naturally occurring proteins, LacI and PurR. Substitution of these residues in LLhP leads to a range of functional effects. Four variants exhibit altered affinity for DNA, with no changes in selectivity or allosteric response. Another two result in proteins that bind operator DNA with very low affinity and no allosteric response, similar to LacI binding nonspecific DNA sequences. Two more substitutions simultaneously diminish affinity, enhance allostery, and profoundly alter DNA ligand selectivity. Thus, positions within the linker can be varied to modulate different aspects of repressor function.

Predicting disulfide bond connectivity in proteins by correlated mutations analysis

MOTIVATION

Prediction of disulfide bond connectivity facilitates structural and functional annotation of proteins. Previous studies suggest that cysteines of a disulfide bond mutate in a correlated manner.

RESULTS

We developed a method that analyzes correlated mutation patterns in multiple sequence alignments in order to predict disulfide bond connectivity. Proteins with known experimental structures and varying numbers of disulfide bonds, and that spanned various evolutionary distances, were aligned. We observed frequent variation of disulfide bond connectivity within members of the same protein families, and it was also observed that in 99% of the cases, cysteine pairs forming non-conserved disulfide bonds mutated in concert. Our data support the notion that substitution of a cysteine in a disulfide bond prompts the substitution of its cysteine partner and that oxidized cysteines appear in pairs. The method we developed predicts disulfide bond connectivity patterns with accuracies of 73, 69 and 61% for proteins with two, three and four disulfide bonds, respectively.

Comparative protein structure modeling by combining multiple templates and optimizing sequence-to-structure alignments

MOTIVATION

Two major bottlenecks in advancing comparative protein structure modeling are the efficient combination of multiple template structures and the generation of a correct input target-template alignment.

RESULTS

A novel method, Multiple Mapping Method with Multiple Templates (M4T) is introduced that implements an algorithm to automatically select and combine Multiple Template structures (MT) and an alignment optimization protocol (Multiple Mapping Method, MMM). The MT module of M4T selects and combines multiple template structures through an iterative clustering approach that takes into account the 'unique' contribution of each template, their sequence similarity among themselves and to the target sequence, and their experimental resolution. MMM is a sequence-to-structure alignment method that optimally combines alternatively aligned regions according to their fit in the structural environment of the template structure. The resulting M4T alignment is used as input to a comparative modeling module. The performance of M4T has been benchmarked on CASP6 comparative modeling target sequences and on a larger independent test set, and showed favorable performance to current state of the art methods.

New structure model for the ATP-binding cassette multidrug transporter LmrA

Multidrug resistance of pathogenic microorganisms and mammalian tumors can be associated with the overexpression of multidrug transporters. These integral membrane proteins are capable of extruding a wide range of structurally unrelated compounds from the cell. Among the different classes of multidrug transporters are the ATP binding cassette (ABC) transporters, which are dependent on the binding and hydrolysis of ATP. In the past five years, many researchers have built homology models of ABC extrusion systems using the atomic coordinates of crystallized MsbA, a lipopolysaccharide transporter in Gram-negative bacteria. Likewise, we have previously used the Vibrio cholera MsbA structure as a template in the modeling of the multidrug transporter LmrA from Lactococcus lactis. In view of the recently discovered inaccuracies in the MsbA structure, we have remodelled LmrA using the atomic coordinates of the MsbA homologue Sav1866 from Staphylococcus aureus. To compare and test our MsbA-based and Sav1866-based LmrA models we performed cysteine cross-linking at three key positions in LmrA. The pattern of cross-linking at these positions was consistent with the overall topology of transmembrane helices in Sav1866, suggesting that its crystal structure might be physiologically relevant. We recently identified E314 as a residue important in proton conduction by LmrA. The predicted location of this residue at the interface between the two half-transporters in the Sav1866-based homodimer, within the inner leaflet of the phospholipid bilayer, provides a new structural basis for the role of E314 in LmrA-mediated transport.

M4T: a comparative protein structure modeling server

Multiple Mapping Method with Multiple Templates (M4T) (http://www.fiserlab.org/servers/m4t) is a fully automated comparative protein structure modeling server. The novelty of M4T resides in two of its major modules, Multiple Templates (MT) and Multiple Mapping Method (MMM). The MT module of M4T selects and optimally combines the sequences of multiple template structures through an iterative clustering approach that takes into account the ‘unique’ contribution of each template, its sequence similarity to other template sequences and to the target sequences, and the quality of its experimental resolution. MMM module is a sequence-to-structure alignment method that is aimed at improving the alignment accuracy, especially at lower sequence identity levels. The current implementation of MMM takes inputs from three profile-to-profile-based alignment methods and iteratively compares and ranks alternatively aligned regions according to their fit in the structural environment of the template structure. The performance of M4T was benchmarked on CASP6 comparative modeling target sequences and on a larger independent test set and showed a favorable performance to current state-of-the-art methods.

The immunoglobulin heavy chain constant region affects kinetic and thermodynamic parameters of antibody variable region interactions with antigen

A central dogma in immunology is that antibody specificity is a function of the variable (V) region. However serological analysis of IgG(1), IgG(2a), and IgG(2b) switch variants of murine monoclonal antibody (mAb) 3E5 IgG(3) with identical V domains revealed apparent specificity differences for Cryptococcus neoformans glucuronoxylomannan (GXM). Kinetic and thermodynamic binding properties of mAbs 3E5 to a 12-mer peptide mimetic of GXM revealed differences in the affinity of these mAbs for a monovalent ligand, a result that implied that the constant (C) region affects the secondary structure of the antigen binding site, thus accounting for variations in specificity. Structural models of mAbs 3E5 suggested that isotype-related differences in binding resulted from amino acid sequence polymorphisms in the C region. This study implies that isotype switching is another mechanism for generating diversity in antigen binding and that isotype restriction of certain antibody responses may reflect structural constraints imposed by C region on V region binding. Furthermore, isotype affected the polyreactivity of V region identical antibodies, implying a role for C region in determining self-reactivity.