Effects of flexibility in coarse-grained models for bovine serum albumin and immunoglobulin G

We construct a coarse-grained, structure-based, low-resolution, 6-bead flexible model of bovine serum albumin (BSA, PDB: 4F5S), which is a popular example of a globular protein in biophysical research. The model is obtained via direct Boltzmann inversion using all-atom simulations of a single molecule, and its particular form is selected from a large pool of 6-bead coarse-grained models using two suitable metrics that quantify the agreement in the distribution of collective coordinates between all-atom and coarse-grained Brownian dynamics simulations of solutions in the dilute limit. For immunoglobulin G (IgG), a similar structure-based 12-bead model has been introduced in the literature [Chaudhri et al., J. Phys. Chem. B 116, 8045 (2012)] and is employed here to compare findings for the compact BSA molecule and the more anisotropic IgG molecule. We define several modified coarse-grained models of BSA and IgG, which differ in their internal constraints and thus account for a variation of flexibility. We study denser solutions of the coarse-grained models with purely repulsive molecules (achievable by suitable salt conditions) and address the effect of packing and flexibility on dynamic and static behavior. Translational and rotational self-diffusivity is enhanced for more elastic models. Finally, we discuss a number of effective sphere sizes for the BSA molecule, which can be defined from its static and dynamic properties. Here, it is found that the effective sphere diameters lie between 4.9 and 6.1 nm, corresponding to a relative spread of about ±10% around a mean of 5.5 nm.

Short-Time Transport Properties of Bidisperse Suspensions of Immunoglobulins and Serum Albumins Consistent with a Colloid Physics Picture

The crowded environment of biological systems such as the interior of living cells is occupied by macromolecules with a broad size distribution. This situation of polydispersity might influence the dependence of the diffusive dynamics of a given tracer macromolecule in a monodisperse solution on its hydrodynamic size and on the volume fraction. The resulting size dependence of diffusive transport crucially influences the function of a living cell. Here, we investigate a simplified model system consisting of two constituents in aqueous solution, namely, of the proteins bovine serum albumin (BSA) and bovine polyclonal gamma-globulin (Ig), systematically depending on the total volume fraction and ratio of these constituents. From high-resolution quasi-elastic neutron spectroscopy, the separate apparent short-time diffusion coefficients for BSA and Ig in the mixture are extracted, which show substantial deviations from the diffusion coefficients measured in monodisperse solutions at the same total volume fraction. These deviations can be modeled quantitatively using results from the short-time rotational and translational diffusion in a two-component hard sphere system with two distinct, effective hydrodynamic radii. Thus, we find that a simple colloid picture well describes short-time diffusion in binary mixtures as a function of the mixing ratio and the total volume fraction. Notably, the self-diffusion of the smaller protein BSA in the mixture is faster than the diffusion in a pure BSA solution, whereas the self-diffusion of Ig in the mixture is slower than in the pure Ig solution.

Structure and diffusive dynamics of aspartate α-decarboxylase (ADC) liganded with D-serine in aqueous solution

Incoherent neutron spectroscopy, in combination with dynamic light scattering, was used to investigate the effect of ligand binding on the center-of-mass self-diffusion and internal diffusive dynamics of Escherichia coli aspartate α-decarboxylase (ADC). The X-ray crystal structure of ADC in complex with the d-serine inhibitor was also determined, and molecular dynamics simulations were used to further probe the structural rearrangements that occur as a result of ligand binding. These experiments reveal that d-serine forms hydrogen bonds with some of the active site residues, that higher order oligomers of the ADC tetramer exist on ns–ms time-scales, and also show that ligand binding both affects the ADC internal diffusive dynamics and appears to further increase the size of the higher order oligomers.

Neutron spectroscopy, dynamic light scattering, X-ray diffraction, and MD-simulations were used to investigate the effect of ligand binding on the structure and diffusive dynamics of Escherichia coli aspartate alpha-decarboxylase.

Notes on Fitting and Analysis Frameworks for QENS Spectra of (Soft) Colloid Suspensions

With continuously improving signal-to-noise ratios, a statistically sound analysis of quasi-elastic neutron scattering (QENS) spectra requires to fit increasingly complex models which poses several challenges. Simultaneous fits of the spectra for all recorded values of the momentum transfer become a standard approach. Spectrometers at spallation sources can have a complicated non-Gaussian resolution function which has to be described most accurately. At the same time, to speed up the fitting, an analytical convolution with this resolution function is of interest. Here, we discuss basic concepts to efficient approaches for fits of QENS spectra based on standard MATLAB and Python fit algorithms. We illustrate the fits with example data from IN16B, BASIS, and BATS.

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.

Molecular Flexibility of Antibodies Preserved Even in the Dense Phase after Macroscopic Phase Separation

Antibody therapies are typically based on high-concentration formulations that need to be administered subcutaneously. These conditions induce several challenges, inter alia a viscosity suitable for injection, sufficient solution stability, and preservation of molecular function. To obtain systematic insights into the molecular factors, we study the dynamics on the molecular level under strongly varying solution conditions. In particular, we use solutions of antibodies with poly(ethylene glycol), in which simple cooling from room temperature to freezing temperatures induces a transition from a well-dispersed solution into a phase-separated and macroscopically arrested system. Using quasi-elastic neutron scattering during in situ cooling ramps and in prethermalized measurements, we observe a strong decrease in antibody diffusion, while internal flexibility persists to a significant degree, thus ensuring the movement necessary for the preservation of molecular function. These results are relevant for a more dynamic understanding of antibodies in high-concentration formulations, which affects the formation of transient clusters governing the solution viscosity.

Zinc determines dynamical properties and aggregation kinetics of human insulin

Protein aggregation is a widespread process leading to deleterious consequences in the organism, with amyloid aggregates being important not only in biology but also for drug design and biomaterial production. Insulin is a protein largely used in diabetes treatment, and its amyloid aggregation is at the basis of the so-called insulin-derived amyloidosis. Here, we uncover the major role of zinc in both insulin dynamics and aggregation kinetics at low pH, in which the formation of different amyloid superstructures (fibrils and spherulites) can be thermally induced. Amyloid aggregation is accompanied by zinc release and the suppression of water-sustained insulin dynamics, as shown by particle-induced x-ray emission and x-ray absorption spectroscopy and by neutron spectroscopy, respectively. Our study shows that zinc binding stabilizes the native form of insulin by facilitating hydration of this hydrophobic protein and suggests that introducing new binding sites for zinc can improve insulin stability and tune its aggregation propensity.

Impact of Sucrose as Osmolyte on Molecular Dynamics of Mouse Acetylcholinesterase

The enzyme model, mouse acetylcholinesterase, which exhibits its active site at the bottom of a narrow gorge, was investigated in the presence of different concentrations of sucrose to shed light on the protein and water dynamics in cholinesterases. The study was conducted by incoherent neutron scattering, giving access to molecular dynamics within the time scale of sub-nano to nanoseconds, in comparison with molecular dynamics simulations. With increasing sucrose concentration, we found non-linear effects, e.g., first a decrease in the dynamics at 5 wt% followed by a gain at 10 wt% sucrose. Direct comparisons with simulations permitted us to understand the following findings: at 5 wt%, sugar molecules interact with the protein surface through water molecules and damp the motions to reduce the overall protein mobility, although the motions inside the gorge are enhanced due to water depletion. When going to 10 wt% of sucrose, some water molecules at the protein surface are replaced by sugar molecules. By penetrating the protein surface, they disrupt some of the intra-protein contacts, and induce new ones, creating new pathways for correlated motions, and therefore, increasing the dynamics. This exhaustive study allowed for an explanation of the detail interactions leading to the observed non-linear behavior.

A neutron scattering perspective on the structure, softness and dynamics of the ligand shell of PbS nanocrystals in solution

Small-angle neutron and X-ray scattering, neutron backscattering and neutron time-of-flight spectroscopy are applied to reveal the structure of the ligand shell, the temperature-dependent diffusion properties and the phonon spectrum of PbS nanocrystals functionalized with oleic acid in deuterated hexane. The nanocrystals decorated with oleic acid as well as the desorbed ligand molecules exhibit simple Brownian diffusion with a Stokes-Einstein temperature-dependence and inhibited freezing. Ligand molecules desorbed from the surface show strong spatial confinement. The phonon spectrum of oleic acid adsorbed to the nanocrystal surface exhibits hybrid modes with a predominant Pb-character. Low-energy surface modes of the NCs are prominent and indicate a large mechanical softness in solution. This work provides comprehensive insights into the ligand-particle interaction of colloidal nanocrystals in solution and highlights its effect on the diffusion and vibrational properties as well as their mechanical softness.

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.

Osmolytes modify protein dynamics and function of tetrameric lactate dehydrogenase upon pressurization

We present a study of the combined effects of natural cosolvents (TMAO, glycine, urea) and pressure on the activity of the tetrameric enzyme lactate dehydrogenase (LDH). To this end, high-pressure stopped-flow methodology in concert with fast UV/Vis spectroscopic detection of product formation was applied. To reveal possible pressure effects on the stability and dynamics of the enzyme, FTIR spectroscopic and neutron scattering measurements were carried out. In neat buffer solution, the catalytic turnover number of the enzyme, k_cat, increases up to 1000 bar, the pressure range where dissociation of the tetrameric species to dimers sets in. Accordingly, we obtain a negative activation volume, ΔV^# = -45.3 mL mol^-1. Further, the enzyme substrate complex has a larger volume compared to the enzyme and substrate in the unbound state. The neutron scattering data show that changes in the fast internal dynamics of the enzyme are not responsible for the increase of k_cat upon compression. Whereas the magnitude of k_cat is similar in the presence of the osmolytes, the pressure of deactivation is modulated by the addition of cosolvents. TMAO and glycine increase the pressure of deactivation, and in accordance with the observed stabilizing effect both cosolvents exhibit against denaturation and/or dissociation of proteins. While urea does not markedly affect the magnitude of the Michaelis constant, K_M, both 1 M TMAO and 1 M glycine exhibit smaller K_M values of about 0.07 mM and 0.05 mM below about 1 kbar. Such positive effect on the substrate affinity could be rationalized by the effect the two cosolutes impose on the thermodynamic activities of the reactants, which reflect changes in water-mediated intermolecular interactions. Our data show that the intracellular milieu, i.e., the solution conditions that have evolved, may be sufficient to maintain enzymatic activity under extreme environmental conditions, including the whole pressure range encountered on Earth.

Protein Short-Time Diffusion in a Naturally Crowded Environment

The interior of living cells is a dense and polydisperse suspension of macromolecules. Such a complex system challenges an understanding in terms of colloidal suspensions. As a fundamental test we employ neutron spectroscopy to measure the diffusion of tracer proteins (immunoglobulins) in a cell-like environment (cell lysate) with explicit control over crowding conditions. In combination with Stokesian dynamics simulation, we address protein diffusion on nanosecond time scales where hydrodynamic interactions dominate over negligible protein collisions. We successfully link the experimental results on these complex, flexible molecules with coarse-grained simulations providing a consistent understanding by colloid theories. Both experiments and simulations show that tracers in polydisperse solutions close to the effective particle radius R_eff = ⟨ R_i³⟩^1/3 diffuse approximately as if the suspension was monodisperse. The simulations further show that macromolecules of sizes R > R_eff ( R < R_eff) are slowed more (less) effectively even at nanosecond time scales, which is highly relevant for a quantitative understanding of cellular processes.

Nanosecond Tracer Diffusion as a Probe of the Solution Structure and Molecular Mobility of Protein Assemblies: The Case of Ovalbumin

Protein diffusion is not only an important process ensuring biological function but can also be used as a probe to obtain information on structural properties of protein assemblies in liquid solutions. Here, we explore the oligomerization state of ovalbumin at high protein concentrations by means of its short-time self-diffusion. We employ high-resolution incoherent quasielastic neutron scattering to access the self-diffusion on nanosecond timescales, on which interparticle contacts are not altered. Our results indicate that ovalbumin in aqueous (D₂O) solutions occurs in increasingly large assemblies of its monomeric subunits with rising protein concentration. It changes from nearly monomeric toward dimeric and ultimately larger than tetrameric complexes. Simultaneously, we access information on the internal molecular mobility of ovalbumin on the nanometer length scale and compare it with results obtained for bovine serum albumin, immunoglobulin, and β-lactoglobulin.

Crowding-Controlled Cluster Size in Concentrated Aqueous Protein Solutions: Structure, Self- and Collective Diffusion

We investigate the concentration-controlled formation of clusters in β-lactoglobulin (BLG) protein solutions combining structural and dynamical scattering techniques. The static structure factor from small-angle X-ray scattering as well as de-Gennes narrowing in the nanosecond diffusion function D(q) from neutron spin echo spectroscopy support a picture of cluster formation. Using neutron backscattering spectroscopy, a monotonous increase of the average hydrodynamic cluster radius is monitored over a broad protein concentration range, corresponding to oligomeric structures of BLG ranging from the native dimers up to roughly four dimers. The results suggest that BLG forms compact clusters that are static on the observation time scale of several nanoseconds. The presented analysis provides a general framework to access the structure and dynamics of macromolecular assemblies in solution.

Salt-Induced Universal Slowing Down of the Short-Time Self-Diffusion of a Globular Protein in Aqueous Solution

The short-time self-diffusion D of the globular model protein bovine serum albumin in aqueous (D2O) solutions has been measured comprehensively as a function of the protein and trivalent salt (YCl3) concentration, noted cp and cs, respectively. We observe that D follows a universal master curve D(cs,cp) = D(cs = 0,cp) g(cs/cp), where D(cs = 0,cp) is the diffusion coefficient in the absence of salt and g(cs/cp) is a scalar function solely depending on the ratio of the salt and protein concentration. This observation is consistent with a universal scaling of the bonding probability in a picture of cluster formation of patchy particles. The finding corroborates the predictive power of the description of proteins as colloids with distinct attractive ion-activated surface patches.