Protein dynamics: hydration and cavities

The temperature-pressure behavior of proteins seems to be unique among the biological macromolecules. Thermodynamic as well as kinetic data show the typical elliptical stability diagram. This may be extended by assuming that the unfolded state gives rise to volume and enthalpy-driven liquid-liquid transitions. A molecular interpretation follows from the temperature and the pressure dependence of the hydration and cavities. We suggest that positron annihilation spectroscopy can provide additional quantitative evidence for the contributions of cavities to the dynamics of proteins. Only mature amyloid fibrils that form from unfolded proteins are very resistant to pressure treatment.

Temperature dependence of protein dynamics simulated with three different water models

The effect of variation of the water model on the temperature dependence of protein and hydration water dynamics is examined by performing molecular dynamics simulations of myoglobin with the TIP3P, TIP4P, and TIP5P water models and the CHARMM protein force field at temperatures between 20 and 300 K. The atomic mean-square displacements, solvent reorientational relaxation times, pair angular correlations between surface water molecules, and time-averaged structures of the protein are all found to be similar, and the protein dynamical transition is described almost indistinguishably for the three water potentials. The results provide evidence that for some purposes changing the water model in protein simulations without a loss of accuracy may be possible.

On the interaction of bovine serum albumin with ionic surfactants: Temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

The interaction of bovine serum albumin (BSA) with the ionic surfactants sodium dodecylsulfate (SDS, anionic), cetyltrimethylammonium chloride (CTAC, cationic) and N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS, zwitterionic) was studied by electron paramagnetic resonance (EPR) spectroscopy of spin label covalently bound to the single free thiol group of the protein. EPR spectra simulation allows to monitor the protein dynamics at the labeling site and to estimate the changes in standard Gibbs free energy, enthalpy and entropy for transferring the nitroxide side chain from the more motionally restricted to the less restricted component. Whereas SDS and CTAC showed similar increases in the dynamics of the protein backbone for all measured concentrations. HPS presented a smaller effect at concentrations above 1.5 mM. At 10 mM of surfactants and 0.15 mM BSA, the standard Gibbs free energy change was consistent with protein backbone conformations more expanded and exposed to the solvent as compared to the native protein, but with a less pronounced effect for HPS. In the presence of the surfactants, the enthalpy change, related to the energy required to dissociate the nitroxide side chain from the protein, was greater, suggesting a lower water activity. The nitroxide side chain also detected a higher viscosity environment in the vicinity of the paramagnetic probe induced by the addition of the surfactants. The results suggest that the surfactant-BSA interaction, at higher surfactant concentration, is affected by the affinities of the surfactant to its own micelles and micelle-like aggregates. Complementary DLS data suggests that the temperature induced changes monitored by the nitroxide probe reflects local changes in the vicinity of the single thiol group of Cys-34 BSA residue. (C) 2011 Elsevier B.V. All rights reserved.

On the interaction of bovine serum albumin with ionic surfactants: Temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Envelope protein dynamics in paramyxovirus entry.

Paramyxoviruses include major pathogens with significant global health and economic impact. This large family of enveloped RNA viruses infects cells by employing two surface glycoproteins that tightly cooperate to fuse their lipid envelopes with the target cell plasma membrane, an attachment and a fusion (F) protein. Membrane fusion is believed to depend on receptor-induced conformational changes within the attachment protein that lead to the activation and subsequent refolding of F. While structural and mechanistic studies have considerably advanced our insight into paramyxovirus cell adhesion and the structural basis of F refolding, how precisely the attachment protein links receptor engagement to F triggering remained poorly understood. Recent reports based on work with several paramyxovirus family members have transformed our understanding of the triggering mechanism of the membrane fusion machinery. Here, we review these recent findings, which (i) offer a broader mechanistic understanding of the paramyxovirus cell entry system, (ii) illuminate key similarities and differences between entry strategies of different paramyxovirus family members, and (iii) suggest new strategies for the development of novel therapeutics.

On the interaction of bovine serum albumin with ionic surfactants: Temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

The interaction of bovine serum albumin (BSA) with the ionic surfactants sodium dodecylsulfate (SDS, anionic), cetyltrimethylammonium chloride (CTAC, cationic) and N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS, zwitterionic) was studied by electron paramagnetic resonance (EPR) spectroscopy of spin label covalently bound to the single free thiol group of the protein. EPR spectra simulation allows to monitor the protein dynamics at the labeling site and to estimate the changes in standard Gibbs free energy, enthalpy and entropy for transferring the nitroxide side chain from the more motionally restricted to the less restricted component. Whereas SDS and CTAC showed similar increases in the dynamics of the protein backbone for all measured concentrations. HPS presented a smaller effect at concentrations above 1.5 mM. At 10 mM of surfactants and 0.15 mM BSA, the standard Gibbs free energy change was consistent with protein backbone conformations more expanded and exposed to the solvent as compared to the native protein, but with a less pronounced effect for HPS. In the presence of the surfactants, the enthalpy change, related to the energy required to dissociate the nitroxide side chain from the protein, was greater, suggesting a lower water activity. The nitroxide side chain also detected a higher viscosity environment in the vicinity of the paramagnetic probe induced by the addition of the surfactants. The results suggest that the surfactant-BSA interaction, at higher surfactant concentration, is affected by the affinities of the surfactant to its own micelles and micelle-like aggregates. Complementary DLS data suggests that the temperature induced changes monitored by the nitroxide probe reflects local changes in the vicinity of the single thiol group of Cys-34 BSA residue. (C) 2011 Elsevier B.V. All rights reserved.

A multidisciplinary approach to the study of protein dynamics and signal transmission

Allostery is a phenomenon of fundamental importance in biology, allowing regulation of function and dynamic adaptability of enzymes and proteins. Despite the allosteric effect was first observed more than a century ago allostery remains a biophysical enigma, defined as the “second secret of life”. The challenge is mainly associated to the rather complex nature of the allosteric mechanisms, which manifests itself as the alteration of the biological function of a protein/enzyme (e.g. ligand/substrate binding at the active site) by binding of “other object” (“allos stereos” in Greek) at a site distant (> 1 nanometer) from the active site, namely the effector site. Thus, at the heart of allostery there is signal propagation from the effector to the active site through a dense protein matrix, with a fundamental challenge being represented by the elucidation of the physico-chemical interactions between amino acid residues allowing communicatio n between the two binding sites, i.e. the “allosteric pathways”. Here, we propose a multidisciplinary approach based on a combination of computational chemistry, involving molecular dynamics simulations of protein motions, (bio)physical analysis of allosteric systems, including multiple sequence alignments of known allosteric systems, and mathematical tools based on graph theory and machine learning that can greatly help understanding the complexity of dynamical interactions involved in the different allosteric systems. The project aims at developing robust and fast tools to identify unknown allosteric pathways. The characterization and predictions of such allosteric spots could elucidate and fully exploit the power of allosteric modulation in enzymes and DNA-protein complexes, with great potential applications in enzyme engineering and drug discovery.

Local and nonlocal integro-differential reaction-diffusion models for protein dynamics in Alzheimer's disease

In this work, integro-differential reaction-diffusion models are presented for the description of the temporal and spatial evolution of the concentrations of Abeta and tau proteins involved in Alzheimer's disease. Initially, a local model is analysed: this is obtained by coupling with an interaction term two heterodimer models, modified by adding diffusion and Holling functional terms of the second type. We then move on to the presentation of three nonlocal models, which differ according to the type of the growth (exponential, logistic or Gompertzian) considered for healthy proteins. In these models integral terms are introduced to consider the interaction between proteins that are located at different spatial points possibly far apart. For each of the models introduced, the determination of equilibrium points with their stability and a study of the clearance inequalities are carried out. In addition, since the integrals introduced imply a spatial nonlocality in the models exhibited, some general features of nonlocal models are presented. Afterwards, with the aim of developing simulations, it is decided to transfer the nonlocal models to a brain graph called connectome. Therefore, after setting out the construction of such a graph, we move on to the description of Laplacian and convolution operations on a graph. Taking advantage of all these elements, we finally move on to the translation of the continuous models described above into discrete models on the connectome. To conclude, the results of some simulations concerning the discrete models just derived are presented.

Structure, Dynamics, and RNA Interaction Analysis of the Human SBDS Protein

Shwachman-Bodian-Diamond syndrome is an autosomal recessive genetic syndrome with pleiotropic phenotypes, including pancreatic deficiencies, bone marrow dysfunctions with increased risk of myelodysplasia or leukemia, and skeletal abnormalities. This syndrome has been associated with mutations in the SBDS gene, which encodes a conserved protein showing orthologs in Archaea and eukaryotes. The Shwachman-Bodian-Diamond syndrome pleiotropic phenotypes may be an indication of different cell type requirements for a fully functional SBDS protein. RNA-binding activity has been predicted for archaeal and yeast SBDS orthologs, with the latter also being implicated in ribosome biogenesis. However, full-length SBDS orthologs function in a species-specific manner, indicating that the knowledge obtained from model systems may be of limited use in understanding major unresolved issues regarding SBDS function, namely, the effect of mutations in human SBDS on its biochemical function and the specificity of RNA interaction. We determined the solution structure and backbone dynamics of the human SBDS protein and describe its RNA binding site using NMR spectroscopy. Similarly to the crystal structures of Archaea, the overall structure of human SBDS comprises three well-folded domains. However, significant conformational exchange was observed in NMR dynamics experiments for the flexible linker between the N-terminal domain and the central domain, and these experiments also reflect the relative motions of the domains. RNA titrations monitored by heteronuclear correlation experiments and chemical shift mapping analysis identified a classic RNA binding site at the N-terminal FYSH (fungal, Yhr087wp, Shwachman) domain that concentrates most of the mutations described for the human SBDS. (C) 2010 Elsevier Ltd. All rights reserved.

The coupling between electron transfer and Protein/Solvent Dynamics in Photosynthetic Reaction Centers: Spectroscopic Studies in Amorphous Matrices

We investigated at the molecular level protein/solvent interactions and their relevance in protein function through the use of amorphous matrices at room temperature. As a model protein, we used the bacterial photosynthetic reaction center (RC) of Rhodobacter sphaeroides, a pigment protein complex which catalyzes the light-induced charge separation initiating the conversion of solar into chemical energy. The thermal fluctuations of the RC and its dielectric conformational relaxation following photoexcitation have been probed by analyzing the recombination kinetics of the primary charge-separated (P+QA-) state, using time resolved optical and EPR spectroscopies. We have shown that the RC dynamics coupled to this electron transfer process can be progressively inhibited at room temperature by decreasing the water content of RC films or of RC-trehalose glassy matrices. Extensive dehydration of the amorphous matrices inhibits RC relaxation and interconversion among conformational substates to an extent comparable to that attained at cryogenic temperatures in water-glycerol samples. An isopiestic method has been developed to finely tune the hydration level of the system. We have combined FTIR spectral analysis of the combination and association bands of residual water with differential light-minus-dark FTIR and high-field EPR spectroscopy to gain information on thermodynamics of water sorption, and on structure/dynamics of the residual water molecules, of protein residues and of RC cofactors. The following main conclusions were reached: (i) the RC dynamics is slaved to that of the hydration shell; (ii) in dehydrated trehalose glasses inhibition of protein dynamics is most likely mediated by residual water molecules simultaneously bound to protein residues and sugar molecules at the protein-matrix interface; (iii) the local environment of cofactors is not involved in the conformational dynamics which stabilizes the P+QA-; (iv) this conformational relaxation appears to be rather delocalized over several aminoacidic residues as well as water molecules weakly hydrogen-bonded to the RC.

High protein flexibility and reduced hydration water dynamics are key pressure adaptive strategies in prokaryotes

Water and protein dynamics on a nanometer scale were measured by quasi-elastic neutron scattering in the piezophile archaeon Thermococcus barophilus and the closely related pressure-sensitive Thermococcus kodakarensis, at 0.1 and 40 MPa. We show that cells of the pressure sensitive organism exhibit higher intrinsic stability. Both the hydration water dynamics and the fast protein and lipid dynamics are reduced under pressure. In contrast, the proteome of T. barophilus is more pressure sensitive than that of T. kodakarensis. The diffusion coefficient of hydration water is reduced, while the fast protein and lipid dynamics are slightly enhanced with increasing pressure. These findings show that the coupling between hydration water and cellular constituents might not be simply a master-slave relationship. We propose that the high flexibility of the T. barophilus proteome associated with reduced hydration water may be the keys to the molecular adaptation of the cells to high hydrostatic pressure.

Optimizing Computational Frameworks to Study the Influence of the Protein Environment on the Individual Site Energies of Chromophores in Photosystem II of Photosynthesis

Photosynthesis is a process in which electromagnetic radiation is converted into chemical energy. Photosystems capture photons with chromophores and transfer their energy to reaction centers using chromophores as a medium. In the reaction center, the excitation energy is used to perform chemical reactions. Knowledge of chromophore site energies is crucial to the understanding of excitation energy transfer pathways in photosystems and the ability to compute the site energies in a fast and accurate manner is mandatory for investigating how protein dynamics ef-fect the site energies and ultimately energy pathways with time. In this work we developed two software frameworks designed to optimize the calculations of chro-mophore site energies within a protein environment. The first is for performing quantum mechanical energy optimizations on molecules and the second is for com-puting site energies of chromophores in a fast and accurate manner using the polar-izability embedding method. The two frameworks allow for the fast and accurate calculation of chromophore site energies within proteins, ultimately allowing for the effect of protein dynamics on energy pathways to be studied. We use these frame-works to compute the site energies of the eight chromophores in the reaction center of photosystem II (PSII) using a 1.9 Å resolution x-ray structure of photosystem II. We compare our results to conflicting experimental data obtained from both isolat-ed intact PSII core preparations and the minimal reaction center preparation of PSII, and find our work more supportive of the former.

Ligand dynamics in myoglobin: calculation of infrared spectra for photodissociated NO

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Combining structure and dynamics: non-denaturing high-pressure effect on lysozyme in solution

Small-angle X-ray scattering (SAXS) and elastic and quasi-elastic neutron scattering techniques were used to investigate the high-pressure-induced changes on interactions, the low-resolution structure and the dynamics of lysozyme in solution. SAXS data, analysed using a global-fit procedure based on a new approach for hydrated protein form factor description, indicate that lysozyme completely maintains its globular structure up to 1500 bar, but significant modi. cations in the protein-protein interaction potential occur at approximately 600-1000 bar. Moreover, the mass density of the protein hydration water shows a clear discontinuity within this pressure range. Neutron scattering experiments indicate that the global and the local lysozyme dynamics change at a similar threshold pressure. A clear evolution of the internal protein dynamics from diffusing to more localized motions has also been probed. Protein structure and dynamics results have then been discussed in the context of protein-water interface and hydration water dynamics. According to SAXS results, the new configuration of water in the first hydration layer induced by pressure is suggested to be at the origin of the observed local mobility changes.