Watching Molecular Nanotubes Self-Assemble in Real Time.

Molecular self-assembly is a fundamental process in nature that can be used to develop novel functional materials for medical and engineering applications. However, their complex mechanisms make the short-lived stages of self-assembly processes extremely hard to reveal. In this article, we track the self-assembly process of a benchmark system, double-walled molecular nanotubes, whose structure is similar to that found in biological and synthetic systems. We selectively dissolved the outer wall of the double-walled system and used the inner wall as a template for the self-reassembly of the outer wall. The reassembly kinetics were followed in real time using a combination of microfluidics, spectroscopy, cryogenic transmission electron microscopy, molecular dynamics simulations, and exciton modeling. We found that the outer wall self-assembles through a transient disordered patchwork structure: first, several patches of different orientations are formed, and only on a longer time scale will the patches interact with each other and assume their final preferred global orientation. The understanding of patch formation and patch reorientation marks a crucial step toward steering self-assembly processes and subsequent material engineering.

Pragmatic Coarse-Graining of Proteins: Models and Applications.

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Modelling structural properties of cyanine dye nanotubes at coarse-grained level.

Self-assembly is a ubiquitous process spanning from biomolecular aggregates to nanomaterials. Even though the resulting aggregates can be studied through experimental techniques, the dynamic pathways of the process and the molecular details of the final structures are not necessarily easy to resolve. Consequently, rational design of self-assembling aggregates and their properties remains extremely challenging. At the same time, modelling the self-assembly with computational methods is not trivial, because its spatio-temporal scales are usually beyond the limits of all-atom based simulations. The use of coarse-grained (CG) models can alleviate this limitation, but usually suffers from the lack of optimised parameters for the molecular constituents. In this work, we describe the procedure of parametrizing a CG Martini model for a cyanine dye (C8S3) that self-assembles into hollow double-walled nanotubes. First, we optimised the model based on quantum mechanics calculations and all-atom reference simulations, in combination with available experimental data. Then, we conducted random self-assembly simulations, and the performance of our model was tested on preformed assemblies. Our simulations provide information on the time-dependent local arrangement of this cyanine dye, when aggregates are being formed. Furthermore, we provide guidelines for designing and optimising parameters for similar self-assembling nanomaterials.

Cryogenic TEM imaging of artificial light harvesting complexes outside equilibrium.

The energy transport in natural light-harvesting complexes can be explored in laboratory conditions via self-assembled supramolecular structures. One such structure arises from the amphiphilic dye C8S3 molecules, which self-assemble in an aqueous medium to a double-wall cylindrical nanotube reminiscent of natural light-harvesting complexes found in green sulphur bacteria. In this paper, we report a way to investigate the structure of inner nanotubes (NTs) alone by dissolving the outer NTs in a microfluidic setting. The resulting thermodynamically unstable system was rapidly frozen, preventing the reassembly of the outer NT from the dissolved molecules, and imaged using cryogenic transmission electron microscopy (cryo-TEM). The experimental cryo-TEM images and the molecular structure were compared by simulating high-resolution TEM images, which were based on the molecular modelling of C8S3 NTs. We found that the inner NT with outer walls removed during the flash-dilution process had a similar size to the parent double-walled NTs. Moreover, no structural inhomogeneity was observed in the inner NT after flash-dilution. This opens up exciting possibilities for functionalisation of inner NTs before the reassembly of the outer NT occurs, which can be broadly extended to modify the intra-architecture of other self-assembled nanostructures.

Comparing Dimerization Free Energies and Binding Modes of Small Aromatic Molecules with Different Force Fields.

Dimerization free energies are fundamental quantities that describe the strength of interaction of different molecules. Obtaining accurate experimental values for small molecules and disentangling the conformations that contribute most to the binding can be extremely difficult, due to the size of the systems and the small energy differences. In many cases, one has to resort to computational methods to calculate such properties. In this work, we used molecular dynamics simulations in conjunction with metadynamics to calculate the free energy of dimerization of small aromatic rings, and compared three models from popular online servers for atomistic force fields, namely G54a7, CHARMM36 and OPLS. We show that, regardless of the force field, the profiles for the dimerization free energy of these compounds are very similar. However, significant care needs to be taken when studying larger molecules, since the deviations from the trends increase with the size of the molecules, resulting in force field dependent preferred stacking modes; for example, in the cases of pyrene and tetracene. Our results provide a useful background study for using topology builders to model systems which rely on stacking of aromatic moieties, and are relevant in areas ranging from drug design to supramolecular assembly.

Martini 3: a general purpose force field for coarse-grained molecular dynamics.

The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.

Molecular versus Excitonic Disorder in Individual Artificial Light-Harvesting Systems.

Natural light-harvesting antennae employ a dense array of chromophores to optimize energy transport via the formation of delocalized excited states (excitons), which are critically sensitive to spatio-energetic variations of the molecular structure. Identifying the origin and impact of such variations is highly desirable for understanding and predicting functional properties yet hard to achieve due to averaging of many overlapping responses from individual systems. Here, we overcome this problem by measuring the heterogeneity of synthetic analogues of natural antennae-self-assembled molecular nanotubes-by two complementary approaches: single-nanotube photoluminescence spectroscopy and ultrafast 2D correlation. We demonstrate remarkable homogeneity of the nanotube ensemble and reveal that ultrafast (∼50 fs) modulation of the exciton frequencies governs spectral broadening. Using multiscale exciton modeling, we show that the dominance of homogeneous broadening at the exciton level results from exchange narrowing of strong static disorder found for individual molecules within the nanotube. The detailed characterization of static and dynamic disorder at the exciton as well as the molecular level presented here opens new avenues in analyzing and predicting dynamic exciton properties, such as excitation energy transport.

Structural characterization of supramolecular hollow nanotubes with atomistic simulations and SAXS.

Self-assembled nanostructures arise when building blocks spontaneously organize into ordered aggregates that exhibit different properties compared to the disorganized monomers. Here, we study an amphiphilic cyanine dye (C8S3) that is known to self-assemble into double-walled, hollow, nanotubes with interesting optical properties. The molecular packing of the dyes inside the nanotubes, however, remains elusive. To reveal the structural features of the C8S3 nanotubes, we performed atomistic Molecular Dynamics simulations of preformed bilayers and nanotubes. We find that different packing arrangements lead to stable structures, in which the tails of the C8S3 molecules are interdigitated. Our results are verified by SAXS experiments. Together our data provide a detailed structural characterization of the C8S3 nanotubes. Furthermore, our approach was able to resolve the ambiguity inherent from cryo-TEM measurements in calculating the wall thickness of similar systems. The insights obtained are expected to be generally useful for understanding and designing other supramolecular assemblies.

Multiscale modeling of molecular structure and optical properties of complex supramolecular aggregates.

Supramolecular aggregates of synthetic dye molecules offer great perspectives to prepare biomimetic functional materials for light-harvesting and energy transport. The design is complicated by the fact that structure-property relationships are hard to establish, because the molecular packing results from a delicate balance of interactions and the excitonic properties that dictate the optics and excited state dynamics, in turn sensitively depend on this packing. Here we show how an iterative multiscale approach combining molecular dynamics and quantum mechanical exciton modeling can be used to obtain accurate insight into the packing of thousands of cyanine dye molecules in a complex double-walled tubular aggregate in close interaction with its solvent environment. Our approach allows us to answer open questions not only on the structure of these prototypical aggregates, but also about their molecular-scale structural and energetic heterogeneity, as well as on the microscopic origin of their photophysical properties. This opens the route to accurate predictions of energy transport and other functional properties.

Fluorinated Alcohols' Effects on Lipid Bilayer Properties.

Fluorinated alcohols (fluoroalcohols) have physicochemical properties that make them excellent solvents of peptides, proteins, and other compounds. Like other alcohols, fluoroalcohols also alter membrane protein function and lipid bilayer properties and stability. Thus, the questions arise: how potent are fluoroalcohols as lipid-bilayer-perturbing compounds, could small residual amounts that remain after adding compounds dissolved in fluoroalcohols alter lipid bilayer properties sufficiently to affect membranes and membrane protein function, and do they behave like other alcohols? To address these questions, we used a gramicidin-based fluorescence assay to determine the bilayer-modifying potency of selected fluoroalcohols: trifluoroethanol (TFE), HFIP, and perfluoro-tert-butanol (PFTB). These fluoroalcohols alter bilayer properties in the low (PFTB) to high (TFE) mM range. Using the same assay, we determined the bilayer partitioning of the alcohols. When referenced to the aqueous concentrations, the fluoroalcohols are more bilayer perturbing than their nonfluorinated counterparts, with the largest fluoroalcohol, PFTB, being the most potent and the smallest, TFE, the least. When referenced to the mole fractions in the membrane, however, the fluoroalcohols have equal or lesser bilayer-perturbing potency than their nonfluorinated counterparts, with TFE being more bilayer perturbing than PFTB. We compared the fluoroalcohols' molecular level bilayer interactions using atomistic molecular dynamics simulations and showed how, at higher concentrations, they can cause bilayer breakdown using absorbance measurements and 31P nuclear magnetic resonance.

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.

Targeted cutaneous delivery of ciclosporin A using micellar nanocarriers and the possible role of inter-cluster regions as molecular transport pathways.

Oral administration of ciclosporin A (CsA) is indicated in the treatment of severe recalcitrant plaque psoriasis. However, CsA is both nephro- and hepatotoxic and its systemic administration also exposes the patient to other severe side effects. Although topical delivery of CsA, targeted directly to psoriatic skin, would offer significant advantages, there are no topical formulations approved for dermatological use. The aim of this work was to formulate CsA loaded polymeric micelles using the biodegradable and biocompatible MPEG-dihexPLA diblock copolymer and to evaluate their potential for delivering the drug selectively into the skin without concomitant transdermal permeation. Micelle formulations were characterised with respect to drug content, size and morphology. Micelle and drug penetration pathways were subsequently visualised with confocal laser scanning microscopy (CLSM) using fluorescein labelled CsA (Fluo-CsA) and Nile-Red (NR) labelled copolymer. Visualisation studies typically use fluorescent dyes as "model drugs"; however, these may have different physicochemical properties to the drug molecule under investigation. Therefore, in this study it was decided to chemically modify CsA and to use this structurally similar fluorescent analogue to visualise molecular distribution and transport pathways. Molecular modelling techniques and experimental determination of log D served as molecular scale and macroscopic methods to compare the lipophilicity of CsA and Fluo-CsA. The spherical, homogeneous and nanometre-scale micelles (with Zav from 25 to 52 nm) increased the aqueous solubility of CsA by 518-fold. Supra-therapeutic amounts of CsA were delivered to human skin (1.4±0.6 µg/cm2, cf. a statistically equivalent 1.1±0.5 µg/cm2 for porcine skin) after application of the formulation with the lowest CsA and copolymer content (1.67±0.03 mg/ml of CsA and 5mg/ml of copolymer) for only 1h without concomitant transdermal permeation. Fluo-CsA was successfully synthesised, characterised and incorporated into fluorescent NR-MPEG-dihexPLA micelles; its conformation was not modified by the addition of fluorescein and its log D, measured from pH4 to 8, was equivalent to that of CsA. Fluo-CsA and NR-MPEG-dihexPLA copolymer were subsequently visualised in skin by CLSM. The images indicated that micelles were preferentially deposited between corneocytes and in the inter-cluster regions (i.e. between the clusters of corneocytes). Fluo-CsA skin penetration was deeper in these structures, suggesting that inter-cluster penetration is probably the preferred transport pathway responsible for the increased cutaneous delivery of CsA.

As good as it gets? Folding molecular dynamics simulations of the LytA choline-binding peptide result to an exceptionally accurate model of the peptide structure.

Folding simulations of a choline-binding peptide derived from the Streptococcus pneumoniae LytA protein converged to a model of the peptide's folded state structure which is in outstanding agreement with the experimentally-determined structures, reaching values for the root mean squared deviation as low as 0.24Å for the peptide's backbone atoms and 0.65Å for all non-hydrogen atoms.