The coexistence region in the Van der Waals fluid and the liquid-liquid phase transitions

Cellular membraneless organelles are thought to be droplets formed within the two-phase region corresponding to proteinaceous systems endowed with the liquid-liquid transition. However, their metastability requires an additional constraint—they arise in a certain region of density and temperature between the spinodal and binodal lines. Here, we consider the well-studied van der Waals fluid as a test model to work out criteria to determine the location of the spinodal line for situations in which the equation of state is not known. Our molecular dynamics studies indicate that this task can be accomplished by considering the specific heat, the surface tension and characteristics of the molecular clusters, such as the number of component chains and radius of gyration

The G\=oMartini approach: Revisiting the concept of contact maps and the modelling of protein complexes

We present a review of a series of contact maps for the determination of native interactions in proteins and nucleic acids based on a distance-threshold. Such contact maps are mostly based on physical and chemical construction, and yet they are sensitive to some parameters (e.g. distances or atomic radii) and can neglect some key interactions. Furthermore, we also comment on a new class of contact maps that only requires geometric arguments. The contact map is a necessary ingredient to build a robust G\=oMartini model for proteins and their complexes in the Martini 3 force field. We present the extension of a popular structure-based G\=o-like approach for the study of protein-sugar complexes, and also limitations of this approach are discussed. The G\=oMartini approach was first introduced by Poma et al. J. Chem. Theory Comput. 2017, 13(3), 1366-1374 in Martini 2 force field and recently, it has gained the status of gold-standard for protein simulation undergoing conformational changes in Martini 3 force field. We discuss several studies that have provided support to this approach in the context of the biophysical community.Comment: 19 pages, 3 figure

Steered molecular dynamics simulations reveal the role of Ca2+ in regulating mechanostability of cellulose-binding proteins

The conversion of cellulosic biomass into biofuels requires degradation of the biomass into fermentable sugars. The most efficient natural cellulase system for carrying out this conversion is an extracellular multi-enzymatic complex named the cellulosome. In addition to temperature and pH stability, mechanical stability is important for functioning of cellulosome domains, and experimental techniques such as Single Molecule Force Spectroscopy (SMFS) have been used to measure the mechanical strength of several cellulosomal proteins. Molecular dynamics computer simulations provide complementary atomic-resolution quantitative maps of domain mechanical stability for identification of experimental leads for protein stabilization. In this study, we used multi-scale steered molecular dynamics computer simulations, benchmarked against new SMFS measurements, to measure the intermolecular contacts that confer high mechanical stability to a family 3 Carbohydrate Binding Module protein (CBM3) derived from the archetypal Clostridium thermocellum cellulosome. Our data predicts that electrostatic interactions in the calcium binding pocket modulate the mechanostability of the cellulose-binding module, which provides an additional design rule for the rational re-engineering of designer cellulosomes for biotechnology. Our data offers new molecular insights into the origins of mechanostability in cellulose binding domains and gives leads for synthesis of more robust cellulose-binding protein modules. On the other hand, simulations predict that insertion of a flexible strand can promote alternative unfolding pathways and dramatically reduce the mechanostability of the carbohydrate binding module, which gives routes to rational design of tailormade fingerprint complexes for force spectroscopy experiments

Polysaccharide–Protein Complexes in a Coarse-Grained Model

We construct two variants of coarse-grained models of three hexaoses: one based on the centers of mass of the monomers and the other associated with the C4 atoms. The latter is found to be better defined and more suitable for studying interactions with proteins described within α-C based models. We determine the corresponding effective stiffness constants through all-atom simulations and two statistical methods. One method is the Boltzmann inversion (BI) and the other, named energy-based (EB), involves direct monitoring of energies as a function of the variables that define the stiffness potentials. The two methods are generally consistent in their account of the stiffness. We find that the elastic constants differ between the hexaoses and are noticeably different from those determined for the crystalline cellulose Iβ. The nonbonded couplings through hydrogen bonds between different sugar molecules are modeled by the Lennard-Jones potentials and are found to be stronger than the hydrogen bonds in proteins. We observe that the EB method agrees with other theoretical and experimental determinations of the nonbonded parameters much better than BI. We then consider the hexaose-Man5B catalytic complexes and determine the contact energies between their the C4−α-C atoms. These interactions are found to be stronger than the proteinic hydrogen bonds: about four times as strong for cellohexaose and two times for mannohexaose. The fluctuational dynamics of the coarse-grained complexes are found to be compatible with previous all-atom studies by Bernardi et al