Validation of helical tilt angles in the solution NMR structure of the Z domain of Staphylococcal protein A by combined analysis of residual dipolar coupling and NOE data

Recent Advances on Affibody- and DARPin-Conjugated Nanomaterials in Cancer Therapy

Affibodies and designed ankyrin repeat proteins (DARPins) are synthetic proteins originally derived from the Staphylococcus aureus virulence factor protein A and the human ankyrin repeat proteins, respectively. The use of these molecules in healthcare has been recently proposed as they are endowed with biochemical and biophysical features heavily demanded to target and fight diseases, as they have a strong binding affinity, solubility, small size, multiple functionalization sites, biocompatibility, and are easy to produce; furthermore, impressive chemical and thermal stability can be achieved. especially when using affibodies. In this sense, several examples reporting on affibodies and DARPins conjugated to nanomaterials have been published, demonstrating their suitability and feasibility in nanomedicine for cancer therapy. This minireview provides a survey of the most recent studies describing affibody- and DARPin-conjugated zero-dimensional nanomaterials, including inorganic, organic, and biological nanoparticles, nanorods, quantum dots, liposomes, and protein- and DNA-based assemblies for targeted cancer therapy in vitro and in vivo.

Implications of the unfolded state in the folding energetics of heterogeneous-backbone protein mimetics

Sequence-encoded folding is the foundation of protein structure and is also possible in synthetic chains of artificial chemical composition. In natural proteins, the characteristics of the unfolded state are as important as those of the folded state in determining folding energetics. While much is known about folded structures adopted by artificial protein-like chains, corresponding information about the unfolded states of these molecules is lacking. Here, we report the consequences of altered backbone composition on the structure, stability, and dynamics of the folded and unfolded states of a compact helix-rich protein. Characterization through a combination of biophysical experiments and atomistic simulation reveals effects of backbone modification that depend on both the type of artificial monomers employed and where they are applied in sequence. In general, introducing artificial connectivity in a way that reinforces characteristics of the unfolded state ensemble of the prototype natural protein minimizes the impact of chemical changes on folded stability. These findings have implications in the design of protein mimetics and provide an atomically detailed picture of the unfolded state of a natural protein and artificial analogues under non-denaturing conditions.

Biophysical experiments and atomistic simulation reveal impacts of protein backbone alteration on the ensemble that defines the unfolded state. These effects have implications on folded stability of protein mimetics.

Integrative Protein Modeling in RosettaNMR from Sparse Paramagnetic Restraints

Computational methods to predict protein structure from nuclear magnetic resonance (NMR) restraints that only require assignment of backbone signals, hold great potential to study larger proteins. Ideally, computational methods designed to work with sparse data need to add atomic detail that is missing in the experimental restraints. We introduce a comprehensive framework into the Rosetta suite that uses NMR restraints derived from paramagnetic labeling. Specifically, RosettaNMR incorporates pseudocontact shifts, residual dipolar couplings, and paramagnetic relaxation enhancements. It continues to use backbone chemical shifts and nuclear Overhauser effect distance restraints. We assess RosettaNMR for protein structure prediction by folding 28 monomeric proteins and 8 homo-oligomeric proteins. Furthermore, the general applicability of RosettaNMR is demonstrated on two protein-protein and three protein-ligand docking examples. Paramagnetic restraints generated more accurate models for 85% of the benchmark proteins and, when combined with chemical shifts, sampled high-accuracy models (≤2Å) in 50% of the cases.

Automatic 13C chemical shift reference correction for unassigned protein NMR spectra

Poor chemical shift referencing, especially for ¹³C in protein Nuclear Magnetic Resonance (NMR) experiments, fundamentally limits and even prevents effective study of biomacromolecules via NMR, including protein structure determination and analysis of protein dynamics. To solve this problem, we constructed a Bayesian probabilistic framework that circumvents the limitations of previous reference correction methods that required protein resonance assignment and/or three-dimensional protein structure. Our algorithm named Bayesian Model Optimized Reference Correction (BaMORC) can detect and correct ¹³C chemical shift referencing errors before the protein resonance assignment step of analysis and without three-dimensional structure. By combining the BaMORC methodology with a new intra-peaklist grouping algorithm, we created a combined method called Unassigned BaMORC that utilizes only unassigned experimental peak lists and the amino acid sequence. Unassigned BaMORC kept all experimental three-dimensional HN(CO)CACB-type peak lists tested within ± 0.4 ppm of the correct ¹³C reference value. On a much larger unassigned chemical shift test set, the base method kept ¹³C chemical shift referencing errors to within ± 0.45 ppm at a 90% confidence interval. With chemical shift assignments, Assigned BaMORC can detect and correct ¹³C chemical shift referencing errors to within ± 0.22 at a 90% confidence interval. Therefore, Unassigned BaMORC can correct ¹³C chemical shift referencing errors when it will have the most impact, right before protein resonance assignment and other downstream analyses are started. After assignment, chemical shift reference correction can be further refined with Assigned BaMORC. These new methods will allow non-NMR experts to detect and correct ¹³C referencing error at critical early data analysis steps, lowering the bar of NMR expertise required for effective protein NMR analysis.

Continuous Interdomain Orientation Distributions Reveal Components of Binding Thermodynamics

The flexibility of biological macromolecules is an important structural determinant of function. Unfortunately, the correlations between different motional modes are poorly captured by discrete ensemble representations. Here, we present new ways to both represent and visualize correlated interdomain motions. Interdomain motions are determined directly from residual dipolar couplings, represented as a continuous conformational distribution, and visualized using the disk-on-sphere representation. Using the disk-on-sphere representation, features of interdomain motions, including correlations, are intuitively visualized. The representation works especially well for multidomain systems with broad conformational distributions.This analysis also can be extended to multiple probability density modes, using a Bingham mixture model. We use this new paradigm to study the interdomain motions of staphylococcal protein A, which is a key virulence factor contributing to the pathogenicity of Staphylococcus aureus. We capture the smooth transitions between important states and demonstrate the utility of continuous distribution functions for computing the reorientational components of binding thermodynamics. Such insights allow for the dissection of the dynamic structural components of functionally important intermolecular interactions.

3-2-1: Structural insights from stepwise shrinkage of a three-helix Fc-binding domain to a single helix

The well-studied B-domain from Staphylococcal protein A is a 59 amino acid three-helix bundle that binds the Fc portion of IgG with a dissociation constant of ~35 nM. The B-domain variant bearing a Gly to Ala mutation (=Z-domain) has been the subject of efforts to minimize a domain's size while retaining its function. We report X-ray crystallographic characterization of three steps in such a process using complexes with Fc: the full three-helix Z-domain, a 34 amino acid two-helix version called Z34C and a 13 amino acid single helix stabilized with an exo-helix tether, called LH1.

Detecting and accounting for multiple sources of positional variance in peak list registration analysis and spin system grouping

Peak lists derived from nuclear magnetic resonance (NMR) spectra are commonly used as input data for a variety of computer assisted and automated analyses. These include automated protein resonance assignment and protein structure calculation software tools. Prior to these analyses, peak lists must be aligned to each other and sets of related peaks must be grouped based on common chemical shift dimensions. Even when programs can perform peak grouping, they require the user to provide uniform match tolerances or use default values. However, peak grouping is further complicated by multiple sources of variance in peak position limiting the effectiveness of grouping methods that utilize uniform match tolerances. In addition, no method currently exists for deriving peak positional variances from single peak lists for grouping peaks into spin systems, i.e. spin system grouping within a single peak list. Therefore, we developed a complementary pair of peak list registration analysis and spin system grouping algorithms designed to overcome these limitations. We have implemented these algorithms into an approach that can identify multiple dimension-specific positional variances that exist in a single peak list and group peaks from a single peak list into spin systems. The resulting software tools generate a variety of useful statistics on both a single peak list and pairwise peak list alignment, especially for quality assessment of peak list datasets. We used a range of low and high quality experimental solution NMR and solid-state NMR peak lists to assess performance of our registration analysis and grouping algorithms. Analyses show that an algorithm using a single iteration and uniform match tolerances approach is only able to recover from 50 to 80% of the spin systems due to the presence of multiple sources of variance. Our algorithm recovers additional spin systems by reevaluating match tolerances in multiple iterations. To facilitate evaluation of the algorithms, we developed a peak list simulator within our nmrstarlib package that generates user-defined assigned peak lists from a given BMRB entry or database of entries. In addition, over 100,000 simulated peak lists with one or two sources of variance were generated to evaluate the performance and robustness of these new registration analysis and peak grouping algorithms.

A Conjugate Based on Anti-HER2 Diaffibody and Auristatin E Targets HER2-Positive Cancer Cells

Antibody-drug conjugates (ADCs) have recently emerged as efficient and selective cancer treatment therapeutics. Currently, alternative forms of drug carriers that can replace monoclonal antibodies are under intensive investigation. Here, a cytotoxic conjugate of an anti-HER2 (Human Epidermal Growth Factor Receptor 2) diaffibody with monomethyl-auristatin E (MMAE) is proposed as a potential anticancer therapeutic. The anti-HER2 diaffibody was based on the Z_HER2:4 affibody amino acid sequence. The anti-HER2 diaffibody has been expressed as a His-tagged protein in E. coli and purified by Ni-nitrilotriacetyl (Ni-NTA) agarose chromatography. The molecule was properly folded, and the high affinity and specificity of its interaction with HER2 was confirmed by surface plasmon resonance (SPR) and flow cytometry, respectively. The (Z_HER2:4)₂DCS-MMAE conjugate was obtained by coupling the maleimide group linked with MMAE to cysteines, which were introduced in a drug conjugation sequence (DCS). Cytotoxicity of the conjugate was evaluated using the 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide MTT assay and the xCELLigence Real-Time Cell Analyzer. Our experiments demonstrated that the conjugate delivered auristatin E specifically to HER2-positive tumor cells, which finally led to their death. These results indicate that the cytotoxic diaffibody conjugate is a highly potent molecule for the treatment of various types of cancer overexpressing HER2 receptors.

Multiscale conformational heterogeneity in staphylococcal protein a: possible determinant of functional plasticity

The Staphylococcus aureus virulence factor staphylococcal protein A (SpA) is a major contributor to bacterial evasion of the host immune system, through high-affinity binding to host proteins such as antibodies. SpA includes five small three-helix-bundle domains (E-D-A-B-C) separated by conserved flexible linkers. Prior attempts to crystallize individual domains in the absence of a binding partner have apparently been unsuccessful. There have also been no previous structures of tandem domains. Here we report the high-resolution crystal structures of a single C domain, and of two B domains connected by the conserved linker. Both structures exhibit extensive multiscale conformational heterogeneity, which required novel modeling protocols. Comparison of domain structures shows that helix1 orientation is especially heterogeneous, coordinated with changes in side chain conformational networks and contacting protein interfaces. This represents the kind of structural plasticity that could enable SpA to bind multiple partners.

An encodable lanthanide binding tag with reduced size and flexibility for measuring residual dipolar couplings and pseudocontact shifts in large proteins

Metal ions serve important roles in structural biology applications from long-range perturbations seen in magnetic resonance experiments to electron-dense signatures in X-ray crystallography data; however, the metal ion must be secured in a molecular framework to achieve the maximum benefit. Polypeptide-based lanthanide-binding tags (LBTs) represent one option that can be directly encoded within a recombinant protein expression construct. However, LBTs often exhibit significant mobility relative to the target molecule. Here we report the characterization of improved LBTs sequences for insertion into a protein loop. These LBTs were inserted to connect two parallel alpha helices of an immunoglobulin G (IgG)-binding Z domain platform. Variants A and B bound Tb(3+) with high affinity (0.70 and 0.13 μM, respectively) and displayed restricted LBT motion. Compared to the parent construct, the metal-bound A experienced a 2.5-fold reduction in tag motion as measured by magnetic field-induced residual dipolar couplings and was further studied in a 72.2 kDa complex with the human IgG1 fragment crystallizable (IgG1 Fc) glycoprotein. The appearance of both pseudo-contact shifts (-0.221 to 0.081 ppm) and residual dipolar couplings (-7.6 to 14.3 Hz) of IgG1 Fc resonances in the IgG1 Fc:(variant A:Tb(3+))2 complex indicated structural restriction of the LBT with respect to the Fc. These studies highlight the applicability of improved LBT sequences with reduced mobility to probe the structure of macromolecular systems.

Structure-based Design of Peptides with High Affinity and Specificity to HER2 Positive Tumors

To identify peptides with high affinity and specificity against human epidermal growth factor receptor 2 (HER2), a series of peptides were designed based on the structure of HER2 and its Z(HER2:342) affibody. By using a combination protocol of molecular dynamics modeling, MM/GBSA binding free energy calculations, and binding free energy decomposition analysis, two novel peptides with 27 residues, pep27 and pep27-24M, were successfully obtained. Immunocytochemistry and flow cytometry analysis verified that both peptides can specifically bind to the extracellular domain of HER2 protein at cellular level. The Surface Plasmon Resonance imaging (SPRi) analysis showed that dissociation constants (K D) of these two peptides were around 300 nmol/L. Furthermore, fluorescence imaging of peptides against nude mice xenografted with SKBR3 cells indicated that both peptides have strong affinity and high specificity to HER2 positive tumors.

Immobilization of immunoglobulin-G-binding domain of Protein A on a gold surface modified with biotin ligase

Protein A from Staphylococcus aureus specifically binds to the Fc region of immunoglobulin G (IgG) and is widely used as a scaffold for the immobilization of IgG antibodies on solid supports. It is known that the oriented immobilization of Protein A on solid supports enhances its antibody-binding capability in comparison with immobilization in a random manner. In the current work, we developed a novel method for the oriented immobilization of the IgG-binding domain of Protein A based on the biotinylation reaction from archaeon Sulfolobus tokodaii. Biotinylation from S. tokodaii has a unique property in that the enzyme, biotin protein ligase (BPL), forms a stable complex with its biotinylated substrate protein, biotin carboxyl carrier protein (BCCP). Here, BCCP was fused to the IgG-binding domain of Protein A, and the resulting fusion protein was immobilized on the BPL-modified gold surface of the sensor chip for quartz crystal microbalance through complexation between BCCP and BPL. The layer of the IgG-binding domain prepared in this way successfully captured the antibody, and the captured antibody retained high antigen-binding capability.

The statistical conformation of a highly flexible protein: small-angle X-ray scattering of S. aureus protein A

Staphylococcal protein A (SpA) is a multidomain protein consisting of five globular IgG binding domains separated by a conserved six- to nine-residue flexible linker. We collected SAXS data on the N-terminal protein-binding half of SpA (SpA-N) and constructs consisting of one to five domain modules in order to determine statistical conformation of this important S. aureus virulence factor. We fit the SAXS data to a scattering function based on a new polymer physics model, which provides an analytical description of the SpA-N statistical conformation. We describe a protocol for systematically determining the appropriate level of modeling to fit a SAXS data set based on goodness of fit and whether the addition of parameters improves it. In the case of SpA-N, the analytical polymer physics description provides a depiction of the statistical conformation of a flexible protein that, while lacking atomistic detail, properly reflects the information content of the data.

Artificial affinity proteins as ligands of immunoglobulins

A number of natural proteins are known to have affinity and specificity for immunoglobulins. Some of them are widely used as reagents for detection or capture applications, such as Protein G and Protein A. However, these natural proteins have a defined spectrum of recognition that may not fit specific needs. With the development of combinatorial protein engineering and selection techniques, it has become possible to design artificial affinity proteins with the desired properties. These proteins, termed alternative scaffold proteins, are most often chosen for their stability, ease of engineering and cost-efficient recombinant production in bacteria. In this review, we focus on alternative scaffold proteins for which immunoglobulin binders have been identified and characterized.

A practical implicit solvent potential for NMR structure calculation

The benefits of protein structure refinement in water are well documented. However, performing structure refinement with explicit atomic representation of the solvent molecules is computationally expensive and impractical for NMR-restrained structure calculations that start from completely extended polypeptide templates. Here we describe a new implicit solvation potential, EEFx (Effective Energy Function for XPLOR-NIH), for NMR-restrained structure calculations of proteins in XPLOR-NIH. The key components of EEFx are an energy term for solvation energy that works together with other nonbonded energy functions, and a dedicated force field for conformational and nonbonded protein interaction parameters. The initial results obtained with EEFx show that significant improvements in structural quality can be obtained. EEFx is computationally efficient and can be used both to fold and refine structures. Overall, EEFx improves the quality of protein conformation and nonbonded atomic interactions. Moreover, such benefits are accompanied by enhanced structural precision and enhanced structural accuracy, reflected in improved agreement with the cross-validated dipolar coupling data. Finally, implementation of EEFx calculations is straightforward and computationally efficient. Overall, EEFx provides a useful method for the practical calculation of experimental protein structures in a physically realistic environment.

Engineered solubility tag for solution NMR of proteins

The low solubility of many proteins hinders large scale expression and purification as well as biophysical measurements. Here, we devised a general strategy to solubilize a protein by conjugating it at a solvent-exposed position to a 6 kDa protein that was re-engineered to be highly soluble. We applied this method to the CARD domain of Apoptosis-associated speck-like protein containing a CARD (ASC), which represents one member of a class of proteins that are notoriously prone to aggregation. Attachment of the tag to a cysteine residue, introduced by site-directed mutagenesis at its self-association interface, improved the solubility of the ASC CARD over 50-fold under physiological conditions. Although it is not possible to use nuclear magnetic resonance (NMR) to obtain a high quality 2D correlation spectrum of the wild type domain under physiological conditions, we demonstrate that NMR relaxation parameters of the solubilized variant are sufficiently improved to facilitate virtually any demanding measurement. The method shown here represents a straightforward approach for dramatically increasing protein solubility, enabled by ease of labeling as well as flexibility in tag placement with minimal perturbation to the target.

Computational protein design with electrostatic focusing: experimental characterization of a conditionally folded helical domain with a reduced amino acid alphabet

Automated methodologies to design synthetic proteins from first principles use energy computations to estimate the ability of the sequences to adopt a targeted structure. This approach is still far from systematically producing native-like sequences, due, most likely, to inaccuracies when modeling the interactions between the protein and its aqueous environment. This is particularly challenging when engineering small protein domains (with less polar pair interactions than with the solvent). We have re-designed a three-helix bundle, domain B, using a fixed backbone and a four amino acid alphabet. We have enlarged the rotamer library with conformers that increase the weight of electrostatic interactions within the design process without altering the energy function used to compute the folding free energy. Our synthetic sequences show less than 15% similarity to any Swissprot sequence. We have characterized our sequences in different solvents using circular dichroism and nuclear magnetic resonance. The targeted structure achieved is dependent on the solvent used. This method can be readily extended to larger domains. Our method will be useful for the engineering of proteins that become active only in a given solvent and for designing proteins in the context of hydrophobic solvents, an important fraction of the situations in the cell.

Design of a modular protein-based MRI contrast agent for targeted application

Magnetic resonance imaging (MRI) offers a non-radioactive alternative for the non-invasive detection of tumours. Low molecular weight MRI contrast agents currently in clinical use suffer either from a lack of specificity for tumour tissue or from low relaxivity and thus low contrast amplification. In this study, we present the newly designed two domain fusion protein Zarvin, which is able to bind to therapeutic IgG antibodies suitable for targeting, while facilitating contrast enhancement through high affinity binding sites for Gd³⁺. We show that the Zarvin fold is stable under serum conditions, specifically targets a cancer cell-line when bound to the Cetuximab IgG, and allows for imaging with high relaxivity, a property that would be advantageous for the detection of small tumours and metastases at 1.5 or 3 T.

Interdomain orientation of cardiac troponin C characterized by paramagnetic relaxation enhancement NMR reveals a compact state

Cardiac troponin C (cTnC) is the calcium binding subunit of the troponin complex that triggers the thin filament response to calcium influx into the sarcomere. cTnC consists of two globular EF-hand domains (termed the N- and C-domains) connected by a flexible linker. While the conformation of each domain of cTnC has been thoroughly characterized through NMR studies involving either the isolated N-domain (N-cTnC) or C-domain (C-cTnC), little attention has been paid to the range of interdomain orientations possible in full-length cTnC that arises as a consequence of the flexibility of the domain linker. Flexibility in the domain linker of cTnC is essential for effective regulatory function of troponin. We have therefore utilized paramagnetic relaxation enhancement (PRE) NMR to assess the interdomain orientation of cTnC. Ensemble fitting of our interdomain PRE measurements reveals that isolated cTnC has considerable interdomain flexibility and preferentially adopts a bent conformation in solution, with a defined range of relative domain orientations.

Lanthanide binding and IgG affinity construct: potential applications in solution NMR, MRI, and luminescence microscopy

Paramagnetic lanthanide ions when bound to proteins offer great potential for structural investigations that utilize solution nuclear magnetic resonance spectroscopy, magnetic resonance imaging, or optical microscopy. However, many proteins do not have native metal ion binding sites and engineering a chimeric protein to bind an ion while retaining affinity for a protein of interest represents a significant challenge. Here we report the characterization of an immunoglobulin G-binding protein redesigned to include a lanthanide binding motif in place of a loop between two helices (Z-L2LBT). It was shown to bind Tb³⁺ with 130 nM affinity. Ions such as Dy³⁺, Yb³⁺, and Ce³⁺ produce paramagnetic effects on NMR spectra and the utility of these effects is illustrated by their use in determining a structural model of the metal-complexed Z-L2LBT protein and a preliminary characterization of the dynamic distribution of IgG Fc glycan positions. Furthermore, this designed protein is demonstrated to be a novel IgG-binding reagent for magnetic resonance imaging (Z-L2LBT:Gd³⁺ complex) and luminescence microscopy (Z-L2LBT: Tb³⁺ complex).

Selective N-terminal fluorescent labeling of proteins using 4-chloro-7-nitrobenzofurazan: a method to distinguish protein N-terminal acetylation

A fluorogenic derivatization method was developed to distinguish the protein N-terminal acetylation status. The unacetylated protein selectively reacted with 4-chloro-7-nitrobenzofurazan (NBD-Cl) at neutral pH to provide high fluorescence. In contrast, the protein with N-terminal acetylation was essentially nonfluorescent under the same conditions despite the presence of many internal lysine residues. Fluorescence of the NBD-labeled protein was very stable, and only micromolar concentrations of proteins were required for reliable detection. This method also provides a general and practical way to quantify proteins when their N-terminal amino group is available.

Engineering of affibody molecules for therapy and diagnostics

Affibody molecules are small and robust non-immunoglobulin affinity ligands capable of binding to a wide range of protein targets. They are selected from combinatorial libraries based on a 58 amino acid, three-alpha-helical Z-domain scaffold. They share no sequence or structural homologies to antibodies and in contrast to antibodies they can be functionally produced both by peptide synthesis and by recombinant expression in Escherichia coli. Protein engineering is used to adapt Affibody molecules binding to a target of interest to the specific demands imposed by the intended application. Obviously, the optimal molecule for molecular imaging will be different from the optimal molecule for therapy. Here, we describe general strategies to optimize Affibody molecules for diagnostic imaging and therapy applications.

Design of an optimized scaffold for affibody molecules

Affibody molecules are non-immunoglobulin-derived affinity proteins based on a three-helical bundle protein domain. Here, we describe the design process of an optimized Affibody molecule scaffold with improved properties and a surface distinctly different from that of the parental scaffold. The improvement was achieved by applying an iterative process of amino acid substitutions in the context of the human epidermal growth factor receptor 2 (HER2)-specific Affibody molecule Z(HER2:342). Replacements in the N-terminal region, loop 1, helix 2 and helix 3 were guided by extensive structural modeling using the available structures of the parent Z domain and Affibody molecules. The effect of several single substitutions was analyzed followed by combination of up to 11 different substitutions. The two amino acid substitutions N23T and S33K accounted for the most dramatic improvements, including increased thermal stability with elevated melting temperatures of up to +12 degrees C. The optimized scaffold contains 11 amino acid substitutions in the nonbinding surface and is characterized by improved thermal and chemical stability, as well as increased hydrophilicity, and enables generation of identical Affibody molecules both by chemical peptide synthesis and by recombinant bacterial expression. A HER2-specific Affibody tracer, [MMA-DOTA-Cys61]-Z(HER2:2891)-Cys (ABY-025), was produced by conjugating MMA-DOTA (maleimide-monoamide-1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid) to the peptide produced either chemically or in Escherichia coli. ABY-025 showed high affinity and specificity for HER2 (equilibrium dissociation constant, K(D), of 76 pM) and detected HER2 in tissue sections of SKOV-3 xenograft and human breast tumors. The HER2-binding capacity was fully retained after three cycles of heating to 90 degrees C followed by cooling to room temperature. Furthermore, the binding surfaces of five Affibody molecules targeting other proteins (tumor necrosis factor alpha, insulin, Taq polymerase, epidermal growth factor receptor or platelet-derived growth factor receptor beta) were grafted onto the optimized scaffold, resulting in molecules with improved thermal stability and a more hydrophilic nonbinding surface.

Synthesis and chemoselective intramolecular crosslinking of a HER2-binding affibody

The human epidermal growth factor receptor HER2 has emerged as an important target for molecular imaging of breast cancer. This article presents the design and synthesis of a HER2-targeting affibody molecule with improved stability and tumor targeting capacity, and with potential use as an imaging agent. The 58 aa three-helix bundle protein was assembled using solid-phase peptide synthesis, and a chemoselective ligation strategy was used to establish an intramolecular thioether bond between the side chain thiol group of a cysteine residue, positioned in the loop between helices I and II, and a chloroacetyl group on the side chain amino group of the C-terminal lysine residue. The tethered protein offered an increased thermal stability, with a melting temperature of 64 degrees C, compared to 54 degrees C for the linear control. The ligation did not have a major influence on the HER2 binding affinity, which was 320 and 380 pM for the crosslinked and linear molecules, respectively. Biodistribution studies were performed both in normal and tumor-bearing mice to evaluate the impact of the crosslinking on the in vivo behavior and on the tumor targeting performance. The distribution pattern was characterized by a low uptake in all organs except kidney, and rapid clearance from blood and normal tissue. Crosslinking of the protein resulted in a significantly increased tumor accumulation, rendering the tethered HER2-binding affibody molecule a valuable lead in the development of superior HER2 imaging agents.

KinImmerse: Macromolecular VR for NMR ensembles

Background

In molecular applications, virtual reality (VR) and immersive virtual environments have generally been used and valued for the visual and interactive experience – to enhance intuition and communicate excitement – rather than as part of the actual research process. In contrast, this work develops a software infrastructure for research use and illustrates such use on a specific case.

Methods

The Syzygy open-source toolkit for VR software was used to write the KinImmerse program, which translates the molecular capabilities of the kinemage graphics format into software for display and manipulation in the DiVE (Duke immersive Virtual Environment) or other VR system. KinImmerse is supported by the flexible display construction and editing features in the KiNG kinemage viewer and it implements new forms of user interaction in the DiVE.

Results

In addition to molecular visualizations and navigation, KinImmerse provides a set of research tools for manipulation, identification, co-centering of multiple models, free-form 3D annotation, and output of results. The molecular research test case analyzes the local neighborhood around an individual atom within an ensemble of nuclear magnetic resonance (NMR) models, enabling immersive visual comparison of the local conformation with the local NMR experimental data, including target curves for residual dipolar couplings (RDCs).

Conclusion

The promise of KinImmerse for production-level molecular research in the DiVE is shown by the locally co-centered RDC visualization developed there, which gave new insights now being pursued in wider data analysis.

Folding processes of the B domain of protein A to the native state observed in all-atom ab initio folding simulations

Reaching the native states of small proteins, a necessary step towards a comprehensive understanding of the folding mechanisms, has remained a tremendous challenge to ab initio protein folding simulations despite the extensive effort. In this work, the folding process of the B domain of protein A (BdpA) has been simulated by both conventional and replica exchange molecular dynamics using AMBER FF03 all-atom force field. Started from an extended chain, a total of 40 conventional (each to 1.0 micros) and two sets of replica exchange (each to 200.0 ns per replica) molecular dynamics simulations were performed with different generalized-Born solvation models and temperature control schemes. The improvements in both the force field and solvent model allowed successful simulations of the folding process to the native state as demonstrated by the 0.80 A C(alpha) root mean square deviation (RMSD) of the best folded structure. The most populated conformation was the native folded structure with a high population. This was a significant improvement over the 2.8 A C(alpha) RMSD of the best nativelike structures from previous ab initio folding studies on BdpA. To the best of our knowledge, our results demonstrate, for the first time, that ab initio simulations can reach the native state of BdpA. Consistent with experimental observations, including Phi-value analyses, formation of helix II/III hairpin was a crucial step that provides a template upon which helix I could form and the folding process could complete. Early formation of helix III was observed which is consistent with the experimental results of higher residual helical content of isolated helix III among the three helices. The calculated temperature-dependent profile and the melting temperature were in close agreement with the experimental results. The simulations further revealed that phenylalanine 31 may play critical to achieve the correct packing of the three helices which is consistent with the experimental observation. In addition to the mechanistic studies, an ab initio structure prediction was also conducted based on both the physical energy and a statistical potential. Based on the lowest physical energy, the predicted structure was 2.0 A C(alpha) RMSD away from the experimentally determined structure.

NMR solution structure and backbone dynamics of domain III of the E protein of tick-borne Langat flavivirus suggests a potential site for molecular recognition

Flaviviruses cause many human diseases, including dengue fever, yellow fever, West Nile viral encephalitis, and hemorrhagic fevers, and are transmitted to their vertebrate hosts by infected mosquitoes and ticks. Domain III of the envelope protein (E-D3) is considered to be the primary viral determinant involved in the virus-host-cell receptor interaction, and thus represents an excellent target for antiviral drug development. Langat (LGT) virus is a naturally attenuated BSL-2 TBE virus and is a model for the pathogenic BSL-3 and BSL-4 viruses in the serogroup. We have determined the solution structure of LGT-E-D3 using heteronuclear NMR spectroscopy. The backbone dynamics of LGT-E-D3 have been investigated using 15N relaxation measurements. A detailed analysis of the solution structure and dynamics of LGT-E-D3 suggests potential residues that could form a surface for molecular recognition, and thereby represent a target site for antiviral therapeutics design.

Thermodynamics of Folding and Binding in an Affibody:Affibody Complex

Affibody binding proteins are selected from phage-displayed libraries of variants of the 58 residue Z domain. Z(Taq) is an affibody originally selected as a binder to Taq DNA polymerase. The anti-Z(Taq) affibody was selected as a binder to Z(Taq) and the Z(Taq):anti-Z(Taq) complex is formed with a dissociation constant K(d)=0.1 microM. We have determined the structure of the Z(Taq):anti-Z(Taq) complex as well as the free state structures of Z(Taq) and anti-Z(Taq) using NMR. Here we complement the structural data with thermodynamic studies of Z(Taq) and anti-Z(Taq) folding and complex formation. Both affibody proteins show cooperative two-state thermal denaturation at melting temperatures T(M) approximately 56 degrees C. Z(Taq):anti-Z(Taq) complex formation at 25 degrees C in 50 mM NaCl and 20 mM phosphate buffer (pH 6.4) is enthalpy driven with DeltaH degrees (bind) = -9.0 (+/-0.1) kcal mol(-1)(.) The heat capacity change DeltaC(P) degrees (,bind)=-0.43 (+/-0.01) kcal mol(-1) K(-1) is in accordance with the predominantly non-polar character of the binding surface, as judged from calculations based on changes in accessible surface areas. A further dissection of the small binding entropy at 25 degrees C (-TDeltaS degrees (bind) = -0.6 (+/-0.1) kcal mol(-1)) suggests that a favourable desolvation of non-polar surface is almost completely balanced by unfavourable conformational entropy changes and loss of rotational and translational entropy. Such effects can therefore be limiting for strong binding also when interacting protein components are stable and homogeneously folded. The combined structure and thermodynamics data suggest that protein properties are not likely to be a serious limitation for the development of engineered binding proteins based on the Z domain.

Structural basis for molecular recognition in an affibody:affibody complex

Affibody molecules constitute a class of engineered binding proteins based on the 58-residue three-helix bundle Z domain derived from staphylococcal protein A (SPA). Affibody proteins are selected as binders to target proteins by phage display of combinatorial libraries in which typically 13 side-chains on the surface of helices 1 and 2 in the Z domain have been randomized. The Z(Taq):anti-Z(Taq) affibody-affibody complex, consisting of Z(Taq), originally selected as a binder to Taq DNA polymerase, and anti-Z(Taq), selected as binder to Z(Taq), is formed with a dissociation constant K(d) approximately 100 nM. We have determined high-precision solution structures of free Z(Taq) and anti-Z(Taq), and the Z(Taq):anti-Z(Taq) complex under identical experimental conditions (25 degrees C in 50 mM NaCl with 20 mM potassium phosphate buffer at pH 6.4). The complex is formed with helices 1 and 2 of anti-Z(Taq) in perpendicular contact with helices 1 and 2 of Z(Taq). The interaction surface is large ( approximately 1670 A(2)) and unusually non-polar (70 %) compared to other protein-protein complexes. It involves all varied residues on anti-Z(Taq), most corresponding (Taq DNA polymerase binding) side-chains on Z(Taq), and several additional side-chain and backbone contacts. Other notable features include a substantial rearrangement (induced fit) of aromatic side-chains in Z(Taq) upon binding, a close contact between glycine residues in the two subunits that might involve aliphatic glycine Halpha to backbone carbonyl hydrogen bonds, and four hydrogen bonds made by the two guanidinium N(eta)H(2) groups of an arginine side-chain. Comparisons of the present structure with other data for affibody proteins and the Z domain suggest that intrinsic binding properties of the originating SPA surface might be inherited by the affibody binders. A thermodynamic characterization of Z(Taq) and anti-Z(Taq) is presented in an accompanying paper.

Rapid NMR data collection

Rapid data collection is an area of intense research in biomolecular NMR spectroscopy, in particular for high-throughput structure determination in structural genomics. NMR data acquisition and processing protocols for rapidly obtaining high-dimensional spectral information aim at avoiding sampling limited data collection and are reviewed here with emphasis on G-matrix Fourier transform NMR spectroscopy.

Determination of protein backbone structures from residual dipolar couplings

There are a number of circumstances in which a focus on determination of the backbone structure of a protein, as opposed to a complete all-atom structure, may be appropriate. This is particularly the case for structures determined as a part of a structural genomics initiative in which computational modeling of many sequentially related structures from the backbone of a single family representative is anticipated. It is, however, also the case when the backbone may be a stepping-stone to more targeted studies of ligand interaction or protein-protein interaction. Here an NMR protocol is described that can produce a backbone structure of a protein without the need for extensive experiments directed at side chain resonance assignment or the collection of structural information on side chains. The procedure relies primarily on orientational constraints from residual dipolar couplings as opposed to distance constraints from NOEs. Procedures for sample preparation, data acquisition, and data analysis are described, along with examples from application to small target proteins of a structural genomics project.