Fiber formation seen through the high-resolution computational microscope.

Supramolecular fibers draw widespread attention due to their role in biological systems and ability to form complex materials exhibiting rich and dynamic behavior. Although the information about the supramolecular structure is encoded in their molecular blocks, a complete understanding of how this information translates into the final structure requires detailed insights into the energy landscape of the process and the possible routes across this landscape. Here, we study the formation of 1,3,5-cyclohexanetricarboxamide fibers by a Markov state model of molecular dynamics simulations with the polarizable CHARMM Drude force-field.

Mechanistic studies of mycobacterial glycolipid biosynthesis by the mannosyltransferase PimE.

Tuberculosis (TB), a leading cause of death among infectious diseases globally, is caused by Mycobacterium tuberculosis (Mtb). The pathogenicity of Mtb is largely attributed to its complex cell envelope, which includes a class of glycolipids called phosphatidyl-myo-inositol mannosides (PIMs). These glycolipids maintain the integrity of the cell envelope, regulate permeability, and mediate host-pathogen interactions.

Structural insights into terminal arabinosylation of mycobacterial cell wall arabinan.

The global challenge of tuberculosis, caused by Mycobacterium tuberculosis (Mtb), is compounded by the emergence of drug-resistant strains. A critical factor in Mtb's pathogenicity is its intricate cell envelope, which acts as a formidable barrier against immune defences and pharmacological interventions. Central to this envelope are arabinogalactan (AG) and lipoarabinomannan (LAM), two complex polysaccharides containing arabinan domains essential for maintaining cell wall structure and function.

OLIVES: A Go̅-like Model for Stabilizing Protein Structure via Hydrogen Bonding Native Contacts in the Martini 3 Coarse-Grained Force Field.

Coarse-grained molecular dynamics simulations enable the modeling of increasingly complex systems at millisecond timescales. The transferable coarse-grained force field Martini 3 has shown great promise in modeling a wide range of biochemical processes, yet folded proteins in Martini 3 are not stable without the application of external bias potentials, such as elastic networks or Go̅-like models. We herein develop an algorithm, called OLIVES, which identifies native contacts with hydrogen bond capabilities in coarse-grained proteins and use it to implement a novel Go̅-like model for Martini 3.

Shocker─A Molecular Dynamics Protocol and Tool for Accelerating and Analyzing the Effects of Osmotic Shocks.

The process of osmosis, a fundamental phenomenon in life, drives water through a semipermeable membrane in response to a solute concentration gradient across this membrane. In vitro, osmotic shocks are often used to drive shape changes in lipid vesicles, for instance, to study fission events in the context of artificial cells. While experimental techniques provide a macroscopic picture of large-scale membrane remodeling processes, molecular dynamics (MD) simulations are a powerful tool to study membrane deformations at the molecular level.

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.

Potential binding modes of the gut bacterial metabolite, 5-hydroxyindole, to the intestinal L-type calcium channels and its impact on the microbiota in rats.

Intestinal microbiota and microbiota-derived metabolites play a key role in regulating the host physiology. Recently, we have identified a gut-bacterial metabolite, namely 5-hydroxyindole, as a potent stimulant of intestinal motility via its modulation of L-type voltage-gated calcium channels located on the intestinal smooth muscle cells. Dysregulation of L-type voltage-gated calcium channels is associated with various gastrointestinal motility disorders, including constipation, making L-type voltage-gated calcium channels an important target for drug development.

Synthetic Membrane Shaper for Controlled Liposome Deformation.

Shape defines the structure and function of cellular membranes. In cell division, the cell membrane deforms into a "dumbbell" shape, while organelles such as the autophagosome exhibit "stomatocyte" shapes. Bottom-up in vitro reconstitution of protein machineries that stabilize or resolve the membrane necks in such deformed liposome structures is of considerable interest to characterize their function.

Martini 3 Coarse-Grained Force Field for Carbohydrates.

The Martini 3 force field is a full reparametrization of the Martini coarse-grained model for biomolecular simulations. Due to the improved interaction balance, it allows for a more accurate description of condensed phase systems. In the present work, we develop a consistent strategy to parametrize carbohydrate molecules accurately within the framework of Martini 3.

Sequence-complementarity dependent co-assembly of phosphodiester-linked aromatic donor-acceptor trimers.

We have created sequenced phosphoester-linked trimers of aromatic donor/acceptors which participate in charge-transfer interactions. Each sequence displays characteristic self-assembly, and complementary sequences interact with each other to produce new nanostructures and thermochromism. This paves the way towards new functional nanomaterials which make bio-analogous use of sequence to tune structure.

Two cooperative binding sites sensitize PI(4,5)P recognition by the tubby domain.

Phosphoinositides (PIs) are lipid signaling molecules that operate by recruiting proteins to cellular membranes via PI recognition domains. The dominant PI of the plasma membrane is phosphatidylinositol 4,5-bisphosphate [PI(4,5)P]. One of only two PI(4,5)P recognition domains characterized in detail is the tubby domain.

Prokaryotic Na/H Exchangers-Transport Mechanism and Essential Residues.

Na/H exchangers are essential for Na and pH homeostasis in all organisms. Human Na/H exchangers are of high medical interest, and insights into their structure and function are aided by the investigation of prokaryotic homologues. Most prokaryotic Na/H exchangers belong to either the Cation/Proton Antiporter (CPA) superfamily, the Ion Transport (IT) superfamily, or the Na-translocating Mrp transporter superfamily.

Hypothesis-Driven, Structure-Based Design in Photopharmacology: The Case of eDHFR Inhibitors.

Photopharmacology uses light to regulate the biological activity of drugs. This precise control is obtained through the incorporation of molecular photoswitches into bioactive molecules. A major challenge for photopharmacology is the rational design of photoswitchable drugs that show light-induced activation.

Polyply; a python suite for facilitating simulations of macromolecules and nanomaterials.

Molecular dynamics simulations play an increasingly important role in the rational design of (nano)-materials and in the study of biomacromolecules. However, generating input files and realistic starting coordinates for these simulations is a major bottleneck, especially for high throughput protocols and for complex multi-component systems. To eliminate this bottleneck, we present the polyply software suite that provides 1) a multi-scale graph matching algorithm designed to generate parameters quickly and for arbitrarily complex polymeric topologies, and 2) a generic multi-scale random walk protocol capable of setting up complex systems efficiently and independent of the target force-field or model resolution.

Protein dynamics and lipid affinity of monomeric, zeaxanthin-binding LHCII in thylakoid membranes.

The xanthophyll cycle in the antenna of photosynthetic organisms under light stress is one of the most well-known processes in photosynthesis, but its role is not well understood. In the xanthophyll cycle, violaxanthin (Vio) is reversibly transformed to zeaxanthin (Zea) that occupies Vio binding sites of light-harvesting antenna proteins. Higher monomer/trimer ratios of the most abundant light-harvesting protein, the light-harvesting complex II (LHCII), usually occur in Zea accumulating membranes and have been observed in plants after prolonged illumination and during high-light acclimation.

Complex nanoemulsion for vitamin delivery: droplet organization and interaction with skin membranes.

Lipid nanoemulsions are promising nanomaterials for drug delivery applications in food, pharmaceutical and cosmetic industries. Despite the noteworthy commercial interest, little is known about their supramolecular organization, especially about how such multicomponent formulations interact with cell membranes. In the present work, coarse-grained molecular dynamics simulations have been employed to study the self-assembly of a 15-component lipid nanoemulsion droplet containing vitamins A and E for skin delivery.

How the Choice of Force-Field Affects the Stability and Self-Assembly Process of Supramolecular CTA Fibers.

In recent years, computational methods have become an essential element of studies focusing on the self-assembly process. Although they provide unique insights, they face challenges, from which two are the most often mentioned in the literature: the temporal and spatial scale of the self-assembly. A less often mentioned issue, but not less important, is the choice of the force-field.

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.

Competing Roles of Ca and Nonmuscle Myosin IIA on the Dynamics of the Metastasis-Associated Protein S100A4.

The calcium-binding protein S100A4 plays an important role in a wide range of biological processes such as cell motility, invasion, angiogenesis, survival, differentiation, contractility, and tumor metastasis and interacts with a range of partners. To understand the functional roles and interplay of S100A4 binding partners such as Ca and nonmuscle myosin IIA (NMIIA), we used molecular dynamics simulations to investigate apo S100A4 and four holo S100A4 structures: S100A4 bound to Ca, S100A4 bound to NMIIA, S100A4 bound to Ca and NMIIA, and a mutated S100A4 bound to Ca and NMIIA. Our results show that two competing factors, namely, Ca-induced activation and NMIIA-induced inhibition, modulate the dynamics of S100A4 in a competitive manner.

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.

Computational Approaches to Explore Bacterial Toxin Entry into the Host Cell.

Many bacteria secrete toxic protein complexes that modify and disrupt essential processes in the infected cell that can lead to cell death. To conduct their action, these toxins often need to cross the cell membrane and reach a specific substrate inside the cell. The investigation of these protein complexes is essential not only for understanding their biological functions but also for the rational design of targeted drug delivery vehicles that must navigate across the cell membrane to deliver their therapeutic payload.

Ceramide structure dictates glycosphingolipid nanodomain assembly and function.

Gangliosides in the outer leaflet of the plasma membrane of eukaryotic cells are essential for many cellular functions and pathogenic interactions. How gangliosides are dynamically organized and how they respond to ligand binding is poorly understood. Using fluorescence anisotropy imaging of synthetic, fluorescently labeled GM1 gangliosides incorporated into the plasma membrane of living cells, we found that GM1 with a fully saturated C16:0 acyl chain, but not with unsaturated C16:1 acyl chain, is actively clustered into nanodomains, which depends on membrane cholesterol, phosphatidylserine and actin.

Creasing of flexible membranes at vanishing tension.

The properties of freestanding tensionless interfaces and membranes at low bending rigidity κ are dominated by strong fluctuations and self-avoidance and are thus outside the range of standard perturbative analysis. We analyze this regime by a simple discretized, self-avoiding membrane model on a frame subject to periodic boundary conditions by use of Monte Carlo simulation and dynamically triangulated surface techniques. We find that at low bending rigidities, the membrane properties fall into three regimes: Below the collapse transition κ_{BP} it is subject to branched polymer instability where the framed surface is not defined, in a range below a threshold rigidity κ_{c} the conformational correlation function are characterized by power-law behavior with a continuously varying exponent α, 2<α≤4 and above κ_{c}, α=4 characteristic for linearized bending excitations.