Beyond Poisson-Boltzmann: Numerical sampling of charge density fluctuations

We present a method aimed at sampling charge density fluctuations in Coulomb systems. The derivation follows from a functional integral representation of the partition function in terms of charge density fluctuations. Starting from the mean-field solution given by the Poisson-Boltzmann equation, an original approach is proposed to numerically sample fluctuations around it, through the propagation of a Langevin like stochastic partial differential equation (SPDE). The diffusion tensor of the SPDE can be chosen so as to avoid the numerical complexity linked to long-range Coulomb interactions, effectively rendering the theory completely local. A finite-volume implementation of the SPDE is described, and the approach is illustrated with preliminary results on the study of a system made of two like-charge ions immersed in a bath of counter-ions

Intermolecular correlations are necessary to explain diffuse scattering from protein crystals

Conformational changes drive protein function, including catalysis, allostery, and signaling. X-ray diffuse scattering from protein crystals has frequently been cited as a probe of these correlated motions, with significant potential to advance our understanding of biological dynamics. However, recent work challenged this prevailing view, suggesting instead that diffuse scattering primarily originates from rigid body motions and could therefore be applied to improve structure determination. To investigate the nature of the disorder giving rise to diffuse scattering, and thus the potential applications of this signal, a diverse repertoire of disorder models was assessed for its ability to reproduce the diffuse signal reconstructed from three protein crystals. This comparison revealed that multiple models of intramolecular conformational dynamics, including ensemble models inferred from the Bragg data, could not explain the signal. Models of rigid body or short-range liquid-like motions, in which dynamics are confined to the biological unit, showed modest agreement with the diffuse maps, but were unable to reproduce experimental features indicative of long-range correlations. Extending a model of liquid-like motions to include disorder across neighboring proteins in the crystal significantly improved agreement with all three systems and highlighted the contribution of intermolecular correlations to the observed signal. These findings anticipate a need to account for intermolecular disorder in order to advance the interpretation of diffuse scattering to either extract biological motions or aid structural inference.Comment: 12 pages, 5 figures (not including Supplementary Information

Identifying Interpretable Visual Features in Artificial and Biological Neural Systems

Single neurons in neural networks are often interpretable in that they represent individual, intuitively meaningful features. However, many neurons exhibit $\textit{mixed selectivity}$, i.e., they represent multiple unrelated features. A recent hypothesis proposes that features in deep networks may be represented in $\textit{superposition}$, i.e., on non-orthogonal axes by multiple neurons, since the number of possible interpretable features in natural data is generally larger than the number of neurons in a given network. Accordingly, we should be able to find meaningful directions in activation space that are not aligned with individual neurons. Here, we propose (1) an automated method for quantifying visual interpretability that is validated against a large database of human psychophysics judgments of neuron interpretability, and (2) an approach for finding meaningful directions in network activation space. We leverage these methods to discover directions in convolutional neural networks that are more intuitively meaningful than individual neurons, as we confirm and investigate in a series of analyses. Moreover, we apply the same method to three recent datasets of visual neural responses in the brain and find that our conclusions largely transfer to real neural data, suggesting that superposition might be deployed by the brain. This also provides a link with disentanglement and raises fundamental questions about robust, efficient and factorized representations in both artificial and biological neural systems

Heterogeneous reconstruction of deformable atomic models in Cryo-EM

Cryogenic electron microscopy (cryo-EM) provides a unique opportunity to study the structural heterogeneity of biomolecules. Being able to explain this heterogeneity with atomic models would help our understanding of their functional mechanisms but the size and ruggedness of the structural space (the space of atomic 3D cartesian coordinates) presents an immense challenge. Here, we describe a heterogeneous reconstruction method based on an atomistic representation whose deformation is reduced to a handful of collective motions through normal mode analysis. Our implementation uses an autoencoder. The encoder jointly estimates the amplitude of motion along the normal modes and the 2D shift between the center of the image and the center of the molecule . The physics-based decoder aggregates a representation of the heterogeneity readily interpretable at the atomic level. We illustrate our method on 3 synthetic datasets corresponding to different distributions along a simulated trajectory of adenylate kinase transitioning from its open to its closed structures. We show for each distribution that our approach is able to recapitulate the intermediate atomic models with atomic-level accuracy.Comment: 8 pages, 1 figur

Reproducibility of protein x-ray diffuse scattering and potential utility for modeling atomic displacement parameters

Protein structure and dynamics can be probed using x-ray crystallography. Whereas the Bragg peaks are only sensitive to the average unit-cell electron density, the signal between the Bragg peaks—diffuse scattering—is sensitive to spatial correlations in electron-density variations. Although diffuse scattering contains valuable information about protein dynamics, the diffuse signal is more difficult to isolate from the background compared to the Bragg signal, and the reproducibility of diffuse signal is not yet well understood. We present a systematic study of the reproducibility of diffuse scattering from isocyanide hydratase in three different protein forms. Both replicate diffuse datasets and datasets obtained from different mutants were similar in pairwise comparisons (Pearson correlation coefficient \u3e/= 0.8). The data were processed in a manner inspired by previously published methods using custom software with modular design, enabling us to perform an analysis of various data processing choices to determine how to obtain the highest quality data as assessed using unbiased measures of symmetry and reproducibility. The diffuse data were then used to characterize atomic mobility using a liquid-like motions (LLM) model. This characterization was able to discriminate between distinct anisotropic atomic displacement parameter (ADP) models arising from different anisotropic scaling choices that agreed comparably with the Bragg data. Our results emphasize the importance of data reproducibility as a model-free measure of diffuse data quality, illustrate the ability of LLM analysis of diffuse scattering to select among alternative ADP models, and offer insights into the design of successful diffuse scattering experiments

La cuisine d’une <i>domus</i> de Javols/<i>Anderitum</i> (Lozère), chef-lieu de cité des Gabales : l’alimentation d’une demeure aisée dans la seconde moitié du II^e s. apr. J.-C.

International audienceA domus has been almost exhaustively excavated within a research program on the site of Javols/Anderitum, the capital of the Gabali in the Roman period. It is equipped, probably around the end of the 1st century AD, with a kitchen which develops until the early 3rd century, with important devices (3 ovens), before being demolished. The excavation, the study of these structures and analysis of the pottery and bioarchaeological remains (anthracology, carpology, archaeozoology, ichtyology) allow an approach of its functioning and part of the diet of a wealthy domus of this Roman town in a middle mountain region.Dans le cadre du programme de recherche concernant le site de Javols/Anderitum, capitale de cité des Gabales à l’époque romaine, une domus à été quasi exhaustivement fouillée. Elle est dotée, sans doute vers la fin du Ier s. apr. J. C., d’une pièce de cuisine qui évolue, avant d’être détruite, jusqu’au début du IIIe s., avec un équipement assez important (3 fours). Sa fouille, l’étude de ses structures puis celle des restes céramiques et bioarchéologiques (anthracologie, carpologie, archéozoologie et ichtyologie) permettent d’approcher son fonctionnement, et une part de l’alimentation d’une domus aisée de cette ville romaine de moyenne montagne

AquaSAXS: a web server for computation and fitting of SAXS profiles with non-uniformally hydrated atomic models

Small Angle X-ray Scattering (SAXS) techniques are becoming more and more useful for structural biologists and biochemists, thanks to better access to dedicated synchrotron beamlines, better detectors and the relative easiness of sample preparation. The ability to compute the theoretical SAXS profile of a given structural model, and to compare this profile with the measured scattering intensity, yields crucial structural informations about the macromolecule under study and/or its complexes in solution. An important contribution to the profile, besides the macromolecule itself and its solvent-excluded volume, is the excess density due to the hydration layer. AquaSAXS takes advantage of recently developed methods, such as AquaSol, that give the equilibrium solvent density map around macromolecules, to compute an accurate SAXS/WAXS profile of a given structure and to compare it to the experimental one. Here, we describe the interface architecture and capabilities of the AquaSAXS web server (http://lorentz.dynstr.pasteur.fr/aquasaxs.php)

Changes in an Enzyme Ensemble During Catalysis Observed by High Resolution XFEL Crystallography

Enzymes populate ensembles of structures with intrinsically different catalytic proficiencies that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography (MISC) at an X-ray free electron laser (XFEL) to observe catalysis in a designed mutant (G150T) isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations and formation of the thioimidate catalytic intermediate selects for catalytically competent substates. A prior proposal for active site cysteine charge-coupled conformational changes in ICH is validated by determining structures of the enzyme over a range of pH values. A combination of large molecular dynamics simulations of the enzyme in crystallo and timeresolved electron density maps shows that ionization of the general acid Asp17 during catalysis causes additional conformational changes that propagate across the dimer interface, connecting the two active sites. These ionization-linked changes in the ICH conformational ensemble permit water to enter the active site in a location that is poised for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics

Reduction of small-angle scattering profiles to finite sets of structural invariants

International audienceThis paper shows how small-angle scattering (SAS) curves can be decomposed in a simple sum using a set of invariant parameters called Kn which are related to the shape of the object of study. These Kn , together with a radius R, give a complete theoretical description of the SAS curve. Adding an overall constant, these parameters are easily fitted against experimental data giving a concise comprehensive description of the data. The pair distance distribution function is also entirely described by this invariant set and the D max parameter can be measured. In addition to the understanding they bring, these invariants can be used to reliably estimate structural moments beyond the radius of gyration, thereby rigorously expanding the actual set of model-free quantities one can extract from experimental SAS data, and possibly paving the way to designing new shape reconstruction strategies