Paramagnetism-Based Refinement Strategy for the Solution Structure of Human α-Parvalbumin

Protein MRI Contrast Agents as an Effective Approach for Precision Molecular Imaging

ABSTRACT

Cancer and other acute and chronic diseases are results of perturbations of common molecular determinants in key biological and signaling processes. Imaging is critical for characterizing dynamic changes in tumors and metastases, the tumor microenvironment, tumor-stroma interactions, and drug targets, at multiscale levels. Magnetic resonance imaging (MRI) has emerged to be a primary imaging modality for both clinical and preclinical applications due to its advantages over other modalities, including sensitivity to soft tissues, nondepth limitations, and the use of nonionizing radiation. However, extending the application of MRI to achieve both qualitative and quantitative precise molecular imaging with the capability to quantify molecular biomarkers for early detection, staging, and monitoring therapeutic treatment requires the capacity to overcome several major challenges including the trade-off between metal-binding affinity and relaxivity, which is an issue frequently associated with small chelator contrast agents. In this review, we will introduce the criteria of ideal contrast agents for precision molecular imaging and discuss the relaxivity of current contrast agents with defined first shell coordination water molecules. We will then report our advances in creating a new class of protein-targeted MRI contrast agents (ProCAs) with contributions to relaxivity largely derived from the secondary sphere and correlation time. We will summarize our rationale, design strategy, and approaches to the development and optimization of our pioneering ProCAs with desired high relaxivity, metal stability, and molecular biomarker-targeting capability, for precision MRI. From first generation (ProCA1) to third generation (ProCA32), we have achieved dual high r1 and r2 values that are 6- to 10-fold higher than clinically approved contrast agents at magnetic fields of 1.5 T, and their relaxivity values at high field are also significantly higher, which enables high resolution during small animal imaging. Further engineering of multiple targeting moieties enables ProCA32 agents that have strong biomarker-binding affinity and specificity for an array of key molecular biomarkers associated with various chronic diseases, while maintaining relaxation and exceptional metal-binding and selectivity, serum stability, and resistance to transmetallation, which are critical in mitigating risks associated with metal toxicity. Our leading product ProCA32.collagen has enabled the first early detection of liver metastasis from multiple cancers at early stages by mapping the tumor environment and early stage of fibrosis from liver and lung in vivo, with strong translational potential to extend to precision MRI for preclinical and clinical applications for precision diagnosis and treatment.

Pseudocontact Shifts in Biomolecular NMR Spectroscopy

Paramagnetic centers in biomolecules, such as specific metal ions that are bound to a protein, affect the nuclei in their surrounding in various ways. One of these effects is the pseudocontact shift (PCS), which leads to strong chemical shift perturbations of nuclear spins, with a remarkably long range of 50 Å and beyond. The PCS in solution NMR is an effect originating from the anisotropic part of the dipole-dipole interaction between the magnetic momentum of unpaired electrons and nuclear spins. The PCS contains spatial information that can be exploited in multiple ways to characterize structure, function, and dynamics of biomacromolecules. It can be used to refine structures, magnify effects of dynamics, help resonance assignments, allows for an intermolecular positioning system, and gives structural information in sensitivity-limited situations where all other methods fail. Here, we review applications of the PCS in biomolecular solution NMR spectroscopy, starting from early works on natural metalloproteins, following the development of non-natural tags to chelate and attach lanthanoid ions to any biomolecular target to advanced applications on large biomolecular complexes and inside living cells. We thus hope to not only highlight past applications but also shed light on the tremendous potential the PCS has in structural biology.

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.

Short two-armed lanthanide-binding tags for paramagnetic NMR spectroscopy based on chiral 1,4,7,10-tetrakis(2-hydroxypropyl)-1,4,7,10-tetraazacyclododecane scaffolds

A new pair of enantiomeric two-armed lanthanide-binding tags have been developed for paramagnetic NMR studies of proteins. The tags produce large and significantly different paramagnetic effects to one another when bound to the same tagging site. Additionally, they are less sensitive to sample pH than our previous two-armed tag designs.

Compact, hydrophilic, lanthanide-binding tags for paramagnetic NMR spectroscopy

The design, synthesis and evaluation of four novel lanthanide-binding tags for paramagnetic NMR spectroscopy are reported.

The design, synthesis and evaluation of four novel lanthanide-binding tags for paramagnetic NMR spectroscopy are reported. Each tag is based on the ((2S,2′S,2′′S,2′′′S)-1,1′,1′′,1′′′-(1,4,7,10-tetraazacyclododecane-1,4,7,10-tetrayl)tetrakis(propan-2-ol)) scaffold, featuring small chiral alcohol coordinating pendants to minimise the size and hydrophobic character of each tag. The tags feature different linkers of variable length for conjugation to protein via a single cysteine residue. Each tag's ability to induce pseudocontact shifts (PCS) was assessed on a ubiquitin A28C mutant. Two enantiomeric tags of particular note, C7 and C8, produced significantly larger Δχ-tensors compared to a previously developed tag, C1, attributed to the extremely short linker utilised, limiting the mobility of the bound lanthanide ion. The C7 and C8 tags' capacity to induce PCSs was further demonstrated on GB1 Q32C and 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) S112C/C80A mutants. Whilst factors such as the choice of lanthanide ion, pH and site of conjugation influence the size of the PCSs obtained, the tags represent a significant advance in the field.

Computational design of binding proteins to EGFR domain II

We developed a process to produce novel interactions between two previously unrelated proteins. This process selects protein scaffolds and designs protein interfaces that bind to a surface patch of interest on a target protein. Scaffolds with shapes complementary to the target surface patch were screened using an exhaustive computational search of the human proteome and optimized by directed evolution using phage display. This method was applied to successfully design scaffolds that bind to epidermal growth factor receptor (EGFR) domain II, the interface of EGFR dimerization, with high reactivity toward the target surface patch of EGFR domain II. One potential application of these tailor-made protein interactions is the development of therapeutic agents against specific protein targets.

Protein structure determination from pseudocontact shifts using ROSETTA

Paramagnetic metal ions generate pseudocontact shifts (PCSs) in nuclear magnetic resonance spectra that are manifested as easily measurable changes in chemical shifts. Metals can be incorporated into proteins through metal binding tags, and PCS data constitute powerful long-range restraints on the positions of nuclear spins relative to the coordinate system of the magnetic susceptibility anisotropy tensor (Δχ-tensor) of the metal ion. We show that three-dimensional structures of proteins can reliably be determined using PCS data from a single metal binding site combined with backbone chemical shifts. The program PCS-ROSETTA automatically determines the Δχ-tensor and metal position from the PCS data during the structure calculations, without any prior knowledge of the protein structure. The program can determine structures accurately for proteins of up to 150 residues, offering a powerful new approach to protein structure determination that relies exclusively on readily measurable backbone chemical shifts and easily discriminates between correctly and incorrectly folded conformations.

Distant homologs of anti-apoptotic factor HAX1 encode parvalbumin-like calcium binding proteins

Background

Apoptosis is a highly ordered and orchestrated multiphase process controlled by the numerous cellular and extra-cellular signals, which executes the programmed cell death via release of cytochrome c alterations in calcium signaling, caspase-dependent limited proteolysis and DNA fragmentation. Besides the general modifiers of apoptosis, several tissue-specific regulators of this process were identified including HAX1 (HS-1 associated protein X-1) - an anti-apoptotic factor active in myeloid cells. Although HAX1 was the subject of various experimental studies, the mechanisms of its action and a functional link connected with the regulation of apoptosis still remains highly speculative.

Findings

Here we provide the data which suggests that HAX1 may act as a regulator or as a sensor of calcium. On the basis of iterative similarity searches, we identified a set of distant homologs of HAX1 in insects. The applied fold recognition protocol gives us strong evidence that the distant insects' homologs of HAX1 are novel parvalbumin-like calcium binding proteins. Although the whole three EF-hands fold is not preserved in vertebrate our analysis suggests that there is an existence of a potential single EF-hand calcium binding site in HAX1. The molecular mechanism of its action remains to be identified, but the risen hypothesis easily translates into previously reported lines of various data on the HAX1 biology as well as, provides us a direct link to the regulation of apoptosis. Moreover, we also report that other family of myeloid specific apoptosis regulators - myeloid leukemia factors (MLF1, MLF2) share the homologous C-terminal domain and taxonomic distribution with HAX1.

Conclusions

Performed structural and active sites analyses gave new insights into mechanisms of HAX1 and MLF families in apoptosis process and suggested possible role of HAX1 in calcium-binding, still the analyses require further experimental verification.

Comparative modelling by restraint-based conformational sampling

BACKGROUND

Although comparative modelling is routinely used to produce three-dimensional models of proteins, very few automated approaches are formulated in a way that allows inclusion of restraints derived from experimental data as well as those from the structures of homologues. Furthermore, proteins are usually described as a single conformer, rather than an ensemble that represents the heterogeneity and inaccuracy of experimentally determined protein structures. Here we address these issues by exploring the application of the restraint-based conformational space search engine, RAPPER, which has previously been developed for rebuilding experimentally defined protein structures and for fitting models to electron density derived from X-ray diffraction analyses.

RESULTS

A new application of RAPPER for comparative modelling uses positional restraints and knowledge-based sampling to generate models with accuracies comparable to other leading modelling tools. Knowledge-based predictions are based on geometrical features of the homologous templates and rules concerning main-chain and side-chain conformations. By directly changing the restraints derived from available templates we estimate the accuracy limits of the method in comparative modelling.

CONCLUSION

The application of RAPPER to comparative modelling provides an effective means of exploring the conformational space available to a target sequence. Enhanced methods for generating positional restraints can greatly improve structure prediction. Generation of an ensemble of solutions that are consistent with both target sequence and knowledge derived from the template structures provides a more appropriate representation of a structural prediction than a single model. By formulating homologous structural information as sets of restraints we can begin to consider how comparative models might be used to inform conformer generation from sparse experimental data.

Solution structure of Ca2+-free rat beta-parvalbumin (oncomodulin)

Relative to other parvalbumin isoforms, the mammalian beta-parvalbumin (oncomodulin) displays attenuated divalent ion affinity. High-resolution structural data for the Ca(2+)-bound protein have provided little insight into the physical basis for this behavior, prompting an examination of the unliganded state. This article describes the solution structure and peptide backbone dynamics of Ca(2+)-free rat beta-parvalbumin (beta-PV). Ca(2+) removal evidently provokes significant structural alterations. Interaction between the D helix and the AB domain in the Ca(2+)-bound protein is greatly diminished in the apo-form, permitting the D helix to straighten. There is also a significant reorganization of the hydrophobic core and a concomitant remodeling of the interface between the AB and CD-EF domains. These modifications perturb the orientation of the C and D helices, and the energetic penalty associated with their reversal could contribute to the low-affinity signature of the CD site. By contrast, Ca(2+) removal causes a comparatively minor perturbation of the E and F helices, consistent with the more typical divalent ion affinity observed for the EF site. Ca(2+)-free rat beta-PV retains structural rigidity on the picosecond-nanosecond timescale. At 20 degrees C, the majority of amide vectors show no evidence for motion on timescales above 20 ps, and the average order parameter for the entire molecule is 0.92.

Strategies for measurements of pseudocontact shifts in protein NMR spectroscopy

Paramagnetic metal ions bound to proteins generate a dipolar field that can be accurately probed by pseudocontact shifts (PCS) displayed by the protein's nuclear spins. PCS are highly useful for determining the coordinates of individual spins in the molecule and for rapid structure determinations of entire protein-protein and protein-ligand complexes. However, PCS measurements require reliable resonance assignments for the molecule in its paramagnetic state and in a diamagnetic reference state. This article discusses different approaches for pairwise resonance assignments, with emphasis on a strategy which exploits chemical exchange between the two states.

Paramagnetism-based NMR restraints provide maximum allowed probabilities for the different conformations of partially independent protein domains

An innovative analytical/computational approach is presented to provide maximum allowed probabilities (MAPs) of conformations in protein domains not rigidly connected. The approach is applied to calmodulin and to its adduct with alpha-synuclein. Calmodulin is a protein constituted by two rigid domains, each of them composed by two calcium-binding EF-hand motifs, which in solution are largely free to move with respect to one another. We used the N60D mutant of calmodulin, which had been engineered to selectively bind a paramagnetic lanthanide ion to only one of its four calcium binding sites, specifically in the second EF-hand motif of the N-terminal domain. In this way, pseudocontact shifts (pcs's) and self-orientation residual dipolar couplings (rdc's) measured on the C-terminal domain provide information on its relative mobility with respect to the domain hosting the paramagnetic center. Available NMR data for terbium(III) and thulium(III) calmodulin were supplemented with additional data for dysprosium(III), analogous data were generated for the alpha-synuclein adduct, and the conformations with the largest MAPs were obtained for both systems. The MAP analysis for calmodulin provides further information on the variety of conformations experienced by the system. Such variety is somewhat reduced in the calmodulin-alpha-synuclein adduct, which however still retains high flexibility. The flexibility of the calmodulin-alpha-synuclein adduct is an unexpected result of this research.

Long-Range-Distance NMR Effects in a Protein Labeled with a Lanthanide–DOTA Chelate

A two-thiol reactive lanthanide-DOTA (1,4,7,10-tetraazacyclododecane-N,N',N'',N'''-tetraacetic acid) chelate, CLaNP-3 (CLaNP=caged lanthanide NMR probe), was synthesized for the rigid attachment to cysteine groups on a protein surface, and used to obtain long-range-distance information from the {15N,1H} HSQC spectra of the protein-lanthanide complex. The DOTA ring exhibits several isomers that are in exchange; however, single resonances were observed for most amide groups in the protein, allowing determination of a single, apparent magnetic-susceptibility tensor. Pseudocontact shifts caused by Yb-containing CLaNP-3 were observed for atoms at 15-35 A from the metal. By using Gd-containing CLaNP-3, relaxation effects were observed, allowing distances up to 30 A from the paramagnetic center to be determined accurately. Similar results were obtained with a Gd-DTPA (diethylene-triaminepentaacetic acid) chelate, CLaNP-1, bound in the same bidentate manner to the protein. This study demonstrates that bidentate attachment of a paramagnetic probe enables determination of long-range distances.

Observation of microsecond time-scale protein dynamics in the presence of Ln3+ ions: application to the N-terminal domain of cardiac troponin C

The microsecond time-scale motions in the N-terminal domain of cardiac troponin C (NcTnC) loaded with lanthanide ions have been investigated by means of a (1)H(N) off-resonance spin-lock experiment. The observed relaxation dispersion effects strongly increase along the series of NcTnC samples containing La(3+), Ce(3+), and Pr(3+) ions. This rise in dispersion effects is due to modulation of long-range pseudocontact shifts by micros time-scale dynamics. Specifically, the motion in the coordination sphere of the lanthanide ion (i.e. in the NcTnC EF-hand motif) causes modulation of the paramagnetic susceptibility tensor which, in turn, causes modulation of pseudocontact shifts. It is also probable that opening/closing dynamics, previously identified in Ca(2+)-NcTnC, contributes to some of the observed dispersions. On the other hand, it is unlikely that monomer-dimer exchange in the solution of NcTnC is directly responsible for the dispersion effects. Finally, on-off exchange of the lanthanide ion does not seem to play any significant role. The amplification of dispersion effects by Ln(3+) ions is a potentially useful tool for studies of micros-ms motions in proteins. This approach makes it possible to observe the dispersions even when the local environment of the reporting spin does not change. This happens, for example, when the motion involves a 'rigid' structural unit such as individual alpha-helix. Even more significantly, the dispersions based on pseudocontact shifts offer better chances for structural characterization of the dynamic species. This method can be generalized for a large class of applications via the use of specially designed lanthanide-binding tags.

Assignment of paramagnetic (15)N-HSQC spectra by heteronuclear exchange spectroscopy

Paramagnetic metal ions in proteins provide a rich source of structural information, but the resonance assignments required to extract the information can be challenging. Here we demonstrate that paramagnetically shifted (15)N-HSQC cross-peaks can be assigned using N(Z)-exchange spectroscopy under conditions in which the paramagnetic form of the protein is in dynamic equilibrium with its diamagnetic form. Even slow exchange of specifically bound metal ions may be detected within the long lifetime of (15)N longitudinal magnetization of large proteins at high magnetic fields. Alternatively, the exchange can be accelerated using an excess of metal ions. In the resulting exchange spectra, paramagnetic (15)N resonances become visible for residues that are not directly observed in a conventional (15)N-HSQC spectrum due to paramagnetic (1)H(N) broadening. The experiments are illustrated by the 30 kDa lanthanide-binding epsilon186/theta complex of DNA polymerase III in the presence of sub-stoichiometric amounts of Dy(3+) or a mixture of Dy(3+) and La(3+).

"Four-dimensional" protein structures: examples from metalloproteins

The fact that an object, for example, a protein, possesses a three-dimensional structure seems an obvious concept. However, when the object is flexible, the concept is less obvious. Growing experimental data over several decades show that proteins are not rigid objects, but they may sample more or less wide ranges of different conformations. To stress this concept, we propose to call the range of sampled conformations the "fourth dimension" of the protein structure. Nuclear magnetic resonance is a precious technique to define this fourth dimension. Examples of conformational heterogeneity taken from the realm of metalloproteins and their functional implications are discussed.

Structural basis for the observed differential magnetic anisotropic tensorial values in calcium binding proteins

Lanthanide ions (Ln(3+)), which have ionic radii similar to those of Ca(2+), can displace the latter in a calcium binding protein, without affecting its tertiary structure. The paramagnetic Ln(3+) possesses large anisotropic magnetic susceptibilities and produce pseudocontact shifts (PCSs), which have r(-3) dependence. The PCS can be seen for spins as far as 45 A from the paramagnetic ion. They aid in structure refinement of proteins by providing long-range distance constraints. Besides, they can be used to determine the interdomain orientation in multidomain proteins. This is particularly important in the context of a calcium binding protein from Entamoeba histolytica (EhCaBP), which consists of two globular domains connected by a flexible linker region containing 8 residues. As a first step to obtain the interdomain orientation in EhCaBP, a suite of 2D and 3D heteronuclear experiments were recorded on EhCaBP by displacing calcium with Ce(3+), Ho(3+), Er(3+), Tm(3+), Dy(3+), and Yb(3+) ions in separate experiments, and the PCS of (1)H(N) and (15)N spins were measured. Such data have been used in the refinement of the individual domain structures of the protein in parallel with the calculation of the respective magnetic anisotropy tensorial values, which differ substantially (2.1-2.8 times) from what is found in other Ca(2+) binding loops. This study provides a structural basis for such variations in the magnetic anisotropy tensorial values.

EF-hand protein dynamics and evolution of calcium signal transduction: an NMR view

Calcium signaling, one of the most widespread signaling mechanisms in cells, is generally carried out by EF-hand proteins, characterized by a helix-loop-helix motif paired in functional domains. EF-hand proteins may be viewed as molecular switches activated by calcium concentration transients. The EF-hand structural database has grown to a point where meaningful inferences on the functional conformational rearrangements upon calcium binding can be made by comparing a fair number of pairs of end points, i.e., the structures of the apo and calcium-bound forms. More compact descriptors of the movement associated with calcium binding, in terms of principal component analysis of the six interhelical angles, have also become available. Dynamic information obtained by NMR, also with the aid of calcium substitution with paramagnetic lanthanides, is shedding light on the intrinsic amplitude of the conformational degrees of freedom sampled by the various members of the EF-hand superfamily, as well as on the time scales of the motions. Particularly, NMR of lanthanide derivatives helps in capturing long time scale motions. Both static and dynamic pictures reveal a large variety of behaviors. It is increasingly recognized that the EF-hand machinery has differentiated its behavior during evolution in several ways, e.g., by modifying one of the loops, by undergoing a further duplication after the initial motif duplication that originated the functional domain, or by acquiring the ability to dimerize.

Efficient chi-tensor determination and NH assignment of paramagnetic proteins

Anisotropic magnetic susceptibility tensors chi of paramagnetic metal ions are manifested in pseudocontact shifts, residual dipolar couplings, and other paramagnetic observables that present valuable long-range information for structure determinations of protein-ligand complexes. A program was developed for automatic determination of the chi-tensor anisotropy parameters and amide resonance assignments in proteins labeled with paramagnetic metal ions. The program requires knowledge of the three-dimensional structure of the protein, the backbone resonance assignments of the diamagnetic protein, and a pair of 2D 15N-HSQC or 3D HNCO spectra recorded with and without paramagnetic metal ion. It allows the determination of reliable chi-tensor anisotropy parameters from 2D spectra of uniformly 15N-labeled proteins of fairly high molecular weight. Examples are shown for the 185-residue N-terminal domain of the subunit epsilon from E. coli DNA polymerase III in complex with the subunit theta and La3+ in its diamagnetic and Dy3+, Tb3+, and Er3+ in its paramagnetic form.

NMR Spectroscopy of Paramagnetic Metalloproteins

This article deals with the solution structure determination of paramagnetic metalloproteins by NMR spectroscopy. These proteins were believed not to be suitable for NMR investigations for structure determination until a decade ago, but eventually novel experiments and software protocols were developed, with the aim of making the approach suitable for the goal and as user-friendly and safe as possible. In the article, we also give hints for the optimization of experiments with respect to each particular metal ion, with the aim of also providing a handy tool for nonspecialists. Finally, a section is dedicated to the significant progress made on 13C direct detection, which reduces the negative effects of paramagnetism and may constitute a new chapter in the whole field of NMR spectroscopy.

Experimentally exploring the conformational space sampled by domain reorientation in calmodulin

The conformational space sampled by the two-domain protein calmodulin has been explored by an approach based on four sets of NMR observables obtained on Tb(3+)- and Tm(3+)-substituted proteins. The observables are the pseudocontact shifts and residual dipolar couplings of the C-terminal domain when lanthanide substitution is at the N-terminal domain. Each set of observables provides independent information on the conformations experienced by the molecule. It is found that not all sterically allowed conformations are equally populated. Taking the N-terminal domain as the reference, the C-terminal domain preferentially resides in a region of space inscribed in a wide elliptical cone. The axis of the cone is tilted by approximately 30 degrees with respect to the direction of the N-terminal part of the interdomain helix, which is known to have a flexible central part in solution. The C-terminal domain also undergoes rotation about the axis defined by the C-terminal part of the interdomain helix. Neither the extended helix conformation initially observed in the solid state for free calcium calmodulin nor the closed conformation(s) adopted by calcium calmodulin either alone or in its adduct(s) with target peptide(s) is among the most preferred ones. These findings are unique, both in terms of structural information obtained on a biomolecule that samples multiple conformations and in terms of the approach developed to achieve the results. The same approach is in principle applicable to other multidomain proteins, as well as to multiple interaction modes between two macromolecular partners.