Search

Complex coacervates are liquid–liquid phase separated systems, typically containing oppositely charged polyelectrolytes. They are widely studied for their functional properties as well as their potential involvement in cellular compartmentalization as biomolecular condensates. Diffusion and partitioning of solutes into a coacervate phase are important to address because their highly dynamic nature is one of their most important functional characteristics in real-world systems, but are difficult to study experimentally or even theoretically without an explicit representation of every molecule in the system. Here, we present an explicit-solvent, molecular dynamics coarse-grain model of complex coacervates, based on the Martini 3.0 force field. We demonstrate the accuracy of the model by reproducing the salt dependent coacervation of poly-lysine and poly-glutamate systems, and show the potential of the model by simulating the partitioning of ions and small nucleotides between the condensate and surrounding solvent phase. Our model paves the way for simulating coacervates and biomolecular condensates in a wide range of conditions, with near-atomic resolution., Martini 3 force field can capture the experimental trends of complex coacervates and can be extended to gain physical insight on the mechanisms that drive the formation of LLPS.

Unraveling how chemistry can give rise to biology is one of the greatest challenges of contemporary science. Achieving life-like properties in chemical systems is therefore a popular topic of research. Synthetic chemical systems are usually deterministic: the outcome is determined by the experimental conditions. In contrast, many phenomena that occur in nature are not deterministic but caused by random fluctuations (stochastic). Here, we report on how, from a mixture of two synthetic molecules, two different self-replicators emerge in a stochastic fashion. Under the same experimental conditions, the two self-replicators are formed in various ratios over several repeats of the experiment. We show that this variation is caused by a stochastic nucleation process and that this stochasticity is more pronounced close to a phase boundary. While stochastic nucleation processes are common in crystal growth and chiral symmetry breaking, it is unprecedented for systems of synthetic self-replicators.

Self-assembly features prominently in fields ranging from materials science to biophysical chemistry. Assembly pathways, often passing through transient intermediates, can control the outcome of assembly processes. Yet, the mechanisms of self-assembly remain largely obscure due to a lack of experimental tools for probing these pathways at the molecular level. Here, the self-assembly of self-replicators into fibers is visualized in real-time by high-speed atomic force microscopy (HS-AFM). Fiber growth requires the conversion of precursor molecules into six-membered macrocycles, which constitute the fibers. HS-AFM experiments, supported by molecular dynamics simulations, revealed that aggregates of precursor molecules accumulate at the sides of the fibers, which then diffuse to the fiber ends where growth takes place. This mechanism of precursor reservoir formation, followed by one-dimensional diffusion, which guides the precursor molecules to the sites of growth, reduces the entropic penalty associated with colocalizing precursors and growth sites and constitutes a new mechanism for supramolecular polymerization.

The possibility of using differential pre-heating prior to supramolecular gelation to control the balance between hydrogen-bonding and aromatic stacking interactions in supramolecular gels and obtain consequent systematic regulation of structure and properties is demonstrated. Using a model aromatic peptide amphiphile, Fmoc-tyrosyl-leucine (Fmoc-YL) and a combination of fluorescence, infrared, circular dichroism and NMR spectroscopy, it is shown that the balance of these interactions can be adjusted by temporary exposure to elevated temperatures in the range 313-365 K, followed by supramolecular locking in the gel state by cooling to room temperature. Distinct regimes can be identified regarding the balance between H-bonding and aromatic stacking interactions, with a transition point at 333 K. Consequently, gels can be obtained with customizable properties, including supramolecular chirality and gel stiffness. The differential supramolecular structures also result in changes in proteolytic stability, highlighting the possibility of obtaining a range of supramolecular architectures from a single molecular structure by simply controlling the pre-assembly temperature.

The self-assembly of bio-inspired supramolecular polymers can be unravelled using molecular dynamics simulations combined with experiments., In bionanotechnology, the field of creating functional materials consisting of bio-inspired molecules, the function and shape of a nanostructure only appear through the assembly of many small molecules together. The large number of building blocks required to define a nanostructure combined with the many degrees of freedom in packing small molecules has long precluded molecular simulations, but recent advances in computational hardware as well as software have made classical simulations available to this strongly expanding field. Here, we review the state of the art in simulations of self-assembling bio-inspired supramolecular systems. We will first discuss progress in force fields, simulation protocols and enhanced sampling techniques using recent examples. Secondly, we will focus on efforts to enable the comparison of experimentally accessible observables and computational results. Experimental quantities that can be measured by microscopy, spectroscopy and scattering can be linked to simulation output either directly or indirectly, via quantum mechanical or semi-empirical techniques. Overall, we aim to provide an overview of the various computational approaches to understand not only the molecular architecture of nanostructures, but also the mechanism of their formation.

Self-replication at the molecular level is often seen as essential to the early origins of life. Recently a mechanism of self-replication has been discovered in which replicator self-assembly drives the process. We have studied one of the examples of such self-assembling self-replicating molecules to a high level of structural detail using a combination of computational and spectroscopic techniques. Molecular Dynamics simulations of self-assembled stacks of peptide-derived replicators provide insights into the structural characteristics of the system and serve as the basis for semiempirical calculations of the UV-vis, circular dichroism (CD) and infrared (IR) absorption spectra that reflect the chiral organization and peptide secondary structure of the stacks. Two proposed structural models are tested by comparing calculated spectra to experimental data from electron microscopy, CD and IR spectroscopy, resulting in a better insight into the specific supramolecular interactions that lead to self-replication. Specifically, we find a cooperative self-assembly process in which β-sheet formation leads to well-organized structures, while also the aromatic core of the macrocycles plays an important role in the stability of the resulting fibers.

Designing molecular disorder Melanins are a group of natural pigments that are the primary factor affecting skin color. Lampel et al. examined a family of melanin-inspired materials based on tripeptides containing tyrosine as precursors for polymeric pigments. They found that the supramolecular organization of the tripeptide assembly is the most important factor for the enzymatic oxidation, with the position of the tyrosine residue playing a dominant role. Thus, simply juggling the order of the peptides allowed tuning of the optical and electrical properties of the resulting polymers. Science , this issue p. 1064

Plastoquinone (PLQ) acts as an electron carrier between photosystem II (PSII) and the cytochrome b6f complex. To understand how PLQ enters and leaves PSII, here we show results of coarse grained molecular dynamics simulations of PSII embedded in the thylakoid membrane, covering a total simulation time of more than 0.5 ms. The long time scale allows the observation of many spontaneous entries of PLQ into PSII, and the unbinding of plastoquinol (PLQol) from the complex. In addition to the two known channels, we observe a third channel for PLQ/PLQol diffusion between the thylakoid membrane and the PLQ binding sites. Our simulations point to a promiscuous diffusion mechanism in which all three channels function as entry and exit channels. The exchange cavity serves as a PLQ reservoir. Our simulations provide a direct view on the exchange of electron carriers, a key step of the photosynthesis machinery., Plastoquinone (PLQ) shuttles electrons between photosystem II (PSII) and cytochrome b6f. Here the authors perform molecular dynamics simulations and propose that PLQ enters the exchange cavity of PSII by a promiscuous diffusion mechanism whereby three different channels each act as entry and exit points.

Stem cells are known to differentiate in response to the chemical and mechanical properties of the substrates on which they are cultured. Thus, supramolecular biomaterials with tunable properties are well suited for the study of stem cell differentiation. In this report, we exploited this phenomenon by combining stem cell differentiation in hydrogels with variable stiffness and metabolomics analysis to identify specific bioactive lipids that are uniquely used up during differentiation. To achieve this, we cultured perivascular stem cells on supramolecular peptide gels of different stiffness, and metabolite depletion followed. On soft (1 kPa), stiff (13 kPa), and rigid (32 kPa) gels, we observed neuronal, chondrogenic, and osteogenic differentiation, respectively, showing that these stem cells undergo stiffness-directed fate selection. By analyzing concentration variances of >600 metabolites during differentiation on the stiff and rigid gels (and focusing on chondrogenesis and osteogenesis as regenerative targets, respectively), we identified that specific lipids (lysophosphatidic acid and cholesterol sulfate, respectively), were significantly depleted. We propose that these metabolites are therefore involved in the differentiation process. In order to unequivocally demonstrate that the lipid metabolites that we identified play key roles in driving differentiation, we subsequently demonstrated that these individual lipids can, when fed to standard stem cell cultures, induce differentiation toward chondrocyte and osteoblast phenotypes. Our concept exploits the design of supramolecular biomaterials as a strategy for discovering cell-directing bioactive metabolites of therapeutic relevance.

The thylakoid membrane has a unique lipid composition, consisting mostly of galactolipids. These thylakoid lipids have important roles in photosynthesis. Here, we investigate to what extent these lipids bind specifically to the Photosystem II complex. To this end, we performed coarse-grain MD simulations of the Photosystem II complex embedded in a thylakoid membrane with realistic composition. Based on >85 μs simulation time, we find that monogalactosyldiacylglycerol and sulfoquinovosyldiacylglycerol lipids are enriched in the annular shell around the protein, and form distinct binding sites. From the analysis of residue contacts, we conclude that electrostatic interactions play an important role in stabilizing these binding sites. Furthermore, we find that chlorophyll a has a prevalent role in the coordination of the lipids. In addition, we observe lipids to diffuse in and out of the plastoquinone exchange cavities, allowing exchange of cocrystallized lipids with the bulk membrane and suggesting a more open nature of the plastoquinone exchange cavity. Together, our data provide a wealth of information on protein-lipid interactions for a key protein in photosynthesis.

Ultrasound, i.e. high frequency oscillating pressure waves, is commonly used to overcome kinetic barriers associated with dissolution, assembly and gelation.We demonstrate that ultrasound energy may also be used to achieve transient reorganization of supramolecular nanostructures, which revert back to the original state when sound is switched off. Aromatic peptide amphiphiles, Fmoc-FL and -YL were used to study the transient acoustic response. These systems showed temporary supramolecular transitions that were sequence dependent. The changes observed were due to an altered balance between H-bonding and p-stacking, giving rise in changes in chiral organisation of peptide building blocks. Transient reconfiguration was visualized by TEM and changes in supramolecular interactions characterized by fluorescence, FT-IR and CD. Remarkably, significant differences are observed when compared to thermal heating, which shows relates to the oscillating and directional characteristics of ultrasound when delivering heat to a system.

In this article, we demonstrate the non-aqueous self-assembly of a low-molecular-mass organic gelator based on an electroactive p-type tetrathiafulvalene (TTF)-dipeptide bioconjugate. We show that a TTF moiety appended with diphenylalanine amide derivative (TTF-FF-NH2) self-assembles into one-dimensional nanofibers that further lead to the formation of self-supporting organogels in chloroform and ethyl acetate. Upon doping of the gels with electron acceptors (TCNQ/iodine vapor), stable two-component charge transfer gels are produced in chloroform and ethyl acetate. These gels are characterized by various spectroscopy (UV-vis-NIR, FTIR and CD), microscopy (AFM and TEM), rheology and cyclic voltammetry techniques. Furthermore, conductivity measurements performed on TTF-FF-NH2 xerogel nanofiber networks formed between gold electrodes on a glass surface indicate that these nanofibers show a remarkable enhancement in the conductivity after doping with TCNQ.

Bovine α-lactalbumin (aLA) and oleate (OA) form a complex that has been intensively studied for its tumoricidal activity. Small-angle X-ray scattering (SAXS) has revealed that this complex consists of a lipid core surrounded by partially unfolded protein. We call this type of complex a liprotide. Little is known of the molecular interactions between OA and aLA, and no technique has so far provided any high-resolution structure of a liprotide. Here we have used coarse-grained (CG) molecular dynamics (MD) simulations, isothermal titration calorimetry (ITC) and SAXS to investigate the interactions between aLA and OA during the process of liprotide formation. With ITC we found that the strongest enthalpic interactions occurred at a molar ratio of 12.0±1.4:1 OA/aLA. Liprotides formed between OA and aLA at several OA/aLA ratios in silico were stable both in CG and in all-atom simulations. From the simulated structures we calculated SAXS spectra that show good agreement with experimentally measured patterns of matching liprotides. The simulations showed that aLA assumes a molten globular (MG) state, exposing several hydrophobic patches involved in interactions with OA. Initial binding of aLA to OA occurs in an area of aLA in which a high amount of positive charge is located, and only later do hydrophobic interactions become important. The results reveal how unfolding of aLA to expose hydrophobic residues is important for complex formation between aLA and OA. Our findings suggest a general mechanism for liprotide formation and might explain the ability of a large number of proteins to form liprotides with OA.

Photosystem II (PSII) is one of the key protein complexes in photosynthesis. We introduce a coarse grained model of PSII and present the analysis of 60 µs molecular dynamics simulations of PSII in both monomeric and dimeric form, embedded in a thylakoid membrane model that reflects its native lipid composition. We describe in detail the setup of the protein complex and the many natural cofactors, and characterize their mobility. Overall we find that the protein subunits and cofactors are more flexible towards the periphery of the complex, as well as near the PLQ exchange cavity and at the dimer interface. Of all cofactors, β-carotenes show the highest mobility. Some of the β-carotenes diffuse in and out of the protein complex via the thylakoid membrane. In contrast to the PSII dimer, the monomeric form adopts a tilted conformation in the membrane, with strong interactions between the soluble PsbO subunit and the glycolipid headgroups. Interestingly, the tilted conformation causes buckling of the membrane. Together, our results provide an unprecedented view of PSII dynamics on a microsecond time scale. Our data may be used as basis for the interpretation of experimental data as well as for theoretical models describing exciton energy transfer.

We demonstrate an in-situ ultrasonic approach to influence self-assembly across the supramolecular to micron length scales, showing enhancement of supramolecular interactions, chirality and orientation, which depends on the peptide sequence and solvent environment. This is the first successful demonstration of using oscillating pressure waves to generate anisotropic organo- and hydro- gels consisting of oriented tripeptides structures.

Biomimetic compounds based upon the active subsite of the [FeFe]-hydrogenase enzyme system have been the focus of much attention as catalysts for hydrogen production: a clean energy vector. Until recently, use of hydrogenase subsite systems for light-driven hydrogen production has typically required the involvement of a photosensitizer, but the molecule [(μ-pdt)(μ-H)Fe2(CO)4(dppv)]+, (1; dppv = cis-1,2-C2H2(PPh2)2; pdt = 1,3-propanedithiolate) has been reported to catalyze the evolution of hydrogen gas under sensitizer-free conditions. Establishing the molecular mechanism that leads to photohydrogen production by 1 is thus an important step that may enable further development of this family of molecules as solar fuel platforms. Here, we report ultrafast UVpump–IRprobe spectroscopy of 1 at three different excitation wavelengths and in a range of solvents, including under the conditions required for H2 production. Combining spectroscopic measurements of the photochemistry and vibrational relaxation dynamics of 1 with ground-state density functional theory (DFT) calculations shows that, irrespective of experimental conditions, near-instantaneous carbonyl ligand loss is the main photochemical channel. No evidence for a long-lived excited electronic state was found. These results provide the first time-resolved data for the photochemistry of 1 and offer an alternative interpretation of the underlying mechanism of light-driven hydrogen generation.

Supramolecular structures were produced by in situ enzymatic condensation of Fmoc–Phe–(4-X), where X denotes electron withdrawing or donating groups, with Phe–NH2. The relative contribution of π-stacking and H-bonding interactions can be regulated by the nature of X, resulting in tuneable nanoscale morphologies.

The results of infrared spectroscopic investigations into the band assignments, vibrational relaxation, and solvation dynamics of the common anti-tuberculosis treatment Isoniazid (INH) are reported. INH is known to inhibit InhA, a 2-trans-enoyl-acyl carrier protein reductase enzyme responsible for the maintenance of cell walls in Mycobacterium tuberculosis but as new drug-resistant strains of the bacterium appear, next-generation therapeutics will be essential to combat the rise of the disease. Small molecules such as INH offer the potential for use as a biomolecular marker through which ultrafast multidimensional spectroscopies can probe drug binding and so inform design strategies but a complete characterization of the spectroscopy and dynamics of INH in solution is required to inform such activity. Infrared absorption spectroscopy, in combination with density functional theory calculations, is used to assign the vibrational modes of INH in the 1400-1700 cm(-1) region of the infrared spectrum while ultrafast multidimensional spectroscopy measurements determine the vibrational relaxation dynamics and the effects of solvation via spectral diffusion of the carbonyl stretching vibrational mode. These results are discussed in the context of previous linear spectroscopy studies on solid-phase INH and its usefulness as a biomolecular probe.