Stochastic Emergence of Two Distinct Self-Replicators from a Dynamic Combinatorial Library

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.

Coacervate formation studied by explicit solvent coarse-grain molecular dynamics with the Martini model

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.

Caught in the Act: Mechanistic Insight into Supramolecular Polymerization-Driven Self-Replication from Real-Time Visualization

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.

Supramolecular Nucleoside-Based Gel: Molecular Dynamics Simulation and Characterization of Its Nanoarchitecture and Self-Assembly Mechanism

Among the diversity of existing supramolecular hydrogels, nucleic acid-based hydrogels are of particular interest for potential drug delivery and tissue engineering applications because of their inherent biocompatibility. Hydrogel performance is directly related to the nanostructure and the self-assembly mechanism of the material, an aspect that is not well-understood for nucleic acid-based hydrogels in general and has not yet been explored for cytosine-based hydrogels in particular. Herein, we use a broad range of experimental characterization techniques along with molecular dynamics (MD) simulation to demonstrate the complementarity and applicability of both approaches for nucleic acid-based gelators in general and propose the self-assembly mechanism for a novel supramolecular gelator, N4-octanoyl-2'-deoxycytidine. The experimental data and the MD simulation are in complete agreement with each other and demonstrate the formation of a hydrophobic core within the fibrillar structures of these mainly water-containing materials. The characterization of the distinct duality of environments in this cytidine-based gel will form the basis for further encapsulation of both small hydrophobic drugs and biopharmaceuticals (proteins and nucleic acids) for drug delivery and tissue engineering applications.

Molecular simulations of self-assembling bio-inspired supramolecular systems and their connection to 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.

Polymeric peptide pigments with sequence-encoded properties

Melanins are a family of heterogeneous polymeric pigments that provide ultraviolet (UV) light protection, structural support, coloration, and free radical scavenging. Formed by oxidative oligomerization of catecholic small molecules, the physical properties of melanins are influenced by covalent and noncovalent disorder. We report the use of tyrosine-containing tripeptides as tunable precursors for polymeric pigments. In these structures, phenols are presented in a (supra-)molecular context dictated by the positions of the amino acids in the peptide sequence. Oxidative polymerization can be tuned in a sequence-dependent manner, resulting in peptide sequence-encoded properties such as UV absorbance, morphology, coloration, and electrochemical properties over a considerable range. Short peptides have low barriers to application and can be easily scaled, suggesting near-term applications in cosmetics and biomedicine.

Prediction of Thylakoid Lipid Binding Sites on Photosystem II

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.

Structural and Spectroscopic Properties of Assemblies of Self-Replicating Peptide Macrocycles

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.

Exchange pathways of plastoquinone and plastoquinol in the photosystem II complex

Plastoquinone (PLQ) acts as an electron carrier between photosystem II (PSII) and the cytochrome b₆f 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.

Molecular Dynamics of Photosystem II Embedded in the Thylakoid Membrane

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 toward 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 with 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.

Alignment of nanostructured tripeptide gels by directional ultrasonication

We demonstrate an in situ ultrasonic approach to generate anisotropic organo- and hydrogels consisting of oriented tripeptides self-assembled structures.

MMP-9 triggered micelle-to-fibre transitions for slow release of doxorubicin

Phenylacetyl-peptide amphiphiles were designed, which upon cleavage by a disease-associated enzyme reconfigure from micellar aggregates to fibres. Upon this morphological change, a doxorubicin payload could be retained in the fibres formed, which makes them valuable carriers for localised formation of nanofibre depots for slow release of hydrophobic anticancer drugs.

Conducting nanofibers and organogels derived from the self-assembly of tetrathiafulvalene-appended dipeptides

We demonstrate the nonaqueous 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.

Stable emulsions formed by self-assembly of interfacial networks of dipeptide derivatives

We demonstrate the use of dipeptide amphiphiles that, by hand shaking of a biphasic solvent system for a few seconds, form emulsions that remain stable for months through the formation of nanofibrous networks at the organic/aqueous interface. Unlike absorption of traditional surfactants, the interfacial networks form by self-assembly through π-stacking interactions and hydrogen bonding. Altering the dipeptide sequence has a dramatic effect on the properties of the emulsions formed, illustrating the possibility of tuning emulsion properties by chemical design. The systems provide superior long-term stability toward temperature and salts compared to with sodium dodecyl sulfate (SDS) and can be enzymatically disassembled causing on-demand demulsification under mild conditions. The interfacial networks facilitate highly tunable and stable encapsulation and compartmentalization with potential applications in cosmetics, therapeutics, and food industry.

Differential self-assembly and tunable emission of aromatic peptide bola-amphiphiles containing perylene bisimide in polar solvents including water

We demonstrate the self-assembly of bola-amphiphile-type conjugates of dipeptides and perylene bisimide (PBI) in water and other polar solvents. Depending on the nature of the peptide used (glycine-tyrosine, GY, or glycine-aspartic acid, GD), the balance between H-bonding and aromatic stacking can be tailored. In aqueous buffer, PBI-[GY]2 forms chiral nanofibers, resulting in the formation of a hydrogel, while for PBI-[GD]2 achiral spherical aggregates are formed, demonstrating that the peptide sequence has a profound effect on the structure formed. In water and a range of other polar solvents, self-assembly of these two PBI-peptides conjugates results in different nanostructures with highly tunable fluorescence performance depending on the peptide sequence employed, e.g., fluorescent emission and quantum yield. Organogels are formed for the PBI-[GD]2 derivative in DMF and DMSO while PBI-[GY]2 gels in DMF. To the best of our knowledge, this is the first successful strategy for using short peptides, specifically, their sequence/structure relationships, to manipulate the PBI nanostructure and consequent optical properties. The combination of controlled self-assembly, varied optical properties, and formation of aqueous and organic gel-phase materials may facilitate the design of devices for various applications related to light harvesting and sensing.

Aromatic peptide amphiphiles: significance of the Fmoc moiety

Aromatic peptide amphiphile hydrogelators commonly utilise the fluorenyl-9-methoxycarbonyl moiety as an N-terminal capping group. Material properties and spectroscopic techniques show the influence of alternative linkers between the fluorenyl moiety and the peptide. This study establishes whether methoxycarbonyl is an optimal or mainly convenient linker, for this class of self-assembling systems.

Biocatalytically triggered co-assembly of two-component core/shell nanofibers

For the development of applications and novel uses for peptide nanostructures, robust routes for their surface functionalization, that ideally do not interfere with their self-assembly properties, are required. Many existing methods rely on covalent functionalization, where building blocks are appended with functional groups, either pre- or post-assembly. A facile supramolecular approach is demonstrated for the formation of functionalized nanofibers by combining the advantages of biocatalytic self-assembly and surfactant/gelator co-assembly. This is achieved by enzymatically triggered reconfiguration of free flowing micellar aggregates of pre-gelators and functional surfactants to form nanofibers that incorporate and display the surfactants' functionality at the surface. Furthermore, by varying enzyme concentration, the gel stiffness and supramolecular organization of building blocks can be varied.

Insights into the coassembly of hydrogelators and surfactants based on aromatic peptide amphiphiles

The coassembly of small molecules is a useful means of increasing the complexity and functionality of their resultant supramolecular constructs in a modular fashion. In this study, we explore the assembly and coassembly of serine surfactants and tyrosine-leucine hydrogelators, capped at the N-termini with either fluorenyl-9-methoxycarbonyl (Fmoc) or pyrene. These systems all exhibit self-assembly behavior, which is influenced by aromatic stacking interactions, while the hydrogelators also exhibit β-sheet-type arrangements, which reinforce their supramolecular structures. We provide evidence for three distinct supramolecular coassembly models; cooperative, disruptive, and orthogonal. The coassembly mode adopted depends on whether the individual constituents (I) are sufficiently different, such that effective segregation and orthogonal assembly occurs; (II) adhere to a communal mode of self-assembly; or (III) act to compromise the assembly of one another via incorporation and disruption. We find that a greater scope for controllable coassembly exists within orthogonal systems; which show minimal relative changes in the native gelator's supramolecular structure by Fourier transform infrared spectroscopy (FTIR), circular dichroism (CD), and fluorescence spectroscopy. This is indicative of the segregation of orthogonal coassembly constituents into distinct domains, where surfactant chemical functionality is presented at the surface of the gelator's supramolecular fibers. Overall, this work provides new insights into the design of modular coassembly systems, which have the potential to augment the chemical and physical properties of existing gelator systems.

Tuneable Fmoc–Phe–(4-X)–Phe–NH2 nanostructures by variable electronic substitution

We show enzymatic introduction of non-natural amino acids with different electronic substituents with a dramatic influence on morphology in self-assembled nanostructures.

Infrared spectroscopy of nicotinamide adenine dinucleotides in one and two dimensions

The development of multidimensional spectroscopic tools capable of resolving site-specific information about proteins and enzymes in the solution phase is an important aid to our understanding of biomolecular mechanisms, structure, and dynamics. Nicotinamide adenine dinucleotide (NAD) is a common biological substrate and so offers significant potential as an intrinsic vibrational probe of protein-ligand interactions but its complex molecular structure and incompletely characterized infrared spectrum currently limit its usefulness. Here, we report the FTIR spectroscopy of the oxidized and reduced forms of NAD at a range of pD values that relate to the "folded" and "unfolded" forms of the molecules that exist in solution. Comparisons with structural analogs and the use of density functional theory simulations provide a full assignment of the observed modes and their complex pD dependencies. Finally, ultrafast two-dimensional infrared spectra of the oxidized and reduced forms of NAD are reported and their usefulness as biomolecular probes is discussed.

Assessing the utility of infrared spectroscopy as a structural diagnostic tool for β-sheets in self-assembling aromatic peptide amphiphiles

β-Sheets are a commonly found structural motif in self-assembling aromatic peptide amphiphiles, and their characteristic "amide I" infrared (IR) absorption bands are routinely used to support the formation of supramolecular structure. In this paper, we assess the utility of IR spectroscopy as a structural diagnostic tool for this class of self-assembling systems. Using 9-fluorene-methyloxycarbonyl dialanine (Fmoc-AA) and the analogous 9-fluorene-methylcarbonyl dialanine (Fmc-AA) as examples, we show that the origin of the band around 1680-1695 cm(-1) in Fourier transform infrared (FTIR) spectra, which was previously assigned to an antiparallel β-sheet conformation, is in fact absorption of the stacked carbamate group in Fmoc-peptides. IR spectra from (13)C-labeled samples support our conclusions. In addition, DFT frequency calculations on small stacks of aromatic peptides help to rationalize these results in terms of the individual vibrational modes.

Photodissociation of van der Waals clusters of isoprene with oxygen, C5H8-O2, in the wavelength range 213-277 nm

The speed and angular distribution of O atoms arising from the photofragmentation of C(5)H(8)-O(2), the isoprene-oxygen van der Waals complex, in the wavelength region of 213-277 nm has been studied with the use of a two-color dissociation-probe method and the velocity map imaging technique. Dramatic enhancement in the O atoms photo-generation cross section in comparison with the photodissociation of individual O(2) molecules has been observed. Velocity map images of these "enhanced" O atoms consisted of five channels, different in their kinetic energy, angular distribution, and wavelength dependence. Three channels are deduced to be due to the one-quantum excitation of the C(5)H(8)-O(2) complex into the perturbed Herzberg III state ((3)Δ(u)) of O(2). This excitation results in the prompt dissociation of the complex giving rise to products C(5)H(8)+O+O when the energy of exciting quantum is higher than the complex photodissociation threshold, which is found to be 41740 ± 200 cm(-1) (239.6±1.2 nm). This last threshold corresponds to the photodissociation giving rise to an unexcited isoprene molecule. The second channel, with threshold shifted to the blue by 1480 ± 280 cm(-1), corresponds to dissociation with formation of rovibrationally excited isoprene. A third channel was observed at wavelengths up to 243 nm with excitation below the upper photodissociation threshold. This channel is attributed to dissociation with the formation of a bound O atom C(5)H(8)-O(2) + hv → C(5)H(8)-O(2)((3)Δ(u)) → C(5)H(8)O + O and/or to dissociation of O(2) with borrowing of the lacking energy from incompletely cooled complex internal degrees of freedom C(5)H(8)*-O(2) + hv → C(5)H(8)*-O(2)((3)Δ(u)) → C(5)H(8) + O + O. The kinetic energy of the O atoms arising in two other observed channels corresponds to O atoms produced by photodissociation of molecular oxygen in the excited a (1)Δ(g) and b (1)Σ(g)(+) singlet states as the precursors. This indicates the formation of singlet oxygen O(2)(a (1)Δ(g)) and O(2)(b (1)Σ(g)(+)) after excitation of the C(5)H(8)-O(2) complex. Cooperative excitation of the complex with a simultaneous change of the spin of both partners (1)X-(3)O(2) + hν → (3)X-(1)O(2) → (3)X + (1)O(2) is suggested as a source of singlet oxygen O(2)(a (1)Δ(g)) and O(2)(b (1)Σ(g)(+)). This cooperative excitation is in agreement with little or no vibrational excitation of O(2)(a (1)Δ(g)), produced from the C(5)H(8)-O(2) complex as studied in the current paper as well as from the C(3)H(6)-O(2) and CH(3)I-O(2) complexes reported in our previous paper [Baklanov et al., J. Chem. Phys. 126, 124316 (2007)]. The formation of O(2)(a (1)Δ(g)) from C(5)H(8)-O(2) was observed at λ(pump) = 213-277 nm with the yield going down towards the long wavelength edge of this interval. This spectral profile is interpreted as the red-side wing of the band of a cooperative transition (1)X-(3)O(2) + hν → (3)X(T(2))-(1)O(2)(a (1)Δ(g)) in the C(5)H(8)-O(2) complex.

Virtual Screening for Dipeptide Aggregation: Toward Predictive Tools for Peptide Self-Assembly

Several short peptide sequences are known to self-assemble into supramolecular nanostructures with interesting properties. In this study, coarse-grained molecular dynamics is employed to rapidly screen all 400 dipeptide combinations and predict their ability to aggregate as a potential precursor to their self-assembly. The simulation protocol and scoring method proposed allows a rapid determination of whether a given peptide sequence is likely to aggregate (an indicator for the ability to self-assemble) under aqueous conditions. Systems that show strong aggregation tendencies in the initial screening are selected for longer simulations, which result in good agreement with the known self-assembly or aggregation of dipeptides reported in the literature. Our extended simulations of the diphenylalanine system show that the coarse-grain model is able to reproduce salient features of nanoscale systems and provide insight into the self-assembly process for this system.