Identification of GB3 as a Novel Biomarker of Tumor-Derived Vasculature in Neuroblastoma Using a Stiffness-Based Model

Neuroblastoma (NB) is a childhood cancer in sympathetic nervous system cells. NB exhibits cellular heterogeneity, with adrenergic and mesenchymal states displaying distinct tumorigenic potentials. NB is highly vascularized, and blood vessels can form through various mechanisms, including endothelial transdifferentiation, leading to the development of tumor-derived endothelial cells (TECs) associated with chemoresistance. We lack specific biomarkers for TECs. Therefore, identifying new TEC biomarkers is vital for effective NB therapies. A stiffness-based platform simulating human arterial and venous stiffness was developed to study NB TECs in vitro. Adrenergic cells cultured on arterial-like stiffness transdifferentiated into TECs, while mesenchymal state cells did not. The TECs derived from adrenergic cells served as a model to explore new biomarkers, with a particular focus on GB3, a glycosphingolipid receptor implicated in angiogenesis, metastasis, and drug resistance. Notably, the TECs unequivocally expressed GB3, validating its novelty as a marker. To explore targeted therapeutic interventions, nanoparticles functionalized with the non-toxic subunit B of the Shiga toxin were generated, because they demonstrated a robust affinity for GB3-positive cells. Our results demonstrate the value of the stiffness-based platform as a predictive tool for assessing NB aggressiveness, the discovery of new biomarkers, and the evaluation of the effectiveness of targeted therapeutic strategies.

Correlated motions of C'-N and Cα-Cβ pairs in protonated and per-deuterated GB3

We investigated correlated µs-ms time scale motions of neighboring ¹³C'-¹⁵N and ¹³C_α-¹³C_β nuclei in both protonated and perdeuterated samples of GB3. The techniques employed, NMR relaxation due to cross-correlated chemical shift modulations, specifically target concerted changes in the isotropic chemical shifts of the two nuclei associated with spatial fluctuations. Field-dependence of the relaxation rates permits identification of the parameters defining the chemical exchange rate constant under the assumption of a two-site exchange. The time scale of motions falls into the intermediate to fast regime (with respect to the chemical shift time scale, 100-400 s^-1 range) for the ¹³C'-¹⁵N pairs and into the slow to intermediate regime for the ¹³C_α-¹³C_β pairs (about 150 s^-1). Comparison of the results obtained for protonated and deuterated GB3 suggests that deuteration has a tendency to reduce these slow scale correlated motions, especially for the ¹³C_α-¹³C_β pairs.

A Liquid Chromatography-Quadrupole-Time-of-Flight Mass Spectrometric Assay for the Quantification of Fabry Disease Biomarker Globotriaosylceramide (GB3) in Fabry Model Mouse

Fabry disease is a rare lysosomal storage disorder resulting from the lack of α-Gal A gene activity. Globotriaosylceramide (GB3, ceramide trihexoside) is a novel endogenous biomarker which predicts the incidence of Fabry disease. At the early stage efficacy/biomarker study, a rapid method to determine this biomarker in plasma and in all relevant tissues related to this disease simultaneously is required. However, the limited sample volume, as well as the various levels of GB3 in different matrices makes the GB3 quantitation very challenging. Hereby we developed a rapid method to identify GB3 in mouse plasma and various tissues. Preliminary stability tests were also performed in three different conditions: short-term, freeze-thaw, long-term. The calibration curve was well fitted over the concentration range of 0.042–10 μg/mL for GB3 in plasma and 0.082–20 μg/g for GB3 in various tissues. This method was successfully applied for the comparison of GB3 levels in Fabry model mice (B6;129-Gla^tm1Kul/J), which has not been performed previously to the best of our knowledge.

A Force Balanced Fragmentation Method for ab Initio Molecular Dynamic Simulation of Protein

A force balanced generalized molecular fractionation with conjugate caps (FB-GMFCC) method is proposed for ab initio molecular dynamic simulation of proteins. In this approach, the energy of the protein is computed by a linear combination of the QM energies of individual residues and molecular fragments that account for the two-body interaction of hydrogen bond between backbone peptides. The atomic forces on the caped H atoms were corrected to conserve the total force of the protein. Using this approach, ab initio molecular dynamic simulation of an Ace-(ALA)₉-NME linear peptide showed the conservation of the total energy of the system throughout the simulation. Further a more robust 110 ps ab initio molecular dynamic simulation was performed for a protein with 56 residues and 862 atoms in explicit water. Compared with the classical force field, the ab initio molecular dynamic simulations gave better description of the geometry of peptide bonds. Although further development is still needed, the current approach is highly efficient, trivially parallel, and can be applied to ab initio molecular dynamic simulation study of large proteins.

Bidirectional band-selective magnetization transfer along the protein backbone doubles the information content of solid-state NMR correlation experiments

Resonance assignment is the first stage towards solving the structure of a protein. This is normally achieved by the employment of separate inter and intra residue experiments. By utilising the mixed rotation and rotary recoupling (MIRROR) condition it is possible to double the information content through the efficient bidirectional transfer of magnetization from the CO to its adjacent Cα and the Cα of the subsequent amino acid. We have incorporated this into a 3D experiment, a 3D-MIRROR-NCOCA, where correlations present in the 3D spectrum permit the sequential assignment of the protein backbone from a single experiment as we have demonstrated on a microcrystalline preparation of GB3. Furthermore, the low-power requirements of the MIRROR recoupling sequence facilitate the development of a low-power 3D-NCOCA experiment. This has enabled us to realise significant reductions in acquisition times, allowing the acquisition of a single 3D-NCOCA spectrum suitable for a full backbone resonance assignment of GB3 in less than 24 h.

Detection of Correlated Protein Backbone and Side-Chain Angle Fluctuations

NMR methods for the characterization of local protein motions have attained a high level of sophistication. Measurement of the synchronization between those motions, however, poses a serious challenge. Such correlated motions are one of the underlying mechanisms for the propagation of local changes to remote sites and as such for information transfer. Here, we demonstrate the experimental detection of the synchronization of motion over an intermediate range. To that purpose, we designed pulse sequences for the measurement of cross-correlated relaxation between the backbone H^N -N and side-chain H^β -C^β dipoles in Ile, Thr, and Val in the protein GB3. These bonds are related through two and three intervening dihedral angles. We show that the correlated motions inherent in a structural ensemble obtained from a large and diverse array of NMR probes are in excellent agreement with our measurements.

Cross-correlated relaxation rates between protein backbone H-X dipolar interactions

The relaxation interference between dipole-dipole interactions of two separate spin pairs carries structural and dynamics information. In particular, when compared to individual dynamic behavior of those spin pairs, such cross-correlated relaxation (CCR) rates report on the correlation between the spin pairs. We have recently mapped out correlated motion along the backbone of the protein GB3, using CCR rates among and between consecutive H^N-N and H^α-C^α dipole-dipole interactions. Here, we provide a detailed account of the measurement of the four types of CCR rates. All rates were obtained from at least two different pulse sequences, of which the yet unpublished ones are presented. Detailed comparisons between the different methods and corrections for unwanted pathways demonstrate that the averaged CCR rates are highly accurate and precise with errors of 1.5-3% of the entire value ranges.

Quantifying protein dynamics in the ps-ns time regime by NMR relaxation

Both ¹⁵N chemical shift anisotropy (CSA) and sufficiently rapid exchange linebroadening transitions exhibit relaxation contributions that are proportional to the square of the magnetic field. Deconvoluting these contributions is further complicated by residue-dependent variations in protein amide ¹⁵N CSA values which have proven difficult to accurately measure. Exploiting recently reported improvements for the implementation of T₁ and T_1ρ experiments, field strength-dependent studies have been carried out on the B3 domain of protein G (GB3) as well as on the immunophilin FKBP12 and a H87V variant of that protein in which the major conformational exchange linebroadening transition is suppressed. By applying a zero frequency spectral density rescaling analysis to the relaxation data collected at magnetic fields from 500 to 900 MHz ¹H, differential residue-specific ¹⁵N CSA values have been obtained for GB3 which correlate with those derived from solid state and liquid crystalline NMR measurements to a level similar to the correlation among those previously reported studies. Application of this analysis protocol to FKBP12 demonstrated an efficient quantitation of both weak exchange linebroadening contributions and differential residue-specific ¹⁵N CSA values. Experimental access to such differential residue-specific ¹⁵N CSA values should significantly facilitate more accurate comparisons with molecular dynamics simulations of protein motion that occurs within the timeframe of global molecular tumbling.

ARTSY-J: Convenient and precise measurement of (3)JHNHα couplings in medium-size proteins from TROSY-HSQC spectra

A new and convenient method, named ARTSY-J, is introduced that permits extraction of the (3)JHNHα couplings in proteins from the relative intensities in a pair of (15)N-(1)H TROSY-HSQC spectra. The pulse scheme includes (3)JHNHα dephasing of the narrower TROSY (1)H(N)-{(15)N} doublet component during a delay, integrated into the regular two-dimensional TROSY-HSQC pulse scheme, and compares the obtained intensity with a reference spectrum where (3)JHNHα dephasing is suppressed. The effect of passive (1)H(α) spin flips downscales the apparent (3)JHNHα coupling by a uniform factor that depends approximately linearly on both the duration of the (3)JHNHα dephasing delay and the (1)H-(1)H cross relaxation rate. Using such a correction factor, which accounts for the effects of both inhomogeneity of the radiofrequency field and (1)H(α) spin flips, agreement between prior and newly measured values for the small model protein GB3 is better than 0.3Hz. Measurement for the HIV-1 protease homodimer (22kDa) yields (3)JHNHα values that agree to better than 0.7Hz with predictions made on the basis of a previously parameterized Karplus equation. Although for Gly residues the two individual (3)JHNHα couplings cannot be extracted from a single set of ARTSY-J spectra, the measurement provides valuable ϕ angle information.