Search

Developing therapeutic antibodies is laborious and costly. Here we report a method for antibody discovery that leverages the Illumina HiSeq platform to, within 3 days, screen in the order of 10 8 antibody-antigen interactions. The method, which we named 'deep screening', involves the clustering and sequencing of antibody libraries, the conversion of the DNA clusters into complementary RNA clusters covalently linked to the instrument's flow-cell surface on the same location, the in situ translation of the clusters into antibodies tethered via ribosome display, and their screening via fluorescently labelled antigens. By using deep screening, we discovered low-nanomolar nanobodies to a model antigen using 4 × 10 6 unique variants from yeast-display-enriched libraries, and high-picomolar single-chain antibody fragment leads for human interleukin-7 directly from unselected synthetic repertoires. We also leveraged deep screening of a library of 2.4 × 10 5 sequences of the third complementarity-determining region of the heavy chain of an anti-human epidermal growth factor receptor 2 (HER2) antibody as input for a large language model that generated new single-chain antibody fragment sequences with higher affinity for HER2 than those in the original library., (© 2023. The Author(s).)

Optimum protein function and biochemical activity critically depends on water availability because solvent thermodynamics drive protein folding and macromolecular interactions 1 . Reciprocally, macromolecules restrict the movement of 'structured' water molecules within their hydration layers, reducing the available 'free' bulk solvent and therefore the total thermodynamic potential energy of water, or water potential. Here, within concentrated macromolecular solutions such as the cytosol, we found that modest changes in temperature greatly affect the water potential, and are counteracted by opposing changes in osmotic strength. This duality of temperature and osmotic strength enables simple manipulations of solvent thermodynamics to prevent cell death after extreme cold or heat shock. Physiologically, cells must sustain their activity against fluctuating temperature, pressure and osmotic strength, which impact water availability within seconds. Yet, established mechanisms of water homeostasis act over much slower timescales 2,3 ; we therefore postulated the existence of a rapid compensatory response. We find that this function is performed by water potential-driven changes in macromolecular assembly, particularly biomolecular condensation of intrinsically disordered proteins. The formation and dissolution of biomolecular condensates liberates and captures free water, respectively, quickly counteracting thermal or osmotic perturbations of water potential, which is consequently robustly buffered in the cytoplasm. Our results indicate that biomolecular condensation constitutes an intrinsic biophysical feedback response that rapidly compensates for intracellular osmotic and thermal fluctuations. We suggest that preserving water availability within the concentrated cytosol is an overlooked evolutionary driver of protein (dis)order and function., (© 2023. The Author(s).)

Poor understanding of intracellular delivery and targeting hinders development of nucleic acid-based therapeutics transported by nanoparticles. Utilizing a siRNA-targeting and small molecule profiling approach with advanced imaging and machine learning biological insights is generated into the mechanism of lipid nanoparticle (MC3-LNP) delivery of mRNA. This workflow is termed Advanced Cellular and Endocytic profiling for Intracellular Delivery (ACE-ID). A cell-based imaging assay and perturbation of 178 targets relevant to intracellular trafficking is used to identify corresponding effects on functional mRNA delivery. Targets improving delivery are analyzed by extracting data-rich phenotypic fingerprints from images using advanced image analysis algorithms. Machine learning is used to determine key features correlating with enhanced delivery, identifying fluid-phase endocytosis as a productive cellular entry route. With this new knowledge, MC3-LNP is re-engineered to target macropinocytosis, and this significantly improves mRNA delivery in vitro and in vivo. The ACE-ID approach can be broadly applicable for optimizing nanomedicine-based intracellular delivery systems and has the potential to accelerate the development of delivery systems for nucleic acid-based therapeutics., (© 2023 The Authors. Small Methods published by Wiley-VCH GmbH.)

The spontaneous emergence of function from diverse RNA sequence pools is widely considered an important transition in the origin of life. Here we show that diverse sequence pools are not a prerequisite for the emergence of function. Starting five independent selection experiments each from a single RNA seed sequence - comprising a central homopolymeric poly-A (or poly-U) segment flanked by different conserved primer binding sites - we observe transformation (continuous drift) of the seeds into low diversity sequence pools by mutation, truncation and recombination without ever reaching that of a random pool even after 24 rounds. Upon continuous error prone replication and selection for ATP binding we isolate specific ATP- or GTP-binding aptamers with low micromolar affinities. Our results have implications for early RNA evolution in the light of the high mutation rates associated with both non-enzymatic and enzymatic prebiotic RNA replication., Competing Interests: Conflict of Interests The authors declare no competing interest.

Angiotensinogen fine-tunes the tightly controlled activity of the renin-angiotensin system by modulating the release of angiotensin peptides that control blood pressure. One mechanism by which this modulation is achieved is via angiotensinogen's Cys18-Cys138 disulfide bond that acts as a redox switch. Molecular dynamics simulations of each redox state of angiotensinogen reveal subtle dynamic differences between the reduced and oxidised forms, particularly at the N-terminus. Surface plasmon resonance data demonstrate that the two redox forms of angiotensinogen display different binding kinetics to an immobilised anti-angiotensinogen monoclonal antibody. Mass spectrometry mapped the epitope for the antibody to the N-terminal region of angiotensinogen. We therefore provide evidence that the different redox forms of angiotensinogen can be detected by an antibody-based detection method., (© 2021 The Author(s). Published by Portland Press Limited on behalf of the Biochemical Society.)

The fibronectin type III (FN3) monobody domain is a promising non-antibody scaffold, which features a less complex architecture than an antibody while maintaining analogous binding loops. We previously developed FN3Con, a hyperstable monobody derivative with diagnostic and therapeutic potential. Prestabilization of the scaffold mitigates the stability-function trade-off commonly associated with evolving a protein domain toward biological activity. Here, we aimed to examine if the FN3Con monobody could take on antibody-like binding to therapeutic targets, while retaining its extreme stability. We targeted the first of the Adnectin derivative of monobodies to reach clinical trials, which was engineered by directed evolution for binding to the therapeutic target VEGFR2; however, this function was gained at the expense of large losses in thermostability and increased oligomerization. In order to mitigate these losses, we grafted the binding loops from Adnectin-anti-VEGFR2 (CT-322) onto the prestabilized FN3Con scaffold to produce a domain that successfully bound with high affinity to the therapeutic target VEGFR2. This FN3Con-anti-VEGFR2 construct also maintains high thermostability, including remarkable long-term stability, retaining binding activity after 2 years of storage at 36 °C. Further investigations into buffer excipients doubled the presence of monomeric monobody in accelerated stability trials. These data suggest that loop grafting onto a prestabilized scaffold is a viable strategy for the development of monobody domains with desirable biophysical characteristics and that FN3Con is therefore well-suited to applications such as the evolution of multiple paratopes or shelf-stable diagnostics and therapeutics., Competing Interests: Conflict of interest The authors declare no competing financial interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Niemann-Pick type C disease is a lysosomal storage disease affecting primarily the nervous system that results in premature death. Here we present the first report and investigation of Niemann-Pick type C disease in Australian Angus/Angus-cross calves. After a preliminary diagnosis of Niemann-Pick type C, samples from two affected calves and two obligate carriers were analysed using single nucleotide polymorphism genotyping and homozygosity mapping, and NPC1 was considered as a positional candidate gene. A likely causal missense variant on chromosome 24 in the NPC1 gene (NM_174758.2:c.2969C>G) was identified by Sanger sequencing of cDNA. SIFT analysis, protein alignment and protein modelling predicted the variant to be deleterious to protein function. Segregation of the variant with disease was confirmed in two additional affected calves and two obligate carrier dams. Genotyping of 403 animals from the original herd identified an estimated allele frequency of 3.5%. The Niemann-Pick type C phenotype was additionally confirmed via biochemical analysis of Lysotracker Green, cholesterol, sphingosine and glycosphingolipids in fibroblast cell cultures originating from two affected calves. The identification of a novel missense variant for Niemann-Pick type C disease in Angus/Angus-cross cattle will enable improved breeding and management of this disease in at-risk populations. The results from this study offer a unique opportunity to further the knowledge of human Niemann-Pick type C disease through the potential availability of a bovine model of disease., Competing Interests: The authors have declared that no competing interests exist.

Conformational diversity and self-cross-reactivity of antigens have been correlated with evasion from neutralizing antibody responses. We utilized single cell B cell sequencing, biolayer interferometry and X-ray crystallography to trace mutation selection pathways where the antibody response must resolve cross-reactivity between foreign and self-proteins bearing near-identical contact surfaces, but differing in conformational flexibility. Recurring antibody mutation trajectories mediate long-range rearrangements of framework (FW) and complementarity determining regions (CDRs) that increase binding site conformational diversity. These antibody mutations decrease affinity for self-antigen 19-fold and increase foreign affinity 67-fold, to yield a more than 1,250-fold increase in binding discrimination. These results demonstrate how conformational diversity in antigen and antibody does not act as a barrier, as previously suggested, but rather facilitates high affinity and high discrimination between foreign and self., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

The ability of reverse transcriptases (RTs) to synthesize a complementary DNA from natural RNA and a range of unnatural xeno nucleic acid (XNA) template chemistries, underpins key methods in molecular and synthetic genetics. However, RTs have proven challenging to discover and engineer, in particular for the more divergent XNA chemistries. Here we describe a general strategy for the directed evolution of RT function for any template chemistry called compartmentalized bead labelling and demonstrate it by the directed evolution of efficient RTs for 2'-O-methyl RNA and hexitol nucleic acids and the discovery of RTs for the orphan XNA chemistries D-altritol nucleic acid and 2'-methoxyethyl RNA, for which previously no RTs existed. Finally, we describe the engineering of XNA RTs with active exonucleolytic proofreading as well as the directed evolution of RNA RTs with very high complementary DNA synthesis fidelities, even in the absence of proofreading.

Resisting and tolerating microbes are alternative strategies to survive infection, but little is known about the evolutionary mechanisms controlling this balance. Here genomic analyses of anatomically modern humans, extinct Denisovan hominins and mice revealed a TNFAIP3 allelic series with alterations in the encoded immune response inhibitor A20. Each TNFAIP3 allele encoded substitutions at non-catalytic residues of the ubiquitin protease OTU domain that diminished IκB kinase-dependent phosphorylation and activation of A20. Two TNFAIP3 alleles encoding A20 proteins with partial phosphorylation deficits seemed to be beneficial by increasing immunity without causing spontaneous inflammatory disease: A20 T108A;I207L, originating in Denisovans and introgressed in modern humans throughout Oceania, and A20 I325N, from an N-ethyl-N-nitrosourea (ENU)-mutagenized mouse strain. By contrast, a rare human TNFAIP3 allele encoding an A20 protein with 95% loss of phosphorylation, C243Y, caused spontaneous inflammatory disease in humans and mice. Analysis of the partial-phosphorylation A20 I325N allele in mice revealed diminished tolerance of bacterial lipopolysaccharide and poxvirus inoculation as tradeoffs for enhanced immunity.

The physicochemical properties of nucleic acids are dominated by their highly charged phosphodiester backbone chemistry. This polyelectrolyte structure decouples information content (base sequence) from bulk properties, such as solubility, and has been proposed as a defining trait of all informational polymers. However, this conjecture has not been tested experimentally. Here, we describe the encoded synthesis of a genetic polymer with an uncharged backbone chemistry: alkyl phosphonate nucleic acids (phNAs) in which the canonical, negatively charged phosphodiester is replaced by an uncharged P-alkyl phosphonodiester backbone. Using synthetic chemistry and polymerase engineering, we describe the enzymatic, DNA-templated synthesis of P-methyl and P-ethyl phNAs, and the directed evolution of specific streptavidin-binding phNA aptamer ligands directly from random-sequence mixed P-methyl/P-ethyl phNA repertoires. Our results establish an example of the DNA-templated enzymatic synthesis and evolution of an uncharged genetic polymer and provide a foundational methodology for their exploration as a source of novel functional molecules.

Serine proteinase inhibitors (serpins), typically fold to a metastable native state and undergo a major conformational change in order to inhibit target proteases. However, conformational lability of the native serpin fold renders them susceptible to misfolding and aggregation, and underlies misfolding diseases such as α 1 -antitrypsin deficiency. Serpin specificity towards its protease target is dictated by its flexible and solvent exposed reactive centre loop (RCL), which forms the initial interaction with the target protease during inhibition. Previous studies have attempted to alter the specificity by mutating the RCL to that of a target serpin, but the rules governing specificity are not understood well enough yet to enable specificity to be engineered at will. In this paper, we use conserpin, a synthetic, thermostable serpin, as a model protein with which to investigate the determinants of serpin specificity by engineering its RCL. Replacing the RCL sequence with that from α1-antitrypsin fails to restore specificity against trypsin or human neutrophil elastase. Structural determination of the RCL-engineered conserpin and molecular dynamics simulations indicate that, although the RCL sequence may partially dictate specificity, local electrostatics and RCL dynamics may dictate the rate of insertion during protease inhibition, and thus whether it behaves as an inhibitor or a substrate. Engineering serpin specificity is therefore substantially more complex than solely manipulating the RCL sequence, and will require a more thorough understanding of how conformational dynamics achieves the delicate balance between stability, folding and function required by the exquisite serpin mechanism of action.

The M1 metallo-aminopeptidase from Plasmodium falciparum, PfA-M1, is an attractive drug target for the design of new antimalarials. Bestatin, a broad-spectrum metalloprotease inhibitor, is a moderate inhibitor of PfA-M1, and has been used to provide structure-activity relationships to inform drug design. The crystal structure of PfA-M1 with bestatin bound within its active site has been determined; however, dynamics of the inhibitor and the association or dissociation pathway have yet to be characterized. Here we present an all-atom molecular dynamics study where we have generated a hidden Markov state model from 2.3 μs of molecular dynamics simulation. Our hidden Markov state model identifies five macrostates that clearly show the events involved in bestatin dissociation from the PfA-M1 active site. The results show for the first time that bestatin can escape the substrate specificity pockets of the enzyme, primarily due to weak interactions within the pockets. Our approach identifies relevant conformational sampling of the inhibitor inside the enzyme and the protein dynamics that could be exploited to produce potent and selective inhibitors that can differentiate between similar members of the M1 aminopeptidase superfamily., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Recombination, the exchange of information between different genetic polymer strands, is of fundamental importance in biology for genome maintenance and genetic diversification and is mediated by dedicated recombinase enzymes. Here, we describe an innate capacity for non-enzymatic recombination (and ligation) in random-sequence genetic oligomer pools. Specifically, we examine random and semi-random eicosamer (N 20 ) pools of RNA, DNA and the unnatural genetic polymers ANA (arabino-), HNA (hexitol-) and AtNA (altritol-nucleic acids). While DNA, ANA and HNA pools proved inert, RNA (and to a lesser extent AtNA) pools displayed diverse modes of spontaneous intermolecular recombination, connecting recombination mechanistically to the vicinal ring cis-diol configuration shared by RNA and AtNA. Thus, the chemical constitution that renders both susceptible to hydrolysis emerges as the fundamental determinant of an innate capacity for recombination, which is shown to promote a concomitant increase in compositional, informational and structural pool complexity and hence evolutionary potential., Competing Interests: HM, AT, BP, AL, GH, MA, PH, PH No competing interests declared, (© 2018, Mutschler et al.)

Canonical mechanisms of protein evolution include the duplication and diversification of pre-existing folds through genetic alterations that include point mutations, insertions, deletions, and copy number amplifications, as well as post-translational modifications that modify processes such as folding efficiency and cellular localization. Following a survey of the human mutation database, we have identified an additional mechanism that we term "structural capacitance," which results in the de novo generation of microstructure in previously disordered regions. We suggest that the potential for structural capacitance confers select proteins with the capacity to evolve over rapid timescales, facilitating saltatory evolution as opposed to gradualistic canonical Darwinian mechanisms. Our results implicate the elements of protein microstructure generated by this distinct mechanism in the pathogenesis of a wide variety of human diseases. The benefits of rapidly furnishing the potential for evolutionary change conferred by structural capacitance are consequently counterbalanced by this accompanying risk. The phenomenon of structural capacitance has implications ranging from the ancestral diversification of protein folds to the engineering of synthetic proteins with enhanced evolvability., (Crown Copyright © 2018. Published by Elsevier Ltd. All rights reserved.)

The M1 and M17 aminopeptidases are metallo-exopeptidases that rely on the presence of divalent cations, usually zinc, in their active site for proteolytic activity. They are from separate protease superfamilies, however, members often have overlapping substrate specificity. Inhibitors of one or both enzymes can be used to modulate hypertension, reduce proliferation of certain types of cancers and control malaria parasites. Current inhibitors act to chelate the zinc ions in the active site, locking the enzymes in an inactive transition state. We were interested in using a computational approach to understand the structure and dynamics of the M1 and M17 aminopeptidases, however, the presence of the essential metal ions in the proteases presents a challenge to classical molecular dynamics (MD) simulation. The zinc amber force field does not contain applicable descriptions of the zinc coordination environment present in either of these two protease families. To provide tools for the study of these two enzymes, we have used the metal centre parameter builder to generate new hybrid bonded/nonbonded force field (FF) parameters to correctly describe the active site architecture for each enzyme. The new parameters were evaluated by fitting the normal mode frequencies of molecular mechanics to the quantum mechanics frequencies and validated by performing short MD simulations. The new FF parameters now enable more accurate and reliable MD simulations for any member of the M1 or M17 aminopeptidase superfamilies.

Antibodies have the specificity to differentiate foreign antigens that mimic self antigens, but it remains unclear how such specificity is acquired. In a mouse model, we generated B cells displaying an antibody that cross-reacts with two related protein antigens expressed on self versus foreign cells. B cell anergy was imposed by self antigen but reversed upon challenge with high-density foreign antigen, leading to germinal center recruitment and antibody gene hypermutation. Single-cell analysis detected rapid selection for mutations that decrease self affinity and slower selection for epistatic mutations that specifically increase foreign affinity. Crystal structures revealed that these mutations exploited subtle topological differences to achieve 5000-fold preferential binding to foreign over self epitopes. Resolution of antigenic mimicry drove the optimal affinity maturation trajectory, highlighting the value of retaining self-reactive clones as substrates for protective antibody responses., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Chemokines and their receptors collectively orchestrate the trafficking of leukocytes in normal immune function and inflammatory diseases. Different chemokines can induce distinct responses at the same receptor. In comparison to monocyte chemoattractant protein-1 (MCP-1; also known as CCL2), the chemokines MCP-2 (CCL8) and MCP-3 (CCL7) are partial agonists of their shared receptor CCR2, a key regulator of the trafficking of monocytes and macrophages that contribute to the pathology of atherosclerosis, obesity, and type 2 diabetes. Through experiments with chimeras of MCP-1 and MCP-3, we identified the chemokine amino-terminal region as being the primary determinant of both the binding and signaling selectivity of these two chemokines at CCR2. Analysis of CCR2 mutants showed that the chemokine amino terminus interacts with the major subpocket in the transmembrane helical bundle of CCR2, which is distinct from the interactions of some other chemokines with the minor subpockets of their receptors. These results suggest the major subpocket as a target for the development of small-molecule inhibitors of CCR2., (Copyright © 2017, American Association for the Advancement of Science.)

Ancestral protein reconstruction allows the resurrection and characterization of ancient proteins based on computational analyses of sequences of modern-day proteins. Unfortunately, many protein families are highly divergent and not suitable for sequence-based reconstruction approaches. This limitation is exemplified by the antigen receptors of jawed vertebrates (B- and T-cell receptors), heterodimers formed by pairs of Ig domains. These receptors are believed to have evolved from an extinct homodimeric ancestor through a process of gene duplication and diversification; however molecular evidence has so far remained elusive. Here, we use a structural approach and laboratory evolution to reconstruct such molecules and characterize their interaction with antigen. High-resolution crystal structures of reconstructed homodimeric receptors in complex with hen-egg white lysozyme demonstrate how nanomolar affinity binding of asymmetrical antigen is enabled through selective recruitment and structural plasticity within the receptor-binding site. Our results provide structural evidence in support of long-held theories concerning the evolution of antigen receptors, and provide a blueprint for the experimental reconstruction of protein ancestry in the absence of phylogenetic evidence.

The favorable biophysical attributes of non-antibody scaffolds make them attractive alternatives to monoclonal antibodies. However, due to the well-known stability-function trade-off, these gains tend to be marginal after functional selection. A notable example is the fibronectin Type III (FN3) domain, FNfn10, which has been previously evolved to bind lysozyme with 1 pM affinity (FNfn10-α-lys), but suffers from poor thermodynamic and kinetic stability. To explore this stability-function compromise further, we grafted the lysozyme-binding loops from FNfn10-α-lys onto our previously engineered, ultra-stable FN3 scaffold, FN3con. The resulting variant (FN3con-α-lys) bound lysozyme with a markedly reduced affinity, but retained high levels of thermal stability. The crystal structure of FNfn10-α-lys in complex with lysozyme revealed unanticipated interactions at the protein-protein interface involving framework residues of FNfn10-α-lys, thus explaining the failure to transfer binding via loop grafting. Utilizing this structural information, we redesigned FN3con-α-lys and restored picomolar binding affinity to lysozyme, while maintaining thermodynamic stability (with a thermal melting temperature 2-fold higher than that of FNfn10-α-lys). FN3con therefore provides an exceptional window of stability to tolerate deleterious mutations, resulting in a substantial advantage for functional design. This study emphasizes the utility of consensus design for the generation of highly stable scaffolds for downstream protein engineering studies., (© The Author 2016. Published by Oxford University Press.)

Enzymes must be ordered to allow the stabilization of transition states by their active sites, yet dynamic enough to adopt alternative conformations suited to other steps in their catalytic cycles. The biophysical principles that determine how specific protein dynamics evolve and how remote mutations affect catalytic activity are poorly understood. Here we examine a 'molecular fossil record' that was recently obtained during the laboratory evolution of a phosphotriesterase from Pseudomonas diminuta to an arylesterase. Analysis of the structures and dynamics of nine protein variants along this trajectory, and three rationally designed variants, reveals cycles of structural destabilization and repair, evolutionary pressure to 'freeze out' unproductive motions and sampling of distinct conformations with specific catalytic properties in bi-functional intermediates. This work establishes that changes to the conformational landscapes of proteins are an essential aspect of molecular evolution and that change in function can be achieved through enrichment of preexisting conformational sub-states.

The kallikrein-related peptidase (KLK) family of proteases is involved in many aspects of human health and disease. One member of this family, KLK4, has been implicated in cancer development and metastasis. Understanding mechanisms of inactivation are critical to developing selective KLK4 inhibitors. We have determined the X-ray crystal structures of KLK4 in complex with both sunflower trypsin inhibitor-1 (SFTI-1) and a rationally designed SFTI-1 derivative to atomic (~1 Å) resolution, as well as with bound nickel. These structures offer a structural rationalization for the potency and selectivity of these inhibitors, and together with MD simulation and computational analysis, reveal a dynamic pathway between the metal binding exosite and the active site, providing key details of a previously proposed allosteric mode of inhibition. Collectively, this work provides insight into both direct and indirect mechanisms of inhibition for KLK4 that have broad implications for the enzymology of the serine protease superfamily, and may potentially be exploited for the design of therapeutic inhibitors.

The rugged folding landscapes of functional proteins puts them at risk of misfolding and aggregation. Serine protease inhibitors, or serpins, are paradigms for this delicate balance between function and misfolding. Serpins exist in a metastable state that undergoes a major conformational change in order to inhibit proteases. However, conformational labiality of the native serpin fold renders them susceptible to misfolding, which underlies misfolding diseases such as α 1 -antitrypsin deficiency. To investigate how serpins balance function and folding, we used consensus design to create conserpin, a synthetic serpin that folds reversibly, is functional, thermostable, and polymerization resistant. Characterization of its structure, folding and dynamics suggest that consensus design has remodeled the folding landscape to reconcile competing requirements for stability and function. This approach may offer general benefits for engineering functional proteins that have risky folding landscapes, including the removal of aggregation-prone intermediates, and modifying scaffolds for use as protein therapeutics.

A popular and successful strategy in semi-rational design of protein stability is the use of evolutionary information encapsulated in homologous protein sequences. Consensus design is based on the hypothesis that at a given position, the respective consensus amino acid contributes more than average to the stability of the protein than non-conserved amino acids. Here, we review the consensus design approach, its theoretical underpinnings, successes, limitations and challenges, as well as providing a detailed guide to its application in protein engineering., (© The Author 2016. Published by Oxford University Press.)

Coiled-coils refer to a bundle of helices coiled together like strands of a rope. It has been estimated that nearly 3% of protein-encoding regions of genes harbour coiled-coil domains (CCDs). Experimental studies have confirmed that CCDs play a fundamental role in subcellular infrastructure and controlling trafficking of eukaryotic cells. Given the importance of coiled-coils, multiple bioinformatics tools have been developed to facilitate the systematic and high-throughput prediction of CCDs in proteins. In this article, we review and compare 12 sequence-based bioinformatics approaches and tools for coiled-coil prediction. These approaches can be categorized into two classes: coiled-coil detection and coiled-coil oligomeric state prediction. We evaluated and compared these methods in terms of their input/output, algorithm, prediction performance, validation methods and software utility. All the independent testing data sets are available at http://lightning.med.monash.edu/coiledcoil/. In addition, we conducted a case study of nine human polyglutamine (PolyQ) disease-related proteins and predicted CCDs and oligomeric states using various predictors. Prediction results for CCDs were highly variable among different predictors. Only two peptides from two proteins were confirmed to be CCDs by majority voting. Both domains were predicted to form dimeric coiled-coils using oligomeric state prediction. We anticipate that this comprehensive analysis will be an insightful resource for structural biologists with limited prior experience in bioinformatics tools, and for bioinformaticians who are interested in designing novel approaches for coiled-coil and its oligomeric state prediction., (© The Author 2015. Published by Oxford University Press. For Permissions, please email: journals.permissions@oup.com.)

Thyroid peroxidase (TPO) catalyses the biosynthesis of thyroid hormones and is a major autoantigen in Hashimoto's disease--the most common organ-specific autoimmune disease. Epitope mapping studies have shown that the autoimmune response to TPO is directed mainly at two surface regions on the molecule: immunodominant regions A and B (IDR-A, and IDR-B). TPO has been a major target for structural studies for over 20 years; however, to date, the structure of TPO remains to be determined. We have used a molecular modelling approach to investigate plausible modes of TPO structure and dimer organisation. Sequence features of the C-terminus are consistent with a coiled-coil dimerization motif that most likely anchors the TPO dimer in the apical membrane of thyroid follicular cells. Two contrasting models of TPO were produced, differing in the orientation and exposure of their active sites relative to the membrane. Both models are equally plausible based upon the known enzymatic function of TPO. The "trans" model places IDR-B on the membrane-facing side of the myeloperoxidase (MPO)-like domain, potentially hindering access of autoantibodies, necessitating considerable conformational change, and perhaps even dissociation of the dimer into monomers. IDR-A spans MPO- and CCP-like domains and is relatively fragmented compared to IDR-B, therefore most likely requiring domain rearrangements in order to coalesce into one compact epitope. Less epitope fragmentation and higher solvent accessibility of the "cis" model favours it slightly over the "trans" model. Here, IDR-B clusters towards the surface of the MPO-like domain facing the thyroid follicular lumen preventing steric hindrance of autoantibodies. However, conformational rearrangements may still be necessary to allow full engagement with autoantibodies, with IDR-B on both models being close to the dimer interface. Taken together, the modelling highlights the need to consider the oligomeric state of TPO, its conformational properties, and its proximity to the membrane, when interpreting epitope-mapping data.

We have previously reported a phage display method for the identification of protein domains on a genome-wide scale (shotgun proteolysis). Here we present the solution structure of a fragment of the Escherichia coli membrane protein yrfF, as identified by shotgun proteolysis, and determined by NMR spectroscopy. Despite the absence of computational predictions, the fragment formed a well-defined beta-barrel structure, distantly falling within the OB-fold classification. Our results highlight the potential of high-throughput experimental approaches for the identification of protein domains for structural studies., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Consensus protein design is a rapid and reliable technique for the improvement of protein stability, which relies on the use of homologous protein sequences. To enhance the stability of a fibronectin type III (FN3) domain, consensus design was employed using an alignment of 2123 sequences. The resulting FN3 domain, FN3con, has unprecedented stability, with a melting temperature >100°C, a ΔG(D-N) of 15.5 kcal mol(-1) and a greatly reduced unfolding rate compared with wild-type. To determine the underlying molecular basis for stability, an X-ray crystal structure of FN3con was determined to 2.0 Å and compared with other FN3 domains of varying stabilities. The structure of FN3con reveals significantly increased salt bridge interactions that are cooperatively networked, and a highly optimized hydrophobic core. Molecular dynamics simulations of FN3con and comparison structures show the cooperative power of electrostatic and hydrophobic networks in improving FN3con stability. Taken together, our data reveal that FN3con stability does not result from a single mechanism, but rather the combination of several features and the removal of non-conserved, unfavorable interactions. The large number of sequences employed in this study has most likely enhanced the robustness of the consensus design, which is now possible due to the increased sequence availability in the post-genomic era. These studies increase our knowledge of the molecular mechanisms that govern stability and demonstrate the rising potential for enhancing stability via the consensus method., (© The Author 2015. Published by Oxford University Press.)

The human neuroendocrine enzyme glutamate decarboxylase (GAD) catalyses the synthesis of the inhibitory neurotransmitter gamma-aminobutyric acid (GABA) using pyridoxal 5'-phosphate as a cofactor. GAD exists as two isoforms named according to their respective molecular weights: GAD65 and GAD67. Although cytosolic GAD67 is typically saturated with the cofactor (holoGAD67) and constitutively active to produce basal levels of GABA, the membrane-associated GAD65 exists mainly as the inactive apo form. GAD65, but not GAD67, is a prevalent autoantigen, with autoantibodies to GAD65 being detected at high frequency in patients with autoimmune (type 1) diabetes and certain other autoimmune disorders. The significance of GAD65 autoinactivation into the apo form for regulation of neurotransmitter levels and autoantibody reactivity is not understood. We have used computational and experimental approaches to decipher the nature of the holo → apo conversion in GAD65 and thus, its mechanism of autoinactivation. Molecular dynamics simulations of GAD65 reveal coupling between the C-terminal domain, catalytic loop, and pyridoxal 5'-phosphate-binding domain that drives structural rearrangement, dimer opening, and autoinactivation, consistent with limited proteolysis fragmentation patterns. Together with small-angle X-ray scattering and fluorescence spectroscopy data, our findings are consistent with apoGAD65 existing as an ensemble of conformations. Antibody-binding kinetics suggest a mechanism of mutually induced conformational changes, implicating the flexibility of apoGAD65 in its autoantigenicity. Although conformational diversity may provide a mechanism for cofactor-controlled regulation of neurotransmitter biosynthesis, it may also come at a cost of insufficient development of immune self-tolerance that favors the production of GAD65 autoantibodies.

Background: Thyroid peroxidase (TPO) is a dimeric membrane-bound enzyme of thyroid follicular cells, responsible for thyroid hormone biosynthesis. TPO is also a common target antigen in autoimmune thyroid disease (AITD). With two active sites, TPO is an unusual enzyme, and thus there is much interest in understanding its structure and role in AITD. Homology modeling has shown TPO to be composed of different structural modules, as well as a propeptide sequence. During the course of studies to obtain homogeneous preparations of recombinant TPO for structural studies, we investigated the role of the large propeptide sequence in TPO., Methods: An engineered recombinant human TPO preparation expressed in Chinese hamster ovary (CHO) cells lacking the propeptide (TPOΔpro; amino acid residues 21-108) was characterized and its properties compared to wild-type TPO. Plasma membrane localization was determined by cell surface protein biotinylation, and biochemical studies were performed to evaluate enzymatic activity and the effect of deglycosylation. Immunological investigations using autoantibodies from AITD patients and other epitope-specific antibodies that recognize conformational determinants on TPO were evaluated for binding to TPOΔpro by flow cytometry, immunocytochemistry, and capture enzyme-linked immunosorbent assay. Molecular modeling and dynamics simulation of TPOΔpro comprising a dimer of myeloperoxidase-like domains was performed in order to investigate the impact of propeptide removal and the role of glycosylation., Results: The TPOΔpro was expressed on the cell surface at comparable levels to wild-type TPO. The TPOΔpro was enzymatically active and recognized by patients' autoantibodies and a panel of epitope-specific antibodies, confirming structural integrity of the two major conformational determinants recognized by autoantibodies. Faithful intracellular trafficking and N-glycosylation of TPOΔpro was also maintained. Molecular modeling and dynamics simulations were consistent with these observations., Conclusions: Our results point to a redundant role for the propeptide sequence in TPO. The successful expression of TPOΔpro in a membrane-anchored, enzymatically active form that is insensitive to intramolecular proteolysis, and importantly is recognized by patients' autoantibodies, is a key advance for purification of substantial quantities of homogeneous preparation of TPO for crystallization, structural, and immunological studies.

In structural biology, management of a large number of Protein Data Bank (PDB) files and raw X-ray diffraction images often presents a major organizational problem. Existing software packages that manipulate these file types were not designed for these kinds of file-management tasks. This is typically encountered when browsing through a folder of hundreds of X-ray images, with the aim of rapidly inspecting the diffraction quality of a data set. To solve this problem, a useful functionality of the Macintosh operating system (OSX) has been exploited that allows custom visualization plugins to be attached to certain file types. Software plugins have been developed for diffraction images and PDB files, which in many scenarios can save considerable time and effort. The direct visualization of diffraction images and PDB structures in the file browser can be used to identify key files of interest simply by scrolling through a list of files.

Dihydrodipicolinate synthase (DHDPS) is an essential enzyme involved in the lysine biosynthesis pathway. DHDPS from E. coli is a homotetramer consisting of a 'dimer of dimers', with the catalytic residues found at the tight-dimer interface. Crystallographic and biophysical evidence suggest that the dimers associate to stabilise the active site configuration, and mutation of a central dimer-dimer interface residue destabilises the tetramer, thus increasing the flexibility and reducing catalytic efficiency and substrate specificity. This has led to the hypothesis that the tetramer evolved to optimise the dynamics within the tight-dimer. In order to gain insights into DHDPS flexibility and its relationship to quaternary structure and function, we performed comparative Molecular Dynamics simulation studies of native tetrameric and dimeric forms of DHDPS from E. coli and also the native dimeric form from methicillin-resistant Staphylococcus aureus (MRSA). These reveal a striking contrast between the dynamics of tetrameric and dimeric forms. Whereas the E. coli DHDPS tetramer is relatively rigid, both the E. coli and MRSA DHDPS dimers display high flexibility, resulting in monomer reorientation within the dimer and increased flexibility at the tight-dimer interface. The mutant E. coli DHDPS dimer exhibits disorder within its active site with deformation of critical catalytic residues and removal of key hydrogen bonds that render it inactive, whereas the similarly flexible MRSA DHDPS dimer maintains its catalytic geometry and is thus fully functional. Our data support the hypothesis that in both bacterial species optimal activity is achieved by fine tuning protein dynamics in different ways: E. coli DHDPS buttresses together two dimers, whereas MRSA dampens the motion using an extended tight-dimer interface.

The crystal structures of unliganded and liganded pMHC molecules provide a structural basis for TCR recognition yet they represent 'snapshots' and offer limited insight into dynamics that may be important for interaction and T cell activation. MHC molecules HLA-B*3501 and HLA-B*3508 both bind a 13 mer viral peptide (LPEP) yet only HLA-B*3508-LPEP induces a CTL response characterised by the dominant TCR clonetype SB27. HLA-B*3508-LPEP forms a tight and long-lived complex with SB27, but the relatively weak interaction between HLA-B*3501-LPEP and SB27 fails to trigger an immune response. HLA-B*3501 and HLA-B*3508 differ by only one amino acid (L/R156) located on α2-helix, but this does not alter the MHC or peptide structure nor does this polymorphic residue interact with the peptide or SB27. In the absence of a structural rationalisation for the differences in TCR engagement we performed a molecular dynamics study of both pMHC complexes and HLA-B*3508-LPEP in complex with SB27. This reveals that the high flexibility of the peptide in HLA-B*3501 compared to HLA-B*3508, which was not apparent in the crystal structure alone, may have an under-appreciated role in SB27 recognition. The TCR pivots atop peptide residues 6-9 and makes transient MHC contacts that extend those observed in the crystal structure. Thus MD offers an insight into 'scanning' mechanism of SB27 that extends the role of the germline encoded CDR2α and CDR2β loops. Our data are consistent with the vast body of experimental observations for the pMHC-LPEP-SB27 interaction and provide additional insights not accessible using crystallography.