key: cord-0768045-maa8c4a4
authors: Zhang, Yuan; Zheng, Nan; Zhong, Yang
title: Computational characterization and design of SARS coronavirus receptor recognition and antibody neutralization
date: 2007-02-17
journal: Comput Biol Chem
DOI: 10.1016/j.compbiolchem.2007.02.005
sha: d2727abafe0e60a9f3d42f5d07d19f832636e47b
doc_id: 768045
cord_uid: maa8c4a4

The sequential determination of crystal structures of the SARS coronavirus spike receptor-binding domain (RBD) in complex with its cellular receptor or neutralizing antibody opened a door for the design and development of antiviral competitive inhibitors. Based on those complex structures, we conduct computational characterization and design of RBD-mediated receptor recognition and antibody neutralization. The comparisons between computational predictions and experimental evidences validate our structural bioinformatics protocols. And the calculations predict a number of single substitutions on RBD, receptor or antibody that could remarkably elevate the binding affinities of those complexes. It is reasonable to anticipate our structure-based computation-derived hypotheses could be informative to the future biochemical and immunological tests.

As an envelope glycoprotein, the spike protein of severe acute respiratory syndrome coronavirus (SARS-CoV) plays a key role in the viral entry and neutralization (Bartlam et al., 2005; Denison, 2004; Lau and Peiris, 2005; Xu and Gao, 2004; Zhu, 2004) . This structural protein consists of two functional regions: the outer globular S1 region responsible for the initial attachment to cellular receptor and the inner stalk S2 region contributing to the subsequent fusion between viral envelope and cellular membrane (Beniac et al., 2006; Hofmann and Pohlmann, 2004; Lin et al., 2005; Xiao and Dimitrov, 2004) . A membrane-associated zinc metallopeptidase, angiotensinconverting enzyme 2 (ACE2), has been identified as the functional receptor for SARS-CoV (Li et al., 2003) . And a soluble form of ACE2 could block the association of S1 region with the permissive Vero E6 cells (Li et al., 2003; Moore et al., 2004) . In addition, a 193-amino acid fragment (residues 318-510), located within the S1 region, was demonstrated as an independently folded receptor-binding domain (RBD) capable of attaching ACE2 more efficiently (IC 50 < 10 nM) compared with the full S1 region (IC 50 ≈ 50 nM) . Besides, this RBD was able to elicit highly potent neutralizing antibodies in the immunized animals, which conferred those animals significant protection from the challenge of pathogenic SARS-CoV (Du et al., 2006; He et al., 2004 He et al., , 2005a He et al., ,b, 2006a Zakhartchouk et al., 2006; Zhao et al., 2006) . Moreover, a human monoclonal antibody 80R, isolated from a nonimmune human antibody library, was shown to potently neutralize SARS-CoV through targeting the RBD and blocking receptor recognition (Sui et al., 2004) . The epitope mapping illustrated a 180-amino acid conformationally sensitive fragment (residues 324-503) within the RBD was the neutralizing epitope of 80R (Sui et al., 2005) . Furthermore, another human monoclonal antibody m396 also exhibited potent neutralization of SARS-CoV by competition with ACE2 for binding to RBD (Prabakaran et al., 2006) . Together those data suggest the receptor association process of SARS-CoV is an attractive opportunity for therapeutic intervention (De Clercq, 2006; He and Jiang, 2005; Hofmann and Pohlmann, 2004; Jiang et al., 2005; Kuhn et al., 2004; Yeung et al., 2006) . The peptide or peptidomimetic antagonist leads, including the SARS-CoV spike RBD, the soluble form of ACE2 and the neutralizing antibodies 80R plus m396, should be able to potently abolish viral attachment to host cells. In this study, we conducted structural bioinformatics analyses on the crystal structures of the SARS-CoV RBD complexed with functional receptor or neutralizing antibody (Hwang et al., 2006; Li et al., 2005a Li et al., ,2006 Prabakaran et al., 2006) to predict single substitutions on spike RBD, receptor or antibodies possibly causing remarkable elevation in the binding affinities of complexes for the design and development of anti-SARS agents.

Three coordinates files were retrieved from the Protein Data Bank (PDB) (Berman et al., 2000) . One file is the ACE2-bound RBD (PDB code: 2AJF) (Li et al., 2005a) , while the others are the RBD complexed with 80R (PDB code: 2GHW) (Hwang et al., 2006) or m396 (PDB code: 2DD8) (Prabakaran et al., 2006) . Both the first and second files harbor a pair of sister complexes. And in the third file, either the heavy chain or the light chain of m396 makes its own contacts with the RBD. Thus, a total of six complex structures (AE/BF for ACE2-RBD, AB/CD for 80R-RBD and HS/LS for m396-RBD) are subject to computational simulations, respectively. Firstly, the program FoldX , based on an empirical effective energy function, was employed for calculation of the binding free energy values of wild type complexes. Then, a computational alanine scanning on the protein-protein interfaces was performed for evaluation of energetical contribution from single binding sites to the complex formation. Those positions yielding a calculated increase in association energy of more than 1 kcal/mol on alanine substitution were defined as energetic hot (important) spots according to previous criteria (Guerois et al., 2002; Guerois and Serrano, 2000; Kiel and Serrano, 2006; Kiel et al., 2004 Kiel et al., , 2005 . The next step was to redesign the interactions between RBD and its binding partners through the software DeepView (Arnold et al., 2006; Guex and Peitsch, 1997) . Each of the binding sites on the RBD, receptor or antibody was saturated with virtual substitutions, i.e., replaced with all the 20 natural amino acid residues except the original one. Finally, the reconstructed models were feed to the program FoldX to compute their binding energies. Here, only the variants rewarded a value of at least 1 kcal/mol lower than that of the wild type were taken into consideration.

The calculated binding energy values and hot spots of the wild type complexes are shown in Table 1 . The complexes ACE2-RBD (AE/BF) and 80R-RBD (AB/CD) show a close correlation between their interaction energies (−15.78/−14.82 kcal/mol versus −18.36/−16.53 kcal/mol) and buried surface area (1700Å 2 versus 2200Å 2 ), gap volume (7000Å 3 versus 4000Å 3 ), or binding affinity (1.70 nM versus 1.59 nM) (Hwang et al., 2006) . Those obvious associations indicate that the higher geometric complementarity, corresponding to the larger buried surface area and the smaller gap volume, offers the complex 80R-RBD rather than the complex ACE2-RBD the lower interaction energy and consequently the stronger binding affinity. Similarly, the correlation of binding energy with buried surface area is also found for the complex m396-RBD in which the heavy chain and the light chain contribute 66% and 34% to the total buried surface (Prabakaran et al., 2006) . And the RBD association energy of the former chain (−8.64 kcal/mol) is remarkably lower than that of the latter (−5.42 kcal/mol). The perfect agreements of computational predictions with structural observations or biochemical evidences strongly suggest the reliability of our protocols.

As to the hot spots of complexes, the consistency between computational predictions and experimental evidences is clearly detected for RBD and receptor. In ACE2-RBD complexes, three receptor residues (GLU37, ASP38 and TYR41 on the chain A of complex AE or GLU37, TYR41 and LYS353 on the chain B of complex BF) form one hot spot cluster interacting with another hot cluster formed by five or six RBD residues (ARG426, TYR436, TYR475, TYR484 and TYR491 on the chain E of complex AE and the chain F of complex BF, whereas ASN473 only on the chain F). The interactions between the two hot clusters make the major contribution to the binding free energy of ACE2-RBD complexes. Notably, our predictions are in agreement with previous experimental alanine mutagenesis, which identified two hot spots on RBD (ARG426 and ASN473) (Chakraborti et al., 2005) and another two on receptor (TYR41 and LYS353) . In addition, computational alanine scanning on the sister complexes AB and CD successfully identify a RBD hot spot (ASP480) revealed in mutational binding analyses (Sui et al., 2005) . In sharp contrast to ACE2, the antibody 80R possesses four hot residues (TYR102, ASN164, ARG223 and TRP226) being scattered on the binding surface rather than centralized into a cluster. The difference in the number and distribution of hot spots might account for the large gap between the interaction energies of 80R-RBD (−18.3627 and −16.5309 kcal/mol) and those of ACE2-RBD (−15.7765 and The complexes AE and BF with the chains A, B for ACE2 and the chains E, F for RBD; the complexes AB and CD with the chains A, C for RBD and the chains B, D for 80R; the complexes HS and LS with the chains H, L and S for heavy and light chains of m396 plus RBD. for association with RBD, or the higher spike-binding affinity of 80R compared to that of receptor. Finally, only two neighboring hot spots (TRP91 and ASP92) are found on the light chain of m396 while none on the heavy chain. Thus, an interesting discovery is the fact that among the five or six ACE2-binding hot spots of RBD, three (TYR436, TYR484 and TYR491) are simultaneously 80R-neutralizing hot spots whereas only one (TYR491) is important for m396 neutralization. This finding indicates that 80R might have the greater potential than m396 for inhibition of spike-mediated infection. In summary, the consistency of calculations with experiments mentioned above further validates our approaches to characterize protein-protein interactions. The predicted replacements on spike RBD, cellular receptor or neutralizing antibody with significant increase in binding affinity are listed in Table 2 . The comparisons between virtual mutants derived from sister complexes of ACE2-RBD or 80R-RBD consistently identify a number of substitutions worth of biochemical and immunological experimental tests. For instance, recent experimental evidences revealed the great potential of ACE2 in the protection of several animal models from SARS-CoV-induced lung injury or severe acute lung failure Kuba et al., 2005 Kuba et al., , 2006 . Simultaneously, the crystal structures of the native and inhibitor-bound forms of ACE2 Towler et al., 2004; Turner et al., 2004) successfully laid a solid foundation for the discovery of novel small-molecule inhibitors of its enzymatic activity or spikemediated virus entry by chemical genetics (Huentelman et al., 2004; Kao et al., 2004) and the identification of its crucial activesite residues by site-directed mutagenesis (Guy et al., 2005a,b) . Very recently, a modest anti-SARS activity (IC 50 ≈ 0.1 mM) was observed for an ACE2-derived peptide containing two segments of receptor (residues 22-44 and 351-357) linked by glycine (Han et al., 2006) . It should be pointed out that both the experimentally confirmed hot spots (TYR41 and LYS353) and the predicted sites for replacements (THR27, LYS31 and HIS34) are nested in those two segments. Similarly, a small peptide derived from spike protein (residues 483-493) also block viral receptor recognition with IC 50 of 6.99 nM (Ho et al., 2006) . And our calculated two hot spots (TYR484 and TYR491) in combination with two target positions (TYR484 and GLN492) are located in this short fragment, too. Consequently, it is reasonable to anticipate that our blueprint could effectively increase the binding affinity of the two novel peptides to disrupt SARS-CoV infection.

The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling

Structural insights into SARS coronavirus proteins

Architecture of the SARS coronavirus prefusion spike

The protein data bank

The SARS coronavirus S glycoprotein receptor binding domain: fine mapping and functional characterization

Potential antivirals and antiviral strategies against SARS coronavirus infections

Severe acute respiratory syndrome coronavirus pathogenesis, disease and vaccines: an update

Receptorbinding domain of SARS-CoV spike protein induces long-term protective immunity in an animal model

The SH3-fold family: experimental evidence and prediction of variations in the folding pathways

Predicting changes in the stability of proteins and protein complexes: a study of more than 1000 mutations

SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling

Identification of critical active-site residues in angiotensin-converting enzyme-2 (ACE2) by site-directed mutagenesis

Membrane-associated zinc peptidase families: comparing ACE and ACE2

Identification of critical determinants on ACE2 for SARS-CoV entry and development of a potent entry inhibitor

Vaccine design for severe acute respiratory syndrome coronavirus

Receptor-binding domain of SARS-CoV spike protein induces highly potent neutralizing antibodies: implication for developing subunit vaccine

Receptor-binding domain of severe acute respiratory syndrome coronavirus spike protein contains multiple conformation-dependent epitopes that induce highly potent neutralizing antibodies

Identification of a critical neutralization determinant of severe acute respiratory syndrome (SARS)-associated coronavirus: importance for designing SARS vaccines

Identification and characterization of novel neutralizing epitopes in the receptor-binding domain of SARS-CoV spike protein: revealing the critical antigenic determinants in inactivated SARS-CoV vaccine

Antigenic and immunogenic characterization of recombinant baculovirus-expressed severe acute respiratory syndrome coronavirus spike protein: implication for vaccine design

A single amino acid substitution (R441A) in the receptor-binding domain of SARS coronavirus spike protein disrupts the antigenic structure and binding activity

Crossneutralization of human and palm civet severe acute respiratory syndrome coronaviruses by antibodies targeting the receptor-binding domain of spike protein

Design and biological activities of novel inhibitory peptides for SARS-CoV spike protein and angiotensin-converting enzyme 2 interaction

Cellular entry of the SARS coronavirus

Structure-based discovery of a novel angiotensin-converting enzyme 2 inhibitor

Structural basis of neutralization by a human anti-severe acute respiratory syndrome spike protein antibody, 80R

Angiotensinconverting enzyme 2 protects from severe acute lung failure

SARS vaccine development

Identification of novel small-molecule inhibitors of severe acute respiratory syndromeassociated coronavirus by chemical genetics

The ubiquitin domain superfold: structure-based sequence alignments and characterization of binding epitopes

A detailed thermodynamic analysis of ras/effector complex interfaces

Recognizing and defining true Ras binding domains. II. In silico prediction based on homology modelling and energy calculations

A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury

Angiotensin-converting enzyme 2 in lung diseases

Angiotensin-converting enzyme 2: a functional receptor for SARS coronavirus

Pathogenesis of severe acute respiratory syndrome

Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus

Structure of SARS coronavirus spike receptor-binding domain complexed with receptor

Receptor and viral determinants of SARS-coronavirus adaptation to human ACE2

Conformational states of the severe acute respiratory syndrome coronavirus spike protein ectodomain

Surface ultrastructure of SARS coronavirus revealed by atomic force microscopy

Retroviruses pseudotyped with the severe acute respiratory syndrome coronavirus spike protein efficiently infect cells expressing angiotensin-converting enzyme 2

Structure of severe acute respiratory syndrome coronavirus receptor-binding domain complexed with neutralizing antibody

The FoldX web server: an online force field 10

Potent neutralization of severe acute respiratory syndrome (SARS) coronavirus by a human mAb to S1 protein that blocks receptor association

Evaluation of human monoclonal antibody 80R for immunoprophylaxis of severe acute respiratory syndrome by an animal study, epitope mapping, and analysis of spike variants

ACE2 X-ray structures reveal a large hinge-bending motion important for inhibitor binding and catalysis

ACE2: from vasopeptidase to SARS virus receptor

A 193-amino acid fragment of the SARS coronavirus S protein efficiently binds angiotensinconverting enzyme 2

Immunological responses against SARS-coronavirus infection in humans

Severe acute respiratory syndrome coronavirus entry into host cells: opportunities for therapeutic intervention

Immunogenicity of a receptor-binding domain of SARS coronavirus spike protein in mice: Implications for a subunit vaccine

A study on antigenicity and receptor-binding ability of fragment 450-650 of the spike protein of SARS coronavirus