Lifetimes of intermediates in the β-sheet to α-helix transition of β-lactoglobulin by using a diffusional IR mixer

The extremely slow α-helix/β-sheet transition of proteins is a crucial step in amylogenic diseases and represents an internal rearrangement of local contacts in an already folded protein. These internal structural rearrangements within an already folded protein are a critical aspect of biological action and are a product of conformational flow along unknown metastable local minima of the energy landscape of the compact protein. We use a diffusional IR mixer with time-resolved Fourier transform IR spectroscopy capable of 400-μs time resolution to show that the trifluoroethanol driven β-sheet to α-helix transition of β-lactoglobulin proceeds via a compact β-sheet intermediate with a lifetime of 7 ms, small compared with the overall folding time of β-lactoglobulin.

Barriers to protein folding: formation of buried polar interactions is a slow step in acquisition of structure.

In the MYL mutant of the Arc repressor dimer, sets of partially buried salt-bridge and hydrogen-bond interactions mediated by Arg-31, Glu-36, and Arg-40 in each subunit are replaced by hydrophobic interactions between Met-31, Tyr-36, and Leu-40. The MYL refolding/dimerization reaction differs from that of wild type in being 10- to 1250-fold faster, having an earlier transition state, and depending upon viscosity but not ionic strength. Formation of the wild-type salt bridges in a hydrophobic environment clearly imposes a kinetic barrier to folding, which can be lowered by high salt concentrations. The changes in the position of the transition state and viscosity dependence can be explained if denatured monomers interact to form a partially folded dimeric intermediate, which then continues folding to form the native dimer. The second step is postulated to be rate limiting for wild type. Replacing the salt bridge with hydrophobic interactions lowers this barrier for MYL. This makes the first kinetic barrier rate limiting for MYL refolding and creates a downhill free-energy landscape in which most molecules which reach the intermediate state continue to form native dimers.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

Foldons, protein structural modules, and exons.

Foldons, which are kinetically competent, quasi-independently folding units of a protein, may be defined using energy landscape analysis. Foldons can be identified by maxima in a scan of the ratio of a contiguous segment's energetic stability gap to the energy variance of that segment's molten globule states, reflecting the requirement of minimal frustration. The predicted foldons are compared with the exons and structural modules for 16 of the 30 proteins studied. Statistical analysis indicates a strong correlation between the energetically determined foldons and Go's geometrically defined structural modules, but there are marked sequence-dependent effects. There is only a weak correlation of foldons to exons. For gammaII-crystallin, myoglobin, barnase, alpha-lactalbumin, and cytochrome c the foldons and some noncontiguous clusters of foldons compare well with intermediates observed in experiment.

A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding.

We propose a general mean field model of ligand-protein interactions to determine the thermodynamic equilibrium of a system at finite temperature. The method is employed in structural assessments of two human immuno-deficiency virus type 1 protease complexes where the gross effects of protein flexibility are incorporated by utilizing a data base of crystal structures. Analysis of the energy spectra for these complexes has revealed that structural and thermo-dynamic aspects of molecular recognition can be rationalized on the basis of the extent of frustration in the binding energy landscape. In particular, the relationship between receptor-specific binding of these ligands to human immunodeficiency virus type 1 protease and a minimal frustration principle is analyzed.

Rapid local synchronization of action potentials: toward computation with coupled integrate-and-fire neurons.

The collective behavior of interconnected spiking nerve cells is investigated. It is shown that a variety of model systems exhibit the same short-time behavior and rapidly converge to (approximately) periodic firing patterns with locally synchronized action potentials. The dynamics of one model can be described by a downhill motion on an abstract energy landscape. Since an energy landscape makes it possible to understand and program computation done by an attractor network, the results will extend our understanding of collective computation from models based on a firing-rate description to biologically more realistic systems with integrate-and-fire neurons.

Error-correcting codes on a Bethe-like lattice

We analyse Gallager codes by employing a simple mean-field approximation that distorts the model geometry and preserves important interactions between sites. The method naturally recovers the probability propagation decoding algorithm as a minimization of a proper free-energy. We find a thermodynamical phase transition that coincides with information theoretical upper-bounds and explain the practical code performance in terms of the free-energy landscape.

Bright and dark vector rogue waves

For an Erbium-doped mode locked fibre laser, we demonstrate experimentally a new type of vector rogue waves (RWs) emergence of which is caused by the coherent coupling of the orthogonal states of polarisation (SOPs). Unlike weak interaction between neighbouring dissipative solitons for the soliton rain, this creates a new type of the energy landscape where the interaction of the orthogonal SOPs leads to polarisation trapping or escapes from the trapping triggered by polarisation instabilities and so results in the pulse dynamics satisfying criteria of the 'dark' and 'bright' RWs. The obtained results, apart from the fundamental interest, can provide a base for development of the rogue waves mitigation techniques in the context of the applications in photonics and beyond.

Numerical detection of the Gardner transition in a mean-field glass former.

Recent theoretical advances predict the existence, deep into the glass phase, of a novel phase transition, the so-called Gardner transition. This transition is associated with the emergence of a complex free energy landscape composed of many marginally stable sub-basins within a glass metabasin. In this study, we explore several methods to detect numerically the Gardner transition in a simple structural glass former, the infinite-range Mari-Kurchan model. The transition point is robustly located from three independent approaches: (i) the divergence of the characteristic relaxation time, (ii) the divergence of the caging susceptibility, and (iii) the abnormal tail in the probability distribution function of cage order parameters. We show that the numerical results are fully consistent with the theoretical expectation. The methods we propose may also be generalized to more realistic numerical models as well as to experimental systems.

Topics in equilibrium and nonequilibrium thermodynamics: computing crystalline free energies and engineering Maxwell’s demon.

This dissertation covers two separate topics in statistical physics. The first part of the dissertation focuses on computational methods of obtaining the free energies (or partition functions) of crystalline solids. We describe a method to compute the Helmholtz free energy of a crystalline solid by direct evaluation of the partition function. In the many-dimensional conformation space of all possible arrangements of N particles inside a periodic box, the energy landscape consists of localized islands corresponding to different solid phases. Calculating the partition function for a specific phase involves integrating over the corresponding island. Introducing a natural order parameter that quantifies the net displacement of particles from lattices sites, we write the partition function in terms of a one-dimensional integral along the order parameter, and evaluate this integral using umbrella sampling. We validate the method by computing free energies of both face-centered cubic (FCC) and hexagonal close-packed (HCP) hard sphere crystals with a precision of $10^{-5}k_BT$ per particle. In developing the numerical method, we find several scaling properties of crystalline solids in the thermodynamic limit. Using these scaling properties, we derive an explicit asymptotic formula for the free energy per particle in the thermodynamic limit. In addition, we describe several changes of coordinates that can be used to separate internal degrees of freedom from external, translational degrees of freedom. The second part of the dissertation focuses on engineering idealized physical devices that work as Maxwell's demon. We describe two autonomous mechanical devices that extract energy from a single heat bath and convert it into work, while writing information onto memory registers. Additionally, both devices can operate as Landauer's eraser, namely they can erase information from a memory register, while energy is dissipated into the heat bath. The phase diagrams and the efficiencies of the two models are solved and analyzed. These two models provide concrete physical illustrations of the thermodynamic consequences of information processing.

Absorption and luminescent excitation spectra of F and M centers in Kcl and NaF

Luminescent excitation spectra were measured for the F and M centers in KCl; in particular, for the F band, M band, and the M2 transition. In all 3 cases, the spectra were nearly double-Gaussian in shape, and the efficiency for luminescence was nearly independent of the wavelength of the exciting light. A comparison of the absorption spectrum with the excitation spectrum of the F-band region of crystals with M centers present and oriented provided further evidence for the existence of the M2 transition of van Doorn and Haven and of Okamoto, and against the energy transfer theory of Lambe and Compton. The efficiency for luminescence of the M center upon M-band excitation was equal to the efficiency for F centers in pulse-annealed crystals of low F-center concentrations. The ratio of the efficiencies of the Ml to M2 transitions was 1.2 ± .25. The oscillator strengths of 3 of the M-center transitions in KCl relative to the oscillator strength for the F center were found to be in better agreement with the results reported by Okamoto, than with the results reported by Delbecq. The polarization of luminescence of M centers in KCl was measured at right angles to the exciting light, and was found to agree with the predictions of the van Doorn-Haven model of the M center. In NaF crystals having no absorption bands to the red side of the M band, the absorption and excitation spectra of the M band were accurately double-Gaussian over a wide range of wavelengths; the efficiency of luminescence of the M center was independent of the wavelength of the exciting light in that range; and the polarization of luminescence upon M-band excitation agreed well with the calculations based on the van DoornHaven model of the M center, In crystals in which the F band was bleached sufficiently to make it smaller in absorption height than the M band, several new color centers appeared on the red side of the M band, in contrast to the results reported by Blum; in these crystals, the polarization of luminescence of the M center upon M-band excitation disagreed strongly with theory, even though the absorptions for the new color centers were small compared to the M-band absorption.

Global change and climate-driven invasion of the Pacific oyster (Crassostrea gigas) along European coasts: a bioenergetics modelling approach

Aim The spread of non-indigenous species in marine ecosystems world-wide is one of today's most serious environmental concerns. Using mechanistic modelling, we investigated how global change relates to the invasion of European coasts by a non-native marine invertebrate, the Pacific oyster Crassostrea gigas. Location Bourgneuf Bay on the French Atlantic coast was considered as the northern boundary of C. gigas expansion at the time of its introduction to Europe in the 1970s. From this latitudinal reference, variations in the spatial distribution of the C. gigas reproductive niche were analysed along the north-western European coast from Gibraltar to Norway. Methods The effects of environmental variations on C. gigas physiology and phenology were studied using a bioenergetics model based on Dynamic Energy Budget theory. The model was forced with environmental time series including in situ phytoplankton data, and satellite data of sea surface temperature and suspended particulate matter concentration. Results Simulation outputs were successfully validated against in situ oyster growth data. In Bourgneuf Bay, the rise in seawater temperature and phytoplankton concentration has increased C. gigas reproductive effort and led to precocious spawning periods since the 1960s. At the European scale, seawater temperature increase caused a drastic northward shift (1400 km within 30 years) in the C. gigas reproductive niche and optimal thermal conditions for early life stage development. Main conclusions We demonstrated that the poleward expansion of the invasive species C. gigas is related to global warming and increase in phytoplankton abundance. The combination of mechanistic bioenergetics modelling with in situ and satellite environmental data is a valuable framework for ecosystem studies. It offers a generic approach to analyse historical geographical shifts and to predict the biogeographical changes expected to occur in a climate-changing world.

‘Mo Rùn am Fearann’ – ‘My Love is the Land’: Gaelic landscapes of the 18th and 19th centuries

The period of the 18th and 19th centuries was one of great change in the Highlands and Islands of Scotland. Improvement, the Jacobite Rebellions, and the Clearances transformed its communities and landscapes. These events have rightly been a focus of research. However, archaeological approaches have often sought simply to illustrate these processes, rather than create new narratives about life in the past. The resulting picture of the period can over-emphasise economic change whilst failing to reflect the richness and variety of everyday life. This thesis aims to suggest a new approach to the place and period, one which addresses matters often ignored in previous work. Whilst it has an archaeological sensibility, it draws on ideas from outside archaeology, such as landscape theory and on Gaelic oral tradition, an underused resource, to create a novel and broad-based approach to the period. An important part of the method is a synchronic approach that seeks to reconstruct the experience of the landscape at very particular times, engaging fully with the everyday experience of landscape rather than grand historical narratives. Two Hebridean case studies are utilised: Hiort (St Kilda) and Loch Aoineart, South Uist. Thematic discussions drawn from these landscapes are intended as critical assessments of the efficacy of the approach, as well as new narratives about life in the past in themselves. The thesis concludes by comparing the two case studies, reflecting on the merits of the approach, discussing recurrent themes in the work, and considering its wider context and implications. It is concluded that taking a novel approach to the case study landscapes can create narratives that often contrast or expand upon those produced by previous scholars, allow for a more detailed consideration of everyday life in the period, and open up new areas for archaeological enquiry. The extensive and critical use of evidence from Gaelic oral tradition is highlighted as crucial in understanding life and society in the period. The thesis questions the utility of grand historical narratives as a framework for archaeological study of post-medieval Gaeldom and suggests that our understanding of the past is best served by approaching the evidence in ways which allows for many different voices and stories from the past to emerge.

Ação de anestésicos gerais em canais iônicos

Dissertação (mestrado)—Universidade de Brasília, Instituto de Ciências Biológicas, Departamento de Biologia Molecular, 2016.