http://orcid.org/0000-0002-9456-3693
Figures
Abstract
A trimeric glycoprotein complex on the surface of human cytomegalovirus (HCMV) binds to platelet-derived growth factor (PDGF) receptor α (PDGFRα) to mediate host cell recognition and fusion of the viral and cellular membranes. Soluble PDGFRα potently neutralizes HCMV in tissue culture, and its potential use as an antiviral therapeutic has the benefit that any escape mutants will likely be attenuated. However, PDGFRα binds multiple PDGF ligands in the human body as part of developmental programs in embryogenesis and continuing through adulthood. Any therapies with soluble receptor therefore come with serious efficacy and safety concerns, especially for the treatment of congenital HCMV. Soluble virus receptors that are orthogonal to human biology might resolve these concerns. This engineering problem is solved by deep mutational scanning on the D2-D3 domains of PDGFRα to identify variants that maintain interactions with the HCMV glycoprotein trimer in the presence of competing PDGF ligands. Competition by PDGFs is conformation-dependent, whereas HCMV trimer binding is independent of proper D2-D3 conformation, and many mutations at the receptor-PDGF interface are suitable for functionally separating trimer from PDGF interactions. Purified soluble PDGFRα carrying a targeted mutation succeeded in displaying wild type affinity for HCMV trimer with a simultaneous loss of PDGF binding, and neutralizes trimer-only and trimer/pentamer-expressing HCMV strains infecting fibroblasts or epithelial cells. Overall, this work makes important progress in the realization of soluble HCMV receptors for clinical application.
Author summary
Human cytomegalovirus (HCMV) causes severe disease in transplant recipients and immunocompromised patients, and infections in a fetus or neonate are responsible for life-long neurological defects. Cell entry is in part mediated by a trimeric glycoprotein complex on the viral surface, which binds tightly to the host receptor PDGFRα. The soluble extracellular region of PDGFRα can be used as an antiviral agent to potently neutralize the virus in vitro. However, PDGFRα ordinarily binds growth factors in the human body to regulate developmental programs, which will limit the in vivo efficacy and safety of soluble PDGFRα. Using saturation mutagenesis and selections in human cell culture, mutations in PDGFRα are identified that eliminate off-target growth factor interactions while preserving HCMV binding and neutralization.
Citation: Park J, Gill KS, Aghajani AA, Heredia JD, Choi H, Oberstein A, et al. (2020) Engineered receptors for human cytomegalovirus that are orthogonal to normal human biology. PLoS Pathog 16(6): e1008647. https://doi.org/10.1371/journal.ppat.1008647
Editor: Brent Ryckman, University of Montana, UNITED STATES
Received: January 23, 2020; Accepted: May 20, 2020; Published: June 19, 2020
Copyright: © 2020 Park et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Raw and processed deep sequencing data are deposited in NCBI’s Gene Expression Omnibus (GEO) under series accession number GSE138169.
Funding: The development of deep mutagenesis to study virus-receptor interactions was supported by NIH award R01AI129719 to E.P. The development of PDGF signaling assays was supported by a seed grant from Illini 4000 and the Cancer Center at Illinois to E.P. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: E.P. is the inventor on a provisional patent filing by the University of Illinois covering aspects of this work. E.P. is a co-founder of Orthogonal Biologics Inc, which has licensed the intellectual property for preclinical development.
Introduction
Human cytomegalovirus (HCMV; or human herpesvirus 5, HHV-5) is a ubiquitous pathogen that has infected most of the global adult population, and typically causes asymptomatic infections [1,2]. In some cases, individuals will experience mononucleosis-like symptoms, and severe disease can occur in immunocompromised or immunosuppressed adults [3]. HCMV infection may also have an oncomodulatory effect and could be associated with tumor progression, most notably in glioblastoma [4–7], although the mechanistic details remain controversial [8]. However, it is amongst pregnant women and newborns that HCMV has a disproportionate impact on public health [9]. Just over 1 in 200 infants born in the United States are congenitally infected with HCMV, and approximately a fifth of these infants will suffer life-long neurological complications [10,11], including hearing and vision loss, seizures, behavioral disorders, and developmental delays [12–15]. Ideally, infection during pregnancy would be avoided, but this is challenging. The high prevalence of HCMV (especially in families with young children [16,17]), easy transmission through contact with bodily fluids, recurrence of latent endogenous infections [18], reinfection by different exogenous strains [19], and difficulties recognizing asymptomatic infected individuals make controlling virus spread difficult. Furthermore, antivirals have not been approved by the U.S. Food and Drug Administration for routine use in the treatment of congenital HCMV where there are unique concerns regarding toxicity, and drug resistance is widely reported [20]. Faced with this reality, the American College of Obstetricians and Gynecologists does not make any specific recommendations for counseling or treating pregnant women for the prevention of HCMV [21].
It is in this environment that the discovery of an important HCMV receptor, the platelet-derived growth factor receptor α (PDGFRα), has generated excitement towards the development of new therapeutics inhibiting virus-host cell attachment and entry. Two glycoprotein complexes on the HCMV surface drive broad host cell tropism [22]. A pentameric complex of glycoprotein H (gH; UL75), gL (UL115), UL128, UL130 and UL131A mediates entry into epithelial and endothelial cells [23–28], possibly through engagement of neuropilin-2 [29] or the olfactory receptor OR14I1 [30] on the host cell. A trimeric complex of gH, gL and gO (UL74) is sufficient for entry into fibroblasts and is also necessary for entry into pentamer-requiring cell types by promoting membrane fusion [31–33]. Multiple lines of evidence indicate that PDGFRα is the receptor for the trimer. Decreased PDGFRα expression reduces HCMV entry into fibroblasts [34–37]; gO can specifically bind PDGFRα-expressing cells [36]; a quaternary complex of soluble extracellular regions (gH-gL-gO-PDGFRα) can be assembled in vitro and the components bind with high nanomolar affinity [29,34]; HCMV strains lacking the trimer fail to enter PDGFRα-dependent cell types to infect cells [36,38]; and the soluble PDGFRα extracellular domain blocks virus entry [29,36,38]. This has led to the soluble domain of PDGFRα or derivative fragments being explored as antivirals [34,36,38], a promising strategy as any mutations in HCMV to escape receptor-based inhibitors would likely decrease binding to the natural receptor and attenuate virulence. Soluble PDGFRα ectodomain effectively blocks virus entry whether applied pre- or post-cell attachment, due to inhibition of both early attachment and fusion steps [38]. Furthermore, soluble PDGFRα inhibits cell entry when only a fraction of glycoprotein trimer is bound and therefore does not require high saturating concentrations; and soluble PDGFRα blocks entry by multiple HCMV strains that express both trimer and pentamer complexes [38].
While appealing, the use of soluble PDGFRα to treat HCMV in humans comes with risk. Indeed, no protein-based viral receptors are currently approved as antiviral drugs by the FDA [39] (although at least one is in a phase II clinical trial for the treatment of COVID-19 by Apeiron Biologics), in part due to safety concerns over interactions with endogenous ligands. By comparison, there are multiple antibody drugs that are specific for viral targets used in the clinic [39]. PDGFRα is an important receptor for platelet-derived growth factors (PDGFs), and regulates cellular proliferation, differentiation and development of multiple tissues during embryogenesis and onwards through adulthood [40]. The PDGF ligands are small polypeptides that covalently associate in disulfide-bonded homo- and heterodimers, which display differing activities towards PDGFRα and its close relative PDGFRβ [41]. Four ligands interact with PDGFRα with high affinity (PDGF-AA, AB, BB, and CC) [42,43] and bind two receptor chains at opposing ends to form a receptor-(ligand)2-receptor signaling complex [44,45]. PDGFs compete with the HCMV trimer for PDGFRα binding and block infection [35,36], suggesting that PDGF binding to PDGFRα sterically hinders trimer binding or allosterically modulates the trimer binding interface in a non-productive manner. Treatment with soluble PDGFRα will almost certainly disrupt physiological PDGF-signaling, and the ability of recombinant receptor to bind virions may be diminished by competing endogenous PDGFs. Ideally, mutations in the receptor can be identified that disrupt PDGF interactions while retaining the site or sites engaged by HCMV.
A high resolution crystal structure of PDGF-BB-bound PDGFRβ has illuminated atomic details of ligand-receptor interactions in the family, revealing a predominantly hydrophobic interface at the cleft between the second (D2) and third (D3) extracellular domains of the receptor [44]. However, the binding site for gO is ambiguous. Cryo-electron microscopy suffers from conformational diversity and low resolution [34]; PDGFRα peptide fragment analysis failed to discover any one peptide solely responsible for high affinity binding [38]; and while domain deletions have shown PDGFRα-D3 is important [37], such studies have coarse resolution. The absence of detailed structural information on the trimer-PDGFRα complex means the engineering of HCMV-specific receptors orthogonal to human biology (i.e. receptors that bind virus but do not interact with endogenous human ligands) remains an unsolved challenge.
This challenge is solved through the use of mutagenesis and in vitro selection. By tracking the enrichment or depletion of sequence variants in a diverse library using next generation sequencing, the relative phenotypes of thousands of mutations can be simultaneously assessed in a single experiment, referred to as deep mutagenesis [46]. Deep mutational scans have been used to address and engineer specificity in proteins that promiscuously bind multiple ligands [47–49], and methods for deep mutagenesis of membrane proteins expressed in human cells have been optimized [50–52]. By screening through all possible single amino acid substitutions in the D2-D3 domains, we discovered that, unlike PDGF binding, HCMV trimer interactions persist when PDGFRα conformation is disrupted by mutations within the folded core. Particularly relevant to therapeutic design, we also identified PDGFRα surface mutations that retain trimer binding in the presence of high PDGF concentrations. This study brings receptor-based inhibition of HCMV infection closer to realization.
Results
A deep mutational scan of PDGFRα D2-D3 domains for trimer binding in the presence of PDGFs
The glycoprotein trimer from HCMV strain TB40/E (BAC4 clone, [53]) was expressed in Expi293F cells, a suspension derivative of human HEK 293. The transmembrane helix of gH was deleted to produce a soluble complex, and the gO subunit, which directly contacts PDGFRα [34,36,54], was fused at its C-terminus to superfolder GFP (sfGFP; [55]) for fluorescence detection. Medium from gH/gL/gO-sfGFP expressing cells was incubated with Expi293F cells transfected with an episomal plasmid encoding PDGFRα. Trimer binding was proportional to PDGFRα surface expression (determined by flow cytometry using antibody staining of an extracellular c-myc epitope tag at the PDGFRα N-terminus; Fig 1A and 1B), and was inhibited by co-incubation with an equimolar mixture of PDGF-AA, AB, BB, and CC (Fig 1D). Medium from cells expressing gO-sfGFP alone, either the wild-type sequence or carrying the mutation C351S (C343S based on numbering of gO from the TB40/E strain) to remove the exposed cysteine that would ordinarily form a disulfide to gL-C144 [56], failed to display high binding signal to PDGFRα positive cells (S1 Fig). The binding signal is therefore dependent on co-expression with gH and gL, suggestive that it is trimer as opposed to monomeric gO that engages the receptor, although this is not explicitly demonstrated.
(A) Expi293F cells expressing myc-tagged PDGFRα were incubated with dilutions of gH/gL/gO-sfGFP-containing medium, and binding (shown as ΔMean Fluorescence Units, ΔMFU, in which Mean Fluorescence Units of a negative control without trimer are subtracted as background) was measured by flow cytometry. Neat, no dilution. (B and C) Cells were transfected with wild-type (WT) PDGFRα (B) or the D2-D3 SSM library (C) under conditions where typically no more than a single sequence variant is expressed by any given cell. Under these conditions, most cells are negative for PDGFRα expression as determined by anti-myc staining. Cells were incubated with a subsaturating (12-fold dilution) of gH/gL/gO-sfGFP medium and analyzed by flow cytometry. (D) Cells were co-incubated with a 2-fold dilution of gH/gL/gO-sfGFP medium and different concentrations of PDGF (an equimolar mixture of PDGF-AA, AB, BB, and CC). The fluorescent signal of gH/gL/gO-sfGFP binding is compared to a sample without competing PDGF. (E and F) Cells expressing WT PDGFRα (E) or the D2-D3 SSM library (F) were incubated with a 3-fold dilution of gH/gL/gO-sfGFP medium and 100 nM PDGF (25 nM of each subtype). The choice of a PDGF concentration that does not completely inhibit trimer binding ensured (i) that there remained some positive binding signal for drawing appropriate sort gates, and (ii) that mutations with weak effects will be distinguished from the main population. During sorting of the D2-D3 library, the top 15% of PDGFRα-expressing cells that bind trimer were collected (UL74-High sort, green gate). Simultaneously, the bottom 15% of PDGFRα-expressing cells that bind trimer were also collected (UL74-Low sort, blue gate). For cells expressing WT PDGFRα, the equivalent gates correspond to 4% and 23% of PDGFRα-expressing cells, respectively.
A single site-saturation mutagenesis (SSM) library of PDGFRα was constructed encompassing all single amino acid mutations in the D2-D3 domains (a.a. D123-E311), and transfected into Expi293F cells under conditions that typically yield no more than one sequence variant per cell, providing a tight link between genotype and phenotype [51]. When incubated with gH/gL/gO-sfGFP medium at subsaturating dilutions, cells expressing the SSM library were surprisingly indistinguishable from cells expressing wild-type PDGFRα (Fig 1B and 1C). Normally, many mutations adversely impact protein activity through destabilization of folded structure or damage to functional sites, and loss-of-function variants tend to dominate naive libraries prior to any selection. That deleterious mutations are not prevalent in the naive SSM library immediately implied that HCMV trimer binding is resistant to most single non-synonymous mutations in the PDGFRα D2-D3 domains.
By comparison, many cells expressing the SSM library displayed higher trimer binding in the presence of PDGFs than cells expressing wild-type PDGFRα (Fig 1E and 1F). There thus appeared to be many PDGFRα mutations that selectively lost PDGF affinity. To identify these mutations, PDGFRα-expressing cells that had elevated HCMV trimer binding in the presence of competing PDGFs were collected by fluorescence-activated cell sorting (FACS). This is referred to as the UL74-High sort (see green gate in Fig 1F). Within the same experiment, cells expressing PDGFRα but displaying low trimer binding were also collected, referred to as the UL74-Low sort (see blue gate in Fig 1F). The UL74 gene encodes gO [31], and the names assigned to the sorted populations correspond to the raw and processed data files deposited in the GEO database [57]. PDGFRα mutants that lose affinity for trimer, or perhaps have enhanced affinity for competing PDGFs, will be preferentially enriched in the UL74-Low sort. PDGFRα mutants that fail to express will be depleted from both sorted populations. Following Illumina sequencing of the naive plasmid library and transcripts from the sorted populations, the enrichment ratios for all 3,780 substitutions in the D2-D3 domains were calculated to define a local mutational landscape (S2A Fig). Data from independent replicates are highly correlated (S2B–S2G Fig), giving confidence in the data’s accuracy.
PDGF interactions are uniquely sensitive to mutational disruption of folded structure in the D2-D3 domains
Enrichment ratios for nonsynonymous mutations are anticorrelated between the UL74-High and UL74-Low sorts, with few mutations other than premature stop codons being depleted in both selections. Hence substitution mutations in PDGFRα D2-D3 domains suffer little cost to surface localization, even though many mutations will almost certainly destabilize the domain fold. This differs markedly from a prior mutational scan of a multidomain membrane protein expressed in human cells, where nonconservative mutations buried within folded structure prevented protein trafficking to the plasma membrane due to intracellular retention by quality control machinery [50]. When mapping experimental conservation scores to the structure of PDGF-bound PDGFRα (modeled from the crystal structure of PDGF-BB-bound PDGFRβ; PDB ID 3MJG [44]), enriched mutations in the UL74-High sort are found to cluster to buried core positions and the PDGF interface (Fig 2). These include substitutions that create cavities, steric clashes, and/or introduce buried ionizable groups; these mutations are incompatible with biophysical principles governing protein folding and will destabilize folded structure. We hypothesize that PDGFRα mutants with disrupted structure have therefore preferentially lost affinity to PDGFs, thereby reducing competition and enhancing binding of the HCMV trimer. The epitope bound by the HCMV trimer is ambiguous from the data (S3 Fig), but must be at least partially independent of proper D2-D3 conformation, although partial structure may remain. This is consistent with either key contacts to the HCMV trimer residing on linear PDGFRα segments and loops, similar to how many antibodies recognize conformation-independent epitopes, or with favorable interactions to the HCMV trimer being distributed over a broad surface (possibly even beyond the D2-D3 domains) such that localized disruptions are tolerated. Libraries encompassing double or higher order PDGFRα mutations may further resolve details of the binding mechanism. While this illuminates aspects of PDGFRα/trimer recognition, using misfolded, destabilized PDGFRα variants to selectively neutralize HCMV is fraught with risks, including toxicity from non-specific interactions and protein aggregation due to the exposure of hydrophobic residues.
(A) Model of PDGFRα D2-D3 domains bound to PDGF-BB (protruding subunit B1 is pale green, receding subunit B2 is pale blue). Conservation scores from the UL74-High sort are mapped to the PDGFRα surface, showing residues where mutations tend to be depleted in orange, and residues where mutations tend to be enriched in dark blue. (B) A cross section through PDGFRα highlights that enriched mutations (dark blue) are concentrated in the protein interior/core. (C-E) Log2 enrichment ratios for individual mutations from the UL74-High sort are plotted from -3 (orange, depleted) to +3 (dark blue, enriched). Amino acid substitutions are indicated on the vertical axis, * is stop codon. The wild-type amino acid is black. Structural views are shown at right, colored as in panel A. Regions shown are the D2 core (C), PDGF-BB interface (D), and PDGF-BB interface and underlying D3 core (E).
Surface mutations at the PDGF binding site favor PDGFRα specificity towards the HCMV trimer
The PDGFRα mutational landscape for high trimer binding in the presence of competing PDGF reveals key surface residues as hot spots for enriched mutations, especially PDGFRα residues L137, L208 and V242 that contact the protruding PDGF subunit, and Y206 which packs between the protruding and receding PDGF subunits (Fig 2D and 2E). Mutations to these residues will disrupt the PDGFRα-PDGF interface, especially by the addition of polar substitutions that invert the chemical properties of the highly hydrophobic binding site, while HCMV trimer interactions persist. Overall, enrichment of mutations for increased HCMV trimer binding in the presence of competing PDGFs is highly correlated to the mutated residue’s connectivity within the D2-D3 core or across the PDGFRα-PDGF interface (S4 Fig), emphasizing the selective importance of PDGFRα conformation and surface contacts in the D2-D3 cleft for high affinity PDGF interactions.
Five mutations on the PDGFRα surface (L137K, L137Q and Y206S in D2, and V242K and V242T in D3, chosen based on their high enrichment in the UL74-High sort and covering multiple surface positions) were validated as having desirable binding properties by targeted mutagenesis. As predicted from the deep mutational scan, these PDGFRα variants maintain near wild-type levels of binding to soluble HCMV trimer from the TB40/E strain, yet have diminished sensitivity to the addition of competing PDGFs (Fig 3A and 3B), demonstrating successful engineering of specificity towards the viral target. Since the mechanism of trimer binding to receptor has not been resolved at atomic or even residue-level resolution, it remains uncertain whether the viral genome could mutate to distinguish wild-type from engineered PDGFRα. This concern is partially alleviated by the observation that the engineered PDGFRα mutants maintain high binding to soluble HCMV trimer from the Merlin clinical strain [58] (Fig 3C), which has relatively high natural sequence variation compared to TB40/E (79% identity between glycoproteins O). Variation is particularly high in the N-terminus of glycoprotein O where important interactions to PDGFRα reside [54].
(A) Expi293F cells expressing wild-type PDGFRα were co-incubated without (black) or with (red) competing PDGFs, and binding of soluble TB40/E trimer (1/3 dilution of medium containing gH/gL/gO-sfGFP) was assessed by flow cytometry. Cells were gated (grey box) to compare binding (measured as ΔMFU) at equivalent receptor expression levels. (B) Binding of soluble TB40/E trimer to PDGFRα variants. Cells were analyzed as outlined in panel A. Data are mean ± SD, n = 3. (C) Cells expressing PDGFRα variants were incubated with Merlin gH/gL/gO-sfGFP containing medium (1/3 dilution) in the absence (grey) or presence of PDGFs (red). Cells were analyzed as outlined in panel A. Data are mean ± SD, n = 4. (D) Expi293F cells were transfected with gL, gO, and full-length gH from strain TB40/E. Binding of sPDGFR-Fc to HCMV trimer anchored at the cell surface was measured by flow cytometry. Data are mean ± SD, n = 2 (sPDGFRβ-Fc), 3 (WT-Legacy and Y206S), 4 (V242K), or 5 (WT). (E) Binding of sPDGFR-Fc to cells expressing membrane-anchored trimer from HCMV strain Merlin. Data are mean ± SD, n = 2 (sPDGFRβ-Fc), or 3 (WT-Legacy, WT, Y206S, and V242K).
Fusion of soluble proteins to immunoglobulin Fc confers multiple desirable properties for a therapeutic [59]. Multimerized Fc chains impose avidity for enhanced apparent affinity, and Fc moieties are engaged by multiple receptors to enhance serum half-life or evoke effector functions, including complement activation and antibody-dependent cell-mediated cytotoxicity [59]. Soluble PDGFRα fused to the Fc region of human IgG1 was shown to bind HCMV trimer when first identified as its receptor [34], and the neutralization properties of Fc-fused soluble PDGFRα (sPDGFRα-Fc) have been explored since in greater detail [38]. These prior studies used commercially supplied sPDGFRα-Fc featuring a random linker of mixed amino acids that connects to the upper hinge of IgG1, upstream of Cys-220 that would ordinarily form a disulfide to the antibody light chain (S5A Fig). The nature of this Fc fusion may cause manufacturing liabilities if Cys-220 is exposed and free, and we refer to this construct as “legacy” sPDGFRα-Fc. We designed an alternative sPDGFRα-Fc construct that includes the PDGFRα signal peptide (a.a. 1–23), ectodomain (a.a. 24–524), a GGGS linker, and finally human IgG1 Fc beginning at Asp-221 (S5A Fig). Both sPDGFRα-Fc constructs bind membrane-localized HCMV trimer with equal affinity (Fig 3D and 3E), and subsequent studies focused exclusively on our redesigned fusion protein.
Soluble PDGFRα-Fc was expressed in Expi293F culture, and was isolated by affinity capture and size exclusion chromatography (SEC) to very high purity (S5B and S5C Fig). The purified protein is monodisperse by SEC, and the final yield is high (> 25 mg / L of transfected culture without any optimization). Soluble PDGFRα-Fc V242K was stressed for 7 days at 40°C, pH 8.5, to promote asparagine deamidation and protein denaturation [60], yet remained mostly stable based on the SEC profile (S6 Fig). A harsher stress for 14 days at 40°C, pH 5.5, to promote asparagine isomerization and protein denaturation [60] caused sPDGFRα-Fc V242K to form soluble aggregates (S6 Fig). We conclude that the protein has moderate stability and is suitable for further manufacturing development.
Two PDGFRα mutations, Y206S and V242K, were explored as Fc-fused soluble orthogonal receptors. Binding of the soluble receptors was measured towards HCMV trimer expressed at the cell surface, with gL and gO subunits tagged at their C-termini with short peptide epitopes, and gH expressed as full-length protein with its native transmembrane domain. The Y206S mutant had reduced binding to HCMV trimer compared to wild-type receptor, but sPDGFRα-Fc V242K bound HCMV trimer from Merlin and TB40/E strains with equal low nanomolar affinity (Fig 3D and 3E). Aspects of these experiments were replicated with independent protein preparations (S7A Fig). Furthermore, fusions of wild-type and V242K sPDGFRα with the Fc region of IgG3 also bound similarly to trimer from the TB40/E and Merlin strains (S7B and S7C Fig). Compared to IgG1, IgG3 has lower affinity for the neonatal Fc receptor (FcRn), which is associated with reduced serum half-life and placental transfer [61,62]. We speculate that these features may be advantageous in the treatment of pregnant women during acute HCMV infection to achieve a lower dose in the developing fetus, thereby further addressing safety concerns.
To further validate that the engineered receptors have lost growth factor interactions, we assessed their ability to inhibit PDGF signaling. Expi293F cells were transfected with a reporter plasmid encoding wild-type PDGFRα and a destabilized fluorescent protein (EGFP-PEST) [63] under the control of a serum response element (SRE) [64]. Treatment with PDGFs upregulates EGFP-PEST expression, which is inhibited by wild-type sPDGFRα-Fc acting as a decoy receptor (Fig 4). As intended, both sPDGFRα-Fc Y206S and V242K failed to inhibit PDGF signaling (Fig 4B). However, we also noticed that the engineered receptors, in particular sPDGFRα-Fc Y206S, unexpectedly enhanced the response to PDGF-A ligands. We suspect this is due to a generic ‘carrier protein’ effect (PDGF ligands and their receptors are highly hydrophobic [44] and are generally formulated with a carrier such as serum albumin), though we are unable to conclude this with certainty, and non-specific gamma globulin has no such effect (S8 Fig). This observation may require future investigation.
(A) Flow cytometry analysis of mock-transfected Expi293F cells (black), or cells transfected with a plasmid driving myc-tagged PDGFRα expression and carrying a SRE-regulated GFP reporter (red). Signaling activity was quantified by measuring mean GFP fluorescence (vertical axis) of responsive cells in the grey gate. Background PDGFRα signaling is at left, PDGF-BB stimulated activity is center, while sPDGFRα-Fc is shown to inhibit PDGF-BB stimulation at right. (B) Stimulation of transfected Expi293F cells with 5 nM PDGF ligands. PDGF ligands were pre-incubated with 0.2% BSA (negative, blue) or with a 3-fold molar excess of wild-type sPDGFRα-Fc (black), sPDGFRα-Fc Y206S (grey), or sPDGFRα-Fc V242K (white). Data are mean ± SD, n = 3 independent replicates.
Orthogonal PDGFRα-based receptors targeting HCMV trimer potently neutralize virus
The two representative orthogonal receptors, PDGFRα Y206S and V242K, were evaluated for their efficacy to neutralize virus infection. HCMV trimer-only lab strain AD169 [25,65] and trimer/pentamer-expressing clinical strain TB40/E [53] were grown on PDGFRα-positive MRC-5 fibroblasts and PDGFRα-negative ARPE-19 epithelial cells, respectively. Three cell lines were infected by both HCMV strains: MRC-5, ARPE-19, and ARPE-19RA which are stably transfected with PDGFRα to confer susceptibility to AD169 [37].
Soluble PDGFRα-Fc potently neutralized both HCMV strains at nanomolar concentrations (Fig 5). Neutralization by orthogonal sPDGFRα-Fc V242K was mostly indistinguishable from wild-type receptor and nearly complete at 1 to 10 nM (Fig 5), demonstrating that virus binding and neutralization can indeed be successfully separated from growth factor interactions by just a single amino acid substitution. Soluble PDGFRα-Fc Y206S also neutralized but with reduced potency, consistent with the biochemical binding studies. However, neutralization of trimer/pentamer-expressing TB40/E remained incomplete even at high sPDGFRα-Fc concentrations. Hyperimmune globulin (sold under the trade name Cytogam by CSL Behring for prophylactic treatment of transplant recipients) completely neutralizes HCMV, albeit at high concentrations. It is unclear if residual pentamer-mediated infection in the presence of high sPDGFRα-Fc would in any way limit therapeutic or prophylactic efficacy in an in vivo setting. Depletion of HCMV-positive patient sera with either trimer or pentamer severely diminishes neutralization of epithelial cell infection [66], suggesting that targeting both glycoprotein complexes can give synergistic effects and is worth further exploration.
MRC-5 (left), ARPE-19 (middle) and ARPE-19RA (right) cells were infected with HCMV strains AD169 (grown in MRC-5; top) and TB40/E (grown in ARPE-19; bottom). Infection by the trimer-only AD169 strain is restricted to PDGFRα-positive MRC-5 and ARPE-19RA lines. Viral inoculum was pre-incubated with sPDGFRα-Fc WT (black), sPDGFRα-Fc Y206S (grey), sPDGFRα-Fc V242K (white), sPDGFRβ-Fc (at a single concentration of 100 nM; red), or hyperimmune globulin (HIG; orange) prior to adding to cells. Data represent mean ± SD, n = 2 independent titrations.
Discussion
Ideas surrounding the use of soluble virus receptors as decoys to treat or prevent infection have a long history [67–70]. Soluble PDGFRα as a prophylactic or treatment for HCMV is particularly promising due to exceptionally tight affinity for glycoprotein trimer and its ability to neutralize both trimer- and pentamer-mediated host cell entry. Here, we have shown how deep mutagenesis can inform the engineering of PDGFRα variants that maintain tight virus binding but have lost growth factor interactions, and are thereby orthogonal to normal human signaling pathways. This is anticipated to improve both efficacy and safety in vivo.
There are no simple, tractable animal models for investigating HCMV infection. Two highly relevant models are the study of rhCMV in rhesus macaques [71] or the engraftment of human tissues into immune-deficient mice [72–74], but these are non-trivial and require specialized expertise. However, it might be possible to test biocompatibility and pharmacokinetics of the orthogonal receptors independent of infection in small animal models. For instance, murine and human PDGF ligands and receptors share high identity (> 85%, with all key binding site residues of PDGFRα conserved), and one standard assessment of purified human PDGF ligands is to measure mitogenic activity towards mouse fibroblast lines [75,76]. This suggests the human and mouse PDGF ligands and receptors will cross-react, though this has not been quantitatively characterized. Mice, especially strains expressing human FcRn for human-like serum stability of Fc fused protein [77,78], may therefore be a suitable model for assessing safety during future development.
The binding site for HCMV gO on PDGFRα is uncertain from the deep mutational scan, however the data unambiguously shows that binding is independent of properly folded structure in the PDGFRα D2-D3 domains. Deletion studies have show that the D3 domain is essential for trimer-mediated HCMV entry [37], although no specific hot spot residues for mediating physical binding have been defined. It is possible that gO makes critical contacts to domains D4 or D5 that are outside the mutagenized region in the PDGFRα D2-D3 library, or simply that negative effects of single amino acid substitutions in the D2-D3 domains are too small to cause an observable decrease in binding to HCMV trimer by flow cytometry. Determining how gO binds PDGFRα at atomic resolution will be important for speculating on the likelihood of HCMV acquiring mutations to distinguish and escape from orthogonal receptors.
In our own prior studies, deep mutagenesis of human membrane proteins was used to understand how sequence relates to folding based on surface localization, as misfolded sequences were retained intracellularly due to quality control machinery [50,51,79]. However, PDGFRα is an exception, and nearly every missense mutation still reaches the plasma membrane, including those that are incompatible with the native folded structure (for example, mutations that disrupt disulfides, create cavities, impose steric clashes, or introduce buried charges). The structural organization of PDGFRα as a series of loosely connected β-sandwich domains may facilitate surface trafficking of folding mutants, whereas previous studies have focused on membrane proteins with complex tertiary and quaternary structures.
Finally, deep mutagenesis has been extensively applied to viral surface proteins based on virus replication and propagation in tissue culture as the selection, often in the presence of neutralizing antibodies to discover principles of immune recognition and escape [80,81]. However, it is by decoupling the sequence library from a virus genome [52] that we have been able to determine a mutational landscape of a virus receptor. The principles and approach outlined here can be easily adapted for the engineering of other viral receptors for improved efficacy and safety to treat or prevent infections. Indeed, the same methods we optimized for the development of a HCMV antiviral were rapidly repurposed for engineering soluble decoy receptors that bind with high affinity to the spike of SARS coronavirus 2, a zoonotic pathogen that recently spilled over to the human population and has profoundly impacted world events [82]. Receptor-based antiviral drugs may be especially desirable for emerging zoonotic viruses or viruses with high sequence diversity, as they potentially have broader specificity than monoclonal antibodies to target all virus strains that utilize a common receptor.
Materials and methods
Plasmids
Human PDGFRα isoform 1 (GenBank NM_006206.4) was cloned in to the NheI-XhoI sites of pCEP4 (Invitrogen) with a N-terminal HA leader (MKTIIALSYIFCLVFA), myc-tag, linker (GSPGGASG), and followed by the mature polypeptide (a.a. 24–1,089). For measuring signaling activity, a minimal promoter under the control of tandem serum response elements was subcloned from SRE reporter vector_559 (Addgene # 82686) [64], a GFP-PEST fluorescent reporter [63] was inserted at the AscI site, and the entire reporter cassette was inserted at the NruI site of pCEP4-myc-PDGFRα by Gibson assembly. As a negative control, human PDGFRβ isoform 1 (GenBank NM_002609.3) was similarly cloned between the NheI-XhoI sites of pCEP4 with a N-terminal HA leader, myc-tag, and linker that connects to the mature protein (a.a. 33–1,106). No interactions between PDGFRβ and HCMV trimer were observed. Soluble PDGFRα-Fc was cloned in to the NheI-XhoI sites of pcDNA3.1(+) (Invitrogen), and encompassed PDGFRα a.a. 1–524, a short connecting linker as outlined in S5A Fig, and the C-terminus of human IgG1 (GenBank KY432415.1) beginning at either C220 or D221 as described in Results. Alternatively, PDGFRα a.a. 1–524 were fused via linker GGGS to D221-K447 of human IgG3 (GenBank P01860.2; this is an R435 allele with reduced binding to FcRn). Synthetic human codon-optimized gene fragments (Integrated DNA Technologies) for HCMV gO (GenBank ABV71596.1 for TB40/E strain, AJY56739.1 for Merlin strain) were genetically fused at the C-terminus to superfolder GFP [55] (for detection of soluble HCMV trimer binding to PDGFRα-positive cells) or to an 8-histidine tag (for expression of membrane-tethered HCMV trimer), and were ligated in to the NheI-XhoI sites of pcDNA3.1(+). Codon-optimized gene fragments for full-length HCMV gH (GenBank ABV71597.1 for TB40/E strain, YP_081523.1 for Merlin Strain) were cloned in to the NheI-XhoI sites of pcDNA3.1(+) for expression of membrane-tethered HCMV trimer. For production of soluble trimer, the extracellular region of gH with the leader peptide (a.a. 1–717 for TB40/E, a.a. 1–716 for Merlin) was subcloned with a 8-histidine tag. Codon-optimized synthetic genes for gL (GenBank ABV71629.1 for TB40/E strain, YP_081555.1 for Merlin) were cloned with C-terminal FLAG tags in to the NheI-XhoI sites of pcDNA3.1(+). Targeted mutations were made by overlap extension PCR. All plasmids were sequence verified (ACGT, Inc) and are deposited with Addgene.
Cells and viruses
Expi293F cells (ThermoFisher) are a suspension culture derivative of HEK293, and were cultured in Expi293 Expression Medium (ThermoFisher) at 125 rpm, 8% CO2, 37°C. MRC-5 embryonic lung fibroblasts and ARPE-19 adult retinal pigment epithelial cells were from the American Type Culture Collection and were grown at 37°C, 5% CO2, in DMEM supplemented with 10% fetal bovine serum, 1 mM sodium pyruvate, 2 mM glutamax (Gibco), 10 mM Hepes pH 7.4, 0.1 mM MEM Non-Essential Amino Acids (Gibco), 100 units/ml penicillin G, and 100 μg/ml streptomycin sulfate. ARPE-19 ectopically expressing PDGFRα (“ARPE19-RA” cells) were generated by lentiviral transduction with pLV-EF1a-PDGFRα-IRES-PURO, generously provided by Kai Wu (Princeton University). HCMV virus stocks were prepared by electroporating purified bacmid DNA into either MRC-5 fibroblasts (AD169; [65]) or ARPE-19 retinal pigment epithelial cells (TB40/E, clone BAC4; [53]). Stocks were amplified twice by infecting cell monolayers in 15 cm tissue culture plates followed by 850 cm2 roller bottles. Virions were concentrated 100× by centrifugation through a 20% sorbitol cushion, resuspended in phosphate-buffered saline (PBS) containing 1% BSA and 7% sucrose, and frozen at -80°C. Virus stocks were titered on MRC-5 and ARPE-19 cells by Immediate-Early Protein 1 (IE1) fluorescent focus assay [83].
Flow cytometry analysis of soluble HCMV trimer binding to receptor
Soluble HCMV trimer was prepared by transfecting Expi293F cells at a density of 2 × 106 / ml with 400 ng pcDNA3-gO-sfGFP, 400 ng pcDNA3-gH-8his and 400 ng pcDNA3-gL-FLAG per ml culture using Expifectamine (ThermoFisher). Transfection Enhancers (ThermoFisher) were added 18 h later, and 4 days post-transfection the culture was centrifuged (1,000 × g, 10 minutes, followed by a second spin of the supernatant at 21,000 × g, 5 minutes). The medium supernatant was stored at -20°C and used directly in binding experiments. Expi293F cells at 2 × 106 / ml were transfected with plasmids encoding receptors (500 ng DNA per ml of culture) using Expifectamine. 24 h post-transfection the cells were washed with PBS supplemented with 0.2% bovine serum albumin (PBS-BSA), incubated for 30 minutes on ice with anti-myc Alexa 647 (clone 9B11, 1/250 dilution; Cell Signaling Technology) and the indicated dilutions of gH/gL/gO-sfGFP-containing medium, washed twice with PBS-BSA, then analyzed on a BD LSR II flow cytometer. For competition assays, washed cells were pre-incubated with PDGF ligands (R&D Systems) for 2 minutes on ice prior to addition of gH/gL/gO-sfGFP-containing medium.
Deep mutagenesis
The D2-D3 domains (a.a. D123-E311) within plasmid pCEP4-myc-PDGFRα were mutagenized by overlap extension PCR [84] using primers with degenerate NNK codons. The plasmid library was transfected in to Expi293F cells using Expifectamine under conditions previously shown to typically give no more than a single coding variant per cell; 1 ng coding plasmid was diluted with 1,500 ng pCEP4-ΔCMV carrier plasmid per ml of cell culture at 2 × 106 / ml, and the medium was replaced 2 h post-transfection. The cells were collected after 24 h, washed with ice-cold PBS-BSA, and incubated for 2 minutes on ice with PDGF-AA, -AB, -BB, and–CC (25 nM of each after addition of soluble trimer-containing medium) and anti-myc Alexa 647 (clone 9B11, 1/250 dilution; Cell Signaling Technology). Medium from cells expressing gH/gL/gO-sfGFP was then added to a final dilution of 1/3, and the cells were incubated on ice for 20 minutes. Cells were washed twice with PBS-BSA, and sorted on a BD FACS Aria II at the Roy J. Carver Biotechnology Center. The main cell population was gated by forward/side scattering to remove debris and doublets, and propidium iodide (1 μg/ml final) was added to the sample to exclude dead cells. Of the myc-positive (Alexa 647) population, the 15% of cells with the highest and lowest GFP fluorescence were collected (Fig 1F) in tubes coated overnight with fetal bovine serum and containing Expi293 Expression Medium. Total RNA was extracted from the collected cells using a GeneJET RNA purification kit (Thermo Scientific), and cDNA was reverse transcribed with high fidelity Accuscript primed with a gene-specific oligonucleotide (RevTrans_PDGFRA_992R, TCATGCAGGTTGACAGCTTC). The diversified region of PDGFRα was PCR amplified as 3 fragments to provide full coverage of the D2-D3 domains during Illumina sequencing. During PCR, flanking sequences on the primers added adapters to the ends of the products for annealing to Illumina sequencing primers, unique barcoding, and for binding the flow cell. Amplicons were sequenced on an Illumina NovaSeq 6000 using a 2×250 nt paired end protocol. Data were analyzed using Enrich [85], and commands are provided in the GEO deposit. Briefly, the frequencies of PDGFRα variants in the transcripts of the sorted populations were compared to their frequencies in the naive plasmid library to calculate an enrichment ratio.
Production of soluble PDGFRα-Fc
Per ml of 2 × 106 Expi293F cells, 500 ng of pcDNA3-sPDGFRα-Fc (IgG1) and 3 μg of polyethylenimine (MW 25,000; Polysciences) were mixed in 100 μl of OptiMEM (Gibco) and incubated for 20 minutes at room temperature prior to adding to cells. Transfection Enhancers were added 18 h post-transfection, and cells were cultured for six to seven days. Cells were removed by centrifugation at 600 × g for 20 minutes at 4°C. Cell debris and precipitates were removed by centrifugation at 18,000 × g for 25 minutes at 4°C. Supernatant was loaded on to KANEKA KanCapA 3G Affinity sorbent (Pall), and the resin was washed with PBS. Soluble PDGFRα-Fc was eluted with 60 mM sodium acetate pH 3.7, and 1 M Tris pH 8.0 was added to the eluate to neutralize the pH. Eluted sPDGFRα-Fc was concentrated with a centrifugal device (MWCO 100 kDa; Sartorius) and NaCl was added to 150 mM final. The protein was separated on a Superdex 200 Increase 10/300 GL (GE Healthcare Life Sciences) size exclusion chromatography column equilibrated with PBS. Peak fractions were pooled, concentrated, and stored at -80°C after snap freezing in liquid nitrogen. The proteins were quantified by absorbance at 280 nm using calculated molar extinction coefficients for the monomeric mature polypeptides.
Fusions of sPDGFRα to the Fc region of IgG3 were expressed as described above, and the expression medium was dialyzed against water. Cell debris was removed by centrifugation at 18,000 × g for 25 minutes at 4°C. Supernatant was loaded on to protein G HTC agarose beads (GoldBio) equilibrated with PBS. Protein was eluted with 100 mM glycine pH 2.5 and the eluate was neutralized by addition of 1 M Tris pH 9.0. Protein was further purified as described for the IgG1 fusions.
Flow cytometry analysis of PDGFRα-Fc binding to HCMV trimer on the cell surface
400 ng of each of pcDNA3-gH (full-length), pcDNA3-gL-FLAG and pcDNA3-gO-8his were transfected in to Expi293F cells at 2 × 106 / ml using Expifectamine. Transfection Enhancers were added 22 h later, and cells were harvested 46 h post-transfection. Cells were washed with PBS-BSA and incubated with PDGFRα-Fc (purified as described) or PDGFRβ-Fc (R&D Systems) for 40 minutes on ice. Cells were then washed 3 times with PBS-BSA and stained with anti-FLAG M2-Cy3 (Sigma), chicken anti-HIS-FITC (polyclonal, Immunology Consultants Laboratory), and anti-human IgG-APC (clone HP6017, BioLegend) for 30 minutes on ice. Cells were washed three times with PBS-BSA and analyzed by flow cytometry. To compare data across different experiments, the change in ΔMFU for each condition was normalized to the ΔMFU at the maximum concentration of WT sPDGFRα-Fc: Relative binding = (MFUsample−MFUbackground) / (MFUmax WT − MFUbackground).
PDGFRα signaling assay
500 ng of PDGFRα reporter plasmid was transfected into 1 ml Expi293F cells at 2 × 106 / ml using Expifectamine. 7.5 μl of 6 μM sPDGFRα-Fc (concentration based on monomer) and 7.5 μl of 2 μM PDGF were mixed, incubated at room temperature for 40 minutes, and then added to 1 ml cells at transfection. Human gamma globulin was from Jackson Immuno Research Labs, and PDGFs were from R&D Systems. Cells were collected 24 h post-transfection and stained with anti-myc-Alexa 647 to detect PDGFRα expression. GFP-PEST reporter expression was measured by flow-cytometry.
Virus neutralization assays
Neutralization assays were performed similarly to Stegmann et al. [38] by serially diluting Fc-fused PDGFRα proteins or anti-HCMV immunoglobulin (Cytogam, CSL Behring) in PBS (vehicle) and incubating diluted proteins with HCMV virions for 1 h at 37°C. Virion-protein mixtures were adsorbed onto target cells for 2 h at 37°C / 5% CO2, washed once with PBS, and infected cells were allowed to recover for 18 h. The ability of each treatment to neutralize HCMV infection was measured using indirect-immunofluorescence by determining the number of cells expressing viral IE1, using anti-IE1 clone 1B12 [83] (generously provided by Thomas Shenk, Princeton University) and counterstaining nuclei with DAPI. A multiplicity of infection of 1 was used for AD169 infections, and a multiplicity of infection of 0.5 was used for TB40/E infections (based on stock titers acquired on MRC-5 fibroblasts).
Structural modeling
The sequence of human PDGFRα was threaded onto the structure of PDGF-BB-bound PDGFRβ (PDB ID 3MJG [44]), with one of the two receptor chains and the D1 domain removed. The model was minimized with ROSETTA using FastRelax [86]. Connectivity was determined using the AverageDegree filter in ROSETTA [87]. Images were rendered with PyMOL (Schrödinger, LLC).
Supporting information
S1 Fig. Binding of gO-sfGFP to PDGFRα-positive cells is dependent on co-expression of gH and gL.
Expi293F cells were incubated with medium from cells expressing (A) gO-sfGFP, (B) gO(C343S)-sfGFP or (C) gH/gL/gO-sfGFP, washed and analyzed by flow cytometry. Cells were either untransfected (black), expressing PDGFRβ (blue) or expressing PDGFRα (red). High binding signal that was specific for PDGFRα-positive cells was only observed for medium containing all three components of the HCMV trimer.
https://doi.org/10.1371/journal.ppat.1008647.s001
(TIF)
S2 Fig. A deep mutational scan of PDGFRα D2-D3 domains for HCMV trimer interactions in the presence of PDGFs.
(A) Log2 enrichment ratios for individual mutations in the UL74-High sort are plotted from -3 (orange, depleted) to +3 (dark blue, enriched). Amino acid position is on the horizontal axis, and substitutions are on the vertical axis. *, stop codon. (B) Agreement between log2 enrichment ratios from independent replicates of the UL74-High sort (positive selection). R2 values are calculated for nonsynonymous mutations in black. Nonsense mutations are red. (C) Agreement between log2 enrichment ratios from replicates of the UL74-Low sort (negative selection). (D) Log2 enrichment ratios for nonsynonymous mutations (black) are anticorrelated between the negative and positive selections. Nonsense mutations (red) are depleted from both sorts due to lost surface expression. (E-G) High correlation between conservation scores (calculated by averaging the log2 enrichment ratios for all nonsynonymous mutations at a given amino acid position) from independent replicates of the UL74-High (E) and UL74-Low (F) sorts. Conservation scores are anticorrelated between the two sorted populations (G).
https://doi.org/10.1371/journal.ppat.1008647.s002
(TIF)
S3 Fig. There are no hot spots for enriched mutations in the negative selection for loss of HCMV trimer binding.
Log2 enrichment ratios for single amino acid substitutions of PDGFRα are plotted based on their enrichment in the UL74-Low sort, from -3 (orange, depleted) to +3 (dark blue, enriched). Amino acid position is on the horizontal axis, and substitutions are on the vertical axis. *, stop codon. Mutations to critical residues for HCMV trimer binding are anticipated to be enriched (dark blue) in this negative selection. However, there are no unambiguous hot spot regions for enriched mutations. Compare to the positive selection shown in S2A Fig, which uses the same color scale.
https://doi.org/10.1371/journal.ppat.1008647.s003
(TIF)
S4 Fig. PDGFRα mutations that increase HCMV trimer binding in the presence of competing PDGFs are biased to structurally connected residues.
Residue conservation scores in the UL74-High deep mutational scan were calculated by averaging the log2 enrichment ratios for all 20 possible amino acids at each diversified position. PDGFRα residues where mutations tend to increase HCMV trimer binding in the presence of competing PDGFs have higher positive scores. A residue’s conservation score is correlated with its connectivity in the modeled PDGF-bound PDGFRα structure, where connectivity is quantified by the number of neighboring residues within a 12 Å radius. Highly connected residues are either buried in the hydrophobic cores of the D2-D3 domains, or are buried at the PDGF binding interface.
https://doi.org/10.1371/journal.ppat.1008647.s004
(TIF)
S5 Fig. Purification of soluble IgG1 Fc-fused PDGFRα.
(A) The extracellular domain of PDGFRα (Gln24-Glu524; grey) was fused via a short linker (purple) to the Fc region of IgG1 (green). The “Legacy” sequence corresponds to the commercially supplied protein (R&D Systems) used in prior publications. The sequence was redesigned for this study. (B) Coomassie-stained SDS gel (run under denaturing and reducing conditions) of wild-type sPDGFRα-Fc eluted from a protein A column. The monomeric protein MW is predicted to be 82 kD. Additional weight may come from glycosylation and/or anomalous electrophoretic mobility. (C) SEC elution of wild-type (solid black line), Y206S (grey line) and V242K (dashed black line) sPDGFRα-Fc. UV absorbance (y-axis) is scaled.
https://doi.org/10.1371/journal.ppat.1008647.s005
(TIF)
S6 Fig. Chemical stress tests of sPDGFRα-Fc.
(A) The most promising engineered orthogonal receptor, sPDGFRα-Fc V242K, was incubated at 40°C for 7 days in 20 mM Tris pH 8.5 with 10 mM EDTA to promote Asn deamidation, or at 40°C for 14 days in 50 mM sodium acetate pH 5.5 to promote Asn isomerization. The control sample in PBS (pH 7.4) was flash frozen and stored at -80°C until analysis. SDS-polyacrylamide gel electrophoresis with Coomassie blue staining shows chemical instability of sPDGFRα-Fc V242K in the harsher pH 5.5 stress test. (B) Stressed proteins were analyzed by SEC on a Superdex 200 Increase 10/300 GL column with PBS pH 7.4 as the running buffer.
https://doi.org/10.1371/journal.ppat.1008647.s006
(TIF)
S7 Fig. Soluble PDGFRα-Fc V242K binds HCMV trimer with comparable affinity to wild-type sPDGFRα-Fc.
(A) Data presented in Fig 3D was replicated using independent preparations of sPDGFRα WT (solid black line) and V242K (broken black line) fused to the Fc region of IgG1. Binding to Expi293F cells expressing full-length gH, gL and gO from the HCMV TB40/E strain was assessed by flow cytometry. (B and C) Soluble PDGFRα WT (solid black line) and V242K (broken black line) were purified as fusions to the Fc region of IgG3, matching the redesigned linker described in S5A Fig. Binding to trimer from (B) TB40/E and (C) Merlin strains expressed on Expi293F cells was measured by flow cytometry. Data are mean ± SD, n = 3 (sPDGFRα-Fc WT and V242K) or 2 (sPDGFRβ-Fc).
https://doi.org/10.1371/journal.ppat.1008647.s007
(TIF)
S8 Fig. Gamma globulin has no effect on PDGF signaling in culture.
PDGF signaling was assessed in Expi293F cells transiently transfected with a PDGFRα reporter plasmid as outlined in Fig 4A. A 3-fold molar excess of wild-type sPDGFRα-Fc (black) blocks 5 nM PDGF-AB signaling. Orthogonal engineered sPDGFRα-Fc V242K (white) does not block signaling and the mutation shows no interactions with PDGFs in competitive binding experiments (Fig 3). Furthermore, wild-type or mutant sPDGFRα-Fc up to 100 nM shows no binding to PDGFRα-positive cells, and we exclude unanticipated interactions between receptor chains. However, sPDGFRα-Fc V242K does cause an increase in signaling of PDGF-A ligands (AA and AB) for unknown reasons. A non-specific carrier effect in which the orthogonal receptor stabilizes the hydrophobic ligand in solution is suspected, especially because the more hydrophobic variant sPDGFRα-Fc Y206S promotes an even bigger increase in PDGF-A signaling (Fig 4B). Human gamma globulin (violet) at the same concentration has no effect in this assay. Data are mean ± SD, n = 2 independent replicates.
https://doi.org/10.1371/journal.ppat.1008647.s008
(TIF)
Acknowledgments
Barbara Pilas and Barbara Balhan assisted with FACS, and Alvaro Hernandez and Chris Wright helped with deep sequencing (UIUC Roy J. Carver Biotechnology Center).
References
- 1. Cannon MJ, Schmid DS, Hyde TB. Review of cytomegalovirus seroprevalence and demographic characteristics associated with infection. Rev Med Virol. 2010;20: 202–213. pmid:20564615
- 2. Staras SAS, Dollard SC, Radford KW, Flanders WD, Pass RF, Cannon MJ. Seroprevalence of cytomegalovirus infection in the United States, 1988–1994. Clin Infect Dis. 2006;43: 1143–1151. pmid:17029132
- 3. Emery VC. Investigation of CMV disease in immunocompromised patients. J Clin Pathol. 2001;54: 84–88. pmid:11215290
- 4. Cobbs CS, Harkins L, Samanta M, Gillespie GY, Bharara S, King PH, et al. Human cytomegalovirus infection and expression in human malignant glioma. Cancer Res. 2002;62: 3347–3350. pmid:12067971
- 5. Mitchell DA, Xie W, Schmittling R, Learn C, Friedman A, McLendon RE, et al. Sensitive detection of human cytomegalovirus in tumors and peripheral blood of patients diagnosed with glioblastoma. Neuro-oncology. 2008;10: 10–18. pmid:17951512
- 6. Scheurer ME, Bondy ML, Aldape KD, Albrecht T, El-Zein R. Detection of human cytomegalovirus in different histological types of gliomas. Acta Neuropathol. 2008;116: 79–86. pmid:18351367
- 7. Michaelis M, Doerr HW, Cinatl J. The story of human cytomegalovirus and cancer: increasing evidence and open questions. Neoplasia. 2009;11: 1–9. pmid:19107226
- 8. Lawler SE. Cytomegalovirus and glioblastoma; controversies and opportunities. J Neurooncol. 2015;123: 465–471. pmid:25682092
- 9. Cannon MJ, Davis KF. Washing our hands of the congenital cytomegalovirus disease epidemic. BMC Public Health. 4 ed. 2005;5: 70. pmid:15967030
- 10. Dollard SC, Grosse SD, Ross DS. New estimates of the prevalence of neurological and sensory sequelae and mortality associated with congenital cytomegalovirus infection. Rev Med Virol. 2007;17: 355–363. pmid:17542052
- 11. Kenneson A, Cannon MJ. Review and meta-analysis of the epidemiology of congenital cytomegalovirus (CMV) infection. Rev Med Virol. 2007;17: 253–276. pmid:17579921
- 12. Grosse SD, Ross DS, Dollard SC. Congenital cytomegalovirus (CMV) infection as a cause of permanent bilateral hearing loss: a quantitative assessment. J Clin Virol. 2008;41: 57–62. pmid:17959414
- 13. Suzuki Y, Toribe Y, Mogami Y, Yanagihara K, Nishikawa M. Epilepsy in patients with congenital cytomegalovirus infection. Brain Dev. 2008;30: 420–424. pmid:18215482
- 14. Zhang X-W, Li F, Yu X-W, Shi X-W, Shi J, Zhang J-P. Physical and intellectual development in children with asymptomatic congenital cytomegalovirus infection: a longitudinal cohort study in Qinba mountain area, China. J Clin Virol. 2007;40: 180–185. pmid:17919973
- 15. Gentile I, Zappulo E, Riccio MP, Binda S, Bubba L, Pellegrinelli L, et al. Prevalence of Congenital Cytomegalovirus Infection Assessed Through Viral Genome Detection in Dried Blood Spots in Children with Autism Spectrum Disorders. In Vivo. 2017;31: 467–473. pmid:28438881
- 16. Staras SAS, Flanders WD, Dollard SC, Pass RF, McGowan JE, Cannon MJ. Cytomegalovirus seroprevalence and childhood sources of infection: A population-based study among pre-adolescents in the United States. J Clin Virol. 2008;43: 266–271. pmid:18778968
- 17. Lanzieri TM, Kruszon-Moran D, Amin MM, Bialek SR, Cannon MJ, Carroll MD, et al. Seroprevalence of cytomegalovirus among children 1 to 5 years of age in the United States from the National Health and Nutrition Examination Survey of 2011 to 2012. Hodinka RL, editor. Clin Vaccine Immunol. American Society for Microbiology; 2015;22: 245–247. pmid:25520150
- 18. Dupont L, Reeves MB. Cytomegalovirus latency and reactivation: recent insights into an age old problem. Rev Med Virol. 2016;26: 75–89. pmid:26572645
- 19. Ross SA, Arora N, Novak Z, Fowler KB, Britt WJ, Boppana SB. Cytomegalovirus reinfections in healthy seroimmune women. J Infect Dis. 2010;201: 386–389. pmid:20039807
- 20. Lurain NS, Chou S. Antiviral drug resistance of human cytomegalovirus. Clin Microbiol Rev. American Society for Microbiology Journals; 2010;23: 689–712. pmid:20930070
- 21.
American College of Obstetricians and Gynecologists. Practice bulletin no. 151: Cytomegalovirus, parvovirus B19, varicella zoster, and toxoplasmosis in pregnancy. Obstetrics and gynecology. 2015. pp. 1510–1525. https://doi.org/10.1097/01.AOG.0000466430.19823.53
- 22. Gerna G, Kabanova A, Lilleri D. Human Cytomegalovirus Cell Tropism and Host Cell Receptors. Vaccines (Basel). 2019;7: 70. pmid:31336680
- 23. Hahn G, Revello MG, Patrone M, Percivalle E, Campanini G, Sarasini A, et al. Human cytomegalovirus UL131-128 genes are indispensable for virus growth in endothelial cells and virus transfer to leukocytes. J Virol. American Society for Microbiology Journals; 2004;78: 10023–10033.
- 24. Wang D, Shenk T. Human cytomegalovirus virion protein complex required for epithelial and endothelial cell tropism. Proc Natl Acad Sci USA. National Academy of Sciences; 2005;102: 18153–18158. pmid:16319222
- 25. Wang D, Shenk T. Human cytomegalovirus UL131 open reading frame is required for epithelial cell tropism. J Virol. American Society for Microbiology Journals; 2005;79: 10330–10338. pmid:16051825
- 26. Ryckman BJ, Chase MC, Johnson DC. HCMV gH/gL/UL128-131 interferes with virus entry into epithelial cells: evidence for cell type-specific receptors. Proc Natl Acad Sci USA. National Academy of Sciences; 2008;105: 14118–14123. pmid:18768787
- 27. Ryckman BJ, Rainish BL, Chase MC, Borton JA, Nelson JA, Jarvis MA, et al. Characterization of the human cytomegalovirus gH/gL/UL128-131 complex that mediates entry into epithelial and endothelial cells. J Virol. American Society for Microbiology Journals; 2008;82: 60–70. pmid:17942555
- 28. Adler B, Scrivano L, Ruzcics Z, Rupp B, Sinzger C, Koszinowski U. Role of human cytomegalovirus UL131A in cell type-specific virus entry and release. J Gen Virol. Microbiology Society; 2006;87: 2451–2460. pmid:16894182
- 29. Martinez-Martin N, Marcandalli J, Huang CS, Arthur CP, Perotti M, Foglierini M, et al. An Unbiased Screen for Human Cytomegalovirus Identifies Neuropilin-2 as a Central Viral Receptor. Cell. 2018;174: 1158–1171.e19. pmid:30057110
- 30. E X, Meraner P, Lu P, Perreira JM, Aker AM, McDougall WM, et al. OR14I1 is a receptor for the human cytomegalovirus pentameric complex and defines viral epithelial cell tropism. Proc Natl Acad Sci USA. National Academy of Sciences; 2019;116: 7043–7052. pmid:30894498
- 31. Huber MT, Compton T. The human cytomegalovirus UL74 gene encodes the third component of the glycoprotein H-glycoprotein L-containing envelope complex. J Virol. American Society for Microbiology (ASM); 1998;72: 8191–8197. pmid:9733861
- 32. Huber MT, Compton T. Intracellular formation and processing of the heterotrimeric gH-gL-gO (gCIII) glycoprotein envelope complex of human cytomegalovirus. J Virol. American Society for Microbiology (ASM); 1999;73: 3886–3892. pmid:10196283
- 33. Zhou M, Lanchy J-M, Ryckman BJ. Human Cytomegalovirus gH/gL/gO Promotes the Fusion Step of Entry into All Cell Types, whereas gH/gL/UL128-131 Broadens Virus Tropism through a Distinct Mechanism. Hutt-Fletcher L, editor. J Virol. American Society for Microbiology Journals; 2015;89: 8999–9009. pmid:26085146
- 34. Kabanova A, Marcandalli J, Zhou T, Bianchi S, Baxa U, Tsybovsky Y, et al. Platelet-derived growth factor-α receptor is the cellular receptor for human cytomegalovirus gHgLgO trimer. Nat Microbiol. 2016;1: 16082. pmid:27573107
- 35. Soroceanu L, Akhavan A, Cobbs CS. Platelet-derived growth factor-alpha receptor activation is required for human cytomegalovirus infection. Nature. 2008;455: 391–395. pmid:18701889
- 36. Wu Y, Prager A, Boos S, Resch M, Brizic I, Mach M, et al. Human cytomegalovirus glycoprotein complex gH/gL/gO uses PDGFR-α as a key for entry. Yurochko A, editor. PLoS Pathog. 2017;13: e1006281. pmid:28403202
- 37. Wu K, Oberstein A, Wang W, Shenk T. Role of PDGF receptor-α during human cytomegalovirus entry into fibroblasts. Proc Natl Acad Sci USA. National Academy of Sciences; 2018;115: E9889–E9898. pmid:30275317
- 38. Stegmann C, Hochdorfer D, Lieber D, Subramanian N, Stöhr D, Laib Sampaio K, et al. A derivative of platelet-derived growth factor receptor alpha binds to the trimer of human cytomegalovirus and inhibits entry into fibroblasts and endothelial cells. Yurochko A, editor. PLoS Pathog. 2017;13: e1006273. pmid:28403220
- 39. De Clercq E, Li G. Approved Antiviral Drugs over the Past 50 Years. Clin Microbiol Rev. American Society for Microbiology Journals; 2016;29: 695–747. pmid:27281742
- 40. Andrae J, Gallini R, Betsholtz C. Role of platelet-derived growth factors in physiology and medicine. Genes Dev. Cold Spring Harbor Lab; 2008;22: 1276–1312. pmid:18483217
- 41. Chen P-H, Chen X, He X. Platelet-derived growth factors and their receptors: structural and functional perspectives. Biochim Biophys Acta. 2013;1834: 2176–2186. pmid:23137658
- 42. Fretto LJ, Snape AJ, Tomlinson JE, Seroogy JJ, Wolf DL, LaRochelle WJ, et al. Mechanism of platelet-derived growth factor (PDGF) AA, AB, and BB binding to alpha and beta PDGF receptor. J Biol Chem. 1993;268: 3625–3631. pmid:7679113
- 43. Li X, Pontén A, Aase K, Karlsson L, Abramsson A, Uutela M, et al. PDGF-C is a new protease-activated ligand for the PDGF alpha-receptor. Nat Cell Biol. 2000;2: 302–309. pmid:10806482
- 44. Shim AH-R, Liu H, Focia PJ, Chen X, Lin PC, He X. Structures of a platelet-derived growth factor/propeptide complex and a platelet-derived growth factor/receptor complex. Proc Natl Acad Sci USA. National Academy of Sciences; 2010;107: 11307–11312. pmid:20534510
- 45. Chen P-H, Unger V, He X. Structure of Full-Length Human PDGFRβ Bound to Its Activating Ligand PDGF-B as Determined by Negative-Stain Electron Microscopy. J Mol Biol. 2015;427: 3921–3934. pmid:26463591
- 46. Fowler DM, Fields S. Deep mutational scanning: a new style of protein science. Nat Methods. Nature Publishing Group; 2014;11: 801–807. pmid:25075907
- 47. Berger S, Procko E, Margineantu D, Lee EF, Shen BW, Zelter A, et al. Computationally designed high specificity inhibitors delineate the roles of BCL2 family proteins in cancer. Elife. 2016;5: 1422. pmid:27805565
- 48. Procko E, Berguig GY, Shen BW, Song Y, Frayo S, Convertine AJ, et al. A computationally designed inhibitor of an Epstein-Barr viral Bcl-2 protein induces apoptosis in infected cells. Cell. 2014;157: 1644–1656. pmid:24949974
- 49. Wrenbeck EE, Azouz LR, Whitehead TA. Single-mutation fitness landscapes for an enzyme on multiple substrates reveal specificity is globally encoded. Nat Commun. 2017;8: 15695. pmid:28585537
- 50. Park J, Selvam B, Sanematsu K, Shigemura N, Shukla D, Procko E. Structural architecture of a dimeric class C GPCR based on co-trafficking of sweet taste receptor subunits. Journal of Biological Chemistry. American Society for Biochemistry and Molecular Biology; 2019;294: 4759–4774. pmid:30723160
- 51. Heredia JD, Park J, Brubaker RJ, Szymanski SK, Gill KS, Procko E. Mapping Interaction Sites on Human Chemokine Receptors by Deep Mutational Scanning. J Immunol. American Association of Immunologists; 2018;200: ji1800343–3839. pmid:29678950
- 52. Heredia JD, Park J, Choi H, Gill KS, Procko E. Conformational Engineering of HIV-1 Env Based on Mutational Tolerance in the CD4 and PG16 Bound States. Kirchhoff F, editor. J Virol. American Society for Microbiology Journals; 2019;93: e00219–19. pmid:30894475
- 53. Sinzger C, Hahn G, Digel M, Katona R, Sampaio KL, Messerle M, et al. Cloning and sequencing of a highly productive, endotheliotropic virus strain derived from human cytomegalovirus TB40/E. J Gen Virol. Microbiology Society; 2008;89: 359–368. pmid:18198366
- 54. Stegmann C, Rothemund F, Laib Sampaio K, Adler B, Sinzger C. The N Terminus of Human Cytomegalovirus Glycoprotein O Is Important for Binding to the Cellular Receptor PDGFRα. Sandri-Goldin RM, editor. J Virol. American Society for Microbiology Journals; 2019;93: 93. pmid:30894468
- 55. Pédelacq J-D, Cabantous S, Tran T, Terwilliger TC, Waldo GS. Engineering and characterization of a superfolder green fluorescent protein. Nat Biotechnol. 2006;24: 79–88. pmid:16369541
- 56. Ciferri C, Chandramouli S, Donnarumma D, Nikitin PA, Cianfrocco MA, Gerrein R, et al. Structural and biochemical studies of HCMV gH/gL/gO and Pentamer reveal mutually exclusive cell entry complexes. Proc Natl Acad Sci USA. National Academy of Sciences; 2015;112: 1767–1772. pmid:25624487
- 57. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. Oxford University Press; 2002;30: 207–210. pmid:11752295
- 58. Stanton RJ, Baluchova K, Dargan DJ, Cunningham C, Sheehy O, Seirafian S, et al. Reconstruction of the complete human cytomegalovirus genome in a BAC reveals RL13 to be a potent inhibitor of replication. J Clin Invest. American Society for Clinical Investigation; 2010;120: 3191–3208. pmid:20679731
- 59. Czajkowsky DM, Hu J, Shao Z, Pleass RJ. Fc-fusion proteins: new developments and future perspectives. EMBO Mol Med. 2012;4: 1015–1028. pmid:22837174
- 60. Lu X, Nobrega RP, Lynaugh H, Jain T, Barlow K, Boland T, et al. Deamidation and isomerization liability analysis of 131 clinical-stage antibodies. MAbs. 2019;11: 45–57. pmid:30526254
- 61. Vidarsson G, Dekkers G, Rispens T. IgG subclasses and allotypes: from structure to effector functions. Front Immunol. Frontiers; 2014;5: 520. pmid:25368619
- 62. Roopenian DC, Christianson GJ, Sproule TJ, Brown AC, Akilesh S, Jung N, et al. The MHC class I-like IgG receptor controls perinatal IgG transport, IgG homeostasis, and fate of IgG-Fc-coupled drugs. J Immunol. American Association of Immunologists; 2003;170: 3528–3533. pmid:12646614
- 63. Li X, Zhao X, Fang Y, Jiang X, Duong T, Fan C, et al. Generation of destabilized green fluorescent protein as a transcription reporter. J Biol Chem. American Society for Biochemistry and Molecular Biology; 1998;273: 34970–34975. pmid:9857028
- 64. Reichhart E, Ingles-Prieto A, Tichy A-M, McKenzie C, Janovjak H. A Phytochrome Sensory Domain Permits Receptor Activation by Red Light. Angew Chem Int Ed Engl. 2016;55: 6339–6342. pmid:27101018
- 65. Yu D, Smith GA, Enquist LW, Shenk T. Construction of a self-excisable bacterial artificial chromosome containing the human cytomegalovirus genome and mutagenesis of the diploid TRL/IRL13 gene. J Virol. American Society for Microbiology Journals; 2002;76: 2316–2328.
- 66. Vanarsdall AL, Chin AL, Liu J, Jardetzky TS, Mudd JO, Orloff SL, et al. HCMV trimer- and pentamer-specific antibodies synergize for virus neutralization but do not correlate with congenital transmission. Proc Natl Acad Sci USA. National Academy of Sciences; 2019;116: 3728–3733. pmid:30733288
- 67. Deen KC, McDougal JS, Inacker R, Folena-Wasserman G, Arthos J, Rosenberg J, et al. A soluble form of CD4 (T4) protein inhibits AIDS virus infection. Nature. 1988;331: 82–84. pmid:3257544
- 68. Fisher RA, Bertonis JM, Meier W, Johnson VA, Costopoulos DS, Liu T, et al. HIV infection is blocked in vitro by recombinant soluble CD4. Nature. 1988;331: 76–78. pmid:2829022
- 69. Marlin SD, Staunton DE, Springer TA, Stratowa C, Sommergruber W, Merluzzi VJ. A soluble form of intercellular adhesion molecule-1 inhibits rhinovirus infection. Nature. 1990;344: 70–72. pmid:1968231
- 70. Norkin LC. Virus receptors: implications for pathogenesis and the design of antiviral agents. Clin Microbiol Rev. American Society for Microbiology (ASM); 1995;8: 293–315. pmid:7621403
- 71. Bialas KM, Tanaka T, Tran D, Varner V, Cisneros De La Rosa E, Chiuppesi F, et al. Maternal CD4+ T cells protect against severe congenital cytomegalovirus disease in a novel nonhuman primate model of placental cytomegalovirus transmission. Proc Natl Acad Sci USA. National Academy of Sciences; 2015;112: 13645–13650. pmid:26483473
- 72. Mocarski ES, Bonyhadi M, Salimi S, McCune JM, Kaneshima H. Human cytomegalovirus in a SCID-hu mouse: thymic epithelial cells are prominent targets of viral replication. Proc Natl Acad Sci USA. National Academy of Sciences; 1993;90: 104–108. pmid:7678330
- 73. Smith MS, Goldman DC, Bailey AS, Pfaffle DL, Kreklywich CN, Spencer DB, et al. Granulocyte-colony stimulating factor reactivates human cytomegalovirus in a latently infected humanized mouse model. Cell Host Microbe. 2010;8: 284–291. pmid:20833379
- 74. Crawford LB, Streblow DN, Hakki M, Nelson JA, Caposio P. Humanized mouse models of human cytomegalovirus infection. Curr Opin Virol. 2015;13: 86–92. pmid:26118890
- 75. Antoniades HN, Scher CD, Stiles CD. Purification of human platelet-derived growth factor. Proc Natl Acad Sci USA. National Academy of Sciences; 1979;76: 1809–1813. pmid:287022
- 76. Raines EW, Ross R. Purification of human platelet-derived growth factor. Meth Enzymol. Elsevier; 1985;109: 749–773. pmid:3990574
- 77. Petkova SB, Akilesh S, Sproule TJ, Christianson GJ, Khabbaz Al H, Brown AC, et al. Enhanced half-life of genetically engineered human IgG1 antibodies in a humanized FcRn mouse model: potential application in humorally mediated autoimmune disease. Int Immunol. 2006;18: 1759–1769. pmid:17077181
- 78. Roopenian DC, Christianson GJ, Proetzel G, Sproule TJ. Human FcRn Transgenic Mice for Pharmacokinetic Evaluation of Therapeutic Antibodies. Methods Mol Biol. New York, NY: Springer New York; 2016;1438: 103–114. pmid:27150086
- 79. McShan AC, Devlin CA, Overall SA, Park J, Toor JS, Moschidi D, et al. Molecular determinants of chaperone interactions on MHC-I for folding and antigen repertoire selection. Proc Natl Acad Sci USA. National Academy of Sciences; 2019;116: 25602–25613. pmid:31796585
- 80. Wu NC, Xie J, Zheng T, Nycholat CM, Grande G, Paulson JC, et al. Diversity of Functionally Permissive Sequences in the Receptor-Binding Site of Influenza Hemagglutinin. Cell Host Microbe. 2017;21: 742–753.e8. pmid:28618270
- 81. Dingens AS, Haddox HK, Overbaugh J, Bloom JD. Comprehensive Mapping of HIV-1 Escape from a Broadly Neutralizing Antibody. Cell Host Microbe. 2017;21: 777–787.e4. pmid:28579254
- 82. Procko E. The sequence of human ACE2 is suboptimal for binding the S spike protein of SARS coronavirus 2. bioRxiv. 2020;: 2020.03.16.994236.
- 83. Zhu H, Shen Y, Shenk T. Human cytomegalovirus IE1 and IE2 proteins block apoptosis. J Virol. American Society for Microbiology (ASM); 1995;69: 7960–7970. pmid:7494309
- 84. Procko E, Hedman R, Hamilton K, Seetharaman J, Fleishman SJ, Su M, et al. Computational design of a protein-based enzyme inhibitor. J Mol Biol. 2013;425: 3563–3575. pmid:23827138
- 85. Fowler DM, Araya CL, Gerard W, Fields S. Enrich: software for analysis of protein function by enrichment and depletion of variants. Bioinformatics. Oxford University Press; 2011;27: 3430–3431. pmid:22006916
- 86. Leaver-Fay A, Tyka M, Lewis SM, Lange OF, Thompson J, Jacak R, et al. ROSETTA3: an object-oriented software suite for the simulation and design of macromolecules. Meth Enzymol. Elsevier; 2011;487: 545–574. pmid:21187238
- 87. Fleishman SJ, Whitehead TA, Strauch E-M, Corn JE, Qin S, Zhou H-X, et al. Community-wide assessment of protein-interface modeling suggests improvements to design methodology. J Mol Biol. 2011;414: 289–302. pmid:22001016