Single Amino Acid Substitution Research Articles

A Single Amino Acid Substitution in Porcine Reproductive and Respiratory Syndrome Virus Glycoprotein 2 Significantly Impairs Its Infectivity in Macrophages.

Porcine reproductive and respiratory syndrome virus (PRRSV) has a restricted tropism for macrophages and CD163 is a key receptor for infection. In this study, the PRRSV strain NCV1 was passaged on MARC-145 cells for 95 passages, and two plaque-clones (C1 and C2) were randomly selected for further analysis. The C1 virus nearly lost the ability to infect porcine alveolar macrophages (PAMs), as well as porcine kidney cells expressing porcine CD163 (PK15-pCD163), while the C2 virus replicates well in these two cell types. Pretreatment of MARC-145 cells with an anti-CD163 antibody nearly blocked C1 virus infection, indicating that the virus still required CD163 to infect cells. The C1 virus carried four unique amino acid substitutions: three in the nonstructural proteins and a K160I in GP2. The introduction of an I160K substitution in GP2 of the C1 virus restored its infectivity in PAMs and PK15-pCD163 cells, while the introduction of a K160I substitution in GP2 of the low-passaged, virulent PRRSV strain NCV13 significantly impaired its infectivity. Importantly, pigs inoculated with the rNCV13-K160I mutant exhibited lower viremia levels and lung lesions than those infected with the parental rNCV13. These results demonstrated that the K160 residue in GP2 is one of the key determinants of PRRSV tropism.

Difference in the Inhibitory Effect of Thiol Compounds and Demetallation Rates from the Zn(II) Active Site of Metallo-β-lactamases (IMP-1 and IMP-6) Associated with a Single Amino Acid Substitution

Gram-negative bacteria producing metallo-β-lactamases (MBLs) have become a considerable threat to public health. MBLs including the IMP, VIM, and NDM types are Zn(II) enzymes that hydrolyze the β-lactam ring present in a broad range of antibiotics, such as N-benzylpenicillin, meropenem, and imipenem. Among IMPs, IMP-1 and IMP-6 differ in a single amino acid substitution at position 262, where serine in IMP-1 is replaced by glycine in IMP-6, conferring a change in substrate specificity. To investigate how this mutation influences enzyme function, we examined lactamase inhibition by thiol compounds. Ethyl 3-mercaptopropionate acted as a competitive inhibitor of IMP-1, but a noncompetitive inhibitor of IMP-6. A comparison of the crystal structures previously reported for IMP-1 (PDB code: 5EV6) and IMP-6 (PDB code: 6LVJ) revealed a hydrogen bond between the side chain of Ser262 and Cys221 in IMP-1 but the absence of hydrogen bond in IMP-6, which affects the Zn2 coordination sphere in its active site. We investigated the demetallation rates of IMP-1 and IMP-6 in the presence of chelating agent ethylenediaminetetraacetic acid (EDTA) and found that the demetallation reactions had fast and slow phases with a first-order rate constant (kfast = 1.76 h-1, kslow = 0.108 h-1 for IMP-1, and kfast = 14.0 h-1 and kslow = 1.66 h-1 for IMP-6). The difference in the flexibility of the Zn2 coordination sphere between IMP-1 and IMP-6 may influence the demetallation rate, the catalytic efficiency against β-lactam antibiotics, and the inhibitory effect of thiol compounds.

304. Viral resistance analysis in the COMET-PEAK study: sotrovimab treatment in participants with mild-to-moderate COVID-19

Abstract Background Sotrovimab is a human monoclonal antibody targeting a conserved region of the SARS-CoV-2 spike (S) protein. COMET-PEAK was a 3-part, phase 2 study that evaluated intravenous (500 mg) and intramuscular (500 mg and 250 mg) administration of sotrovimab in outpatients (n=353) with mild-to-moderate COVID-19. We assessed amino acid substitutions in the SARS-CoV-2 S protein and circulating variants of concern/interest (VOC/VOI) in COMET-PEAK participants (enrolled Feb–July 2021). Methods Mid-turbinate (Part A) or nasopharyngeal (Part B/C) samples were obtained from all participants at Baseline and Post-Baseline visits. Next generation sequencing of the SARS-CoV-2 S gene was conducted using Illumina MiSeq with a ≥5% allelic frequency cut-off for samples with a viral load above 3.0 log10 copies/mL. Baseline, post-baseline and treatment-emergent (TE) substitutions were assessed, and prevalence of VOC/VOI was evaluated. Phenotypic analyses of epitope substitutions were conducted using a pseudotyped virus assay. Results In total, 282/353 participants had sequencing results for ≥1 visit (253 baseline, 248 post-baseline), and 219 had paired baseline and post-baseline sequences; among these, 149 (68%) had TE substitutions in the S protein (26 [12%] in the epitope). E340K was the predominant TE substitution in the epitope (15 [7%]). Across all arms 92/245 (38%) experienced virologic rebound, 8 of whom (Part A: 2; Part B: 2; Part C: 4) had TE substitutions in the epitope; none had evidence of clinical progression to severe disease. Prevalence of VOC/VOI or single amino acid substitutions of concern was 94% (266/282); the most frequent were Alpha (Part A: 8/16 [50%]; Part B: 75/128 [59%]) and Delta (Part C: 99/122 [81%]). Of 7 participants with evidence of clinical progression, none had S protein substitutions in the epitope and all had VOC/VOI (3 Alpha, 3 Delta, 1 Gamma). Sotrovimab effectively neutralized most epitope substitutions tested in vitro; P337L and E340A/K/V conferred significantly reduced susceptibility. Conclusion There was no evidence that sotrovimab epitope substitutions were associated with clinical progression or virologic rebound. These data are consistent with those from the COMET-ICE study. Funding Vir Biotechnology/GSK (NCT04779879). Disclosures Phillip J Yates, PhD, GlaxoSmithKline: Employee Jennifer Moore, MD, GlaxoSmithKline: Employee Jennifer Han, MD, GlaxoSmithKline: Employee Andrew Skingsley, MD, GlaxoSmithKline: Employee|GlaxoSmithKline: Stocks/Bonds Gretja Schnell, PhD, Vir Biotechnology: Employee|Vir Biotechnology: Stocks/Bonds Andrea L. Cathcart, PhD, Vir Biotechnology: Employee|Vir Biotechnology: Stocks/Bonds Melissa Aldinger, PharmD, Vir Biotechnology: Employee|Vir Biotechnology: Stocks/Bonds Amanda Peppercorn, MD, GlaxoSmithKline: Employee|GlaxoSmithKline: Stocks/Bonds Jill Walker, PhD, GlaxoSmithKline: Employee.

Single amino acid residue mediates reciprocal specificity in two mosquito odorant receptors.

The southern house mosquito, Culex quinquefasciatus, utilizes two odorant receptors, CquiOR10 and CquiOR2, narrowly tuned to oviposition attractants and well conserved among mosquito species. They detect skatole and indole, respectively, with reciprocal specificity. We swapped the transmembrane (TM) domains of CquiOR10 and CquiOR2 and identified TM2 as a specificity determinant. With additional mutations, we showed that CquiOR10A73L behaved like CquiOR2. Conversely, CquiOR2L74A recapitulated CquiOR10 specificity. Next, we generated structural models of CquiOR10 and CquiOR10A73L using RoseTTAFold and AlphaFold and docked skatole and indole using RosettaLigand. These modeling studies suggested space-filling constraints around A73. Consistent with this hypothesis, CquiOR10 mutants with a bulkier residue (Ile, Val) were insensitive to skatole and indole, whereas CquiOR10A73G retained the specificity to skatole and showed a more robust response than the wildtype receptor CquiOR10. On the other hand, Leu to Gly mutation of the indole receptor CquiOR2 reverted the specificity to skatole. Lastly, CquiOR10A73L, CquiOR2, and CquiOR2L74I were insensitive to 3-ethylindole, whereas CquiOR2L74A and CquiOR2L74G gained activity. Additionally, CquiOR10A73G gave more robust responses to 3-ethylindole than CquiOR10. Thus, we suggest the specificity of these receptors is mediated by a single amino acid substitution, leading to finely tuned volumetric space to accommodate specific oviposition attractants.

Genome-wide unraveling SNP pairwise epistatic effects associated with sheep body weight

Epistatic effects are an important part of the genetic effect of complex traits in livestock. In this study, we used 218 synthetic ewes from the Xinjiang Academy of Agricultural Reclamation in China to identify interacting paired with genome-wide single nucleotide polymorphisms (SNPs) associated with birth weight, weaning weight, and one-yearling weight. We detected 2 and 66 SNP-SNP interactions of sheep birth weight and weaning weight, respectively. No significant epistatic interaction of one-year-old body weight was detected. The genetic interaction of sheep body weight is dynamic and time-dependent. Most significant interactions of weaning body weight contributed 1% or higher. In the weaning weight trait, 66 significant SNP pairs consisted of 98 single SNPs covering 23 chromosomes, 5 of which were nonsynonymous SNPs (nsSNPs), resulting in single amino acid substitution. We found that genes that interact with transcription factors (TFs) are target genes for the corresponding TFs. Four epitron networks affecting weaning weight, including subnetworks of HIVEP3 and BACH2 transcription factors, constructed using significant SNP pairs, were also analyzed and annotated. These results suggest that transcription factors may play an important role in explaining epistatic effects. It provides a new idea to study the genetic mechanism of weight developing.

A fast-growing dengue virus mutant reveals a dual role of STING in response to infection.

The four dengue viruses (DENVs) have evolved multiple mechanisms to ensure its survival. Among these mechanisms is the ability to regulate its replication rate, which may contribute to avoiding premature immune activation that limit infection dissemination: DENVs associated with dengue epidemics have shown slower replication rate than pre-epidemic strains. Correspondingly, wild-type DENVs replicate more slowly than their clinically attenuated derivatives. To understand how DENVs 'make haste slowly', we generated and screened for DENV2 mutants with accelerated replication that also induced high type-I interferon (IFN) expression in infected cells. We chanced upon a single NS2B-I114T amino acid substitution, in an otherwise highly conserved amino acid residue. Accelerated DENV2 replication damaged host DNA as mutant infection was dependent on host DNA damage repair factors, namely RAD21, EID3 and NEK5. DNA damage induced cGAS/STING signalling and activated early type-I IFN response that inhibited infection dissemination. Unexpectedly, STING activation also supported mutant DENV replication in infected cells through STING-induced autophagy. Our findings thus show that DENV NS2B has multi-faceted role in controlling DENV replication rate and immune evasion and suggest that the dual role of STING in supporting virus replication within infected cells but inhibiting infection dissemination could be particularly advantageous for live attenuated vaccine development.

Prophage-encoded small protein YqaH counteracts the activities of the replication initiator DnaA in Bacillus subtilis.

Bacterial genomes harbour cryptic prophages that are mostly transcriptionally silent with many unannotated genes. Still, cryptic prophages may contribute to their host fitness and phenotypes. In Bacillus subtilis, the yqaF-yqaN operon belongs to the prophage element skin, and is tightly repressed by the Xre-like repressor SknR. This operon contains several small ORFs (smORFs) potentially encoding small-sized proteins. The smORF-encoded peptide YqaH was previously reported to bind to the replication initiator DnaA. Here, using a yeast two-hybrid assay, we found that YqaH binds to the DNA binding domain IV of DnaA and interacts with Spo0A, a master regulator of sporulation. We isolated single amino acid substitutions in YqaH that abolished the interaction with DnaA but not with Spo0A. Then, using a plasmid-based inducible system to overexpress yqaH WT and mutant derivatives, we studied in B. subtilis the phenotypes associated with the specific loss-of-interaction with DnaA (DnaA_LOI). We found that expression of yqaH carrying DnaA_LOI mutations abolished the deleterious effects of yqaH WT expression on chromosome segregation, replication initiation and DnaA-regulated transcription. When YqaH was induced after vegetative growth, DnaA_LOI mutations abolished the drastic effects of YqaH WT on sporulation and biofilm formation. Thus, YqaH inhibits replication, sporulation and biofilm formation mainly by antagonizing DnaA in a manner that is independent of the cell cycle checkpoint Sda.

Accurate protein stability predictions from homology models

Calculating changes in protein stability (ΔΔG) has been shown to be central for predicting the consequences of single amino acid substitutions in protein engineering as well as interpretation of genomic variants for disease risk. Structure-based calculations are considered most accurate, however the tools used to calculate ΔΔGs have been developed on experimentally resolved structures. Extending those calculations to homology models based on related proteins would greatly extend their applicability as large parts of e.g. the human proteome are not structurally resolved. In this study we aim to investigate the accuracy of ΔΔG values predicted on homology models compared to crystal structures. Specifically, we identified four proteins with a large number of experimentally tested ΔΔGs and templates for homology modeling across a broad range of sequence identities, and selected three methods for ΔΔG calculations to test. We find that ΔΔG-values predicted from homology models compare equally well to experimental ΔΔGs as those predicted on experimentally established crystal structures, as long as the sequence identity of the model template to the target protein is at least 40%. In particular, the Rosetta cartesian_ddg protocol is robust against the small perturbations in the structure which homology modeling introduces. In an independent assessment, we observe a similar trend when using ΔΔGs to categorize variants as low or wild-type-like abundance. Overall, our results show that stability calculations performed on homology models can substitute for those on crystal structures with acceptable accuracy as long as the model is built on a template with sequence identity of at least 40% to the target protein.

Engineered Single Amino Acid Substitutions Protect Hematopoietic Stem and Progenitor Cells from CD123 Targeted Immunotherapy

Targeted eradication of transformed or otherwise dysregulated cells using monoclonal antibodies (mAb), T cell engagers (TCE) or chimeric antigen receptor (CAR) cells is very effective for hematologic diseases. CD123, the interleukin-3 (IL-3) receptor alpha-chain, is highly expressed in various hematological malignancies, including acute myeloid leukemia (AML) and blastic plasmacytoid dendritic cell neoplasm (BPDCN). However, shared expression of CD123 on hematopoietic stem and progenitor cells (HSPCs) bears the risk for extensive myelotoxicity upon targeted depletion. To enable therapeutic targeting of non-essential proteins post HSC transplantation (HSCT), it was proposed to engineer HSPCs that lack the target protein. Here, we show that expression and function of a genetically modified CD123 receptor can be preserved in genome engineered therapeutic cells despite full shielding from targeted therapy. We used the known co-crystal structure of CD123 bound to the mAb CSL362 (talacotuzumab) for computer-aided protein engineering. In silico, 28 candidate single amino acid (aa) substitutions at 3 key residues were designed to interfere with CSL362 binding while potentially preserving its function. Flow cytometry confirmed no binding of a CSL362 biosimilar to 13 variants and weak binding to 12 variants while 3 variants showed similar binding as wildtype (wt) CD123. Protein expression of all variants was comparable to wt CD123. Using real-time, label-free protein interaction measurements, we could not detect interactions between the non-binding variant receptors and the CSL362 biosimilar, even at high protein concentrations. Despite this, the non-binding variants preserved the normal binding of IL-3. Protein characterization of selected non- and weak binding variants demonstrated normal biophysical properties of all but one variant. All 13 non-binding variants abolished antibody dependent cellular cytotoxicity (ADCC) of the CSL362 IgG1 biosimilar in a surrogate reporter assay. In addition, a CSL362-derived TCE resulted in T cell activation and target cell killing of wt and strong binding CD123 variants while non-binding variants neither activated the T cells nor resulted in cell killing. Finally, the majority of non-binding variants (12/13) were completely protected from CSL362-CAR T cell cytotoxicity. To assess the function of the variant receptors we engineered the top 2 variants into the IL-3 dependent cell line TF-1 using CRISPR/Cas9-mediated homology directed repair (HDR). Engineered variant TF-1 cells demonstrated IL-3 dose-dependent growth comparable to their unedited counterparts. The addition of the CSL362 biosimilar inhibited IL-3 dependent cell growth in unedited TF-1 cells but not cells expressing variant receptors. Similarly, CRISPR/Cas9-engineered human CD34+ HSPC knock-in (KI) cells expressing the selected variants (Fig. 1a) were shielded from TCE killing using HLA-matched T cells in vitro (Fig. 1b). In addition, a set of in vitro experiments performed with CD34+ HSPC KI cells suggested normal differentiation potential and intact IL-3 signaling as measured by pSTAT5. In summary, our data demonstrate that single aa substitutions in CD123 introduced to HSPCs by CRISPR/Cas9-mediated HDR enable selective shielding from various immunotherapies, while preserving biophysical and functional properties. Variant IL-3 receptor engineered HSPC's and paired immunotherapies could enable novel and targeted therapeutic interventions post HSCT to either treat minimal residual disease or be used as a salvage therapy. Such approaches could expand the therapeutic window and broaden indications and applications for targeted immunotherapies for various malignant and benign conditions. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

Mutations in the α4-α5 allosteric lobe of RAS do not significantly impair RAS signaling or self-association

Mutations in one of the three RAS genes (HRAS, KRAS, and NRAS) are present in nearly 20% of all human cancers. These mutations shift RAS to the GTP-loaded active state due to impairment in the intrinsic GTPase activity and disruption of GAP-mediated GTP hydrolysis, resulting in constitutive activation of effectors such as RAF. Because activation of RAF involves dimerization, RAS dimerization has been proposed as an important step in RAS-mediated activation of effectors. The α4-α5 allosteric lobe of RAS has been proposed as a RAS dimerization interface. Indeed, the NS1 monobody, which binds the α4-α5 region within the RAS G domain, inhibits RAS-dependent signaling and transformation as well as RAS nanoclustering at the plasma membrane. Although these results are consistent with a model in which the G domain dimerizes through the α4-α5 region, the isolated G domain of RAS lacks intrinsic dimerization capacity. Furthermore, prior studies analyzing α4-α5 point mutations have reported mixed effects on RAS function. Here, we evaluated the activity of a panel of single amino acid substitutions in the α4-α5 region implicated in RAS dimerization. We found that these proposed "dimerization-disrupting" mutations do not significantly impair self-association, signaling, or transformation of oncogenic RAS. These results are consistent with a model in which activated RAS protomers cluster in close proximity to promote the dimerization of their associated effector proteins (e.g., RAF) without physically associating into dimers mediated by specific molecular interactions. Our findings suggest the need for a nonconventional approach to developing therapeutics targeting the α4-α5 region.

Cryo-EM structures of cancer-specific helical and kinase domain mutations of PI3Kα

Phosphoinositide 3-kinases (PI3Ks) are a family of lipid kinases that perform multiple and important cellular functions. The protein investigated here belongs to class IA of the PI3Ks; it is a dimer consisting of a catalytic subunit, p110α, and a regulatory subunit, p85α, and is referred to as PI3Kα. The catalytic subunit p110α is frequently mutated in cancer. The mutations induce a gain of function and constitute a driving force in cancer development. About 80% of these mutations lead to single-amino-acid substitutions in one of three sites of p110α: two in the helical domain of the protein (E542K and E545K) and one at the C-terminus of the kinase domain (H1047R). Here, we report the cryo-electron microscopy structures of these mutants in complex with the p110α-specific inhibitor BYL-719. The H1047R mutant rotates its sidechain to a new position and weakens the kα11 activation loop interaction, thereby reducing the inhibitory effect of p85α on p110α. E542K and E545K completely abolish the tight interaction between the helical domain of p110α and the N-terminal SH2 domain of p85α and lead to the disruption of all p85α binding and a dramatic increase in flexibility of the adaptor-binding domain (ABD) in p110α. Yet, the dimerization of PI3Kα is preserved through the ABD-p85α interaction. The local and global structural features induced by these mutations provide molecular insights into the activation of PI3Kα, deepen our understanding of the oncogenic mechanism of this important signaling molecule, and may facilitate the development of mutant-specific inhibitors.

Deletion of the col-26 Transcription Factor Gene and a Point Mutation in the exo-1 F-Box Protein Gene Confer Sorbose Resistance in Neurospora crassa.

L-Sorbose induces hyperbranching of hyphae, which results in colonial growth in Neurospora crassa. The sor-4 gene, which encodes a glucose sensor that acts in carbon catabolite repression (CCR), has been identified as a sorbose resistance gene. In this study, we found that the deletion mutant of col-26, which encodes an AmyR-like transcription factor that acts in CCR, displayed sorbose resistance. In contrast, the deletion mutants of other CCR genes, such as a hexokinase (hxk-2), an AMP-activated S/T protein kinase (prk-10), and a transcription factor (cre-1), showed no sorbose resistance. Double mutant analysis revealed that the deletion of hxk-2, prk-10, and cre-1 did not affect the sorbose resistance of the col-26 mutant. Genes for a glucoamylase (gla-1), an invertase (inv), and glucose transporters (glt-1 and hgt-1) were highly expressed in the cre-1 mutant, even in glucose-rich conditions, but this upregulation was suppressed in the Δcre-1;Δcol-26a double-deletion mutant. Furthermore, we found that a dgr-2(L1)a mutant with a single amino-acid substitution, S11L, in the F-box protein EXO-1 displayed sorbose resistance, unlike the deletion mutants of exo-1, suggesting that the function of EXO-1 is crucial for the resistance. Our data strongly suggest that CCR directly participates in sorbose resistance, and that COL-26 and EXO-1 play important roles in regulating the amylase and glucose transporter genes during CCR.

Increasing protein stability by inferring substitution effects from high-throughput experiments

We apply a computational model, global multi-mutant analysis (GMMA), to inform on effects of most amino acid substitutions from a randomly mutated gene library. Using a high mutation frequency, the method can determine mutations that increase the stability of even very stable proteins for which conventional selection systems have reached their limit. As a demonstration of this, we screened a mutant library of a highly stable and computationally redesigned model protein using an invivo genetic sensor for folding and assigned a stability effect to 374 of 912 possible single amino acid substitutions. Combining the top 9 substitutions increased the unfolding energy 47 to 69kJ/mol in a single engineering step. Crystal structures of stabilized variants showed small perturbations in helices 1 and 2, which rendered them closer in structure to the redesign template. This case study illustrates the capability of the method, which is applicable to any screen for protein function.

Dimers of D76N-β2-microglobulin display potent antiamyloid aggregation activity

Self-association of WT β2-microglobulin (WT-β2m) into amyloid fibrils is associated with the disorder dialysis related amyloidosis. In the familial variant D76N-β2m, the single amino acid substitution enhances the aggregation propensity of the protein dramatically and gives rise to a disorder that is independent of renal dysfunction. Numerous biophysical and structural studies on WT- and D76N-β2m have been performed in order to better understand the structure and dynamics of the native proteins and their different potentials to aggregate into amyloid. However, the structural properties of transient D76N-β2m oligomers and their role(s) in assembly remained uncharted. Here, we have utilized NMR methods, combined with photo-induced crosslinking, to detect, trap, and structurally characterize transient dimers of D76N-β2m. We show that the crosslinked D76N-β2m dimers have different structures from those previously characterized for the on-pathway dimers of ΔN6-β2m and are unable to assemble into amyloid. Instead, the crosslinked D76N-β2m dimers are potent inhibitors of amyloid formation, preventing primary nucleation and elongation/secondary nucleation when added in substoichiometric amounts with D76N-β2m monomers. The results highlight the specificity of early protein-protein interactions in amyloid formation and show how mapping these interfaces can inform new strategies to inhibit amyloid assembly.

SARS-CoV-2 variant evasion of monoclonal antibodies based on in vitro studies.

Monoclonal antibodies (mAbs) offer a treatment option for individuals with severe COVID-19 and are especially important in high-risk individuals where vaccination is not an option. Given the importance of understanding the evolution of resistance to mAbs by SARS-CoV-2, we reviewed the available in vitro neutralization data for mAbs against live variants and viral constructs containing spike mutations of interest. Unfortunately, evasion of mAb-induced protection is being reported with new SARS-CoV-2 variants. The magnitude of neutralization reduction varied greatly among mAb-variant pairs. For example, sotrovimab retained its neutralization capacity against Omicron BA.1 but showed reduced efficacy against BA.2, BA.4 and BA.5, and BA.2.12.1. At present, only bebtelovimab has been reported to retain its efficacy against all SARS-CoV-2 variants considered here. Resistance to mAb neutralization was dominated by the action of epitope single amino acid substitutions in the spike protein. Although not all observed epitope mutations result in increased mAb evasion, amino acid substitutions at non-epitope positions and combinations of mutations also contribute to evasion of neutralization. This Review highlights the implications for the rational design of viral genomic surveillance and factors to consider for the development of novel mAb therapies.

A single amino acid substitution in the VP2 protein of Indian foot-and-mouth disease virus serotype O vaccine strain confers thermostability and protective immunity in cattle.

Foot-and-mouth disease (FMD) is a significant threat to animal health globally. Prophylactic vaccination using inactivated FMD virus (FMDV) antigen is being practised for the control in endemic countries. A major limitation of the current vaccine is its susceptibility to high environmental temperature causing loss of immunogenicity, thus necessitating the cold chain for maintenance of its efficacy. Hence, the FMD vaccine with thermostable virus particles will be highly useful in sustaining the integrity of whole virus particle (146S)during storage at 4°C. In this study, 12 recombinant mutants of Indian vaccine strain of FMDV serotype O (O/IND/R2/1975) were generated through reverse genetics approach and evaluated for thermostability. One of the mutant viruses, VP2_Y98F was more thermostable than other mutants and the parent FMDV. The oil-adjuvanted vaccine formulated with the inactivated VP2_Y98F mutant FMDV was stable up to 8 months when stored at 4°C and induced protective antibody response till dpv 180after primary vaccination. It is concluded that the VP2_Y98F mutant FMDV was thermostable and has the potential to replace the parent vaccine strain.

Deep mutational scanning and massively parallel kinetics of plasminogen activator inhibitor-1 functional stability to probe its latency transition

Plasminogen activator inhibitor-1 (PAI-1), a member of the serine protease inhibitor superfamily of proteins, is unique among serine protease inhibitors for exhibiting a spontaneous conformational change to a latent or inactive state. The functional half-life for this transition at physiologic temperature and pH is ∼1 to 2h. To better understand the molecular mechanisms underlying this transition, we now report on the analysis of a comprehensive PAI-1 variant library expressed on filamentous phage and selected for functional stability after 48h at 37 °C. Of the 7201 possible single amino acid substitutions in PAI-1, we identified 439 that increased the functional stability of PAI-1 beyond that of the WT protein. We also found 1549 single amino acid substitutions that retained inhibitory activity toward the canonical target protease of PAI-1 (urokinase-like plasminogen activator), whereas exhibiting functional stability less than or equal to that of WT PAI-1. Missense mutations that increase PAI-1 functional stability are concentrated in highly flexible regions within the PAI-1 structure. Finally, we developed a method for simultaneously measuring the functional half-lives of hundreds of PAI-1 variants in a multiplexed, massively parallel manner, quantifying the functional half-lives for 697 single missense variants of PAI-1 by this approach. Overall, these findings provide novel insight into the mechanisms underlying the latency transition of PAI-1 and provide a database for interpreting human PAI-1 genetic variants.

Comprehensive In Silico Functional Prediction Analysis of CDKL5 by Single Amino Acid Substitution in the Catalytic Domain.

Cyclin-dependent kinase-like 5 (CDKL5) is a serine/threonine protein kinase whose pathological mutations cause CDKL5 deficiency disorder. Most missense mutations are concentrated in the catalytic domain. Therefore, anticipating whether mutations in this region affect CDKL5 function is informative for clinical diagnosis. This study comprehensively predicted the pathogenicity of all 5700 missense substitutions in the catalytic domain of CDKL5 using in silico analysis and evaluating their accuracy. Each missense substitution was evaluated as “pathogenic” or “benign”. In silico tools PolyPhen-2 HumDiv mode/HumVar mode, PROVEAN, and SIFT were selected individually or in combination with one another to determine their performance using 36 previously reported mutations as a reference. Substitutions predicted as pathogenic were over 88.0% accurate using each of the three tools. The best performance score (accuracy, 97.2%; sensitivity, 100%; specificity, 66.7%; and Matthew’s correlation coefficient (MCC), 0.804) was achieved by combining PolyPhen-2 HumDiv, PolyPhen-2 HumVar, and PROVEAN. This provided comprehensive information that could accurately predict the pathogenicity of the disease, which might be used as an aid for clinical diagnosis.

Real-time dynamic single-molecule protein sequencing on an integrated semiconductor device.

Studies of the proteome would benefit greatly from methods to directly sequence and digitally quantify proteins and detect posttranslational modifications with single-molecule sensitivity. Here, we demonstrate single-molecule protein sequencing using a dynamic approach in which single peptides are probed in real time by a mixture of dye-labeled N-terminal amino acid recognizers and simultaneously cleaved by aminopeptidases. We annotate amino acids and identify the peptide sequence by measuring fluorescence intensity, lifetime, and binding kinetics on an integrated semiconductor chip. Our results demonstrate the kinetic principles that allow recognizers to identify multiple amino acids in an information-rich manner that enables discrimination of single amino acid substitutions and posttranslational modifications. With further development, we anticipate that this approach will offer a sensitive, scalable, and accessible platform for single-molecule proteomic studies and applications.

Loss of stability and unfolding cooperativity in hPGK1 upon gradual structural perturbation of its N-terminal domain hydrophobic core

Phosphoglycerate kinase has been a model for the stability, folding cooperativity and catalysis of a two-domain protein. The human isoform 1 (hPGK1) is associated with cancer development and rare genetic diseases that affect several of its features. To investigate how mutations affect hPGK1 folding landscape and interaction networks, we have introduced mutations at a buried site in the N-terminal domain (F25 mutants) that either created cavities (F25L, F25V, F25A), enhanced conformational entropy (F25G) or introduced structural strain (F25W) and evaluated their effects using biophysical experimental and theoretical methods. All F25 mutants folded well, but showed reduced unfolding cooperativity, kinetic stability and altered activation energetics according to the results from thermal and chemical denaturation analyses. These alterations correlated well with the structural perturbation caused by mutations in the N-terminal domain and the destabilization caused in the interdomain interface as revealed by H/D exchange under native conditions. Importantly, experimental and theoretical analyses showed that these effects are significant even when the perturbation is mild and local. Our approach will be useful to establish the molecular basis of hPGK1 genotype–phenotype correlations due to phosphorylation events and single amino acid substitutions associated with disease.