Thermal fluctuations of haemoglobin from different species: adaptation to temperature via conformational dynamics

Human Hemoglobin and Antipsychotics Clozapine, Ziprasidone and Sertindole: Friends or Foes?

Packed with hemoglobin, an essential protein for oxygen transport, human erythrocytes are a suitable model system for testing the pleiotropic effects of lipophilic drugs. Our study investigated the interaction between antipsychotic drugs clozapine, ziprasidone, sertindole, and human hemoglobin under simulated physiological conditions. Analysis of protein fluorescence quenching at different temperatures and data obtained from the van't Hoff diagram and molecular docking indicate that the interactions are static and that the tetrameric human hemoglobin has one binding site for all drugs in the central cavity near αβ interfaces and is dominantly mediated through hydrophobic forces. The association constants were lower-moderate strength (~10⁴ M^-1), the highest observed for clozapine (2.2 × 10⁴ M^-1 at 25 °C). The clozapine binding showed "friendly" effects: increased α-helical content, a higher melting point, and protein protection from free radical-mediated oxidation. On the other hand, bound ziprasidone and sertindole had a slightly pro-oxidative effect, increasing ferrihemoglobin content, a possible "foe". Since the interaction of proteins with drugs plays a vital role in their pharmacokinetic and pharmacodynamic properties, the physiological significance of the obtained findings is briefly discussed.

A novel single sensor hemoglobin domain from the thermophilic cyanobacteria Thermosynechococcus elongatus BP-1 exhibits higher pH but lower thermal stability compared to globins from mesophilic organisms

Thermosynechococcus elongatus-BP1 belongs to the class of photoautotrophic cyanobacterial organisms. The presence of chlorophyll a, carotenoids, and phycocyanobilin are the characteristics that categorize T. elongatus as a photosynthetic organism. Here, we report the structural and spectroscopic characteristics of novel hemoglobin (Hb) Synel Hb from T.elongatus, synonymous with Thermosynechococcus vestitus BP-1. The X-ray crystal structure (2.15 Å) of Synel Hb suggests the presence of a globin domain with a pre-A helix similar to the sensor domain (S) family of Hbs. The rich hydrophobic core accommodates heme in a penta-coordinated state and readily binds an extraneous ligand(imidazole). The absorption and circular dichroic spectral analysis of Synel Hb reiteratedthat the heme is in Fe^III+ state with a predominantly α-helical structure similar to myoglobin. Synel Hb displays higher resistance to structural perturbations induced via external stresses like pH and guanidium hydrochloride, which is comparable to Synechocystis Hb. However, Synel Hb exhibited lower thermal stability compared to mesophilic hemoglobins. Overall, the data is suggestive of the structural sturdiness of Synel Hb, which probably corroborates its origin in extreme thermophilic conditions. The stable globin provides scope for further investigation and may lead to new insights with scope for engineering stability in hemoglobin-based oxygen carriers.

Design and simulation of a caprylic acid enzymatically modified phosphatidylcholine micelle using a coarse-grained molecular dynamics simulations approach

Computationally simulated micelle models provide useful structural information on the molecular and biological sciences. One strategy to study the self-aggregation process of surfactant molecules that make up a micelle is through molecular dynamics (MD) simulations. In this study, a theoretical approach with a coarse-grained MD simulation (CG-MD) was employed to evaluate the critical micellar concentration (CMC), the micellization process, building a tridimensional (3D) model system of a micelle using data from the experimentally enzymatically modified phospholipids (PL) by phospholipase A1 (PA1). This required enzymatic interesterification of soybean phosphatidylcholine (PC) with caprylic acid, along with purification and characterization by chromatographic techniques to measure the esterified fatty acids and the corresponding PL composition. The number of molecules used in the CG-MD simulation system was determined from the experimental CMC data which was 0.025%. The molecular composition of the system is: 1 C 18:2, 2 C 8:0/8:0, 3 C 8:0/18:3n-9, 4 C 8:0/18:0, 5 C8:0/18:2n-6, 6 C8:0/18:1n-9, and 7 C 8:0/16:0. According to our theoretical results, the micelle model is structurally stable with an average Rg of 3.64 ± 0.10 Å, and might have an elliptical form with a radius of 24.6 Å. Regarding CMC value there was a relationship between the experimental data of the modified PLs and the theoretical analysis by GC-MD, which suggest that the enzymatic modification of PLs does not affect their self-aggregation properties. Finally, the micellar system obtained in the current research can be used as a simple and useful model to design optimal biocompatible nanoemulsions as possible vehicles for bioactive small molecules.Communicated by Ramaswamy H. Sarma.

Modeling the Dynamics of Protein–Protein Interfaces, How and Why?

Protein–protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein–protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein–protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein–protein interface.

Regional thermal hyperemia in the human leg: Evidence of the importance of thermosensitive mechanisms in the control of the peripheral circulation

Hyperthermia is thought to increase limb blood flow through the activation of thermosensitive mechanisms within the limb vasculature, but the precise vascular locus in which hyperthermia modulates perfusion remains elusive. We tested the hypothesis that local temperature-sensitive mechanisms alter limb hemodynamics by regulating microvascular blood flow. Temperature and oxygenation profiles and leg hemodynamics of the common (CFA), superficial (SFA) and profunda (PFA) femoral arteries, and popliteal artery (POA) of the experimental and control legs were measured in healthy participants during: (1) 3 h of whole leg heating (WLH) followed by 3 h of recovery (n = 9); (2) 1 h of upper leg heating (ULH) followed by 30 min of cooling and 1 h ULH bout (n = 8); and (3) 1 h of lower leg heating (LLH) (n = 8). WLH increased experimental leg temperature by 4.2 ± 1.2ºC and blood flow in CFA, SFA, PFA, and POA by ≥3-fold, while the core temperature essentially remained stable. Upper and lower leg blood flow increased exponentially in response to leg temperature and then declined during recovery. ULH and LLH similarly increased the corresponding segmental leg temperature, blood flow, and tissue oxygenation without affecting these responses in the non-heated leg segment, or perfusion pressure and conduit artery diameter across all vessels. Findings demonstrate that whole leg hyperthermia induces profound and sustained elevations in upper and lower limb blood flow and that segmental hyperthermia matches the regional thermal hyperemia without causing thermal or hemodynamic alterations in the non-heated limb segment. These observations support the notion that heat-activated thermosensitive mechanisms in microcirculation regulate limb tissue perfusion during hyperthermia.

Photo-Switching of Protein Dynamical Collectivity

We examine changes in the picosecond structural dynamics with irreversible photobleaching of red fluorescent proteins (RFP) mCherry, mOrange2 and TagRFP-T. Measurements of the protein dynamical transition using terahertz time-domain spectroscopy show in all cases an increase in the turn-on temperature in the bleached state. The result is surprising given that there is little change in the protein surface, and thus, the solvent dynamics held responsible for the transition should not change. A spectral analysis of the measurements guided by quasiharmonic calculations of the protein absorbance reveals that indeed the solvent dynamical turn-on temperature is independent of the thermal stability/photostate however the protein dynamical turn-on temperature shifts to higher temperatures. This is the first demonstration of switching the protein dynamical turn-on temperature with protein functional state. The observed shift in protein dynamical turn-on temperature relative to the solvent indicates an increase in the required mobile waters necessary for the protein picosecond motions, that is, these motions are more collective. Melting-point measurements reveal that the photobleached state is more thermally stable, and structural analysis of related RFP’s shows that there is an increase in internal water channels as well as a more uniform atomic root mean squared displacement. These observations are consistent with previous suggestions that water channels form with extended light excitation providing O2 access to the chromophore and subsequent fluorescence loss. We report that these same channels increase internal coupling enhancing thermal stability and collectivity of the picosecond protein motions. The terahertz spectroscopic characterization of the protein and solvent dynamical onsets can be applied generally to measure changes in collectivity of protein motions.

Effect of red blood cell shape changes on haemoglobin interactions and dynamics: a neutron scattering study

By using a combination of experimental neutron scattering techniques, it is possible to obtain a statistical perspective on red blood cell (RBC) shape in suspensions, and the inter-relationship with protein interactions and dynamics inside the confinement of the cell membrane. In this study, we examined the ultrastructure of RBC and protein-protein interactions of haemoglobin (Hb) in them using ultra-small-angle neutron scattering and small-angle neutron scattering (SANS). In addition, we used the neutron backscattering method to access Hb motion on the ns time scale and Å length scale. Quasi-elastic neutron scattering (QENS) experiments were performed to measure diffusive motion of Hb in RBCs and in an RBC lysate. By using QENS, we probed both internal Hb dynamics and global protein diffusion, on the accessible time scale and length scale by QENS. Shape changes of RBCs and variation of intracellular Hb concentration were induced by addition of the Na⁺-selective ionophore monensin and the K⁺-selective one, valinomycin. The experimental SANS and QENS results are discussed within the framework of crowded protein solutions, where free motion of Hb is obstructed by mutual interactions.

The molecular dynamics of bacterial spore and the role of calcium dipicolinate in core properties at the sub-nanosecond time-scale

Bacterial spores are among the most resistant forms of life on Earth. Their exceptional resistance properties rely on various strategies, among them the core singular structure, organization and hydration. By using elastic incoherent neutron scattering, we probed the dynamics of Bacillus subtilis spores to determine whether core macromolecular motions at the sub-nanosecond timescale could also contribute to their resistance to physical stresses. In addition, in order to better specify the role of the various spore components, we used different mutants lacking essential structure such as the coat (PS4150 mutant), or the calcium dipicolinic acid complex (CaDPA) located in the core (FB122 mutant). PS4150 allows to better probe the core’s dynamics, as proteins of the coat represent an important part of spore proteins, and FB122 gives information about the role of the large CaDPA depot for the mobility of core’s components. We show that core’s macromolecular mobility is not particularly constrained at the sub-nanosecond timescale in spite of its low water content as some dynamical characteristics as force constants are very close to those of vegetative bacteria such as Escherichia coli or to those of fully hydrated proteins. Although the force constants of the coatless mutant are similar to the wild-type’s ones, it has lower mean square displacements (MSDs) at high Q showing that core macromolecules are somewhat more constrained than the rest of spore components. However, no behavior reflecting the glassy state regularly evoked in the literature could be drawn from our data. As hydration and macromolecules’ mobility are highly correlated, the previous assumption, that core low water content might explain spores’ exceptional resistance properties seems unlikely. Thus, we confirm recent theories, suggesting that core water is mostly as free as bulk water and proteins/macromolecules are fully hydrated. The germination of spores leads to a much less stable system with a force constant of 0.1 N/m and MSDs ~2.5 times higher at low Q than in the dormant state. DPA has also an influence on core mobility with a slightly lower force constant for the DPA-less mutant than for the wild-type, and MSDs that are ~ 1.8 times higher on average than for the wild-type at low Q. At high Q, germinated and DPA-less spores were very similar to the wild-type ones, showing that DPA and core compact structure might influence large amplitude motions rather than local dynamics of macromolecules.

Ternary Complex Formation and Photoactivation of a Photoenzyme Results in Altered Protein Dynamics

The interplay between protein dynamics and catalysis remains a fundamental question in enzymology. We here investigate the ns-timescale dynamics of a light-dependent NADPH:protochlorophyllide oxidoreductase (LPOR), a photoenzyme crucial for chlorophyll synthesis. LPORs catalyze the light-triggered trans addition of a hydride and a proton across the C17═C18 double bond of the chlorophyll precursor protochlorophyllide (Pchlide). Because of the lack of an LPOR structure, the global structural and dynamic consequences of LPOR/Pchlide/NADPH ternary complex formation remain elusive. Moreover, photoactivation of LPORs by low-light preillumination is controversially discussed as unequivocal proof for this phenomenon is lacking. By employing quasielastic neutron spectroscopy (QENS), we show that the formation of the ternary holoprotein complex as well as photoactivation lead to progressive rigidification of the protein. These findings are supported by thermostability measurements, which reveal different melting behavior and thermostabilities for the apo- and holoprotein ternary complexes. Molecular dynamics simulations in good agreement with the experimental QENS results suggest that the increased flexibility observed for the apoprotein stems from structural fluctuations of the NADPH and Pchlide substrate binding sites of the enzyme. On the basis of our results, in conjunction with activity and stability measurements, we provide independent proof for LPOR photoactivation, defined as a process that modifies the protein structure and dynamics, resulting in an increased substrate turnover. Our findings advance the structural and dynamic understanding of LPORs and provide a first link between protein dynamics and catalysis for this enzyme class.

Temperature dependence of nanoscale dynamic processes measured in living millipedes by high resolution inelastic and elastic neutron scattering

We have used high energy-resolution neutron scattering to probe nanoscale dynamic processes in living millipedes (Narceus americanus). We have measured the temperature dependence of the intensity of scattered neutrons that do not exchange energy with the living samples on the 1.5 ns time scale, thereby excluding the signal from the highly mobile intra- and extra-cellular bulk-like aqueous constituents in the sample. This measured “elastic” scattering intensity exhibits a non-monotonic temperature dependence, with a noticeable systematic decrease detected between 295 and 303 K on warming up from 283 to 310 K. This decrease demonstrates an excellent inverse correlation with the non-monotonic, as a function of temperature, increase in the slow diffusivity previously observed in planarian flatworms and housefly larvae. This correlation suggests the existence of a biological mechanism, possibly common between different classes (Insects and Myriapods) and even phyla (Arthropods and Platyhelminthes), that dampens the slow nanoscopic dynamics in ectothermic organisms in response to the temperature of the environment exceeding the physiologically optimal range.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Dynamics of a family of cyan fluorescent proteins probed by incoherent neutron scattering

Cyan fluorescent proteins (CFPs) are variants of green fluorescent proteins in which the central tyrosine of the chromophore has been replaced by a tryptophan. The increased bulk of the chromophore within a compact protein and the change in the positioning of atoms capable of hydrogen bonding have made it difficult to optimize their fluorescence properties, which took approximately 15 years between the availability of the first useable CFP, enhanced cyan fluorescent protein (ECFP), and that of a variant with almost perfect fluorescence efficiency, mTurquoise2. To understand the molecular bases of the progressive improvement in between these two CFPs, we have studied by incoherent neutron scattering the dynamics of five different variants exhibiting progressively increased fluorescence efficiency along the evolution pathway. Our results correlate well with the analysis of the previously determined X-ray crystallographic structures, which show an increase in flexibility between ECFP and the second variant, Cerulean, which is then hindered in the three later variants, SCFP3A (Super Cyan Fluorescent Protein 3A), mTurquoise and mTurquoise2. This confirms that increasing the rigidity of the direct environment of the fluorescent chromophore is not the sole parameter leading to brighter fluorescent proteins and that increased flexibility in some cases may be helpful.

Mechanisms for the control of local tissue blood flow during thermal interventions: influence of temperature-dependent ATP release from human blood and endothelial cells

New Findings

What is the central question of this study?

Skin and muscle blood flow increases with heating and decreases with cooling, but the temperature‐sensitive mechanisms underlying these responses are not fully elucidated.

What is the main finding and its importance?

We found that local tissue hyperaemia was related to elevations in ATP release from erythrocytes. Increasing intravascular ATP augmented skin and tissue perfusion to levels equal or above thermal hyperaemia. ATP release from isolated erythrocytes was altered by heating and cooling. Our findings suggest that erythrocytes are involved in thermal regulation of blood flow via modulation of ATP release.

Local tissue perfusion changes with alterations in temperature during heating and cooling, but the thermosensitivity of the vascular ATP signalling mechanisms for control of blood flow during thermal interventions remains unknown. Here, we tested the hypotheses that the release of the vasodilator mediator ATP from human erythrocytes, but not from endothelial cells or other blood constituents, is sensitive to both increases and reductions in temperature and that increasing intravascular ATP availability with ATP infusion would potentiate thermal hyperaemia in limb tissues. We first measured blood temperature, brachial artery blood flow and plasma [ATP] during passive arm heating and cooling in healthy men and found that they increased by 3.0 ± 1.2°C, 105 ± 25 ml min⁻¹ °C⁻¹ and twofold, respectively, (all P < 0.05) with heating, but decreased or remained unchanged with cooling. In additional men, infusion of ATP into the brachial artery increased skin and deep tissue perfusion to levels equal or above thermal hyperaemia. In isolated erythrocyte samples exposed to different temperatures, ATP release increased 1.9‐fold from 33 to 39°C (P < 0.05) and declined by ∼50% at 20°C (P < 0.05), but no changes were observed in cultured human endothelial cells, plasma or serum samples. In conclusion, increases in plasma [ATP] and skin and deep tissue perfusion with limb heating are associated with elevations in ATP release from erythrocytes, but not from endothelial cells or other blood constituents. Erythrocyte ATP release is also sensitive to temperature reductions, suggesting that erythrocytes may function as thermal sensors and ATP signalling generators for control of tissue perfusion during thermal interventions.

The effect of the pathological V72I, D109N and T190M missense mutations on the molecular structure of α-dystroglycan

Dystroglycan (DG) is a highly glycosylated protein complex that links the cytoskeleton with the extracellular matrix, mediating fundamental physiological functions such as mechanical stability of tissues, matrix organization and cell polarity. A crucial role in the glycosylation of the DG α subunit is played by its own N-terminal region that is required by the glycosyltransferase LARGE. Alteration in this O-glycosylation deeply impairs the high affinity binding to other extracellular matrix proteins such as laminins. Recently, three missense mutations in the gene encoding DG, mapped in the α-DG N-terminal region, were found to be responsible for hypoglycosylated states, causing congenital diseases of different severity referred as primary dystroglycanopaties.To gain insight on the molecular basis of these disorders, we investigated the crystallographic and solution structures of these pathological point mutants, namely V72I, D109N and T190M. Small Angle X-ray Scattering analysis reveals that these mutations affect the structures in solution, altering the distribution between compact and more elongated conformations. These results, supported by biochemical and biophysical assays, point to an altered structural flexibility of the mutant α-DG N-terminal region that may have repercussions on its interaction with LARGE and/or other DG-modifying enzymes, eventually reducing their catalytic efficiency.

Determinants of neuroglobin plasticity highlighted by joint coarse-grained simulations and high pressure crystallography

Investigating the effect of pressure sheds light on the dynamics and plasticity of proteins, intrinsically correlated to functional efficiency. Here we detail the structural response to pressure of neuroglobin (Ngb), a hexacoordinate globin likely to be involved in neuroprotection. In murine Ngb, reversible coordination is achieved by repositioning the heme more deeply into a large internal cavity, the “heme sliding mechanism”. Combining high pressure crystallography and coarse-grain simulations on wild type Ngb as well as two mutants, one (V101F) with unaffected and another (F106W) with decreased affinity for CO, we show that Ngb hinges around a rigid mechanical nucleus of five hydrophobic residues (V68, I72, V109, L113, Y137) during its conformational transition induced by gaseous ligand, that the intrinsic flexibility of the F-G loop appears essential to drive the heme sliding mechanism, and that residue Val 101 may act as a sensor of the interaction disruption between the heme and the distal histidine.

The fluctuating ribosome: thermal molecular dynamics characterized by neutron scattering

Conformational changes associated with ribosome function have been identified by X-ray crystallography and cryo-electron microscopy. These methods, however, inform poorly on timescales. Neutron scattering is well adapted for direct measurements of thermal molecular dynamics, the ‘lubricant’ for the conformational fluctuations required for biological activity. The method was applied to compare water dynamics and conformational fluctuations in the 30 S and 50 S ribosomal subunits from Haloarcula marismortui, under high salt, stable conditions. Similar free and hydration water diffusion parameters are found for both subunits. With respect to the 50 S subunit, the 30 S is characterized by a softer force constant and larger mean square displacements (MSD), which would facilitate conformational adjustments required for messenger and transfer RNA binding. It has been shown previously that systems from mesophiles and extremophiles are adapted to have similar MSD under their respective physiological conditions. This suggests that the results presented are not specific to halophiles in high salt but a general property of ribosome dynamics under corresponding, active conditions. The current study opens new perspectives for neutron scattering characterization of component functional molecular dynamics within the ribosome.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.

Bridging Enzymatic Structure Function via Mechanics: A Coarse-Grain Approach

Flexibility is a central aspect of protein function, and ligand binding in enzymes involves a wide range of structural changes, ranging from large-scale domain movements to small loop or side-chain rearrangements. In order to understand how the mechanical properties of enzymes, and the mechanical variations that are induced by ligand binding, relate to enzymatic activity, we carried out coarse-grain Brownian dynamics simulations on a set of enzymes whose structures in the unbound and ligand-bound forms are available in the Protein Data Bank. Our results show that enzymes are remarkably heterogeneous objects from a mechanical point of view and that the local rigidity of individual residues is tightly connected to their part in the protein's overall structure and function. The systematic comparison of the rigidity of enzymes in their unbound and bound forms highlights the fact that small conformational changes can induce large mechanical effects, leading to either more or less flexibility depending on the enzyme's architecture and the location of its ligand-biding site. These mechanical variations target a limited number of specific residues that occupy key locations for enzymatic activity, and our approach thus offers a mean to detect perturbation-sensitive sites in enzymes, where the addition or removal of a few interactions will lead to important changes in the proteins internal dynamics.

Molecular adaptation and salt stress response of Halobacterium salinarum cells revealed by neutron spectroscopy

Halobacterium salinarum is an extreme halophile archaeon with an absolute requirement for a multimolar salt environment. It accumulates molar concentrations of KCl in the cytosol to counterbalance the external osmotic pressure imposed by the molar NaCl. As a consequence, cytosolic proteins are permanently exposed to low water activity and highly ionic conditions. In non-adapted systems, such conditions would promote protein aggregation, precipitation, and denaturation. In contrast, in vitro studies showed that proteins from extreme halophilic cells are themselves obligate halophiles. In this paper, adaptation via dynamics to low-salt stress in H. salinarum cells was measured by neutron scattering experiments coupled with microbiological characterization. The molecular dynamic properties of a proteome represent a good indicator for environmental adaptation and the neutron/microbiology approach has been shown to be well tailored to characterize these modifications. In their natural setting, halophilic organisms often have to face important variations in environmental salt concentration. The results showed deleterious effects already occur in the H. salinarum proteome, even when the external salt concentration is still relatively high, suggesting the onset of survival mechanisms quite early when the environmental salt concentration decreases.

Analysis of biosurfaces by neutron reflectometry: from simple to complex interfaces

Because of its high sensitivity for light elements and the scattering contrast manipulation via isotopic substitutions, neutron reflectometry (NR) is an excellent tool for studying the structure of soft-condensed material. These materials include model biophysical systems as well as in situ living tissue at the solid-liquid interface. The penetrability of neutrons makes NR suitable for probing thin films with thicknesses of 5-5000 Å at various buried, for example, solid-liquid, interfaces [J. Daillant and A. Gibaud, Lect. Notes Phys. 770, 133 (2009); G. Fragneto-Cusani, J. Phys.: Condens. Matter 13, 4973 (2001); J. Penfold, Curr. Opin. Colloid Interface Sci. 7, 139 (2002)]. Over the past two decades, NR has evolved to become a key tool in the characterization of biological and biomimetic thin films. In the current report, the authors would like to highlight some of our recent accomplishments in utilizing NR to study highly complex systems, including in-situ experiments. Such studies will result in a much better understanding of complex biological problems, have significant medical impact by suggesting innovative treatment, and advance the development of highly functionalized biomimetic materials.

Correlation between supercoiling and conformational motions of the bacterial flagellar filament

The bacterial flagellar filament is a very large macromolecular assembly of a single protein, flagellin. Various supercoiled states of the filament exist, which are formed by two structurally different conformations of flagellin in different ratios. We investigated the correlation between supercoiling of the protofilaments and molecular dynamics in the flagellar filament using quasielastic and elastic incoherent neutron scattering on the picosecond and nanosecond timescales. Thermal fluctuations in the straight L- and R-type filaments were measured and compared to the resting state of the wild-type filament. Amplitudes of motion on the picosecond timescale were found to be similar in the different conformational states. Mean-square displacements and protein resilience on the 0.1 ns timescale demonstrate that the L-type state is more flexible and less resilient than the R-type, whereas the wild-type state lies in between. Our results provide strong support that supercoiling of the protofilaments in the flagellar filament is determined by the strength of molecular forces in and between the flagellin subunits.

Structure and dynamics of G-protein coupled receptors

G-protein coupled receptors (GPCRs) are seven helical transmembrane proteins that mediate cell-to-cell communication. They also form the largest superfamily of drug targets. Hence detailed studies of the three dimensional structure and dynamics are critical to understanding the functional role of GPCRs in signal transduction pathways, and for drug design. In this chapter we compare the features of the crystal structures of various biogenic amine receptors, such as β1 and β2 adrenergic receptors, dopamine D3 receptor, M2 and M3 muscarinic acetylcholine receptors. This analysis revealed that conserved residues are located facing the inside of the transmembrane domain in these GPCRs improving the efficiency of packing of these structures. The NMR structure of the chemokine receptor CXCR1 without any ligand bound, shows significant dynamics of the transmembrane domain, especially the helical kink angle on the transmembrane helix6. The activation mechanism of the β2-adrenergic receptor has been studied using multiscale computational methods. The results of these studies showed that the receptor without any ligand bound, samples conformations that resemble some of the structural characteristics of the active state of the receptor. Ligand binding stabilizes some of the conformations already sampled by the apo receptor. This was later observed in the NMR study of the dynamics of human β2-adrenergic receptor. The dynamic nature of GPCRs leads to a challenge in obtaining purified receptors for biophysical studies. Deriving thermostable mutants of GPCRs has been a successful strategy to reduce the conformational heterogeneity and stabilize the receptors. This has lead to several crystal structures of GPCRs. However, the cause of how these mutations lead to thermostability is not clear. Computational studies are beginning to shed some insight into the possible structural basis for the thermostability. Molecular Dynamics simulations studying the conformational ensemble of thermostable mutants have shown that the stability could arise from both enthalpic and entropic factors. There are regions of high stress in the wild type GPCR that gets relieved upon mutation conferring thermostability.

Thermostabilization of the β1-adrenergic receptor correlates with increased entropy of the inactive state

The dynamic nature of GPCRs is a major hurdle in their purification and crystallization. Thermostabilization can facilitate GPCR structure determination, as has been shown by the structure of the thermostabilized β1-adrenergic receptor (β1AR) mutant, m23-β1AR, which has been thermostabilized in the inactive state. However, it is unclear from the structure how the six thermostabilizing mutations in m23-β1AR affect receptor dynamics. We have used molecular dynamics simulations in explicit solvent to compare the conformational ensembles for both wild type β1AR (wt-β1AR) and m23-β1AR. Thermostabilization results in an increase in the number of accessible microscopic conformational states within the inactive state ensemble, effectively increasing the side chain entropy of the inactive state at room temperature, while suppressing large-scale main chain conformational changes that lead to activation. We identified several diverse mechanisms of thermostabilization upon mutation. These include decrease of long-range correlated movement between residues in the G-protein coupling site to the extracellular region (Y227A(5.58), F338M(7.48)), formation of new hydrogen bonds (R68S), and reduction of local stress (Y227(5.58), F327(7.37), and F338(7.48)). This study provides insights into microscopic mechanisms underlying thermostability that leads to an understanding of the effect of these mutations on the structure of the receptor.