Molecular Characterization of Adeno-Associated Viruses Infecting Children

Predicted VP1 capsid amino acid alignment of AAV2 and novel human AAVs. Diagram shows sequence alignment using the CLUSTAL W program. Black boxes designate amino acid substitutions compared to the AAV2 sequence. The locations of previously identified HVR regions (11) are labeled (HVR 1 to 12), as is an additional region (HVR 2′) that possesses several substitutions. Several HRV regions (5 to 7, 9, and 10) are colored to facilitate visualization of these regions onto the known atomic structure of AAV2, while invariant HVRs are labeled with black boxes. The locations of R585S and R588T are starred, and arrows denote the approximate locations of nested primers used to amplify the 255-bp HVR3 fragment.

Phylogenetic analysis of VP1 capsid nucleotide sequences. A neighbor-joining program with a Kimura two-parameter setting was used to derive phylogenetic distances based on 2,200 bp of VP1 sequence. Recently described AAV clade nomenclature (12) was adopted and organized by vertical brackets. The human isolates identified herein are designated in teal type. Due to space restrictions, only a few representative isolates from clades A, D, and E are shown. Sequence isolates are labeled with reference to the source species (bb, baboon; ch, chimpanzee; cy, cynomolgus macaque; hu, human; rh, rhesus macaque). Clade B sequences possessing R585 and R588 amino acids and predicted to bind HSPG efficiently are labeled in red type. The scale for genetic distance is indicated in the bottom left corner.

S17 sequence homology comparison with AAV2 and AAV3. Simplot analyses of similarity percentages of S17 VP1 versus AAV 2 (red) and AAV3 (blue) are shown. Data were plotted within a sliding window of 200 bp, centered on the position plotted, with a step size between data points of 20 bp. Positions containing gaps were excluded from the comparison. The bar on the top shows the predicted composition of the S17 capsid gene. The corresponding positions of the HVRs are labeled as magenta boxes (HVR 2′ in gray).

Surface diagrams of AAV2 trimer atomic models. (A) Electrostatic surface potential of the VP3 AAV2 trimer viewed down the threefold axis (yellow triangle) calculated with GRASP (23) running from negative (red) to positive (blue). Labeled arrows indicate the positions of residues implicated in HSPG binding. (B) Predicted electrostatic surface potential of AAV2 VP3 trimer as a result of R585S and R588T substitutions. Amino acid substitutions were modeled using energy minimization simulations with Quanta (Accelrys, San Diego, CA) prior to generating the electrostatic potential map in GRASP. The surface electrostatic potential scale is the same as depicted in panel A. Highlighted regions denote predicted HSPG coreceptor engagement domains in the VP3 trimer.

Ribbon diagrams of atomic models of AAV2 VP3 trimers showing the location of predicted amino acid substitutions in the human AAV isolates. (A) Ribbon drawing viewed down the threefold axis of symmetry of the AAV2 VP3 trimer. C_α backbones for the three VP3 monomers are rendered as teal ribbons. Predicted locations of the observed amino acid substitutions present within the eight AAV2-like sequences are color coded to reflect HVR location (HVR 5 to 7, 9, and 10) within the primary sequence (Fig. 2). White space-filling amino acid substitutions mapped outside the known HVRs. (B) Side view of the predicted location of the observed amino acid substitution demonstrating surface display (right side). (C) Superimposition of observed S17 amino acid substitutions relative to the AAV2 VP3 trimer atomic structure viewed down the threefold axis. (D) Side view of the predicted location of the observed amino acid substitutions in isolate S17 (surface display oriented on right side). Images were generated in NAMD/VMD (UIUC Theoretical Biophysics Group) and rendered using Raster3D.