Dynamics of human acetylcholinesterase bound to non-covalent and covalent inhibitors shedding light on changes to the water network structure

Dynamic interchange between two protonation states is characteristic of active sites of cholinesterases

Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.

Role of fast dynamics in the complexation of G-quadruplexes with small molecules

G-quadruplexes (G4s) formed by the human telomeric sequence AG₃ (TTAG₃)₃ (Tel22) play a key role in cancer and aging. We combined elastic incoherent neutron scattering (EINS) and quasielastic incoherent neutron scattering (QENS) to characterize the internal dynamics of Tel22 G4s and to assess how it is affected by complexation with two standard ligands, Berberine and BRACO19. We show that the interaction with the two ligands induces an increase of the overall mobility of Tel22 as quantified by the mean squared displacements (MSD) of hydrogen atoms. At the same time, the complexes display a lower stiffness than G4 alone. Two different types of motion characterize the G4 nanosecond timescale dynamics. Upon complexation, an increasing fraction of G4 atomic groups participate in this fast dynamics, along with an increase in the relevant characteristic length scales. We suggest that the entropic contribution to the conformational free energy of these motions might be crucial for the complexation mechanisms.

Structural and dynamic effects of paraoxon binding to human acetylcholinesterase by X-ray crystallography and inelastic neutron scattering

Organophosphorus (OP) compounds, including nerve agents and some pesticides, covalently bind to the catalytic serine of human acetylcholinesterase (hAChE), thereby inhibiting acetylcholine hydrolysis necessary for efficient neurotransmission. Oxime antidotes can reactivate the OP-conjugated hAChE, but reactivation efficiency can be low for pesticides, such as paraoxon (POX). Understanding structural and dynamic determinants of OP inhibition and reactivation can provide insights to design improved reactivators. Here, X-ray structures of hAChE with unaged POX, with POX and oximes MMB4 and RS170B, and with MMB4 are reported. A significant conformational distortion of the acyl loop was observed upon POX binding, being partially restored to the native conformation by oximes. Neutron vibrational spectroscopy combined with molecular dynamics simulations showed that picosecond vibrational dynamics of the acyl loop soften in the ∼20-50 cm-1 frequency range. The acyl loop structural perturbations may be correlated with its picosecond vibrational dynamics to yield more comprehensive template for structure-based reactivator design.

Broad‐Spectrum Antidote Discovery by Untangling the Reactivation Mechanism of Nerve‐Agent‐Inhibited Acetylcholinesterase

Reactivators are vital for the treatment of organophosphorus nerve agent (OPNA) intoxication but new alternatives are needed due to their limited clinical applicability. The toxicity of OPNAs stems from covalent inhibition of the essential enzyme acetylcholinesterase (AChE), which reactivators relieve via a chemical reaction with the inactivated enzyme. Here, we present new strategies and tools for developing reactivators. We discover suitable inhibitor scaffolds by using an activity‐independent competition assay to study non‐covalent interactions with OPNA‐AChEs and transform these inhibitors into broad‐spectrum reactivators. Moreover, we identify determinants of reactivation efficiency by analysing reactivation and pre‐reactivation kinetics together with structural data. Our results show that new OPNA reactivators can be discovered rationally by exploiting detailed knowledge of the reactivation mechanism of OPNA‐inhibited AChE.

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm^-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol^-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

Molecular Dynamics of Lysozyme Amyloid Polymorphs Studied by Incoherent Neutron Scattering

Lysozyme amyloidosis is a hereditary disease, which is characterized by the deposition of lysozyme amyloid fibrils in various internal organs. It is known that lysozyme fibrils show polymorphism and that polymorphs formed at near-neutral pH have the ability to promote more monomer binding than those formed at acidic pH, indicating that only specific polymorphs become dominant species in a given environment. This is likely due to the polymorph-specific configurational diffusion. Understanding the possible differences in dynamical behavior between the polymorphs is thus crucial to deepen our knowledge of amyloid polymorphism and eventually elucidate the molecular mechanism of lysozyme amyloidosis. In this study, molecular dynamics at sub-nanosecond timescale of two kinds of polymorphic fibrils of hen egg white lysozyme, which has long been used as a model of human lysozyme, formed at pH 2.7 (LP27) and pH 6.0 (LP60) was investigated using elastic incoherent neutron scattering (EINS) and quasi-elastic neutron scattering (QENS). Analysis of the EINS data showed that whereas the mean square displacement of atomic motions is similar for both LP27 and LP60, LP60 contains a larger fraction of atoms moving with larger amplitudes than LP27, indicating that the dynamical difference between the two polymorphs lies not in the averaged amplitude, but in the distribution of the amplitudes. Furthermore, analysis of the QENS data showed that the jump diffusion coefficient of atoms is larger for LP60, suggesting that the atoms of LP60 undergo faster diffusive motions than those of LP27. This study thus characterizes the dynamics of the two lysozyme polymorphs and reveals that the molecular dynamics of LP60 is enhanced compared with that of LP27. The higher molecular flexibility of the polymorph would permit to adjust its conformation more quickly than its counterpart, facilitating monomer binding.

Dynamics of Apolipoprotein B-100 in Interaction with Detergent Probed by Incoherent Neutron Scattering

Apolipoprotein B-100 (apo B-100) is the protein moiety of both low- and very-low-density lipoproteins, whose role is crucial to cholesterol and triglyceride transport. Aiming at the molecular dynamics' details of apo B-100, scarcely studied, we performed elastic and quasi-elastic incoherent neutron scattering (EINS, QENS) experiments combining different instruments and time scales. Similar to classical membrane proteins, the solubilization results in remaining detergent, here Nonidet P-40 (NP40). Therefore, we propose a framework for QENS studies of protein-detergent complexes, with the introduction of a combined model, including the experimental apo B-100/NP40 ratio. Relying on the simultaneous analysis of all QENS amplitudes, this approach is sensitive enough to separate both contributions. Its application identified two points: (i) apo B-100 slow dynamics and (ii) the acceleration of NP40 dynamics in the presence of apo B-100. Direct translation of the exposed methodology now makes the investigation of more membrane proteins by neutron spectroscopy achievable.

A Comprehensive Review of Cholinesterase Modeling and Simulation

The present article reviews published efforts to study acetylcholinesterase and butyrylcholinesterase structure and function using computer-based modeling and simulation techniques. Structures and models of both enzymes from various organisms, including rays, mice, and humans, are discussed to highlight key structural similarities in the active site gorges of the two enzymes, such as flexibility, binding site location, and function, as well as differences, such as gorge volume and binding site residue composition. Catalytic studies are also described, with an emphasis on the mechanism of acetylcholine hydrolysis by each enzyme and novel mutants that increase catalytic efficiency. The inhibitory activities of myriad compounds have been computationally assessed, primarily through Monte Carlo-based docking calculations and molecular dynamics simulations. Pharmaceutical compounds examined herein include FDA-approved therapeutics and their derivatives, as well as several other prescription drug derivatives. Cholinesterase interactions with both narcotics and organophosphate compounds are discussed, with the latter focusing primarily on molecular recognition studies of potential therapeutic value and on improving our understanding of the reactivation of cholinesterases that are bound to toxins. This review also explores the inhibitory properties of several other organic and biological moieties, as well as advancements in virtual screening methodologies with respect to these enzymes.

Substrate Binding Stiffens Aspartate Aminotransferase by Altering the Enzyme Picosecond Vibrational Dynamics

Protein dynamics on various time scales from femtoseconds to milliseconds impacts biological function by driving proteins to conformations conducive to ligand binding and creating functional states in enzyme catalysis. Neutron vibrational spectroscopy carried out by measuring inelastic neutron scattering from protein molecules in combination with molecular simulations has the unique ability of detecting and visualizing changes in the picosecond protein vibrational dynamics due to ligand binding. Here we present neutron vibrational spectra of a homodimeric pyridoxal 5'-phosphate-dependent enzyme, aspartate aminotransferase, obtained from the open internal aldimine and closed external aldimine conformational states. We observe that in the external aldimine state the protein structure stiffens relative to the internal aldimine state, indicating rigidified vibrational dynamics on the picosecond time scale in the low-frequency regime of 5-50 cm^-1. Our molecular dynamics simulations indicate substantial changes in the picosecond dynamics of the enzyme secondary structure elements upon substrate binding, with the largest contributions from just two helices and the β-sheet.

Quasielastic neutron scattering studies on couplings of protein and water dynamics in hydrated elastin

Performing quasielastic neutron scattering measurements and analyzing both elastic and quasielasic contributions, we study protein and water dynamics of hydrated elastin. At low temperatures, hydration-independent methyl group rotation dominates the findings. It is characterized by a Gaussian distribution of activation energies centered at about E_m = 0.17 eV. At ∼195 K, coupled protein-water motion sets in. The hydration water shows diffusive motion, which is described by a Gaussian distribution of activation energies with E_m = 0.57 eV. This Arrhenius behavior of water diffusion is consistent with previous results for water reorientation, but at variance with a fragile-to-strong crossover at ∼225 K. The hydration-related elastin backbone motion is localized and can be attributed to the cage rattling motion. We speculate that its onset at ∼195 K is related to a secondary glass transition, which occurs when a β relaxation of the protein has a correlation time of τ_β ∼ 100 s. Moreover, we show that its temperature-dependent amplitude has a crossover at the regular glass transition T_g = 320 K of hydrated elastin, where the α relaxation of the protein obeys τ_α ∼ 100 s. By contrast, we do not observe a protein dynamical transition when water dynamics enters the experimental time window at ∼240 K.

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.

Differences between calcium rich and depleted alpha-lactalbumin investigated by molecular dynamics simulations and incoherent neutron scattering

We present a study comparing atomic motional amplitudes in calcium rich and depleted alpha-lactalbumin. The investigations were performed by elastic incoherent neutron scattering (EINS) and molecular dynamics (MD) simulations. As the variations were expected to be very small, three different hydration levels and timescales (instrumental resolutions) were measured. In addition, we used two models to extract the mean square displacements (MSDs) from the EINS data, one taking into account the motional heterogeneity of the MSD. At a timescale of several nanoseconds, small differences in the amplitudes between the calcium enriched and depleted alpha-lactalbumin are visible, whereas at lower timescales no changes can be concluded within the statistics. The results are compared to MD simulations at 280 and 300 K by extracting the MSDs of the trajectories in two separate ways: first by direct calculation, and second by a virtual neutron experiment using the same models as for the experimental data. We show that the simulated data give qualitatively similar results as the experimental data but quantitatively there are differences. Furthermore, the distribution of the MSDs in the simulations suggests that the inclusion of heterogeneity is reasonable for alpha-lactalbumin, but a bi-or trimodal approach may be sufficient.

Frequency domain modeling of quasielastic neutron scattering from hydrated protein powders: Application to free and inhibited human acetylcholinesterase

This article reports on a frequency domain analysis of quasielastic neutron scattering spectra from free and Huperzine-A-inhibited human acetylcholinesterase, extending a recent time domain analysis of the same experimental data [M. Saouessi et al., J. Chem. Phys. 150, 161104 (2019)]. An important technical point here is the construction of a semianalytical model for the resolution-broadened dynamic structure factor that can be fitted to the experimental spectra. We find comparable parameters as in our previous study and demonstrate that our model is sensitive to subpercent changes in the experimental data, which are caused by reversible binding of the inhibitor Huperzine A.

In vitro and in silico evaluation of fucosterol from Sargassum horridum as potential human acetylcholinesterase inhibitor

The fucosterol has been reported numerous biological activities. In this study, the activity in vitro of the fucosterol from Sargassum horridum as potential human acetylcholinesterase inhibitor was evaluated. The structural identification was obtained by nuclear magnetic resonance (NMR) spectroscopy and based on experimental data, we combined docking and molecular dynamics simulations coupled to the molecular-mechanics-generalized-born-surface-area approach to evaluating the structural and energetic basis for the molecular recognition of fucosterol and neostigmine at the binding site of acetylcholinesterase (AChE). In addition, the Lineweaver-Burk plot showed the nature of a non-competitive inhibition. The maximum velocity (Vmax) and the constant of Michaelis-Menten (Km) estimated for fucosterol (0.006 µM) were 0.015 1/Vo (ΔA/h and 6.399 1/[ACh] mM-1, respectively. While, for neostigmine (0.14 µM), the Vmax was 0.022 1/Vo (ΔA/h) and Km of 6.726 1/[ACh] mM-1, these results showed a more effective inhibition by fucosterol respect to neostigmine. Structural analysis revealed that neostigmine reaches the AChE binding site reported elsewhere, whereas fucosterol can act as a no-competitive and competitive acetylcholinesterase inhibitor, in agree with kinetic enzymatic experiments. Binding free energy calculations revealed that fucosterol reaches the acetylcholinesterase binding site with higher affinity than neostigmine, which is according to experimental results. Whereas the per-residue decomposition free energy analysis let us identify crucial residues involved in the molecular recognition of ligands by AChE. Results corroborate the ability of theoretical methods to provide crucial information at the atomic level about energetic and structural differences in the binding interaction and affinity from fucosterol with AChE. Communicated by Ramaswamy H. Sarma.

Asymptotic analysis of quasielastic neutron scattering data from human acetylcholinesterase reveals subtle dynamical changes upon ligand binding

In this paper, we show that subtle changes in the internal dynamics of human acetylcholinesterase upon ligand binding can be extracted from quasielastic neutron scattering data by employing a nonexponential relaxation model for the intermediate scattering function. The relaxation is here described by a stretched Mittag-Leffler function, which exhibits slow power law decay for long times. Our analysis reveals that binding of a Huperzine A ligand increases the atomic motional amplitudes of the enzyme and slightly slows down its internal diffusive motions. This result is interpreted within an energy landscape picture for the motion of the hydrogen atoms.

Dynamics Properties of Photosynthetic Microorganisms Probed by Incoherent Neutron Scattering

Studies on the dynamical properties of photosynthetic membranes of land plants and purple bacteria have been previously performed by neutron spectroscopy, revealing a tight coupling between specific photochemical reactions and macromolecular dynamics. Here, we probed the intrinsic dynamics of biotechnologically useful mutants of the green alga Chlamydomonas reinhardtii by incoherent neutron scattering coupled with prompt chlorophyll fluorescence experiments. We brought to light that single amino acid replacements in the plastoquinone (PQ)-binding niche of the photosystem II D1 protein impair electron transport (ET) efficiency between quinones and confer increased flexibility to the host membranes, expanding to the entire cells. Hence, a more flexible environment in the PQ-binding niche has been associated to a less efficient ET. A similar function/dynamics relationship was also demonstrated in Rhodobacter sphaeroides reaction centers having inhibited ET, indicating that flexibility at the quinones region plays a crucial role in evolutionarily distant organisms. Instead, a different functional/dynamical correlation was observed in algal mutants hosting a single amino acid replacement residing in a D1 domain far from the PQ-binding niche. Noteworthy, this mutant displayed the highest degree of flexibility, and besides having a nativelike ET efficiency in physiological conditions, it acquired novel, to our knowledge, phenotypic traits enabling it to preserve a high maximal quantum yield of photosystem II photochemistry in extreme habitats. Overall, in the nanosecond timescale, the degree of the observed flexibility is related to the mutation site; in the picosecond timescale, we highlighted the presence of a more pronounced dynamic heterogeneity in all mutants compared to the native cells, which could be related to a marked chemically heterogeneous environment.

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.

Analysis of elastic incoherent neutron scattering data beyond the Gaussian approximation

This work addresses the use of the Gaussian approximation as a common tool to extract atomic motions in proteins from elastic incoherent neutron scattering and whether improvements in data analysis and additional information can be obtained when going beyond that. We measured alpha-lactalbumin with different levels of hydration on three neutron backscattering spectrometers, to be able to resolve a wide temporal and spatial range for dynamics. We demonstrate that the Gaussian approximation gives qualitatively similar results to models that include heterogeneity, if one respects a certain procedure to treat the intercept of the elastic intensities with the momentum transfer axis. However, the inclusion of motional heterogeneity provides better fits to the data. Our analysis suggests an approach of limited heterogeneity, where including only two kinds of motions appears sufficient to obtain more quantitative results for the mean square displacement. Finally, we note that traditional backscattering spectrometers pose a limit on the lowest accessible momentum transfer. We therefore suggest that complementary information about the spatial evolution of the elastic intensity close to zero momentum transfer can be obtained using other neutron methods, in particular, neutron spin-echo together with polarization analysis.

Water structure changes in oxime-mediated reactivation process of phosphorylated human acetylcholinesterase

The role of water in oxime-mediated reactivation of phosphylated cholinesterases (ChEs) has been asked with recurrence. To investigate oximate water structure changes in this reaction, reactivation of paraoxon-inhibited human acetylcholinesterase (AChE) was performed by the oxime asoxime (HI-6) at different pH in the presence and absence of lyotropic salts: a neutral salt (NaCl), a strong chaotropic salt (LiSCN) and strong kosmotropic salts (ammonium sulphate and phosphate HPO42-). At the same time, molecular dynamic (MD) simulations of enzyme reactivation under the same conditions were performed over 100 ns. Reactivation kinetics showed that the low concentration of chaotropic salt up to 75 mM increased the percentage of reactivation of diethylphosphorylated AChE whereas kosmotropic salts lead only to a small decrease in reactivation. This indicates that water-breaker salt induces destructuration of water molecules that are electrostricted around oximate ions. Desolvation of oximate favors nucleophilic attack on the phosphorus atom. Effects observed at high salt concentrations (>100 mM) result either from salting-out of the enzyme by kosmotropic salts (phosphate and ammonium sulphate) or denaturing action of chaotropic LiSCN. MDs simulations of diethylphosphorylated hAChE complex with HI-6 over 100 ns were performed in the presence of 100 mM (NH4)2SO4 and 50 mM LiSCN. In the presence of LiSCN, it was found that protein and water have a higher mobility, i.e. water is less organized, compared with the ammonium sulphate system. LiSCN favors protein solvation (hydrophobic hydration) and breakage of elelectrostricted water molecules around of oximate ion. As a result, more free water molecules participated to reaction steps accompanying oxime-mediated dephosphorylation.

EFFECT OF ORAL MICROBIOTA ON THE DEVELOPMENT OF INFLAMMATION AND SOMATIC DISEASES

The review presents the modern concept of the oral microbioma, its species diversity, the properties of microbial associations and their effect on the immunobiological resistance of the organism, the development of chronic inflammation and the two leading infectious and inflammatory diseases, dental caries and periodontal diseases. Briefly summarized information about the presence and nature of the relationship between the composition of periodontal microorganisms and concomitant somatic diseases: in atherosclerotic plaques in cardiovascular diseases, with pathology of the gastrointestinal tract, respiratory system and other organs and systems.

The effect of pH on the dynamics of natural membranes

Pure phospholipids and membrane fragments from bacterial cells living under various conditions were studied against the influence of the surrounding acidity on the internal dynamics. For that we compared mean square displacements extracted from elastic incoherent neutron scattering data, measured both at low and at neutral pH, of the phospholipids 1,2-dimyristoyl-sn-glycero-3-phosphocholine and of samples from neutralophilic and acidophilic micro-organisms (some being hyperthermophilic and others mesophilic). The lipids showed a slight shift in the phase transition temperature of about 4 degrees under pH variation and became slightly more mobile at lower pH. The membrane fragments not used to extreme acidic conditions were significantly more sensitive to variations in the pH values, whereas the acidophilic and -tolerant samples were much less influenced by this parameter. They presented the higher softness at low pH, which was closer to their native condition. Such findings might be a hint for adaptation mechanisms to different acidity conditions.

Thermodynamics of lipid multi-lamellar vesicles in presence of sterols at high hydrostatic pressure

We compared the effect of cholesterol at different concentration on the phase behaviour of DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) multilamellar vesicles. We used pressure perturbation differential scanning calorimetry (PPC) that studies a system on the whole by giving access to relevant thermodynamic quantities, and elastic incoherent neutron scattering (EINS) that probes local motions of a system at the atomic level by allowing extraction of dynamical parameters. PPC revealed that the volume expansion coefficient of DMPC and DMPC/Cholesterol samples with 13 and 25 mol% cholesterol is a linear function of the heat capacity measured by differential scanning calorimetry. Neutron backscattering spectroscopy showed that the mean square displacements of H atoms do exhibit an increase with temperature and a decrease under increasing pressure. Cholesterol added at concentrations of 25 and 50 mol% led to suppression of the main phase transition. Taking advantage of these results, the present study aims (i) to show that calorimetry and EINS using the Bicout and Zaccai model equally permit to get access to thermodynamic quantities characterizing pure DMPC and DMPC/cholesterol mixtures, thus directly confirming the theoretical method, and (ii) to validate our approach as function of temperature and of pressure, as both are equally important and complementary thermodynamic variables.

Computational Studies on Acetylcholinesterases

Functions of biomolecules, in particular enzymes, are usually modulated by structural fluctuations. This is especially the case in a gated diffusion-controlled reaction catalyzed by an enzyme such as acetylcholinesterase. The catalytic triad of acetylcholinesterase is located at the bottom of a long and narrow gorge, but it catalyzes the extremely rapid hydrolysis of the neurotransmitter, acetylcholine, with a reaction rate close to the diffusion-controlled limit. Computational modeling and simulation have produced considerable advances in exploring the dynamical and conformational properties of biomolecules, not only aiding in interpreting the experimental data, but also providing insights into the internal motions of the biomolecule at the atomic level. Given the remarkably high catalytic efficiency and the importance of acetylcholinesterase in drug development, great efforts have been made to understand the dynamics associated with its functions by use of various computational methods. Here, we present a comprehensive overview of recent computational studies on acetylcholinesterase, expanding our views of the enzyme from a microstate of a single structure to conformational ensembles, strengthening our understanding of the integration of structure, dynamics and function associated with the enzyme, and promoting the structure-based and/or mechanism-based design of new inhibitors for it.

The Dynamical Transition of Lipid Multilamellar Bilayers as a Matter of Cooperativity

The present study is the application of a two-state model formerly developed by Bicout and Zaccai [ Bicout , D. J. and Zaccai , G. Biophys. J. 2001 , 80 ( 3 ), 1115 - 1123 ] to describe the dynamical transition exhibited in the atomic mean square displacements of biological samples in terms of dynamic and thermodynamic parameters. Data were obtained by elastic incoherent neutron scattering on 1,2-dimyristoyl-sn-glycero-3-phosphocholine lipid membranes in various hydration states and on one partially per-deuterated lipid membrane. Fitting the data with the model allowed investigating which parts of lipid molecules were mainly involved in the dynamical transition, heads, tails, or both. Clear differences were found between the fully protonated and partially deuterated membranes. These findings shed light on the question of what is the degree of dynamical cooperativity of the atoms during the transition. Whereas the level of hydration does not significantly affect it, as the dry, the intermediate dry, and fully hydrated membranes all undergo a rather broad transition, the transition of the lipid tails is much sharper and sets in at much lower temperature than that of the heads. Therefore, the dynamical cooperativity appears high among the particles in the tails. Moreover, the transition of the lipid tails has to be completed first before the one of the head groups starts.

Low-density lipoproteins investigated under high hydrostatic pressure by elastic incoherent neutron scattering

Human low-density lipoprotein (LDL) is a highly complex nano-particle built up of various lipid classes and a single large protein moiety (apoB-100) owning essential physiological functions in the human body. Besides its vital role as a supplier of cholesterol and fat for peripheral tissues and cells, it is also a known key player in the formation of atherosclerosis. Due to these important roles in physiology and pathology the elucidation of structural and dynamical details is of great interest. In the current study we drew a broader picture of LDL dynamics using elastic incoherent neutron scattering (EINS) as a function of specified temperature and pressure points. We not only investigated a normolipidemic LDL sample, but also a triglyceride-rich and an oxidized one to mimic pathologic conditions as found under hyperlipidemic conditions or in atherosclerotic plaques, respectively. We could show that pressure has a significant effect on atomic motions in modified forms of LDL, whereas the normolipidemic sample seems to cope much better with high-pressure conditions irrespective of temperature. These findings might be explained by the altered lipid composition, which is either caused through elevated triglyceride content or modifications through lipid peroxidation.

Changes in dynamics of α-chymotrypsin due to covalent inhibitors investigated by elastic incoherent neutron scattering

The dynamics of chymotrypsin increases when bound to two different covalent inhibitors. These effects were analyzed by univariate and multivariate methods.