Search

Human serum albumin binds a wide variety of drugs with different structure and affinity to two main binding sites, drug site 1 (DS1) and drug site 2 (DS2), which partially or totally overlap with fatty acid (FA) sites. Although multiple binding sites are available for endogenous compounds, FAs are the primary physiological ligands of albumin and their competition in the occupancy of DS1 and DS2 affects the binding of exogenous molecules, with a possible impact on drug delivery. In this work, we have investigated the simultaneous binding of oleic acid, warfarin and ibuprofen to albumin using differential scanning calorimetry and fluorescence to evaluate the impact on the conformational stability of the protein. The two drugs are widely used for their anticoagulant (warfarin) and anti-inflammatory (ibuprofen) properties, and can be also considered as site markers to probe DS1 and DS2, respectively. Oleic acid is one of the most important fatty acids from a physiological point of view for its role as a source of energy for cells, and also it binds albumin with the highest association constant. When complexed with oleic acid the calorimetric profile of albumin shows a biphasic trend whose line shape depends on the ligand concentration. The binding capacity of either warfarin or ibuprofen to albumin is modulated by oleate molecules in a concentration-dependent mode being synergic cooperative (warfarin) or competitive-like (ibuprofen). The overall results provide insights on the dynamics of albumin/ligands complex, which in turn may have important pharmacokinetic and pharmacodynamic implications.

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.

Among the prognostic and predictive biomarkers of breast cancer (BC), the role of estrogen receptor (ER)α wild‐type has been acknowledged, although the action of certain ERα splice variants has not been elucidated. Insulin/insulin receptor (IR) axis has also been involved in the progression and metastasis of BC. For instance, hyperinsulinemia, which is often associated with obesity and type 2 diabetes, may be a risk factor for BC. Similarly, an aberrant expression of IR or its hyperactivation may correlate with aggressive BC phenotypes. In the present study, we have shown that a novel naturally immortalized BC cell line (named BCAHC‐1) is characterized by a unique expression of 46 kDa ERα splice variant (ERα46) along with IR. Moreover, we have shown that a multifaceted crosstalk between ERα46 and IR occurs in BCAHC‐1 cells upon estrogen and insulin exposure for growth and pulmonary metastasis. Through high‐throughput RNA sequencing analysis, we have also found that the cytokine interleukin‐11 (IL11) is the main factor linking BCAHC‐1 cells to breast cancer‐associated fibroblasts (CAFs). In particular, we have found that IL11 induced by estrogens and insulin in BCAHC‐1 cells regulates pro‐tumorigenic genes of the “extracellular matrix organization” signaling pathway, such as ICAM‐1 and ITGA5, and promotes both migratory and invasive features in breast CAFs. Overall, our results may open a new scientific avenue to identify additional prognostic and therapeutic targets in BC., –The crosstalk between ERα46 and IR signaling is involved in the growth and pulmonary metastasis of BCAHC‐1 cells.–IL11 is the most induced gene by E2 and insulin stimulation in BCAHC‐1 cells, as assessed by high‐throughput RNA sequencing analysis.–IL11 produced by BCAHC‐1 cells promotes gene expression changes and invasive features in CAFs

Spin labels based on cinobufagin, a specific inhibitor of the Na,K-ATPase, have proved valuable tools to characterize the binding site of cardiotonic steroids (CTSs), which also constitutes the extracellular cation pathway. Because existing literature suggests variations in the physiological responses caused by binding of different CTSs, we extended the original set of spin-labeled inhibitors to the more potent bufalin derivatives. Positioning of the spin labels within the Na,K-ATPase site was defined and visualized by molecular docking. Although the original cinobufagin labels exhibited lower affinity, continuous-wave electron paramagnetic resonance spectra of spin-labeled bufalins and cinobufagins revealed a high degree of pairwise similarity, implying that these two types of CTS bind in the same way. Further analysis of the spectral lineshapes of bound spin labels was performed with emphasis on their structure (PROXYL vs. TEMPO), as well as length and rigidity of the linkers. For comparable structures, the dynamic flexibility increased in parallel with linker length, with the longest linker placing the spin label at the entrance to the binding site. Temperature-related changes in spectral lineshapes indicate that six-membered nitroxide rings undergo boat-chair transitions, showing that the binding-site cross section can accommodate the accompanying changes in methyl-group orientation. D2O-electron spin echo envelope modulation in pulse-electron paramagnetic resonance measurements revealed high water accessibilities and similar polarity profiles for all bound spin labels, implying that the vestibule leading to steroid-binding site and cation-binding sites is relatively wide and water-filled.

16 pags., 8 figs., 4 tabs. -- This article belongs to the Special Issue Personalized Medicine for Multiple Sclerosis, Background: Multiple sclerosis (MS) is frequently characterized by a variety of clinical signs, often exhibiting little specificity. The diagnosis requires a combination of medical observations and instrumental tests, and any support for its objective assessment is helpful. Objective: Herein, we describe the application of thermal liquid biopsy (TLB) of blood plasma samples, a methodology for predicting the occurrence of MS with a noninvasive, quick blood test. Methods: TLB allows one to define an index (TLB score), which provides information about overall real-time alterations in plasma proteome that may be indicative of MS. Results: This pilot study, based on 85 subjects (45 MS patients and 40 controls), showed good performance indexes (sensitivity and specificity both around 70%). The diagnostic methods better discriminate between early stage and low-burden MS patients, and it is not influenced by gender, age, or assumption of therapeutic drugs. TLB is more accurate for patients having low disability level (≤ 3.0, measured by the expanded disability status scale, EDSS) and a relapsing–remitting diagnosis. Conclusion: Our results suggest that TLB can be applied to MS, especially in an initial phase of the disease when diagnosis is difficult and yet more important (in such cases, accuracy of prediction is close to 80%), as well as in personalized patient periodic monitoring. The next step will be determining its utility in differentiating between MS and other disorders, in particular in inflammatory diseases., This research was funded by the Spanish Ministry of Economy and Competitiveness and European ERDF Funds (MCIU/AEI/FEDER, EU) (BFU2016-78232-P to A.V.C.); Projects funded by Instituto de Salud Carlos III and co-funded by European Union (ESF, “Investing in your future”): “PI15/00663 (FIS project to O.A)”, “PI18/00349 (FIS project to O.A.)”, “FI19/00146 (PFIS contract for SHD)”; Diputación General de Aragón (Protein Targets and Bioactive Compounds Group E45_20R to A.V.C. and Digestive Pathology Group B25_20R to O.A.); and the Centro de Investigación Biomédica en Red en Enfermedades Hepáticas y Digestivas (CIBERehd).

ALP thanks Professors Jose Manuel Sanchez-Ruiz and Beatriz Ibarra-Molero (both from the University of Granada) for providing access and advice on their home-built software for electrostatic calculations. BR acknowledges kind hospitality and use of computational resources in the European Magnetic Resonance Center (CERM), Sesto Fiorentino (Florence), Italy. This work was supported by Spanish Ministry of Economy and Competitiveness and European ERDF Funds (MCIU/AEI/FEDER, EU) [RTI2018-097991-BI00 to JLN and RTI2018-096246-B-I00 to ALP]; FEDER/Junta de Andalucia-Consejeria de Transformacion Economica, Industria, Conocimiento y Universidades [Grant P18-RT-2413 to ALP]. Financial support from EU Horizon 2020 project EU FT-ICR MS (731077) as well as institutional (CZ.1.05/1.1.00/02.0109) and MS facility support (LM2018127 CIISB) are gratefully acknowledged. Funding for open charge: Universidad de Granada/CBUA., The phosphomimetic mutation S82D in the cancer-associated, FADdependent human NADP(H):quinone oxidoreductase 1 (hNQO1) causes a decrease in flavin-adenine dinucleotide-binding affinity and intracellular stability. We test in this work whether the evolutionarily recent neutral mutation R80H in the vicinity of S82 may alter the strong functional effects of S82 phosphorylation through electrostatic interactions. We show using biophysical and bioinformatic analyses that the reverse mutation H80R prevents the effects of S82D phosphorylation on hNQO1 by modulating the local stability. Consistently, in rat NQO1 (rNQO1) which contains R80, the effects of phosphorylation were milder, resembling the behaviour found in hNQO1 when this residue was humanized in rNQO1 (by the R80H mutation). Thus, apparently neutral and evolutionarily divergent mutations may determine the functional response of mammalian orthologues towards phosphorylation., Spanish Government, European ERDF Funds (MCIU/AEI/FEDER, EU) RTI2018-097991-BI00 RTI2018-096246-B-I00, FEDER/Junta de Andalucia-Consejeria de Transformacion Economica, Industria, Conocimiento y Universidades P18-RT-2413, EU Horizon 2020 project EU FT-ICR MS 731077, Universidad de Granada/CBUA CZ.1.05/1.1.00/02.0109 LM2018127 CIISB

Spin labels based on cinobufagin, a specific inhibitor of the Na,K-ATPase, have proved valuable tools to characterize the binding site of cardiotonic steroids (CTSs), which also constitutes the extracellular cation pathway. Because existing literature suggests variations in the physiological responses caused by binding of different CTSs, we extended the original set of spin-labeled inhibitors to the more potent bufalin derivatives. Positioning of the spin labels within the Na,K-ATPase site was defined and visualized by molecular docking. Although the original cinobufagin labels exhibited lower affinity, continuous-wave electron paramagnetic resonance spectra of spin-labeled bufalins and cinobufagins revealed a high degree of pairwise similarity, implying that these two types of CTS bind in the same way. Further analysis of the spectral lineshapes of bound spin labels was performed with emphasis on their structure (PROXYL vs. TEMPO), as well as length and rigidity of the linkers. For comparable structures, the dynamic flexibility increased in parallel with linker length, with the longest linker placing the spin label at the entrance to the binding site. Temperature-related changes in spectral lineshapes indicate that six-membered nitroxide rings undergo boat-chair transitions, showing that the binding-site cross section can accommodate the accompanying changes in methyl-group orientation. D

Over a quarter million of protein phosphorylation sites have been identified so far, although the effects of site-specific phosphorylation on protein function and stability, as well as their possible impact in the phenotypic manifestation in genetic diseases are vastly unknown. We investigated here the effects of phosphorylating S82 in human NADP(H):quinone oxidoreductase 1, a representative example of disease-associated flavoprotein in which protein stability is coupled to the intracellular flavin levels. Additionally, the cancer-associated P187S polymorphism causes inactivation and destabilization of the enzyme. By using extensive in vitro and in silico characterization of phosphomimetic S82D mutations, we showed that S82D locally affected the flavin binding site of the wild-type (WT) and P187S proteins thus altering flavin binding affinity, conformational stability and aggregation propensity. Consequently, the phosphomimetic S82D may destabilize the WT protein intracellularly by promoting the formation of the degradation-prone apo-protein. Noteworthy, WT and P187S proteins respond differently to the phosphomimetic mutation in terms of intracellular stability, further supporting differences in molecular recognition of these two variants by the proteasomal degradation pathway. We propose that phosphorylation could have critical consequences on stability and function of human flavoproteins, important for our understanding of genotype-phenotype relationships in their related genetic diseases.

Ibuprofen is a non-steroidal anti-inflammatory drug widely used to treat inflammatory diseases, and for its analgesic and antipyretic activity. Although operating as a protein inhibitor, it is also known to interact with lipid membranes. We combined calorimetry, electron spin resonance, attenuated total reflectance-Fourier transform infrared and molecular docking to characterize the interaction of ibuprofen with dimyristyolphosphatidylcholine (DMPC) bilayers, as a function of temperature and drug concentration. At increasing concentration, ibuprofen first perturbs and then suppresses the DMPC pre-transition, stabilizes the fluid state, and favours gel-fluid phase coexistence. The drug decreases the molecular packing of the polar heads and of the first methylene segments of lipid membranes in the gel phase, whereas it leaves unperturbed the chain flexibility in the liquid-crystalline phase. The action of ibuprofen also leads to a higher degree of hydration of the bilayer polar heads and favours hydrogen bond formation with solvent molecules. The overall results reveal that ibuprofen affects a number of key molecular properties of DMPC bilayers by binding through non-specific interactions at the polar/apolar interface.

Neutral cyclometalated Ir(III) complexes of general formula [Ir(ppy)2(O^O)squi], where ppy = 2-phenylpyridine and (O^O)squi = TBC (tetrabromocatechol) or TCC (tetrachlorocatechol) in their semiquinone (squi) monoanionic redox state, were synthesized by chemically oxidizing the anionic parent complexes NBu4[Ir(ppy)2(O^O)cat], in which (O^O)cat represents the corresponding ancillary dioxolene ligand in its dianionic catecholate (cat) redox state. This chemical oxidation leads to the modification of both the photophysical and the magnetic properties of the complexes. While the NBu4[Ir(ppy)2(O^O)cat] complexes are diamagnetic (D) and yellow-orange solids, the corresponding oxidized complexes [Ir(ppy)2(O^O)squi] display paramagnetic (P) properties and are characterized by a dark-green color. The conversion between the two forms (squi vs. cat) is electrochemically and chemically fully reversible. Indeed, the anionic NBu4[Ir(ppy)2(O^O)cat] complexes are quantitatively restored by chemical reduction of the neutral [Ir(ppy)2(O^O)squi] parents. These complexes therefore represent interesting redox based switches between multi-parameter states since they allow switching from a neutral paramagnetic to an anionic diamagnetic form together with a significant change in chromicity. Taking advantage of the significant color difference between the oxidized and the reduced form, an electrochromic cell was prepared with [Ir(ppy)2(TBC)squi] and its spectroelectrochemical properties are reported.

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.

The focus of this work is pointing out the different behavior of two structurally related Pt

Fully hydrated bilayers of monounsaturated palmitoyloleoylphosphatidylcholine (POPC) and diunsaturated dioleoylphosphatidylcholine (DOPC) lipids have low main phase transition temperatures (271 K for POPC and 253 K for DOPC). Two-pulse echo detected spectra, combined with continuous wave electron paramagnetic resonance spectroscopy, are employed to study the low-temperature lamellar phases of the POPC and DOPC unsaturated bilayers that are usually studied in the fluid state. Phosphatidylcholine spin-labeled at C-5 and C-16 carbon atom positions along the acyl chain were used and the temperature varied over the range 77–270 K. Segmental chain librational oscillations of small amplitude and with correlation time in the subnanosecond to nanosecond range are found in both membranes. The mean-square angular amplitude, 〈α2〉, of librations increases with temperature, is larger close to the bilayer midplane than close to the first acyl chain segments, and is larger in diunsaturated than in monounsaturated bilayers. In the inner hydrocarbon region of both lipid matrices, 〈α2〉 increases first slowly and linearly with temperature and then more rapidly, and a dynamical transition is detected in the range 190–210 K. Compared to dipalmitoylphosphatidylcholine bilayers of fully saturated symmetric chain lipids, the presence of double bonds in the acyl chain enhances the intensity of librational motion which is characterized by larger angular variations at the terminal methyl ends. These findings highlight biophysical properties of unsaturated bilayers in the frozen state, including a detailed characterization of segmental chain dynamics and the evidence of a dynamical transition that appears to be a generic feature in hydrated macromolecular systems. These results can also be relevant in regulating membrane physical properties and function at higher physiological temperatures.

Background The chemical carcinogen 3-methylcholanthrene (3MC) binds to the aryl hydrocarbon receptor (AHR) that regulates the expression of cytochrome P450 (CYP) enzymes as CYP1B1, which is involved in the oncogenic activation of environmental pollutants as well as in the estrogen biosynthesis and metabolism. 3MC was shown to induce estrogenic responses binding to the estrogen receptor (ER) α and stimulating a functional interaction between AHR and ERα. Recently, the G protein estrogen receptor (GPER) has been reported to mediate certain biological responses induced by endogenous estrogens and environmental compounds eliciting an estrogen-like activity. Methods Molecular dynamics and docking simulations were performed to evaluate the potential of 3MC to interact with GPER. SkBr3 breast cancer cells and cancer-associated fibroblasts (CAFs) derived from breast tumor patients were used as model system. Real-time PCR and western blotting analysis were performed in order to evaluate the activation of transduction mediators as well as the mRNA and protein levels of CYP1B1 and cyclin D1. Co-immunoprecipitation studies were performed in order to explore the potential of 3MC to trigger the association of GPER with AHR and EGFR. Luciferase assays were carried out to determine the activity of CYP1B1 promoter deletion constructs upon 3MC exposure, while the nuclear shuttle of AHR induced by 3MC was assessed through confocal microscopy. Cell proliferation stimulated by 3MC was determined as biological counterpart of the aforementioned experimental assays. The statistical analysis was performed by ANOVA. Results We first ascertained by docking simulations the ability of 3MC to interact with GPER. Thereafter, we established that 3MC activates the EGFR/ERK/c-Fos transduction signaling through both AHR and GPER in SkBr3 cells and CAFs. Then, we found that these receptors are involved in the up-regulation of CYP1B1 and cyclin D1 as well as in the stimulation of growth responses induced by 3MC. Conclusions In the present study we have provided novel insights regarding the molecular mechanisms by which 3MC may trigger a physical and functional interaction between AHR and GPER, leading to the stimulation of both SkBr3 breast cancer cells and CAFs. Altogether, our results indicate that 3MC may engage both GPER and AHR transduction pathways toward breast cancer progression. Electronic supplementary material The online version of this article (10.1186/s13046-019-1337-2) contains supplementary material, which is available to authorized users.

Lubrication of soft polydimethylsiloxane (PDMS) elastomer interfaces was studied in aqueous mixtures of porcine gastric mucin (PGM) and branched polyethyleneimine (b-PEI) at neutral pH and various ionic strengths (0.1-1.0 M). While neither PGM nor b-PEI improved lubrication compared to polymer-free buffer solution, their mixtures produced a synergistic lubricating effect by reducing friction coefficients by nearly two orders of magnitude, especially at slow sliding speed in the boundary lubrication regime. An array of spectroscopic studies revealed that small cationic b-PEI molecules were able to strongly bind and penetrate into large anionic PGM molecules, producing an overall contraction of the randomly coiled PGM conformation. The interaction also affected the structure of the folded PGM protein terminals, decreased the surface potential and increased light absorbance in PGM:b-PEI mixtures. Adding an electrolyte (NaCl) weakened the aggregation between PGM and b-PEI, and degraded the lubrication synergy, indicating that electrostatic interactions drive PGM:b-PEI complexation.

Several phenolic compounds bind to proteins and show the ability to interfere with their aggregation process. The impact of the natural polyphenol resveratrol on the stability and heat induced aggregation of human serum albumin (HSA) was investigated by differential scanning calorimetry (DSC), attenuated total reflectance Fourier transform infrared (ATR-FTIR), UV-vis absorbance, ThT fluorescence, atomic force microscopy (AFM) and molecular modeling. The binding of resveratrol to HSA improves the stability of the protein to thermal unfolding, particularly for the energetic domain containing the ligand binding site, as modeled by computational techniques. The thermal unfolding is irreversible and after the melting the protein aggregates, either with or without the ligand. The kinetics of HSA aggregation between 70 and 80°C shows an exponential growth of the absorbance change and it slows down when resveratrol is added. The aggregates have fibril-like morphology and resveratrol attenuates the formation of β-structured species. The overall results suggest that resveratrol stabilizes the protein structure and modulates the formation of fibrils along the initial stage of the HSA aggregation pathway.

Warfarin is a coumarin derivative drug widely used for its anticoagulant properties. The interaction of warfarin with fully hydrated lipid bilayers has been studied by combining differential scanning calorimetry, spectrophotometry, electron spin resonance of chain-labelled lipids and molecular docking. Bilayers formed by lipids with different chemico-physical properties were considered, namely dimyristoyl-phosphatidylcholine (DMPC), dimyristoyl-phosphatidylglycerol (DMPG), and dioleoyltrimethyl-ammoniumpropane (DOTAP). We observed in all cases the binding of warfarin in proximity of the surface of the bilayers, leading to a variety of distinct effects on key molecular properties of the membranes. The drug associates with the lipid bilayers in the deprotonated open chain form, with an association constant similar for DMPC and DMPG (1.27·104 and 2.82·104 M−1, respectively) and lower for DOTAP (0.46·104 M−1). In DMPC bilayers, which are zwitterionic and with saturated symmetrical chains, warfarin at 10 mol% suppresses the pre-transition, slightly stabilizes the fluid state and reduces the cooperativity of the main transition. Moreover, it alters the lateral packing density of the chain segments close to the polar/apolar interface at any temperature through the gel phase. In anionic DMPG bilayers, the drug slightly perturbs the thermotropic phase behavior, and at 10 mol% markedly loosens the compact gel phase packing of the first chain segments. In cationic DOTAP bilayers, possessing unsaturated acyl chains, the drug induces a slightly higher degree of order and motional restriction in the outer hydrocarbon region in the frozen state. In all cases, as a surface adsorbed molecule, warfarin does not affect the segmental chain order and dynamics for temperatures in the fluid phase. The overall results provide an outline of the action of warfarin on membranes formed by lipids of different types.

The focus of this work is pointing out the different behavior of two structurally related Pt(II) complexes, the anionic cyclometalated NBu4[(Bzq)Pt(Thio)], 1 and the neutral [(Phen)Pt(Thio)], 2, (Bzq = benzo[h]quinoline, Phen = 1,10-phenantroline, Thio = 1,2-benzenedithiolate), on the interaction with human serum albumin (HSA), a key drug-delivery protein in the bloodstream. Being very limited the number of anionic Pt(II) complexes reported to date, this is a pioneering example of report on a protein-ligand interaction involving a negatively charged platinum compound. The study was carried out by using fluorescence spectroscopy, differential scanning calorimetry and molecular docking simulations. The results revealed a strong binding affinity between the anionic compound and the protein, whereas a weak/moderate binding interaction was highlighted for the neutral one. Comparative studies with site specific ligands (warfarin and ibuprofen), allowed us to identify the protein binding sites of the two compounds. The work aims to shed light on the relevance of the charge in designing new drugs with a favorable binding affinity for HSA, which strongly contributes to influence their pharmacological and toxicological profile.

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.

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 (tbind1 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.

Electron spin echo envelope modulation (ESEEM) and conventional electron paramagnetic resonance (EPR) of site-specifically spin-labelled phospholipids are used to investigate the effect of ether-linked chains on the water-penetration and polarity profiles, as well as the phase behaviour and chain flexibility profiles, of phospholipid membranes. D2O-ESEEM reveals that water exposure of the terminal methyl groups in the interdigitated phase of dihexadecyl phosphatidylcholine (DHPC) is comparable to that of the methylene groups at the polar head-group end of the chains. Similarly, an uniform transmembrane polarity profile is obtained from the dependence of the outer 14N-hyperfine splitting on the spin-label position along the chain in frozen interdigitated DHPC dispersions. Two-component conventional EPR spectra of spin labels at the terminal methyl end of the chain reveal that the intermediate gel phase above the pretransition of DHPC contains components in which the lipid chains are interdigitated. The polarity and chain-flexibility profiles in the fluid Lα-phase of DHPC with ether-linked chains are shifted outwards, towards the polar-apolar interface, as compared with that of dihexadecanoyl phosphatidylcholine (DPPC) with ester-linked chains. Also, the polarity profile of DHPC is shifted upwards, to higher polarities. These differences reflect those in hydrocarbon thickness and area/lipid molecule reported by x-ray diffraction for the Lα-phases of the two lipids.

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, ⟨α2⟩, 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, ⟨α2⟩ 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 sa...

Ibuprofen is one of the most used anti-inflammatory drugs, and it is transported in the blood by human serum albumin, a major plasmatic protein with a peculiar adaptability in the binding of several different ligands. We have characterized the interaction between albumin and ibuprofen, either in racemic mixture, or in the S(+) and R(-) enantiomeric forms, by using differential scanning calorimetry, attenuated total reflectance Fourier transform infrared spectroscopy, and molecular dynamics simulation. The results show that increasing concentrations of ibuprofen (up to sixfold drug/protein molar ratio) improve the protein resistance to thermal unfolding without altering the secondary structure. Deconvolution of the calorimetric thermal profiles at different albumin/ibuprofen molar ratios demonstrates a selective stability of the protein domains where the binding sites of the drug are localized. At the highest ibuprofen concentration, the melting temperature increased by about 10°C with respect to the drug-free protein, whereas the unfolding enthalpy maintains an almost constant value. Furthermore, the degree of protein stabilization depends upon the chirality of the drug, and the R(-) enantiomer is more effective compared to the S(+) form. The stability is supported by molecular dynamics simulations, showing that ibuprofen maintains a stable coordination in the most favorable binding sites, leading to a more compact protein structure at high temperature. The overall results attest that the binding of ibuprofen determines on albumin a stereoselective and domain-specific stabilization with a predominantly entropic character, contributing to clarify significant aspects of the molecular mechanism of protein/drug interaction.

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.

The interaction of saturated fatty acids of different length (C8:0 to C18:0) with β-lactoglobulin (βLG) was investigated by molecular dynamics simulation and docking approaches. The results show that the presence of such ligands in the hydrophobic central cavity of βLG, known as the protein calyx, determines an enhancement of atomic fluctuations compared with the unliganded form, especially for loops at the entrance of the binding site. Concerted motions are evidenced for protein regions that could favor the binding of ligands. The mechanism of anchoring of fatty acids of different length is similar for the carboxylate head-group, through electrostatic interactions with the side chains of Lys60/Lys69. The key protein residues to secure the hydrocarbon chain are Phe105/Met107, which adapt their conformation upon ligand binding. In particular, Phe105 provides an additional hydrophobic clamp only for the tail of the two fatty acids with the longest chains, palmitic, and stearic acid, which are known to bind βLG with a high affinity. The search of additional external binding sites for fatty acids, distinct from the calyx, was also carried out for palmitic acid. Two external sites with a lower affinity were identified as secondary sites, one consisting in a hydrophobic cavity allowing two distinct binding modes for the fatty acid, and the other corresponding to a surface crevice close to the protein α-helix. The overall results provide a comprehensive picture of the dynamical behavior of βLG in complex with fatty acids, and elucidate the structural basis of the binding of these physiological ligands.

Protein-drug interaction is of prominent interest in determining the pharmacokinetic and pharmacodynamic consequences on drug delivery. Warfarin is a widely used anticoagulant drug in the treatment of venous thrombosis and pulmonary embolism and is carried in the blood almost exclusively by human serum albumin. The effects of the binding of warfarin to the native state of albumin were characterized by UV-vis absorption, conventional and synchronous fluorescence, isothermal titration calorimetry, differential scanning calorimetry and molecular dynamics simulation. The overall results indicate that, under physiological condition, the binding of warfarin in site DS1 of albumin promotes local stabilization with resulting effects on the global protein dynamics. The increase of the protein stability has both an enthalpic and entropic character. Under denaturing condition, the stabilizing effect of warfarin is evidenced by an increase of both the melting temperature and unfolding enthalpy of albumin with the drug/protein molar ratio. More importantly, thermal resistance is increased due to selective effect on the specific protein lobe that includes the main drug binding site. The comparison of the thermal behavior of the protein-warfarin complex with that in the presence of a typical ligand of the other main protein binding site, i.e. drug site DS2, provides key insight on domain-specific stabilization effects on albumin.

The energy landscape of proteins is characterized by a hierarchy of substates, which give rise to conformational heterogeneity at low temperatures. In multiply spin-labeled membranous Na, K-ATPase, this heterogeneous population of conformations is manifest by strong inhomogeneous broadening of the electron paramagnetic resonance (EPR) line shapes and nonexponential spin-echo decays, which undergo a transition to homogeneous broadening and exponential relaxation at higher temperatures (previous study). In this study, we apply these EPR methods to small water-soluble proteins, of the type for which the existence of conformational substates is well established. Both alpha-helical and beta-sheet aqueous proteins that are spin-labeled on a single cysteine residue display spin-echo decays with a single phase-memory time T-2M and conventional EPR line shapes with predominantly homogeneous broadening, over a broad range of temperatures from 77 K to similar to 250 K or higher. Above similar to 200 K, the residual inhomogeneous broadening is reduced almost to zero. In contrast, both the proteins and the spin label alone, when in a glycerol-water mixture below the glass transition, display heterogeneity in spin-echo phase-memory time and a stronger inhomogeneous broadening of the conventional line shapes, similar to multiply spin-labeled membranous Na, K-ATPase below 200 K. Above 200 K (or the glass transition), a single phase-memory time and predominantly homogeneous broadening are found in both spin-label systems. The results are discussed in terms of solvent-mediated protein transitions, the ability of single spin-label sites to detect conformational heterogeneity, and the desirability of exploring multiple sites for proteins with the size and complexity of the Na,K-ATPase.

The spectroscopic, thermal, and functional properties of blue copper proteins can be modulated by mutations in the metal binding loop. Molecular dynamics simulation was used to compare the conformational properties of azurin and two chimeric variants, which were obtained by inserting into the azurin scaffold the copper binding loop of amicyanin and plastocyanin, respectively. Simulations at room temperature show that the proteins retain their overall structure and exhibit concerted motions among specific inner regions, as revealed by principal component analysis. Molecular dynamics at high temperature indicates that the first events in the unfolding pathway are structurally similar in the three proteins and unfolding starts from the region of the alpha-helix that is far from the metal binding loop. The results provide details of the denaturation process that are consistent with experimental data and in close agreement with other computational approaches, suggesting a distinct mechanism of unfolding of azurin and its chimeric variants. Moreover, differences observed in the dynamics of specific regions in the three proteins correlate with their thermal behavior, contributing to the determination of the basic factors that influence the stability.

Resveratrol is a natural polyphenol found in various plants with potential therapeutic activity as anti-oxidant, anti-inflammatory, cardioprotective and anti-tumoral. Lipid membranes are among cellular components that are targets of its action. In this work ESR of chain labeled lipids, calorimetry, X-ray diffraction and molecular docking are used to study the interaction of resveratrol with membrane model systems of dipalmitoylphosphatidylcholine (DPPC) as a function of resveratrol concentration (0–30 mol% of the lipid) and temperature (10–50 °C). Resveratrol incorporated in DPPC bilayers induces considerable motional restriction at the lipid tail termini, removing the gradient of increasing mobility along the chain found in DPPC bilayers in the gel phase. In contrast, it leaves unperturbed the DPPC chain flexibility profile in the liquid-crystalline phase. At low concentration, resveratrol progressively reduces the pre-transition temperature and eliminates the pre-transition for content ≥5 mol%. A reduced cooperativity and a downshift of the main transition temperature are observed, especially at high content. The typical diffraction pattern of DPPC multibilayers in the Lβ′ phase is converted to a lamellar pattern with reduced d-spacing of untilted lipid chain in a hexagonal packing at 30 mol% of resveratrol. Molecular docking indicates that the energetically favoured anchoring site is the polar headgroup region, where resveratrol acts as a spacer. The overall results are consistent with the formation in DPPC of an interdigitated Lβi gel phase induced by 30 mol% resveratrol.

The copper site and overall structures of azurin (AZ) variants in which the amicyanin (AMI) and plastocyanin (PC) metal binding loops have been introduced, AZAMI and AZPC, respectively, are similar to that of AZ, whereas the loop conformations resemble those in the native proteins. To assess the influence of these loop mutations on stability, the thermal unfolding of AZAMI and AZPC has been investigated by differential scanning calorimetry, absorption and fluorescence spectroscopy. The calorimetric profiles of both variants exhibit a complex shape consisting of two endothermic peaks and an exothermic peak. The temperature of the maximum heat of absorption for the single endothermic peak is 82.7 °C for AZ, whereas for AZAMI and AZPC the most intense endothermic peaks are at 74.9 and 68.1 °C comparable to values for AMI and PC, respectively. Denaturation investigated using the temperature dependence of the absorbance at ∼600 nm and Trp emission, also demonstrates decreased stability for both loop mutants. The thermal transition between the native and the denaturated states is irreversible, scan rate dependent and consistent with the two-state irreversible model. The structure of the active-site loop has a dramatic effect on the kinetic stability and the unfolding pathway of cupredoxins.

Two-pulse, echo-detected electron paramagnetic resonance (ED-EPR) spectra and continuous-wave EPR (CW-EPR) spectra were used to investigate the solvent effect on the librational motion of human haemoglobin spin-labelled on cysteine β93 with the nitroxide derivative of maleimide, 6-MSL. Protein samples fully hydrated in phosphate buffer solution (PBS), in a 60% v/v glycerol/water mixture and in the lyophilized form were measured at cryogenic temperature in the frozen state. The protein librational motion was characterized by the amplitude-correlation time product, τ(c), deduced from the ED-EPR spectra. The librational amplitude, τ(c), was determined independently, from the motionally averaged hyperfine splitting in the CW-EPR spectra, and the librational correlation time, τ(c), was derived from the combination of the pulsed and conventional EPR data. Rapid librational motion of small amplitude was detected in all samples. In each case, the librational dynamics was restricted up to 180 K, beyond which it increased steeply for the hydrated protein in PBS and in the presence of glycerol. In contrast, in the dehydrated protein, the librational dynamics was hindered and less dependent on temperature up to ~240 K. In all samples, deviated from small values only for T > 200 K, where a rapid increase of was evident for the hydrated samples, whereas limited temperature variation was shown in the lyophilized samples. The librational correlation time was in the sub-nanosecond regime and weakly dependent on temperature. The results evidence that solvent favours protein dynamics.

Using electron scanning microscopy, we have studied the protein deposit left on silicon and mica substrates by dried droplets of aqueous solutions of bovine beta-lactoglobulin at low concentration and pH = 2-7. We have observed different self-assembled structures: homogeneous layers, hexagonal platelets and flower-shaped patterns laying flat on the surface, and rods formed by columns. Homogeneous layers covered the largest area of the droplet deposit. The other structures were found in small isolated regions, where the protein solution dried in the form of microdroplets. The presence of hexagonal platelets, flower-shaped patterns and columnar rods shows that beta-lactoglobulin self-assembles at the surface in a hexagonal columnar phase, which has never been observed in solution. A comparison with proteins showing similar aggregates suggests that beta-lactoglobulin structures grow from hexagonal germs composed of discotic nanometric building blocks, possibly possessing an octameric structure. We propose that discotic building blocks of beta-lactoglobulin may be produced by the anisotropic interaction with the solid surface.

Membranous Na,K-ATPase from shark salt gland and from pig kidney was spin-labeled on class I -SH groups in the presence of glycerol, or on class II -SH groups in the absence of glycerol. The class I-labeled preparations retain full enzymatic activity, whereas the class II-labeled preparations are at least partially inactivated. This provides an excellent testbed on which to demonstrate how advanced electron paramagnetic resonance (EPR) can provide novel information on specific residues in unique environments in a complex, membrane-bound transport system. The polarity of the environment, and the librational dynamics and conformational exchange, of the spin-labeled groups were studied with pulsed EPR by using electron spin echo envelope modulation (ESEEM) spectroscopy and spin-echo detected (ED) EPR spectroscopy, respectively. 2H-ESEEM spectra of membranes dispersed in D2O reveal that class I groups of the shark enzyme are more exposed to water than are those of the pig enzyme or class II groups of either species, consistent with the more superficial membrane location in the former case. Spin-echo decay curves indicate conformational heterogeneity at low temperatures (150 K), but a more homogeneous conformational state at higher temperatures that is characterized by a single phase-memory T2M relaxation time. Conventional EPR lineshapes also demonstrate conformational microheterogeneity at low temperatures: the inhomogeneously broadened lines narrow progressively with increasing temperature reaching an almost pure Lorentzian line shape at temperatures of ca. 220 K and above. The inhomogeneous broadening at low temperature is well described by a Gaussian distribution of Lorentzian lines. ED spectra as a function of echo-delay time demonstrate the onset of rapid librational motions of appreciable amplitude, and slower conformational exchange, at temperatures above 220 K. These motions could drive transitions between the different conformational substates, which are frozen in at lower temperatures but contribute to the pathways between the principal enzymatic intermediates at higher temperatures.

Molecular dynamics simulation has been carried out for the blue copper protein amicyanin from two different sources, Paracoccus denitrificans and Paraccocus versutus, to investigate the structural and dynamical properties common to the two molecules and to identify prominent features shared with proteins of the same family, the monomeric cupredoxins. The two amicyanins have almost identical secondary and tertiary structure. In the simulation, they differ for the number of hydrogen bonds in the main chain and the conformation of some beta-strands. However, they strictly maintain the arrangement of the portions of the beta-barrel that are conserved in the folding architecture of the blue copper proteins. Paracoccus versutus amicyanin equilibrates more rapidly, shows lower atomic deviation values, and is less rigid with respect to Paracoccus denitrificans amicyanin. Principal component analysis reveals that the conformational subspaces corresponding to eigenvectors with the same index for each of the two molecules are not necessarily equivalent. Nevertheless, a core scaffold with constrained dynamics exist for both amicyanins. In addition, two fairly flexible regions that are located on the opposite side with respect to the interaction sites with the partner molecules in the redox process have been evidenced in the protein structure. This description of amicyanin, with a few mobile regions remote from the active site and a rigid scaffold including most of the protein beta-barrel, has a close similarity with that of azurin and plastocyanin, two other cupredoxins previously investigated in simulation. Furthermore, similarities in the distribution of the atomic fluctuations indicate that amicyanin, azurin, and plastocyanin possess common dynamical features, in spite of differences in their structure. On the basis of these findings, we suggest that topological constraints imposed by the folding in correspondence of protein regions that are the most conserved determine the protein dynamics of the cupredoxin family. The dynamical properties of the cupredoxins might be controlled for functional advantages that include the binding mechanism with the biological partners and the collective inner motions of the protein matrix required for the electron transfer, whereas long-range conformational changes in the redox reaction should be excluded.

The thermal stability of umecyanin, a stellacyanin from horseradish roots, has been investigated by differential scanning calorimetry, optical absorption and fluorescence spectroscopy at neutral and alkaline pH. Above pH 9 the Cu(II) protein experiences a blue shift of the main visible absorption band at approximately 600 nm and changes colour from blue to violet. The thermal transition of the protein is irreversible and occurs between 61.4 and 68.8 degrees C at pH 7.5 and between 50.7 and 57.4 degrees C at pH 9.8. The calorimetric data indicates that at both pH values the thermally induced transition of the protein between the native and denaturated states can be described in terms of the classical Lumry-Eyring unfolding model Native--Unfolded--Final. The analysis of the reversible step in the unfolding pathway demonstrates a significant reduction in conformational stability (DeltaG) of the alkaline form of the protein. Such a reduction is consistent with an enhanced flexibility of UMC at high pH and has mainly entropic character.

The energetics of amyloid fibrillar aggregation of β-lactoglobulin (βLG) following incubation at high temperature and acid pH was studied by differential scanning calorimetry in the presence of Cu(2+) or Fe(3+) cations, and without any metal. Cu(2+) and metal-free protein solutions showed a distinct exothermic response that disappeared almost completely when the Fe(3+) molar concentration was ten times greater than the βLG concentration. Thioflavin T fluorescence studies in solution and atomic force microscopy analysis of the deposit left on flat mica substrates by heat-incubated βLG solutions correlated the absence of exothermic response of Fe(3+)-βLG solutions with a lack of fibril production. In contrast, abundant fibril deposits were observed for Cu(2+)-βLG solutions, with a rich polymorphism of multistrand fibrillar structures. Electron paramagnetic resonance revealed that Fe(3+) permanently binds to βLG in the aggregate state whereas Cu(2+) plays a catalytic role without binding to the protein. We propose that Fe(3+) inhibits fibril production after binding to a key region of the protein sequence, possibly interfering with the nucleation step of the fibrillation process and opening a nonfibrillar aggregation pathway. These findings suggest that transition metal ions can be utilized to effectively modulate protein self-assembly into a variety of structures with distinct morphologies at the nanoscale level.

Recently, developments in time-resolved spin-label electron spin resonance (ESR) spectroscopy have contributed considerably to the study of biomembranes. Two different applications of electron spin echo spectroscopy of spin-labelled phospholipids are reviewed here: (1) the use of partially relaxed echo-detected ESR spectra to study the librational lipid-chain motions in the low-temperature phases of phospholipid bilayers; (2) the use of electron spin echo envelope modulation spectroscopy to determine the penetration of water into phospholipid membranes. Results are described for phosphatidylcholine bilayer membranes, with and without equimolar cholesterol, that are obtained with phosphatidylcholine spin probes site-specifically labelled throughout the sn-2 chain.

The contribution of the copper ion to the stability and to the unfolding pathway of pseudoazurin was investigated by a comparative analysis of the thermal unfolding of the Cu(II)-holo and apo form of the protein. The unfolding has been followed by calorimetry, fluorescence, optical density, and electron paramagnetic resonance (EPR) spectroscopy. The thermal transition of Cu(II)-holo pseudoazurin is irreversible and occurs between 60.0 and 67.3°C, depending on the scan rate and technique used. The denaturation pathway of Cu(II)-holo pseudoazurin can be described by the Lumry–Eyring model: N U F; the protein reversibly goes from the native (N) to the unfolded (U) state, and then irreversibly to the final (F) state. The simulation of the experimental calorimetric profiles, according to this model, allowed us to determine the thermodynamic and kinetic parameters of the two steps. The ΔG value calculated for the Cu(II)-holo pseudoazurin is 39.2 kJ·mol−1 at 25°C. The sequence of events in the denaturation process of Cu(II)-holo pseudoazurin emergence starts with the disruption of the copper site and the hydrophobic core destabilization followed by the global protein unfolding. According to the EPR findings, the native type-1 copper ion shows type-2 copper features after the denaturation. The removal of the copper ion (apo form) significantly reduces the stability of the protein as evidenced by a ΔG value of 16.5 kJ·mol−1 at 25°C. Moreover, the apo Paz unfolding occurs at 41.8°C and is compatible with a two-state reversible process N U.© 2006 Wiley Periodicals, Inc. Biopolymers 83: 487–497, 2006 This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com

Two-pulse, echo-detected (ED) electron paramagnetic resonance (EPR) spectroscopy was used to study the librational motions of spin-labeled lipids in membranes of dipalmitoylphosphatidylcholine + 50 mol % cholesterol. The temperature dependence, over the range 77-240 K, and the dependence on position of spin-labeling in the sn-2 chain (n=5, 7, 10, 12, and 14) of the phospholipid, were characterized in detail. The experimental ED-spectra were corrected for instantaneous spin diffusion arising from static spin-spin interactions, by using spectra recorded at 77 K, where motional contributions are negligible. Simulations according to a model of rapid, small-amplitude librations about an axis whose direction is randomly distributed are able to describe the experimental spectra. Calibrations, in terms of the amplitude-correlation time product, alpha2tauc, were constructed for diagnostic spectral line-height ratios at different echo delay times, and for relaxation spectra obtained from the ratio of ED-spectra recorded at two different echo delays. The librational amplitude, alpha2, was determined for a spin label at the 14-C position of the lipid chain from the partially motionally averaged hyperfine splitting in the conventional EPR spectra. The librational correlation time, tauc, which is deduced from combination of the conventional and ED-EPR results, lies in the subnanosecond regime and depends only weakly on temperature. The temperature dependence of the ED-EPR spectra arises mainly from an increase in librational amplitude with increasing temperature, and position down the lipid chain. A gradual transition takes place at higher temperatures, from a situation in which segmental torsional librations are cumulative, i.e., the contributions of the individual segments add up progressively upon going down the chain, to one of concerted motion only weakly dependent on chain position. Such librational motions are important for glass-like states and are generally relevant to high lipid packing densities, e.g., in cholesterol-containing raft domains and condensed complexes.

High-temperature molecular dynamics simulations were used to gain insight into the early steps in the unfolding pathway of azurin, a blue copper protein with a beta-barrel structure formed by two sheets arranged in a Greek key folding topology. The results reveal that unfolding of the beta-barrel in azurin is associated with dislocation of its unique alpha-helix with respect to the protein scaffold. Exposure of the hydrophobic core to solvent precedes complete disruption of secondary and tertiary structure, and modifications in the region around the active site are directly connected with this event. Denaturation of the protein initiates from the sheet coordinating the copper ion, and the other sheet maintains its topology. Results derived from the simulation were compared with experimental data obtained with different techniques, showing excellent agreement and providing a framework to understand the process of disruption and formation of the beta-barrel in azurin.

Active site modeling in molecular dynamics simulations is investigated for the reduced state of copper azurin. Five simulation runs (5 ns each) were performed at room temperature to study the consequences of a mixed electrostatic/constrained modeling for the coordination between the metal and the polypeptide chain, using for the ligand residues a set of charges that is modified with respect to the apo form of the protein by the presence of the copper ion. The results show that the different charge values do not lead to relevant effects on the geometry of the active site of the protein, as long as bond distance constraints are used for all the five ligand atoms. The distance constraint on the O atom of Gly45 can be removed without altering the active site geometry. The coordination between Cu and the other axial ligand Met121 is outlined as being flexible. Differences are found between the bonds of the copper ion with the two apparently equivalent Nδ1 atoms of His46 and His117. The overall findings are discussed in connection with the issue of determining a model for the active site of azurin suitable to be used in molecular dynamics simulations under unfolding conditions. Figure Model of azurin active site. Copper ligand residues are cut off at Cα position except Gly45, for which the portion of backbone connecting it to His46 is shown. Only polar H atoms are shown. All atoms are in standard colors (Cu in violet), and the five ligands are labeled