Towards hybrid molecular simulations

In many biology, chemistry and physics applications molecular simulations can be used to study material and process properties. The level of detail needed in such simulations depends on the application. In some cases quantum mechanical simulations are indispensable. However, traditional ab-initio methods, usually employing plane waves or a linear combination of atomic orbitals as a basis, are extremely expensive in terms of computational as well as memory requirements. The well-known fact that electronic wave functions vary much more rapidly near the atomic nuclei than in inter-atomic regions calls for a multi-resolution approach, allowing one to use low resolution and to add extra resolution only in those regions where necessary, so limiting the costs. This is provided by an alternative basis formed of wavelets. Using such a wavelet basis, a method has been developed for solving electronic structure problems that has been applied successfully to 2D quantum dots and 3D molecular systems. In other cases, it suffices to use effective potentials to describe the atomic interaction instead of the use of the electronic structure, enabling the simulation of larger systems. Molecular dynamics simulations with such effective potentials have been used for a systematic study of surface wettability influence on particle and heat flow in nanochannels, showing that the effects at the solid-gas interface are crucial for the behavior of the whole nanochannel. Again in other cases even coarse grained models can be used where the average behavior of several atoms is combined into a single particle. Such a model, refraining from as much detail as possible while maintaining realistic behavior, has been developed for lipids and with this model the dynamics of membranes and vesicle formation have been studied in detail. A disadvantage of molecular dynamics simulations with effective potentials is that no reactions are possible. Therefore a new method has been developed, where molecular dynamics is coupled with stochastic reactions. Using this method, both unilamellar and multilamellar vesicle formation, and vesicle growth, bursting, and healing are shown. Still larger systems can be simulated using other methods, like the direct simulation Monte Carlo method. However, as shown for nanochannels, these methods are not always accurate enough. But, exploiting again that the finest level of detail is often only needed in part of the domain, a hybrid method has been developed coupling molecular dynamics, where needed for accuracy, and direct simulation Monte Carlo, where possible in order to speed up the calculation. Further development of such hybrid simulations will further increase molecular simulation’s scientific role

Coarse Grained Molecular Dynamics Simulations of the Fusion of Vesicles Incorporating Water Channels

As the dynamics of the cell membrane and the working mechanisms of proteins cannot be readily asserted at a molecular level, many different hypotheses exist that try to predict and explain these processes, for instance vesicle fusion. Therefore, we use coarse grained molecular dynamics simulations to elucidate the fusion mechanism of vesicles. The implementation of this method with hydrophilic and hydrophobic particles is known for its valid representation of bilayers. With a minimalistic approach, using only 3 atom types, 12 atoms per two-tailed phospholipids and incorporating only a bond potential and Lennard-Jones potential, phospholipid bilayers and vesicles can be simulated exhibiting authentic dynamics. We have simulated the spontaneous full fusion of both tiny (6 nm diameter) and larger (13 nm diameter) vesicles. We showed that, without applying constraints to the vesicles, the initial contact between two fusing vesicles, the stalk, is initiated by a bridging lipid tail that extends from the membrane spontaneously. Subsequently it is observed that the evolution of the stalk can proceed via two pathways, anisotropic and radial expansion, which is in accordance with literature. Contrary to the spherical vesicles of in vitro experiments, the fused vesicles remain tubular since the internal volume of these vesicles is too small compared to their membrane area. While the lipid bilayer has some permeability for water, it is not high enough to allow for the large flux required to equilibrate the vesicle content in the time accessible to our simulations. To increase the membrane permeability, we incorporate proteinaceous water channels, by applying the coarse grained technique to aquaporin. Even though incorporating water channels in the vesicles does significantly increase water permeability, the vesicles do not become spherical. Presumably the lipids have to be redistributed as well

Coarse Grained Molecular Dynamics Simulations of Transmembrane Protein-Lipid Systems

Many biological cellular processes occur at the micro- or millisecond time scale. With traditional all-atom molecular modeling techniques it is difficult to investigate the dynamics of long time scales or large systems, such as protein aggregation or activation. Coarse graining (CG) can be used to reduce the number of degrees of freedom in such a system, and reduce the computational complexity. In this paper the first version of a coarse grained model for transmembrane proteins is presented. This model differs from other coarse grained protein models due to the introduction of a novel angle potential as well as a hydrogen bonding potential. These new potentials are used to stabilize the backbone. The model has been validated by investigating the adaptation of the hydrophobic mismatch induced by the insertion of WALP-peptides into a lipid membrane, showing that the first step in the adaptation is an increase in the membrane thickness, followed by a tilting of the peptide

LiFE Assessment Tool

As part of an ongoing study to construct a molecular Turing machine in which a polymer chain is encoded via allosteric information transfer between macrocyclic complexes, we describe the thermodynamic and kinetic characterization of a multicomponent self-assembled system based on a zinc porphyrin macrocyclic compound, a bidentate ligand (1,4-diazabicyclo[2.2.2]octane, DABCO), and a viologen-substituted polymer guest. Initial addition of DABCO to the porphyrin macrocycle in chloroform solution leads to the formation of a stable 2:1 (porphyrin:DABCO) dimeric complex, even under dilute conditions, by means of strong cooperative interactions involving hydrogen and metal-ligand bonds. Further titration of the porphyrin-DABCO mixtures with the polymer gives rise to a complex array of species in the solution. The system is analyzed in detail by a combination of spectroscopic measurements and computational modeling. Each association constant in the binding scheme and the fraction of each individual complex that is formed in solution are determined precisely using a mass-balance model. Kinetic studies revealed that the rates of the polymer threading and dethreading in and out of the dimeric system are remarkably slow, indicating that the polymer is locked inside the cavity of the stable 2:1 dimeric complex as a result of strong allosteric interactions

Unequal Perylene Diimide Twins in a Quadruple Assembly

Natural light-harvesting (LH) systems can divide identical dyes into unequal aggregate states, thereby achieving intelligent "allocation of labor". From a synthetic point of view, the construction of such kinds of unequal and integrated systems without the help of proteinaceous scaffolding is challenging. Here, we show that four octatetrayne-bridged ortho-perylene diimide (PDI) dyads (POPs) self-assemble into a quadruple assembly (POP)4 both in solution and in the solid state. The two identical PDI units in each POP are compartmentalized into weakly coupled PDIs (P520) and closely stacked PDIs (P550) in (POP)4 . The two extreme pools of PDI chromophores were unambiguously confirmed by single-crystal X-ray crystallography and NMR spectroscopy. To interpret the formation of the discrete quadruple assembly, we also developed a two-step cooperative model. Quantum-chemical calculations indicate the existence of multiple couplings within and across P520 and P550, which can satisfactorily describe the photophysical properties of the unequal quadruple assembly. This finding is expected to help advance the rational design of dye stacks to emulate functions of natural LH systems.</p

Supramolecular copolymers: structure and composition revealed by theoretical modeling

Supramolecular copolymers, non-covalent analogues of synthetic copolymers, constitute a new and promising class of polymers. In contrast to their covalent counterparts, the details of their mechanism of formation, as well as the factors determining their composition and length, are still poorly understood. Here, the supramolecular copolymerization between two slightly structurally different benzene-1,3,5-tricarboxamide (BTA) monomers functionalized with either oligodimethylsiloxane (oDMSi) or alkyl side chains is unraveled by combining experimental and theoretical approaches. By applying the “sergeant-and-soldiers” approach using circular dichroism (CD) experiments, we are able to obtain detailed insights into the structure and composition of these supramolecular copolymers. Moreover, we observe an unexpected chiral induction upon mixing two independently CD-silent solutions of the achiral (soldier) and chiral (sergeant) monomers. We find that the subtle differences in the chemical structure of the two monomers impact their homopolymerization mechanism: whereas alkyl-BTAs cooperatively self-assemble, oDMSi-BTAs self-assemble in an isodesmic manner. The effect of these mechanistic differences in the supramolecular copolymerization process is investigated as a function of the composition of the two monomers and explicitly rationalized by mathematical modeling. The results show that, at low fractions of oDMSi-BTA sergeants (25 mol%), the isodesmic assembly of the increasing amounts of sergeant becomes more dominant, and different species start to coexist in the copolymerization process. The analysis of the experimental data with a newly developed theoretical model allows us to quantify the thermodynamic parameters, the distribution of different species, and the compositions and stack lengths of the formed supramolecular copolymers existing at various feed ratios of the two monomers.This work was financially supported by The Netherlands Organization for Scientific Research (NWO-TOP PUNT Grant 10018944) and the Dutch Ministry of Education, Culture and Science (Gravity Program 024.001.035)