Denisovan, modern human and mouse TNFAIP3 alleles tune A20 phosphorylation and immunity.

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.

Macromolecular condensation buffers intracellular water potential.

Optimum protein function and biochemical activity critically depends on water availability because solvent thermodynamics drive protein folding and macromolecular interactions1. 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 timescales2,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.

Conformational diversity facilitates antibody mutation trajectories and discrimination between foreign and self-antigens.

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.

Mutational and biophysical robustness in a prestabilized monobody.

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.

Molecular basis of a redox switch: molecular dynamics simulations and surface plasmon resonance provide insight into reduced and oxidised angiotensinogen.

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.

Structural reconstruction of protein ancestry.

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.

Critical evaluation of in silico methods for prediction of coiled-coil domains in proteins.

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 role of protein dynamics in the evolution of new enzyme function.

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.

Cofactor-dependent conformational heterogeneity of GAD65 and its role in autoimmunity and neurotransmitter homeostasis.

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.

Rapid discovery of high-affinity antibodies via massively parallel sequencing, ribosome display and affinity screening.

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 108 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 × 106 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 × 105 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.

Epitope flexibility and dynamic footprint revealed by molecular dynamics of a pMHC-TCR complex.

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.

Emergence of ATP- and GTP-Binding Aptamers from Single RNA Sequences by Error-Prone Replication and Selection.

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.

Discovery and evolution of RNA and XNA reverse transcriptase function and fidelity.

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.

Reactive centre loop dynamics and serpin specificity.

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.

A synthetic genetic polymer with an uncharged backbone chemistry based on alkyl phosphonate nucleic acids.

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.

Generation of AMBER force field parameters for zinc centres of M1 and M17 family aminopeptidases.

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.

Mapping the Pathway and Dynamics of Bestatin Inhibition of the Plasmodium falciparum M1 Aminopeptidase PfA-M1.

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.

Random-sequence genetic oligomer pools display an innate potential for ligation and recombination.

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

Structural Capacitance in Protein Evolution and Human Diseases.

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.