Multi-state modeling of antibody-antigen complexes with SAXS profiles and deep-learning models.

Antibodies are an established class of human therapeutics. Epitope characterization is an important part of therapeutic antibody discovery. However, structural characterization of antibody-antigen complexes remains challenging. On the one hand, X-ray crystallography or cryo-electron microscopy provide atomic resolution characterization of the epitope, but the data collection process is typically long and the success rate is low. On the other hand, computational methods for modeling antibody-antigen structures from the individual components frequently suffer from a high false positive rate, rarely resulting in a unique solution. Recent deep learning models for structure prediction are also successful in predicting protein-protein complexes. However, they do not perform well for antibody-antigen complexes. Small Angle X-ray Scattering (SAXS) is a reliable technique for rapid structural characterization of protein samples in solution albeit at low resolution. Here, we present an integrative approach for modeling antigen-antibody complexes using the antibody sequence, antigen structure, and experimentally determined SAXS profiles of the antibody, antigen, and the complex. The method models antibody structures using a novel deep-learning approach, NanoNet. The structures of the antibodies and antigens are represented using multiple 3D conformations to account for compositional and conformational heterogeneity of the protein samples that are used to collect the SAXS data. The complexes are predicted by integrating the SAXS profiles with scoring functions for protein-protein interfaces that are based on statistical potentials and antibody-specific deep-learning models. We validated the method via application to four Fab:EGFR and one Fab:PCSK9 antibody:antigen complexes with experimentally available SAXS datasets. The integrative approach returns accurate predictions (interface RMSD<4Å) in the top five predictions for four out of five complexes (respective interface RMSD values of 1.95, 2.18, 2.66 and 3.87Å), providing support for the utility of such a computational pipeline for epitope characterization during therapeutic antibody discovery.

Pharmacogenomics Clinical Annotation Tool (PharmCAT).

Pharmacogenomics (PGx) decision support and return of results is an active area of precision medicine. One challenge of implementing PGx is extracting genomic variants and assigning haplotypes in order to apply prescribing recommendations and information from the Clinical Pharmacogenetics Implementation Consortium (CPIC), the US Food and Drug Administration (FDA), the Pharmacogenomics Knowledgebase (PharmGKB), etc. Pharmacogenomics Clinical Annotation Tool (PharmCAT) (i) extracts variants specified in guidelines from a genetic data set derived from sequencing or genotyping technologies, (ii) infers haplotypes and diplotypes, and (iii) generates a report containing genotype/diplotype-based annotations and guideline recommendations. We describe PharmCAT and a pilot validation project comparing results for 1000 Genomes Project sequences of Coriell samples with corresponding Genetic Testing Reference Materials Coordination Program (GeT-RM) sample characterization. PharmCAT was highly concordant with the GeT-RM data. PharmCAT is available in GitHub to evaluate, test, and report results back to the community. As precision medicine becomes more prevalent, our ability to consistently, accurately, and clearly define and report PGx annotations and prescribing recommendations is critical.

Prion Protein-Antibody Complexes Characterized by Chromatography-Coupled Small-Angle X-Ray Scattering.

Aberrant self-assembly, induced by structural misfolding of the prion proteins, leads to a number of neurodegenerative disorders. In particular, misfolding of the mostly α-helical cellular prion protein (PrP(C)) into a ß-sheet-rich disease-causing isoform (PrP(Sc)) is the key molecular event in the formation of PrP(Sc) aggregates. The molecular mechanisms underlying the PrP(C)-to-PrP(Sc) conversion and subsequent aggregation remain to be elucidated. However, in persistently prion-infected cell-culture models, it was shown that treatment with monoclonal antibodies against defined regions of the prion protein (PrP) led to the clearing of PrP(Sc) in cultured cells. To gain more insight into this process, we characterized PrP-antibody complexes in solution using a fast protein liquid chromatography coupled with small-angle x-ray scattering (FPLC-SAXS) procedure. High-quality SAXS data were collected for full-length recombinant mouse PrP [denoted recPrP(23-230)] and N-terminally truncated recPrP(89-230), as well as their complexes with each of two Fab fragments (HuM-P and HuM-R1), which recognize N- and C-terminal epitopes of PrP, respectively. In-line measurements by fast protein liquid chromatography coupled with SAXS minimized data artifacts caused by a non-monodispersed sample, allowing structural analysis of PrP alone and in complex with Fab antibodies. The resulting structural models suggest two mechanisms for how these Fabs may prevent the conversion of PrP(C) into PrP(Sc).

The accurate assessment of small-angle X-ray scattering data.

Small-angle X-ray scattering (SAXS) has grown in popularity in recent times with the advent of bright synchrotron X-ray sources, powerful computational resources and algorithms enabling the calculation of increasingly complex models. However, the lack of standardized data-quality metrics presents difficulties for the growing user community in accurately assessing the quality of experimental SAXS data. Here, a series of metrics to quantitatively describe SAXS data in an objective manner using statistical evaluations are defined. These metrics are applied to identify the effects of radiation damage, concentration dependence and interparticle interactions on SAXS data from a set of 27 previously described targets for which high-resolution structures have been determined via X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. The studies show that these metrics are sufficient to characterize SAXS data quality on a small sample set with statistical rigor and sensitivity similar to or better than manual analysis. The development of data-quality analysis strategies such as these initial efforts is needed to enable the accurate and unbiased assessment of SAXS data quality.

Mutations in RASA1 and GDF2 identified in patients with clinical features of hereditary hemorrhagic telangiectasia.

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant vascular disorder caused by mutations in ENG, ACVRL1 and SMAD4, which function in regulating the transforming growth factor beta and bone morphogenetic protein signaling pathways. Symptoms of HHT can be present in individuals who test negative for mutations in these three genes indicating other genes may be involved. In this study, we tested for mutations in two genes, RASA1 and GDF2, which were recently reported to be involved in vascular disorders. To determine whether RASA1 and GDF2 have phenotypic overlap with HHT and should be included in diagnostic testing, we developed a next-generation sequencing assay to detect mutations in 93 unrelated individuals who previously tested negative for mutations in ENG, ACVRL1 and SMAD4, but were clinically suspected to have HHT. Pathogenic mutations in RASA1 were identified in two samples (2.15%) and a variant of unknown significance in GDF2 was detected in one sample. All three individuals experienced epistaxis with dermal lesions described in medical records as telangiectases. These results indicate that the inclusion of RASA1 and GDF2 screening in individuals suspected to have HHT will increase the detection rate and aid clinicians in making an accurate diagnosis.

Nematic director reorientation at solid and liquid interfaces under flow: SAXS studies in a microfluidic device.

In this work we investigate the interplay between flow and boundary condition effects on the orientation field of a thermotropic nematic liquid crystal under flow and confinement in a microfluidic device. Two types of experiments were performed using synchrotron small-angle X-ray-scattering (SAXS). In the first, a nematic liquid crystal flows through a square-channel cross section at varying flow rates, while the nematic director orientation projected onto the velocity/velocity gradient plane is measured using a 2D detector. At moderate-to-high flow rates, the nematic director is predominantly aligned in the flow direction, but with a small tilt angle of ∼±11° in the velocity gradient direction. The director tilt angle is constant throughout most of the channel width but switches sign when crossing the center of the channel, in agreement with the Ericksen-Leslie-Parodi (ELP) theory. At low flow rates, boundary conditions begin to dominate, and a flow profile resembling the escaped radial director configuration is observed, where the director is seen to vary more smoothly from the edges (with homeotropic alignment) to the center of the channel. In the second experiment, hydrodynamic focusing is employed to confine the nematic phase into a sheet of liquid sandwiched between two layers of Triton X-100 aqueous solutions. The average nematic director orientation shifts to some extent from the flow direction toward the liquid boundaries, although it remains unclear if one tilt angle is dominant through most of the nematic sheet (with abrupt jumps near the boundaries) or if the tilt angle varies smoothly between two extreme values (∼90 and 0°). The technique presented here could be applied to perform high-throughput measurements for assessing the influence of different surfactants on the orientation of nematic phases and may lead to further improvements in areas such as boundary lubrication and clarifying the nature of defect structures in LC displays.

Structural basis of the pH-dependent assembly of a botulinum neurotoxin complex.

Botulinum neurotoxins (BoNTs) are among the most poisonous biological substances known. They assemble with non-toxic non-hemagglutinin (NTNHA) protein to form the minimally functional progenitor toxin complexes (M-PTC), which protects BoNT in the gastrointestinal tract and releases it upon entry into the circulation. Here we provide molecular insight into the assembly between BoNT/A and NTNHA-A using small-angle X-ray scattering. We found that the free form BoNT/A maintains a pH-independent conformation with limited domain flexibility. Intriguingly, the free form NTNHA-A adopts pH-dependent conformational changes due to a torsional motion of its C-terminal domain. Once forming a complex at acidic pH, they each adopt a stable conformation that is similar to that observed in the crystal structure of the M-PTC. Our results suggest that assembly of the M-PTC depends on the environmental pH and that the complex form of BoNT/A is induced by interacting with NTNHA-A at acidic pH.

Crystal structure of Streptococcus pyogenes EndoS, an immunomodulatory endoglycosidase specific for human IgG antibodies.

To evade host immune mechanisms, many bacteria secrete immunomodulatory enzymes. Streptococcus pyogenes, one of the most common human pathogens, secretes a large endoglycosidase, EndoS, which removes carbohydrates in a highly specific manner from IgG antibodies. This modification renders antibodies incapable of eliciting host effector functions through either complement or Fc γ receptors, providing the bacteria with a survival advantage. On account of this antibody-specific modifying activity, EndoS is being developed as a promising injectable therapeutic for autoimmune diseases that rely on autoantibodies. Additionally, EndoS is a key enzyme used in the chemoenzymatic synthesis of homogenously glycosylated antibodies with tailored Fc γ receptor-mediated effector functions. Despite the tremendous utility of this enzyme, the molecular basis of EndoS specificity for, and processing of, IgG antibodies has remained poorly understood. Here, we report the X-ray crystal structure of EndoS and provide a model of its encounter complex with its substrate, the IgG1 Fc domain. We show that EndoS is composed of five distinct protein domains, including glycosidase, leucine-rich repeat, hybrid Ig, carbohydrate binding module, and three-helix bundle domains, arranged in a distinctive V-shaped conformation. Our data suggest that the substrate enters the concave interior of the enzyme structure, is held in place by the carbohydrate binding module, and that concerted conformational changes in both enzyme and substrate are required for subsequent antibody deglycosylation. The EndoS structure presented here provides a framework from which novel endoglycosidases could be engineered for additional clinical and biotechnological applications.

Structure of the arginine methyltransferase PRMT5-MEP50 reveals a mechanism for substrate specificity.

The arginine methyltransferase PRMT5-MEP50 is required for embryogenesis and is misregulated in many cancers. PRMT5 targets a wide variety of substrates, including histone proteins involved in specifying an epigenetic code. However, the mechanism by which PRMT5 utilizes MEP50 to discriminate substrates and to specifically methylate target arginines is unclear. To test a model in which MEP50 is critical for substrate recognition and orientation, we determined the crystal structure of Xenopus laevis PRMT5-MEP50 complexed with S-adenosylhomocysteine (SAH). PRMT5-MEP50 forms an unusual tetramer of heterodimers with substantial surface negative charge. MEP50 is required for PRMT5-catalyzed histone H2A and H4 methyltransferase activity and binds substrates independently. The PRMT5 catalytic site is oriented towards the cross-dimer paired MEP50. Histone peptide arrays and solution assays demonstrate that PRMT5-MEP50 activity is inhibited by substrate phosphorylation and enhanced by substrate acetylation. Electron microscopy and reconstruction showed substrate centered on MEP50. These data support a mechanism in which MEP50 binds substrate and stimulates PRMT5 activity modulated by substrate post-translational modifications.

A model for 3-methyladenine recognition by 3-methyladenine DNA glycosylase I (TAG) from Staphylococcus aureus.

The removal of chemically damaged DNA bases such as 3-methyladenine (3-MeA) is an essential process in all living organisms and is catalyzed by the enzyme 3-MeA DNA glycosylase I. A key question is how the enzyme selectively recognizes the alkylated 3-MeA over the much more abundant adenine. The crystal structures of native and Y16F-mutant 3-MeA DNA glycosylase I from Staphylococcus aureus in complex with 3-MeA are reported to 1.8 and 2.2 Å resolution, respectively. Isothermal titration calorimetry shows that protonation of 3-MeA decreases its binding affinity, confirming previous fluorescence studies that show that charge-charge recognition is not critical for the selection of 3-MeA over adenine. It is hypothesized that the hydrogen-bonding pattern of Glu38 and Tyr16 of 3-MeA DNA glycosylase I with a particular tautomer unique to 3-MeA contributes to recognition and selection.

The Protein Structure Initiative Structural Biology Knowledgebase Technology Portal: a structural biology web resource.

The Technology Portal of the Protein Structure Initiative Structural Biology Knowledgebase (PSI SBKB; http://technology.sbkb.org/portal/ ) is a web resource providing information about methods and tools that can be used to relieve bottlenecks in many areas of protein production and structural biology research. Several useful features are available on the web site, including multiple ways to search the database of over 250 technological advances, a link to videos of methods on YouTube, and access to a technology forum where scientists can connect, ask questions, get news, and develop collaborations. The Technology Portal is a component of the PSI SBKB ( http://sbkb.org ), which presents integrated genomic, structural, and functional information for all protein sequence targets selected by the Protein Structure Initiative. Created in collaboration with the Nature Publishing Group, the SBKB offers an array of resources for structural biologists, such as a research library, editorials about new research advances, a featured biological system each month, and a functional sleuth for searching protein structures of unknown function. An overview of the various features and examples of user searches highlight the information, tools, and avenues for scientific interaction available through the Technology Portal.

Ni(II) coordination to mixed sites modulates DNA binding of HpNikR via a long-range effect.

Helicobacter pylori NikR (HpNikR) is a nickel-dependent transcription factor that regulates multiple genes in the H. pylori pathogen. There are conflicting data regarding the locations of the Ni(II) sites and the role of Ni(II) coordination in DNA recognition. Herein, we report crystal structures of (i) the metal-binding domain (MBD) of HpNikR (3.08 Å) and (ii) a mutant, H74A (2.04 Å), designed to disrupt native Ni(II) coordination. In the MBD structure, four nickel ions are coordinated to two different types of nickel sites (4-coordinate, square planar, and 5/6-coordinate, square pyramidal/octahedral). In the H74A structure, all four nickel ions are coordinated to 4-coordinate square-planar sites. DNA-binding studies reveal tighter binding for target DNA sequences for holo-HpNikR compared with the affinities of Ni(II) reconstituted apo-HpNikR and H74A for these same DNA targets, supporting a role for Ni(II) coordination to 5/6 sites in DNA recognition. Small-angle X-ray scattering studies of holo-HpNikR and H74A reveal a high degree of conformational flexibility centered at the DNA-binding domains of H74A, which is consistent with disorder observed in the crystal structure of the protein. A model of DNA recognition by HpNikR is proposed in which Ni(II) coordination to specific sites in the MBD have a long-range effect on the flexibility of the DNA-binding domains and, consequently, the DNA recognition properties.

The Structural Biology Knowledgebase: a portal to protein structures, sequences, functions, and methods.

The Protein Structure Initiative's Structural Biology Knowledgebase (SBKB, URL: http://sbkb.org ) is an open web resource designed to turn the products of the structural genomics and structural biology efforts into knowledge that can be used by the biological community to understand living systems and disease. Here we will present examples on how to use the SBKB to enable biological research. For example, a protein sequence or Protein Data Bank (PDB) structure ID search will provide a list of related protein structures in the PDB, associated biological descriptions (annotations), homology models, structural genomics protein target status, experimental protocols, and the ability to order available DNA clones from the PSI:Biology-Materials Repository. A text search will find publication and technology reports resulting from the PSI's high-throughput research efforts. Web tools that aid in research, including a system that accepts protein structure requests from the community, will also be described. Created in collaboration with the Nature Publishing Group, the Structural Biology Knowledgebase monthly update also provides a research library, editorials about new research advances, news, and an events calendar to present a broader view of structural genomics and structural biology.

Inhibition of the PLP-dependent enzyme serine palmitoyltransferase by cycloserine: evidence for a novel decarboxylative mechanism of inactivation.

Cycloserine (CS, 4-amino-3-isoxazolidone) is a cyclic amino acid mimic that is known to inhibit many essential pyridoxal 5'-phosphate (PLP)-dependent enzymes. Two CS enantiomers are known; D-cycloserine (DCS, also known as Seromycin) is a natural product that is used to treat resistant Mycobacterium tuberculosis infections as well as neurological disorders since it is a potent NMDA receptor agonist, and L-cycloserine (LCS) is a synthetic enantiomer whose usefulness as a drug has been hampered by its inherent toxicity arising through inhibition of sphingolipid metabolism. Previous studies on various PLP-dependent enzymes revealed a common mechanism of inhibition by both enantiomers of CS; the PLP cofactor is disabled by forming a stable 3-hydroxyisoxazole/pyridoxamine 5'-phosphate (PMP) adduct at the active site where the cycloserine ring remains intact. Here we describe a novel mechanism of CS inactivation of the PLP-dependent enzyme serine palmitoyltransferase (SPT) from Sphingomonas paucimobilis. SPT catalyses the condensation of l-serine and palmitoyl-CoA, the first step in the de novo sphingolipid biosynthetic pathway. We have used a range of kinetic, spectroscopic and structural techniques to postulate that both LCS and DCS inactivate SPT by transamination to form a free pyridoxamine 5'-phosphate (PMP) and beta-aminooxyacetaldehyde that remain bound at the active site. We suggest this occurs by ring opening of the cycloserine ring followed by decarboxylation. Enzyme kinetics show that inhibition is reversed by incubation with excess PLP and that LCS is a more effective SPT inhibitor than DCS. UV-visible spectroscopic data, combined with site-directed mutagenesis, suggest that a mobile Arg(378) residue is involved in cycloserine inactivation of SPT.

The Scottish Structural Proteomics Facility: targets, methods and outputs.

The Scottish Structural Proteomics Facility was funded to develop a laboratory scale approach to high throughput structure determination. The effort was successful in that over 40 structures were determined. These structures and the methods harnessed to obtain them are reported here. This report reflects on the value of automation but also on the continued requirement for a high degree of scientific and technical expertise. The efficiency of the process poses challenges to the current paradigm of structural analysis and publication. In the 5 year period we published ten peer-reviewed papers reporting structural data arising from the pipeline. Nevertheless, the number of structures solved exceeded our ability to analyse and publish each new finding. By reporting the experimental details and depositing the structures we hope to maximize the impact of the project by allowing others to follow up the relevant biology.

Extensive DNA mimicry by the ArdA anti-restriction protein and its role in the spread of antibiotic resistance.

The ardA gene, found in many prokaryotes including important pathogenic species, allows associated mobile genetic elements to evade the ubiquitous Type I DNA restriction systems and thereby assist the spread of resistance genes in bacterial populations. As such, ardA contributes to a major healthcare problem. We have solved the structure of the ArdA protein from the conjugative transposon Tn916 and find that it has a novel extremely elongated curved cylindrical structure with defined helical grooves. The high density of aspartate and glutamate residues on the surface follow a helical pattern and the whole protein mimics a 42-base pair stretch of B-form DNA making ArdA by far the largest DNA mimic known. Each monomer of this dimeric structure comprises three alpha-beta domains, each with a different fold. These domains have the same fold as previously determined proteins possessing entirely different functions. This DNA mimicry explains how ArdA can bind and inhibit the Type I restriction enzymes and we demonstrate that 6 different ardA from pathogenic bacteria can function in Escherichia coli hosting a range of different Type I restriction systems.

Structural and functional characterisation of a conserved archaeal RadA paralog with antirecombinase activity.

DNA recombinases (RecA in bacteria, Rad51 in eukarya and RadA in archaea) catalyse strand exchange between homologous DNA molecules, the central reaction of homologous recombination, and are among the most conserved DNA repair proteins known. RecA is the sole protein responsible for this reaction in bacteria, whereas there are several Rad51 paralogs that cooperate to catalyse strand exchange in eukaryotes. All archaea have at least one (and as many as four) RadA paralog, but their function remains unclear. Herein, we show that the three RadA paralogs encoded by the Sulfolobus solfataricus genome are expressed under normal growth conditions and are not UV inducible. We demonstrate that one of these proteins, Sso2452, which is representative of the large archaeal RadC subfamily of archaeal RadA paralogs, functions as an ATPase that binds tightly to single-stranded DNA. However, Sso2452 is not an active recombinase in vitro and inhibits D-loop formation by RadA. We present the high-resolution crystal structure of Sso2452, which reveals key structural differences from the canonical RecA family recombinases that may explain its functional properties. The possible roles of the archaeal RadA paralogs in vivo are discussed.

The external aldimine form of serine palmitoyltransferase: structural, kinetic, and spectroscopic analysis of the wild-type enzyme and HSAN1 mutant mimics.

Sphingolipid biosynthesis begins with the condensation of L-serine and palmitoyl-CoA catalyzed by the PLP-dependent enzyme serine palmitoyltransferase (SPT). Mutations in human SPT cause hereditary sensory autonomic neuropathy type 1, a disease characterized by loss of feeling in extremities and severe pain. The human enzyme is a membrane-bound hetereodimer, and the most common mutations are located in the enzymatically incompetent monomer, suggesting a "dominant" or regulatory effect. The molecular basis of how these mutations perturb SPT activity is subtle and is not simply loss of activity. To further explore the structure and mechanism of SPT, we have studied the homodimeric bacterial enzyme from Sphingomonas paucimobilis. We have analyzed two mutants (N100Y and N100W) engineered to mimic the mutations seen in hereditary sensory autonomic neuropathy type 1 as well as a third mutant N100C designed to mimic the wild-type human SPT. The N100C mutant appears fully active, whereas both N100Y and N100W are significantly compromised. The structures of the holoenzymes reveal differences around the active site and in neighboring secondary structure that transmit across the dimeric interface in both N100Y and N100W. Comparison of the l-Ser external aldimine structures of both native and N100Y reveals significant differences that hinder the movement of a catalytically important Arg(378) residue into the active site. Spectroscopic analysis confirms that both N100Y and N100W mutants subtly affect the chemistry of the PLP. Furthermore, the N100Y and R378A mutants appear less able to stabilize a quinonoid intermediate. These data provide the first experimental insight into how the most common disease-associated mutations of human SPT may lead to perturbation of enzyme activity.

AcsD catalyzes enantioselective citrate desymmetrization in siderophore biosynthesis.

Bacterial pathogens need to scavenge iron from their host for growth and proliferation during infection. They have evolved several strategies to do this, one being the biosynthesis and excretion of small, high-affinity iron chelators known as siderophores. The biosynthesis of siderophores is an important area of study, not only for potential therapeutic intervention but also to illuminate new enzyme chemistries. Two general pathways for siderophore biosynthesis exist: the well-characterized nonribosomal peptide synthetase (NRPS)-dependent pathway and the NRPS-independent siderophore (NIS) pathway, which relies on a different family of sparsely investigated synthetases. Here we report structural and biochemical studies of AcsD from Pectobacterium (formerly Erwinia) chrysanthemi, an NIS synthetase involved in achromobactin biosynthesis. The structures of ATP and citrate complexes provide a mechanistic rationale for stereospecific formation of an enzyme-bound (3R)-citryladenylate, which reacts with L-serine to form a likely achromobactin precursor. AcsD is a unique acyladenylate-forming enzyme with a new fold and chemical catalysis strategy.