Exploring the Ligand Binding and Conformational Dynamics of the Substrate-Binding Domain 1 of the ABC Transporter GlnPQ.

The adenosine triphosphate (ATP)-binding cassette (ABC) importer GlnPQ from Lactococcus lactis has two sequential covalently linked substrate-binding domains (SBDs), which capture the substrates and deliver them to the translocon. The two SBDs differ in their ligand specificities, binding affinities and the distance to the transmembrane domain; interestingly, both SBDs can bind their ligands simultaneously without affecting each other. In this work, we studied the binding of ligands to both SBDs using X-ray crystallography and molecular dynamics simulations. We report three high-resolution structures of SBD1, namely, the wild-type SBD1 with bound asparagine or arginine, and E184D SBD1 with glutamine bound. Molecular dynamics (MD) simulations provide a detailed insight into the dynamics associated with open-closed transitions of the SBDs.

Membrane manipulation by free fatty acids improves microbial plant polyphenol synthesis.

Microbial synthesis of nutraceutically and pharmaceutically interesting plant polyphenols represents a more environmentally friendly alternative to chemical synthesis or plant extraction. However, most polyphenols are cytotoxic for microorganisms as they are believed to negatively affect cell integrity and transport processes. To increase the production performance of engineered cell factories, strategies have to be developed to mitigate these detrimental effects. Here, we examine the accumulation of the stilbenoid resveratrol in the cell membrane and cell wall during its production using Corynebacterium glutamicum and uncover the membrane rigidifying effect of this stilbenoid experimentally and with molecular dynamics simulations. A screen of free fatty acid supplements identifies palmitelaidic acid and linoleic acid as suitable additives to attenuate resveratrol's cytotoxic effects resulting in a three-fold higher product titer. This cost-effective approach to counteract membrane-damaging effects of product accumulation is transferable to the microbial production of other polyphenols and may represent an engineering target for other membrane-active bioproducts.

Seeing the Invisibles: Detection of Peptide Enantiomers, Diastereomers, and Isobaric Ring Formation in Lanthipeptides Using Nanopores.

Mass spectrometry (MS) is widely used in proteomic analysis but cannot differentiate between molecules with the same mass-to-charge ratio. Nanopore technology might provide an alternative method for the rapid and cost-effective analysis and sequencing of proteins. In this study, we demonstrate that nanopore currents can distinguish between diastereomeric and enantiomeric differences in l- and d-peptides, not observed by conventional MS analysis, down to individual d-amino acids in small opioid peptides. Molecular dynamics simulations suggest that similar to chiral chromatography the resolution likely arises from multiple chiral interactions during peptide transport across the nanopore. Additionally, we used nanopore recordings to rapidly assess 4- and 11-amino acid ring formation in lanthipeptides, a process used in the synthesis of pharmaceutical peptides. The cyclization step requires distinguishing between constitutional isomers, which have identical MS signals and typically involve numerous tedious experiments to confirm. Hence, nanopore technology offers new possibilities for the rapid and cost-effective analysis of peptides, including those that cannot be easily differentiated by mass spectrometry.

Expulsion mechanism of the substrate-translocating subunit in ECF transporters.

Energy-coupling factor (ECF)-type transporters mediate the uptake of micronutrients in many bacteria. They consist of a substrate-translocating subunit (S-component) and an ATP-hydrolysing motor (ECF module) Previous data indicate that the S-component topples within the membrane to alternately expose the binding site to either side of the membrane. In many ECF transporters, the substrate-free S-component can be expelled from the ECF module. Here we study this enigmatic expulsion step by cryogenic electron microscopy and reveal that ATP induces a concave-to-convex shape change of two long helices in the motor, thereby destroying the S-component's docking site and allowing for its dissociation. We show that adaptation of the membrane morphology to the conformational state of the motor may favour expulsion of the substrate-free S-component when ATP is bound and docking of the substrate-loaded S-component after hydrolysis. Our work provides a picture of bilayer-assisted chemo-mechanical coupling in the transport cycle of ECF transporters.

Molecular insights into the effects of focused ultrasound mechanotherapy on lipid bilayers: Unlocking the keys to design effective treatments.

Administration of focused ultrasounds (US) represents an attractive complement to classical therapies for a wide range of maladies, from cancer to neurological pathologies, as they are non-invasive, easily targeted, their dosage is easy to control, and they involve low risks. Different mechanisms have been proposed for their activity but the direct effect of their interaction with cell membranes is not well understood at the molecular level. This is in part due to the difficulty of designing experiments able to probe the required spatio-temporal resolutions. Here we use Molecular Dynamics (MD) simulations at two resolution levels and machine learning (ML) classification tools to shed light on the effects that focused US mechanotherapy methods have over a range of lipid bilayers. Our results indicate that the dynamic-structural response of the membrane models to the mechanical perturbations caused by the sound waves strongly depends on the lipid composition. The analyses performed on the MD trajectories contribute to a better understanding of the behavior of lipid membranes, and to open up a path for the rational design of new therapies for the long list of diseases characterized by specific lipid profiles of pathological membrane cells.

Molecular dynamics simulation of an entire cell.

The ultimate microscope, directed at a cell, would reveal the dynamics of all the cell's components with atomic resolution. In contrast to their real-world counterparts, computational microscopes are currently on the brink of meeting this challenge. In this perspective, we show how an integrative approach can be employed to model an entire cell, the minimal cell, JCVI-syn3A, at full complexity. This step opens the way to interrogate the cell's spatio-temporal evolution with molecular dynamics simulations, an approach that can be extended to other cell types in the near future.

Quantification of Protein Glycosylation Using Nanopores.

Although nanopores can be used for single-molecule sequencing of nucleic acids using low-cost portable devices, the characterization of proteins and their modifications has yet to be established. Here, we show that hydrophilic or glycosylated peptides translocate too quickly across FraC nanopores to be recognized. However, high ionic strengths (i.e., 3 M LiCl) and low pH (i.e., pH 3) together with using a nanopore with a phenylalanine at its constriction allows the recognition of hydrophilic peptides, and to distinguish between mono- and diglycosylated peptides. Using these conditions, we devise a nanopore method to detect, characterize, and quantify post-translational modifications in generic proteins, which is one of the pressing challenges in proteomic analysis.

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. We benchmark quality and performance of the approach by creating realistic coordinates for polymer melt simulations, single-stranded as well as circular single-stranded DNA. We further demonstrate the power of our approach by setting up a microphase-separated block copolymer system, and by generating a liquid-liquid phase separated system inside a lipid vesicle.

Structure of the human signal peptidase complex reveals the determinants for signal peptide cleavage.

The signal peptidase complex (SPC) is an essential membrane complex in the endoplasmic reticulum (ER), where it removes signal peptides (SPs) from a large variety of secretory pre-proteins with exquisite specificity. Although the determinants of this process have been established empirically, the molecular details of SP recognition and removal remain elusive. Here, we show that the human SPC exists in two functional paralogs with distinct proteolytic subunits. We determined the atomic structures of both paralogs using electron cryo-microscopy and structural proteomics. The active site is formed by a catalytic triad and abuts the ER membrane, where a transmembrane window collectively formed by all subunits locally thins the bilayer. Molecular dynamics simulations indicate that this unique architecture generates specificity for SPs based on the length of their hydrophobic segments.

Piezo1 Forms Specific, Functionally Important Interactions with Phosphoinositides and Cholesterol.

Touch, hearing, and blood pressure regulation require mechanically gated ion channels that convert mechanical stimuli into electrical currents. One such channel is Piezo1, which plays a key role in the transduction of mechanical stimuli in humans and is implicated in diseases, such as xerocytosis and lymphatic dysplasia. There is building evidence that suggests Piezo1 can be regulated by the membrane environment, with the activity of the channel determined by the local concentration of lipids, such as cholesterol and phosphoinositides. To better understand the interaction of Piezo1 with its environment, we conduct simulations of the protein in a complex mammalian bilayer containing more than 60 different lipid types together with electrophysiology and mutagenesis experiments. We find that the protein alters its local membrane composition, enriching specific lipids and forming essential binding sites for phosphoinositides and cholesterol that are functionally relevant and often related to Piezo1-mediated pathologies. We also identify a number of key structural connections between the propeller and pore domains located close to lipid-binding sites.

Membrane mediated toppling mechanism of the folate energy coupling factor transporter.

Energy coupling factor (ECF) transporters are responsible for the uptake of micronutrients in bacteria and archaea. They consist of an integral membrane unit, the S-component, and a tripartite ECF module. It has been proposed that the S-component mediates the substrate transport by toppling over in the membrane when docking onto an ECF module. Here, we present multi-scale molecular dynamics simulations and in vitro experiments to study the molecular toppling mechanism of the S-component of a folate-specific ECF transporter. Simulations reveal a strong bending of the membrane around the ECF module that provides a driving force for toppling of the S-component. The stability of the toppled state depends on the presence of non-bilayer forming lipids, as confirmed by folate transport activity measurements. Together, our data provide evidence for a lipid-dependent toppling-based mechanism for the folate-specific ECF transporter, a mechanism that might apply to other ECF transporters.

Unidirectional rotating molecular motors dynamically interact with adsorbed proteins to direct the fate of mesenchymal stem cells.

Artificial rotary molecular motors convert energy into controlled motion and drive a system out of equilibrium with molecular precision. The molecular motion is harnessed to mediate the adsorbed protein layer and then ultimately to direct the fate of human bone marrow-derived mesenchymal stem cells (hBM-MSCs). When influenced by the rotary motion of light-driven molecular motors grafted on surfaces, the adsorbed protein layer primes hBM-MSCs to differentiate into osteoblasts, while without rotation, multipotency is better maintained. We have shown that the signaling effects of the molecular motion are mediated by the adsorbed cell-instructing protein layer, influencing the focal adhesion-cytoskeleton actin transduction pathway and regulating the protein and gene expression of hBM-MSCs. This unique molecular-based platform paves the way for implementation of dynamic interfaces for stem cell control and provides an opportunity for novel dynamic biomaterial engineering for clinical applications.

Two distinct anionic phospholipid-dependent events involved in SecA-mediated protein translocation.

Protein translocation across the bacterial cytoplasmic membrane is an essential process catalyzed by the Sec translocase, which in its minimal form consists of the protein-conducting channel SecYEG, and the motor ATPase SecA. SecA binds via its positively charged N-terminus to membranes containing anionic phospholipids, leading to a lipid-bound intermediate. This interaction induces a conformational change in SecA, resulting in a high-affinity association with SecYEG, which initiates protein translocation. Here, we examined the effect of anionic lipids on the SecA-SecYEG interaction in more detail, and discovered a second, yet unknown, anionic lipid-dependent event that stimulates protein translocation. Based on molecular dynamics simulations we identified an anionic lipid-enriched region in vicinity of the lateral gate of SecY. Here, the anionic lipid headgroup accesses the lateral gate, thereby stabilizing the pre-open state of the channel. The simulations suggest flip-flop movement of phospholipid along the lateral gate. Electrostatic contribution of the anionic phospholipids at the lateral gate may directly stabilize positively charged residues of the signal sequence of an incoming preprotein. Such a mechanism allows for the correct positioning of the entrant peptide, thereby providing a long-sought explanation for the role of anionic lipids in signal sequence folding during protein translocation.

Binding of quinazolinones to c-KIT G-quadruplex; an interplay between hydrogen bonding and π-π stacking.

Stabilization of G-quadruplex structures in the c-KIT promoter with the aid of ligands has become an area of great interest in potential cancer therapeutics. Understanding the binding process between ligands and G-quadruplex is essential for a discovery of selective ligands with high binding affinity to G-quadruplex. In the present work, binding mechanisms of 4-quinazolinones to c-KIT G-quadruplex were investigated theoretically by means of molecular dynamics (MD) simulations. To explore the binding affinity of ligands, binding free energy calculations were performed using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. We demonstrate that the key interactions in G-quadruplex-ligand complexes are π-π stacking and hydrogen bond interactions. However, neither of these two interactions alone determines the stability of the G-quadruplex-ligand complexes; rather, it is the result of an intricate interplay between the two. To further examine the nature of the binding, a free energy decomposition analysis at residue level was carried out. The results clearly demonstrate the crucial roles of two hot spot residues (DG4 and DG8) for the binding of ligands to c-KIT G-quadruplex, and highlight the importance of the planar aromatic moiety of ligands in G-quadruplex stabilization via π-π stacking interactions. Our study can assist in the design of new derivatives of 4-quinazolinone with high binding affinity for c-KIT G-quadruplex.

Ceramides bind VDAC2 to trigger mitochondrial apoptosis.

Ceramides draw wide attention as tumor suppressor lipids that act directly on mitochondria to trigger apoptotic cell death. However, molecular details of the underlying mechanism are largely unknown. Using a photoactivatable ceramide probe, we here identify the voltage-dependent anion channels VDAC1 and VDAC2 as mitochondrial ceramide binding proteins. Coarse-grain molecular dynamics simulations reveal that both channels harbor a ceramide binding site on one side of the barrel wall. This site includes a membrane-buried glutamate that mediates direct contact with the ceramide head group. Substitution or chemical modification of this residue abolishes photolabeling of both channels with the ceramide probe. Unlike VDAC1 removal, loss of VDAC2 or replacing its membrane-facing glutamate with glutamine renders human colon cancer cells largely resistant to ceramide-induced apoptosis. Collectively, our data support a role of VDAC2 as direct effector of ceramide-mediated cell death, providing a molecular framework for how ceramides exert their anti-neoplastic activity.

Transferable MARTINI Model of Poly(ethylene Oxide).

Motivated by the deficiencies of the previous MARTINI models of poly(ethylene oxide) (PEO), we present a new model featuring a high degree of transferability. The model is parametrized on (a) a set of 8 free energies of transfer of dimethoxyethane (PEO dimer) from water to solvents of varying polarity; (b) the radius of gyration in water at high dilution; and (c) matching angle and dihedral distributions from atomistic simulations. We demonstrate that our model behaves well in five different areas of application: (1) it produces accurate densities and phase behavior or small PEO oligomers and water mixtures; (2) it yields chain dimensions in good agreement with the experiment in three different solvents (water, diglyme, and benzene) over a broad range of molecular weights (∼1.2 kg/mol to 21 kg/mol); (3) it reproduces qualitatively the structural features of lipid bilayers containing PEGylated lipids in the brush and mushroom regime; (4) it is able to reproduce the phase behavior of several PEO-based nonionic surfactants in water; and (5) it can be combined with the existing MARTINI PS to model PS-PEO block copolymers. Overall, the new PEO model outperforms previous models and features a high degree of transferability.

Energetics Underlying Twist Polymorphisms in Amyloid Fibrils.

Amyloid fibrils are highly ordered protein aggregates associated with more than 40 human diseases. The exact conditions under which the fibrils are grown determine many types of reported fibril polymorphism, including different twist patterns. Twist-based polymorphs display unique mechanical properties in vitro, and the relevance of twist polymorphism in amyloid diseases has been suggested. We present transmission electron microscopy images of Aß42-derived (amyloid ß) fibrils, which are associated with Alzheimer's disease, demonstrating the presence of twist variability even within a single long fibril. To better understand the molecular underpinnings of twist polymorphism, we present a structural and thermodynamics analysis of molecular dynamics simulations of the twisting of ß-sheet protofilaments of a well-characterized cross-ß model: the GNNQQNY peptide from the yeast prion Sup35. The results show that a protofilament model of GNNQQNY is able to adopt twist angles from -11° on the left-hand side to +8° on the right-hand side in response to various external conditions, keeping an unchanged peptide structure. The potential of mean force (PMF) of this cross-ß structure upon twisting revealed that only ∼2kBT per peptide are needed to stabilize a straight conformation with respect to the left-handed free-energy minimum. The PMF also shows that the canonical structural core of ß-sheets, i.e., the hydrogen-bonded backbone ß-strands, favors the straight conformation. However, the concerted effects of the side chains contribute to twisting, which provides a rationale to correlate polypeptide sequence, environmental growth conditions and number of protofilaments in a fibril with twist polymorphisms.

Structural insights into K48-linked ubiquitin chain formation by the Pex4p-Pex22p complex.

Pex4p is a peroxisomal E2 involved in ubiquitinating the conserved cysteine residue of the cycling receptor protein Pex5p. Previously, we demonstrated that Pex4p from the yeast Saccharomyces cerevisiae binds directly to the peroxisomal membrane protein Pex22p and that this interaction is vital for receptor ubiquitination. In addition, Pex22p binding allows Pex4p to specifically produce lysine 48 linked ubiquitin chains in vitro through an unknown mechanism. This activity is likely to play a role in targeting peroxisomal proteins for proteasomal degradation. Here we present the crystal structures of Pex4p alone and in complex with Pex22p from the yeast Hansenula polymorpha. Comparison of the two structures demonstrates significant differences to the active site of Pex4p upon Pex22p binding while molecular dynamics simulations suggest that Pex22p binding facilitates active site remodelling of Pex4p through an allosteric mechanism. Taken together, our data provide insights into how Pex22p binding allows Pex4p to build K48-linked Ub chains.

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.

Ganglioside-Lipid and Ganglioside-Protein Interactions Revealed by Coarse-Grained and Atomistic Molecular Dynamics Simulations.

Gangliosides are glycolipids in which an oligosaccharide headgroup containing one or more sialic acids is connected to a ceramide. Gangliosides reside in the outer leaflet of the plasma membrane and play a crucial role in various physiological processes such as cell signal transduction and neuronal differentiation by modulating structures and functions of membrane proteins. Because the detailed behavior of gangliosides and protein-ganglioside interactions are poorly known, we investigated the interactions between the gangliosides GM1 and GM3 and the proteins aquaporin (AQP1) and WALP23 using equilibrium molecular dynamics simulations and potential of mean force calculations at both coarse-grained (CG) and atomistic levels. In atomistic simulations, on the basis of the GROMOS force field, ganglioside aggregation appears to be a result of the balance between hydrogen bond interactions and steric hindrance of the headgroups. GM3 clusters are slightly larger and more ordered than GM1 clusters due to the smaller headgroup of GM3. The different structures of GM1 and GM3 clusters from atomistic simulations are not observed at the CG level based on the Martini model, implying a difference in driving forces for ganglioside interactions in atomistic and CG simulations. For protein-ganglioside interactions, in the atomistic simulations, GM1 lipids bind to specific sites on the AQP1 surface, whereas they are depleted from WALP23. In the CG simulations, the ganglioside binding sites on the AQP1 surface are similar, but ganglioside aggregation and protein-ganglioside interactions are more prevalent than in the atomistic simulations. Using the polarizable Martini water model, results were closer to the atomistic simulations. Although experimental data for validation is lacking, we proposed modified Martini parameters for gangliosides to more closely mimic the sizes and structures of ganglioside clusters observed at the atomistic level.