Dynamics and binding affinity of spin-labeled stearic acids in β-lactoglobulin: evidences from EPR spectroscopy and molecular dynamics simulation

Librational Dynamics of Spin-Labeled Membranes at Cryogenic Temperatures From Echo-Detected ED-EPR Spectra

Methods of electron spin echo of pulse electron paramagnetic resonance (EPR) spectroscopy are increasingly employed to investigate biophysical properties of nitroxide-labeled biosystems at cryogenic temperatures. Two-pulse echo-detected ED-spectra have proven to be valuable tools to describe the librational dynamics in the low-temperature phases of both lipids and proteins in membranes. The motional parameter, α2τC, given by the product of the mean-square angular amplitude, α2, and the rotational correlation time, τC, of the motion, is readily determined from the nitroxide ED-spectra as well as from the W-relaxation rate curves. An independent evaluation of α2 is obtained from the motionally averaged ¹⁴N-hyperfine splitting separation in the continuous wave cw-EPR spectra. Finally, the rotational correlation time τC can be estimated by combining ED- and cw-EPR data. In this mini-review, results on the librational dynamics in model and natural membranes are illustrated.

Low-temperature librations and dynamical transition in proteins at differing hydration levels

Hydration of water affects the dynamics and in turn the activity of biomacromolecules. We investigated the dependence of the librational oscillations and the dynamical transition on the hydrating conditions of two globular proteins with different structure and size, namely β-lactoglobulin (βLG) and human serum albumin (HSA), by spin-label electron paramagnetic resonance (EPR) in the temperature range of 120–270 K. The proteins were spin-labeled with 5-maleimide spin-label on free cysteins and prepared in the lyophilized state, at low (h = 0.12) and full (h = 2) hydration levels in buffer. The angular amplitudes of librations are small and almost temperature independent for both lyophilized proteins. Therefore, in these samples, the librational dynamics is restricted and the dynamical transition is absent. In the small and compact beta-structured βLG, the angular librational amplitudes increase with temperature and hydrating conditions, whereas hydration-independent librational oscillations whose amplitudes rise with temperature are recorded in the large and flexible alpha-structured HSA. Both βLG and HSA at low and fully hydration levels undergo the dynamical transition at about 230 K. The overall results indicate that protein librational dynamics is activated at the low hydration level h = 0.12 and highlight biophysical properties that are common to other biosamples at cryogenic temperatures.

The armadillo-repeat domain of plakophilin 1 binds the C-terminal sterile alpha motif (SAM) of p73

BACKGROUND

Plakophilin 1 (PKP1) is a component of desmosomes, which are key structural components for cell-cell adhesion, and can also be found in other cell locations. The p53, p63 and p73 proteins belong to the p53 family of transcription factors, playing crucial roles in tumour suppression. The α-splice variant of p73 (p73α) has at its C terminus a sterile alpha motif (SAM); such domain, SAMp73, is involved in the interaction with other macromolecules.

METHODS

We studied the binding of SAMp73 with the armadillo domain of PKP1 (ARM-PKP1) in the absence and the presence of 100 mM NaCl, by using several biophysical techniques, namely fluorescence, far-ultraviolet circular dichroism (CD), nuclear magnetic resonance (NMR), isothermal titration calorimetry (ITC), and molecular docking and simulations.

RESULTS

Association was observed between the two proteins, with a dissociation constant of ~5 μM measured by ITC and fluorescence in the absence of NaCl. The binding region of SAMp73 involved residues of the so-called "middle-loop-end-helix" binding region (i.e., comprising the third helix, together with the C terminus of the second one, and the N-cap of the fourth), as shown by ¹⁵N, ¹H- HSQC-NMR spectra. Molecular modelling provided additional information on the possible structure of the binding complex.

CONCLUSIONS

This newly-observed interaction could have potential therapeutic relevance in the tumour pathways where PKP1 is involved, and under conditions when there is a possible inactivation of p53.

GENERAL SIGNIFICANCE

The discovery of the binding between SAMp73 and ARM-PKP1 suggests a functional role for their interaction, including the possibility that SAMp73 could assist PKP1 in signalling pathways.

Rutin Is a Low Micromolar Inhibitor of SARS-CoV-2 Main Protease 3CLpro: Implications for Drug Design of Quercetin Analogs

The pandemic, due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has stimulated the search for antivirals to tackle COVID-19 infection. Molecules with known pharmacokinetics and already approved for human use have been demonstrated or predicted to be suitable to be used either directly or as a base for a scaffold-based drug design. Among these substances, quercetin is known to be a potent in vitro inhibitor of 3CLpro, the SARS-CoV-2 main protease. However, its low in vivo bioavailability calls for modifications to its molecular structure. In this work, this issue is addressed by using rutin, a natural flavonoid that is the most common glycosylated conjugate of quercetin, as a model. Combining experimental (spectroscopy and calorimetry) and simulation techniques (docking and molecular dynamics simulations), we demonstrate that the sugar adduct does not hamper rutin binding to 3CLpro, and the conjugated compound preserves a high potency (inhibition constant in the low micromolar range, K_i = 11 μM). Although showing a disruption of the pseudo-symmetry in the chemical structure, a larger steric volume and molecular weight, and a higher solubility compared to quercetin, rutin is able to associate in the active site of 3CLpro, interacting with the catalytic dyad (His41/Cys145). The overall results have implications in the drug-design of quercetin analogs, and possibly other antivirals, to target the catalytic site of the SARS-CoV-2 3CLpro.

Effects of Polar Head Nature and Tail Length of Single-Chain Lipids on the Conformational Stability of β-Lactoglobulin

Interaction between β-lactoglobulin and single-chain lipids, differing for either the length of the aliphatic chain or the molecular properties of the headgroup, was investigated at neutral and acidic pH to determine the impact on the thermal stability of the protein. Differential scanning calorimetry results with different fatty acids (from C10:0 to C18:0) show a correlation of both melting temperature and unfolding enthalpy of the protein with the ligand binding affinity, and the maximum effect was found for palmitic acid (PLM). The influence of the lipid polar head was investigated by comparing PLM with lyso-palmitoylphosphatidylcholine (LPC), which possesses the same aliphatic chain. At neutral pH, the stabilizing effect of LPC is less favorable compared to PLM. However, fluorescence results revealed that LPC can bind into the protein calyx even at acidic pH, at variance with fatty acids. Molecular dynamics simulations indicated that this difference is due to the ability of the polar head of LPC to interact with the protein loop that regulates the shift (Tanford transition) between open and closed state of the binding site of β-lactoglobulin. The results provide a rationale for how a ligand has the ability to access the protein active site at acidic conditions by overcoming the Tanford transition, and they demonstrate that β-lactoglobulin can deliver ligands with tailored properties of the polar head in a wide pH range.

Low-Temperature Dynamics of Chain-Labeled Lipids in Ester- and Ether-Linked Phosphatidylcholine Membranes

Continuous wave electron paramagnetic resonance spectroscopy and two-pulse echo detected spectra of chain-labeled lipids are used to study the dynamics of frozen lipid membranes over the temperature range 77-260 K. Bilayers of ester-linked dihexadecanoylphosphatidylcholine (DPPC) with noninterdigitated chains and ether-linked dihexadecyl phosphatidylcholine (DHPC) with interdigitated chains are considered. Rapid stochastic librations of small angular amplitude are found in both lipid matrices. In noninterdigitated DPPC bilayers, the mean-square angular amplitude, [Formula: see text], of the motion increases with temperature and it is larger close to the chain termini than close to the polar/apolar interface. In contrast, in interdigitated DHPC lamellae, [Formula: see text] is small and temperature and label-position independent at low temperature and increases steeply at high temperature. The rotational correlation time, τ_c, of librations lies in the subnanosecond range for DPPC and in the nanosecond range for DHPC. In all membrane samples, the temperature dependence of [Formula: see text] resembles that of the mean-square atomic displacement revealed by neutron scattering and a dynamical transition is detected in the range 210-240 K. The results highlight the librational oscillations and the glass-like behavior in bilayer and interdigitated lipid membranes.

Intrinsically disordered chromatin protein NUPR1 binds to the C-terminal region of Polycomb RING1B

Intrinsically disordered proteins (IDPs) are ubiquitous in eukaryotes, and they are often associated with diseases in humans. The protein NUPR1 is a multifunctional IDP involved in chromatin remodeling and in the development and progression of pancreatic cancer; however, the details of such functions are unknown. Polycomb proteins are involved in specific transcriptional cascades and gene silencing. One of the proteins of the Polycomb complex is the Ring finger protein 1 (RING1). RING1 is related to aggressive tumor features in multiple cancer types. In this work we characterized the interaction between NUPR1 and the paralogue RING1B in vitro, in silico, and in cellulo. The interaction occurred through the C-terminal region of RING1B (C-RING1B), with an affinity in the low micromolar range (∼10 μM). The binding region of NUPR1, mapped by NMR, was a hydrophobic polypeptide patch at the 30s region of its sequence, as pinpointed by computational results and site-directed mutagenesis at Ala33. The association between C-RING1B and wild-type NUPR1 also occurred in cellulo as tested by protein ligation assays; this interaction is inhibited by trifluoperazine, a drug known to hamper binding of wild-type NUPR1 with other proteins. Furthermore, the Thr68Gln and Ala33Gln/Thr68Gln mutants had a reduction in the binding toward C-RING1B as shown by in vitro, in silico, and in cellulo studies. This is an example of a well-folded partner of NUPR1, because its other interacting proteins are also unfolded. We hypothesize that NUPR1 plays an active role in chromatin remodeling and carcinogenesis, together with Polycomb proteins.

Electron spin resonance of spin-labeled lipid assemblies and proteins

Spin-label electron spin resonance (ESR) spectroscopy is a valuable means to study molecular mobility and interactions in biological systems. This paper deals with conventional, continuous wave ESR of nitroxide spin-labels at 9-GHz providing an introduction to the basic principles of the technique and applications to self-assembled lipid aggregates and proteins. Emphasis is given to segmental lipid chain order and rotational dynamics of lipid structures, environmental polarity of membranes and proteins, structure and conformational dynamics of proteins.

Fatty acid binding into the highest affinity site of human serum albumin observed in molecular dynamics simulation

Multiple molecular dynamics simulations were performed to investigate the association of stearic acid into the highest affinity binding site of human serum albumin. All binding events ended with a rapid (<10 ps) lock-in of the fatty acid due to formation of a hydrogen bond with Tyr401. The kinetics and energetics of the penetration process both depended linearly on the positional shift of the fatty acid, with an average insertion time and free energy reduction of, respectively, 32 ± 20 ps and 0.70 ± 0.15 kcal/mol per methylene group absorbed. Binding events of longer duration (tbind>1 ns) were characterized by a slow exploration of the pocket entry and, frequently, of a nearby protein crevice corresponding to a metastable state along the route to the binding site. Taken all together, these findings reconstruct the following pathway for the binding process of stearic acid: (i) contact with the protein surface, possibly facilitated by the presence of an intermediate location, (ii) probing of the site entry, (iii) insertion into the protein, and (iv) lock-in at the final position. This general description may also apply to other long-chain fatty acids binding into any of the high-affinity sites of albumin, or to specific sites of other lipid-binding proteins.

Water penetration profile at the protein-lipid interface in Na,K-ATPase membranes

The affinity of ionized fatty acids for the Na,K-ATPase is used to determine the transmembrane profile of water penetration at the protein-lipid interface. The standardized intensity of the electron spin echo envelope modulation (ESEEM) from (2)H-hyperfine interaction with D2O is determined for stearic acid, n-SASL, spin-labeled systematically at the C-n atoms throughout the chain. In both native Na,K-ATPase membranes from shark salt gland and bilayers of the extracted membrane lipids, the D2O-ESEEM intensities of fully charged n-SASL decrease progressively with position down the fatty acid chain toward the terminal methyl group. Whereas the D2O intensities decrease sharply at the n = 9 position in the lipid bilayers, a much broader transition region in the range n = 6 to 10 is found with Na,K-ATPase membranes. Correction for the bilayer population in the membranes yields the intrinsic D2O-intensity profile at the protein-lipid interface. For positions at either end of the chains, the D2O concentrations at the protein interface are greater than in the lipid bilayer, and the positional profile is much broader. This reveals the higher polarity, and consequently higher intramembrane water concentration, at the protein-lipid interface. In particular, there is a significant water concentration adjacent to the protein at the membrane midplane, unlike the situation in the bilayer regions of this cholesterol-rich membrane. Experiments with protonated fatty acid and phosphatidylcholine spin labels, both of which have a considerably lower affinity for the Na,K-ATPase, confirm these results.

Stability of trans-resveratrol associated with transport proteins

Spectrophotometry and fluorescence combined with docking and molecular dynamics simulations are used to study the effect of the carrier proteins β-lactoglobulin and human serum albumin on the degradative trans-to-cis conversion of resveratrol. The spectroscopic measurements quantify the concentration of resveratrol isoforms after 2 h of irradiation with light at 340 nm, showing that their ratio depends linearly on temperature between 20 and 50 °C and obeys an Arrhenius law with activation energies of photoisomerization of 7.8 and 11.2 kcal/mol for β-lactoglobulin and albumin, respectively, compared to 5.1 kcal/mol in solution. Thus, both proteins protect trans-resveratrol from degradation, with albumin being more effective than β-lactoglobulin. The computational techniques clarify details of the binding of trans-resveratrol to the proteins and show that the stabilizing effect correlates with an increase of the dihedral order parameter of the ligand. These findings suggest that transport proteins are viable carriers to stabilize and deliver resveratrol in vivo in the biologically effective trans form.