Overlay databank unlocks data-driven analyses of biomolecules for all.

Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.

A Perspective on the Prospective Use of AI in Protein Structure Prediction.

AlphaFold2 (AF2) and RoseTTaFold (RF) have revolutionized structural biology, serving as highly reliable and effective methods for predicting protein structures. This article explores their impact and limitations, focusing on their integration into experimental pipelines and their application in diverse protein classes, including membrane proteins, intrinsically disordered proteins (IDPs), and oligomers. In experimental pipelines, AF2 models help X-ray crystallography in resolving the phase problem, while complementarity with mass spectrometry and NMR data enhances structure determination and protein flexibility prediction. Predicting the structure of membrane proteins remains challenging for both AF2 and RF due to difficulties in capturing conformational ensembles and interactions with the membrane. Improvements in incorporating membrane-specific features and predicting the structural effect of mutations are crucial. For intrinsically disordered proteins, AF2's confidence score (pLDDT) serves as a competitive disorder predictor, but integrative approaches including molecular dynamics (MD) simulations or hydrophobic cluster analyses are advocated for accurate dynamics representation. AF2 and RF show promising results for oligomeric models, outperforming traditional docking methods, with AlphaFold-Multimer showing improved performance. However, some caveats remain in particular for membrane proteins. Real-life examples demonstrate AF2's predictive capabilities in unknown protein structures, but models should be evaluated for their agreement with experimental data. Furthermore, AF2 models can be used complementarily with MD simulations. In this Perspective, we propose a "wish list" for improving deep-learning-based protein folding prediction models, including using experimental data as constraints and modifying models with binding partners or post-translational modifications. Additionally, a meta-tool for ranking and suggesting composite models is suggested, driving future advancements in this rapidly evolving field.

Incidence and Predictors of Ventricular Arrhythmias in Transthyretin Amyloid Cardiomyopathy.

BACKGROUND: Wild-type transthyretin amyloid cardiomyopathy (wtATTR-CM) is characterized by heart failure, conduction abnormalities and arrhythmias. The incidence of ventricular arrhythmias, particularly ventricular tachycardias (VTs), in wtATTR-CM is unclear. With the development of targeted therapies and improved overall prognosis, there is an unmet need to identify patients at high risk for VTs who might benefit from ICD therapy. METHODS: Between 2017 and 2022, 72 patients diagnosed with wtATTR-CM were prospectively evaluated for the presence of ventricular arrhythmias using a Holter ECG. VTs were defined as >3 consecutive beats with a heart rate > 100 beats per minute originating from a ventricle. RESULTS: The incidence of VTs was 44% (n = 32/72) in unselected wtATTR-CM patients. Patients with VT showed significantly more severe left ventricular (LV) hypertrophy (septum diameter 21 ± 2.6 vs. 19 ± 3.0 mm, p = 0.006), reduced LV ejection fraction (47 ± 8 vs. 52 ± 8%, p = 0.014) and larger left atria (32 ± 7 vs. 28 ± 6 mm2, p = 0.020), but no differences in cardiac markers such as NTproBNP and troponin. In a multivariable model, LV hypertrophy (LV mass indexed, OR = 1.02 [1.00-1.03], p = 0.031), LV end-diastolic diameter (OR = 0.85 [0.74-0.98], p = 0.021) and LV end-systolic diameter (OR = 1.19 [1.03-1.349], p = 0.092) were predictive for VT occurrence with an area under the receiver operating characteristic of 0.76 [0.65-0.87]. CONCLUSIONS: The incidence of ventricular arrhythmia in wtATTR-CM is high and is associated with an advanced stage of left ventricular disease. Further studies are needed evaluating the role of VTs in predicting sudden cardiac death and the benefit of ICD therapy in wtATTR-CM.

Incomplete spectrum QSM using support information.

Introduction: Reconstructing a bounded object from incomplete k-space data is a well posed problem, and it was recently shown that this incomplete spectrum approach can be used to reconstruct undersampled MRI images with similar quality to compressed sensing approaches. Here, we apply this incomplete spectrum approach to the field-to-source inverse problem encountered in quantitative magnetic susceptibility mapping (QSM). The field-to-source problem is an ill-posed problem because of conical regions in frequency space where the dipole kernel is zero or very small, which leads to the kernel's inverse being ill-defined. These "ill-posed" regions typically lead to streaking artifacts in QSM reconstructions. In contrast to compressed sensing, our approach relies on knowledge of the image-space support, more commonly referred to as the mask, of our object as well as the region in k-space with ill-defined values. In the QSM case, this mask is usually available, as it is required for most QSM background field removal and reconstruction methods. Methods: We tuned the incomplete spectrum method (mask and band-limit) for QSM on a simulated dataset from the most recent QSM challenge and validated the QSM reconstruction results on brain images acquired in five healthy volunteers, comparing incomplete spectrum QSM to current state-of-the art-methods: FANSI, nonlinear dipole inversion, and conventional thresholded k-space division. Results: Without additional regularization, incomplete spectrum QSM performs slightly better than direct QSM reconstruction methods such as thresholded k-space division (PSNR of 39.9 vs. 39.4 of TKD on a simulated dataset) and provides susceptibility values in key iron-rich regions similar or slightly lower than state-of-the-art algorithms, but did not improve the PSNR in comparison to FANSI or nonlinear dipole inversion. With added (ℓ1-wavelet based) regularization the new approach produces results similar to compressed sensing based reconstructions (at sufficiently high levels of regularization). Discussion: Incomplete spectrum QSM provides a new approach to handle the "ill-posed" regions in the frequency-space data input to QSM.

Triglyceride lipolysis triggers liquid crystalline phases in lipid droplets and alters the LD proteome.

Lipid droplets (LDs) are reservoirs for triglycerides (TGs) and sterol-esters (SEs), but how these lipids are organized within LDs and influence their proteome remain unclear. Using in situ cryo-electron tomography, we show that glucose restriction triggers lipid phase transitions within LDs generating liquid crystalline lattices inside them. Mechanistically this requires TG lipolysis, which decreases the LD's TG:SE ratio, promoting SE transition to a liquid crystalline phase. Molecular dynamics simulations reveal TG depletion promotes spontaneous TG and SE demixing in LDs, additionally altering the lipid packing of the PL monolayer surface. Fluorescence imaging and proteomics further reveal that liquid crystalline phases are associated with selective remodeling of the LD proteome. Some canonical LD proteins, including Erg6, relocalize to the ER network, whereas others remain LD-associated. Model peptide LiveDrop also redistributes from LDs to the ER, suggesting liquid crystalline phases influence ER-LD interorganelle transport. Our data suggests glucose restriction drives TG mobilization, which alters the phase properties of LD lipids and selectively remodels the LD proteome.

Anionic lipids induce a fold-unfold transition in the membrane-translocating Engrailed homeodomain.

Homeoprotein transcription factors have the property of interacting with membranes through their DNA-binding homeodomain, which is involved in unconventional internalization and secretion. Both processes depend on membrane-translocating events but their detailed molecular mechanisms are still poorly understood. We have previously characterized the conformational properties of Engrailed 2 homeodomain (EnHD) in aqueous solution and in micelles as membrane-mimetic environments. In the present study, we used small isotropic lipid bicelles as a more relevant membrane-mimetic model to characterize the membrane-bound state of EnHD. We show that lipid bicelles, in contrast to micelles, adequately reproduce the requirement of anionic lipids in the membrane binding and conformational transition of EnHD. The fold-unfold transition of EnHD induced by anionic lipids was characterized by NMR using 1H, 13C, 15N chemical shifts, nuclear Overhauser effects, residual dipolar couplings, intramolecular and intermolecular paramagnetic relaxation enhancements induced by site-directed spin-label or paramagnetic lipid probe, respectively. A global unpacking of EnHD helices is observed leading to a loss of the native fold. However, near-native propensities of EnHD backbone conformation are maintained in membrane environment, including not only the three helices but also the turn connecting helices H2 and H3. NMR and coarse-grained molecular dynamics simulations reveal that the EnHD adopts a shallow insertion in the membrane, with the three helices oriented parallel to the membrane. EnHD explores extended conformations and closed U-shaped conformations, which are stabilized by anionic lipid recruitment.

Structural Bases for the Involvement of Phosphatidylinositol-4,5-bisphosphate in the Internalization of the Cell-Penetrating Peptide Penetratin.

Cell-penetrating peptides cross cell membranes through various parallel internalization pathways. Herein, we analyze the role of the negatively charged lipid phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2) in the internalization of Penetratin. Contributions of both inner leaflet and outer leaflet pools of PI(4,5)P2 were revealed by quantifying the internalization of Penetratin in cells treated with PI(4,5)P2 binders. Studies on model systems showed that Penetratin has a strong affinity for PI(4,5)P2 and interacts selectively with this lipid, even in the presence of other negatively charged lipids, as demonstrated by affinity photo-crosslinking experiments. Differential scanning calorimetry experiments showed that Penetratin induces lateral segregation in PI(4,5)P2-containing liposomes, which was confirmed by coarse-grained molecular dynamics simulations. NMR experiments indicated that Penetratin adopts a stabilized helical conformation in the presence of PI(4,5)P2-containing membranes, with an orientation parallel to the bilayer plane, which was also confirmed by all-atom simulations. NMR and photo-crosslinking experiments also suggest a rather shallow insertion of the peptide in the membrane. Put together, our findings suggest that PI(4,5)P2 is a privileged interaction partner for Penetratin and that it plays an important role in Penetratin internalization.

Specific Game-Based Performance in Elite Male Adolescent Team Handball Players.

PURPOSE: In elite team handball, talent identification and selection of the best young players is a fundamental process in several national federations and clubs; however, literature addressing the specific game-based performance in team handball is almost nonexistent. Consequently, the aim of the study was to assess and compare the team-handball-specific game-based performance of elite male team handball players of different ages. METHODS: Twelve under-23, 10 under-19, 10 under-17, and 10 under-15 elite male players performed the team-handball game-based performance test. During testing, oxygen uptake, heart rate, sprinting time in defense, offense, fast breaks, and fast retreats, as well as ball velocity and jump height in the jump shot, were measured. RESULTS: Significant differences (P < .05) between under-23, under-19, under-17, and under-15 players were found for absolute peak oxygen uptake, defense, offense and fast break time, ball velocity, and jump height in the game-based performance test, as well as in body weight and height. CONCLUSION: The results revealed that with increasing age, elite male team handball players are heavier and taller (body weight and height); faster (team-handball offense, defense, and fast break); jump higher and throw faster (in the team-handball jump shot); and perform better aerobically (absolute peak oxygen uptake). The better performance in the under-23 and under-19 players compared with male adult players competing in a lower National Federation league (not on top-elite level) demonstrates that highly specific game-based physical performance determines the potential for developing young male team handball players for competition at the top level.

How to best estimate the viscosity of lipid bilayers.

The viscosity of lipid bilayers is a property relevant to biological function, as it affects the diffusion of membrane macromolecules. To determine its value, and hence portray the membrane, various literature-reported techniques lead to significantly different results. Herein we compare the results issuing from two widely used techniques to determine the viscosity of membranes: the Fluorescence Lifetime Imaging Microscopy (FLIM), and Fluorescence Recovery After Photobleaching (FRAP). FLIM relates the time of rotation of a molecular rotor inserted into the membrane to the macroscopic viscosity of a fluid. Whereas FRAP measures molecular diffusion coefficients. This approach is based on a hydrodynamic model connecting the mobility of a membrane inclusion to the viscosity of the membrane. We show that: This article emphasizes the pitfalls to be avoided and the rules to be observed in order to obtain a value of the bilayer viscosity that characterizes the bilayer instead of interactions between the bilayer and the embedded probe.

Cooperation of Conical and Polyunsaturated Lipids to Regulate Initiation and Processing of Membrane Fusion.

The shape of lipids has long been suspected to be a critical determinant for the control of membrane fusion. To experimentally test this assertion, we used conical and malleable lipids and measured their influence on the fusion kinetics. We found that, as previously suspected, both types of lipids accelerate fusion. However, the implicated molecular mechanisms are strikingly different. Malleable lipids, with their ability to change shape with low energy cost, favor fusion by decreasing the overall activation energy. On the other hand, conical lipids, with their small polar head relative to the area occupied by the hydrophobic chains, tend to make fusion less energetically advantageous because they tend to migrate towards the most favorable lipid leaflet, hindering fusion pore opening. They could however facilitate fusion by generating hydrophobic defects on the membranes; this is suggested by the similar trend observed between the experimental rate of fusion nucleation and the surface occupied by hydrophobic defects obtained by molecular simulations. The synergy of dual-process, activation energy and nucleation kinetics, could facilitate membrane fusion regulation in vivo.

Inverse Conformational Selection in Lipid-Protein Binding.

Interest in lipid interactions with proteins and other biomolecules is emerging not only in fundamental biochemistry but also in the field of nanobiotechnology where lipids are commonly used, for example, in carriers of mRNA vaccines. The outward-facing components of cellular membranes and lipid nanoparticles, the lipid headgroups, regulate membrane interactions with approaching substances, such as proteins, drugs, RNA, or viruses. Because lipid headgroup conformational ensembles have not been experimentally determined in physiologically relevant conditions, an essential question about their interactions with other biomolecules remains unanswered: Do headgroups exchange between a few rigid structures, or fluctuate freely across a practically continuous spectrum of conformations? Here, we combine solid-state NMR experiments and molecular dynamics simulations from the NMRlipids Project to resolve the conformational ensembles of headgroups of four key lipid types in various biologically relevant conditions. We find that lipid headgroups sample a wide range of overlapping conformations in both neutral and charged cellular membranes, and that differences in the headgroup chemistry manifest only in probability distributions of conformations. Furthermore, the analysis of 894 protein-bound lipid structures from the Protein Data Bank suggests that lipids can bind to proteins in a wide range of conformations, which are not limited by the headgroup chemistry. We propose that lipids can select a suitable headgroup conformation from the wide range available to them to fit the various binding sites in proteins. The proposed inverse conformational selection model will extend also to lipid binding to targets other than proteins, such as drugs, RNA, and viruses.

Molecular Study of Ultrasound-Triggered Release of Fluorescein from Liposomes.

Several investigations have suggested that ultrasound triggers the release of drugs encapsulated into liposomes at acoustic pressures low enough to avoid cavitation or high hyperthermia. However, the mechanism leading to this triggered release as well as the adequate composition of the liposome membrane remains unknown. Here, we investigate the ultrasound-triggered release of fluorescein disodium salt encapsulated into liposomes made of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) or 1,2-distearoylphosphatidyl-ethanolamine (DSPC) lipids with various concentrations of cholesterol (from 0 to 44 mol %). The passive release of encapsulated fluorescein was first characterized. It was observed to be higher when the membrane is in a fluid phase and increased with temperature but decreased upon addition of cholesterol. Next, the release of fluorescein was measured at different acoustic frequencies (0.8, 1.1, and 3.3 MHz) and peak-to-peak pressures (0, 2, 2.5, 5, and 8 MPa). Measurements were performed at temperatures where DOPC and DSPC liposomes were, respectively, in the fluid or gel phase. We found that the release rate did not depend on the ultrasound frequency. For DOPC liposomes, the ultrasound-triggered release of fluorescein decreased with increasing concentration of cholesterol in liposomes, while the behavior was more complex for DSPC liposomes. Overall, the triggered release from DSPC liposomes was up to ten times less than DOPC liposomes. Molecular dynamics simulations performed on a pure DOPC membrane showed that a membrane experiences, under a directional pressure of ±2.4 MPa, various changes in properties such as the area per lipid (APL). An increase in the APL was notably observed when the simulation box was laterally stretched or perpendicularly compressed, which was accompanied by an increase in the number of water molecules crossing the membrane. This suggests that ultrasound most probably enhances the diffusion of encapsulated molecules at small acoustic pressures by increasing the distance between lipids.

On the Effect of the Various Assumptions and Approximations used in Molecular Simulations on the Properties of Bio-Molecular Systems: Overview and Perspective on Issues.

Computer simulations of molecular systems enable structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic or supra-molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio-molecular systems are considered. What the developers of bio-molecular force fields and simulation software can do to facilitate and broaden research involving bio-molecular simulations is also discussed.

Transceive phase corrected 2D contrast source inversion-electrical properties tomography.

PURPOSE: To remove the necessity of the tranceive phase assumption for CSI-EPT and show electrical properties maps reconstructed from measured data obtained using a standard 3T birdcage body coil setup. METHODS: The existing CSI-EPT algorithm is reformulated to use the transceive phase rather than relying on the transceive phase assumption. Furthermore, the radio frequency (RF)-shield is numerically implemented to accurately model the RF fields inside the MRI scanner. We verify that the reformulated two-dimensional (2D) CSI-EPT algorithm can reconstruct electrical properties maps given 2D electromagnetic simulations. Afterward, the algorithm is tested with three-dimensional (3D) FDTD simulations to investigate if the 2D CSI-EPT can retrieve the electrical properties for 3D RF fields. Finally, an MR experiment at 3T with a phantom is performed. RESULTS: From the results of the 2D simulations, it is seen that CSI-EPT can reconstruct the electrical properties using MRI accessible quantities. For 3D simulations, it is observed that the electrical properties are underestimated, nonetheless, CSI-EPT has a lower standard deviation than the standard Helmholtz based methods. Finally, the first CSI-EPT reconstructions based on measured data are presented showing comparable accuracy and precision to reconstructions based on simulated data, and demonstrating the feasibility of CSI-EPT. CONCLUSIONS: The CSI-EPT algorithm was rewritten to use MRI accessible quantities. This allows for CSI-EPT to fully exploit the benefits of the higher static magnetic field strengths with a standard quadrature birdcage coil setup.

CX3CL1 homo-oligomerization drives cell-to-cell adherence.

During inflammatory response, blood leukocytes adhere to the endothelium. This process involves numerous adhesion molecules, including a transmembrane chemokine, CX3CL1, which behaves as a molecular cluster. How this cluster assembles and whether this association has a functional role remain unknown. The analysis of CX3CL1 clusters using native electrophoresis and single molecule fluorescence kinetics shows that CX3CL1 is a homo-oligomer of 3 to 7 monomers. Fluorescence recovery after photobleaching assays reveal that the CX3CL1-transmembrane domain peptide self-associates in both cellular and acellular lipid environments, while its random counterpart (i.e. peptide with the same residues in a different order) does not. This strongly indicates that CX3CL1 oligomerization is driven by its intrinsic properties. According to the molecular modeling, CX3CL1 does not associate in compact bundles but rather with monomers linearly assembled side by side. Finally, the CX3CL1 transmembrane peptide inhibits both the CX3CL1 oligomerization and the adhesive function, while its random counterpart does not. This demonstrates that CX3CL1 oligomerization is mandatory for its adhesive potency. Our results provide a new direction to control CX3CL1-dependent cellular adherence in key immune processes.