CHARMM-GUI Martini Maker for modeling and simulation of complex bacterial membranes with lipopolysaccharides

Creating Coarse-Grained Systems with COBY: Toward Higher Accuracy of Complex Biological Systems

Current trends in molecular modeling are geared toward increasingly realistic representations of the biological environments reflected in larger, more complex systems. The complexity of the system-building procedure is ideally handled by software that converts user-provided descriptors into system coordinates. This, however, is not a trivial task, as building algorithms use simplifications that result in inaccuracies in the system properties. We created COBY, a coarse-grained system builder that can create a large variety of systems in a single command call with an improved accuracy of the complex membrane and solvent building procedures. COBY also contains features for building diverse systems in a single step, and has functionalities aiding force field development. COBY is an open-source software written in Python 3, and the code, documentation, and tutorials are hosted at https://github.com/MikkelDA/COBY.

Increased burden of rare risk variants across gene expression networks predisposes to sporadic Parkinson's disease

Alpha-synuclein (αSyn) is an intrinsically disordered protein that accumulates in the brains of patients with Parkinson's disease (PD). Through a high-throughput screen, we recently identified 38 genes whose knockdown modulates αSyn propagation. Here, we show that, among those, TAX1BP1 regulates how αSyn interacts with lipids, and ADAMTS19 modulates how αSyn phase separates into inclusions, adding to the growing body of evidence implicating those processes in PD. Through RNA sequencing, we identify several genes that are differentially expressed after knockdown of TAX1BP1 or ADAMTS19 and carry an increased frequency of rare risk variants in patients with PD versus healthy controls. Those differentially expressed genes cluster within modules in regions of the brain that develop high degrees of αSyn pathology. We propose a model for the genetic architecture of sporadic PD: increased burden of risk variants across genetic networks dysregulates pathways underlying αSyn homeostasis and leads to pathology and neurodegeneration.

Computational Methods for Modeling Lipid-Mediated Active Pharmaceutical Ingredient Delivery

Lipid-mediated delivery of active pharmaceutical ingredients (API) opened new possibilities in advanced therapies. By encapsulating an API into a lipid nanocarrier (LNC), one can safely deliver APIs not soluble in water, those with otherwise strong adverse effects, or very fragile ones such as nucleic acids. However, for the rational design of LNCs, a detailed understanding of the composition-structure-function relationships is missing. This review presents currently available computational methods for LNC investigation, screening, and design. The state-of-the-art physics-based approaches are described, with the focus on molecular dynamics simulations in all-atom and coarse-grained resolution. Their strengths and weaknesses are discussed, highlighting the aspects necessary for obtaining reliable results in the simulations. Furthermore, a machine learning, i.e., data-based learning, approach to the design of lipid-mediated API delivery is introduced. The data produced by the experimental and theoretical approaches provide valuable insights. Processing these data can help optimize the design of LNCs for better performance. In the final section of this Review, state-of-the-art of computer simulations of LNCs are reviewed, specifically addressing the compatibility of experimental and computational insights.

Systematic Approach to Parametrization of Disaccharides for the Martini 3 Coarse-Grained Force Field

Sugars are ubiquitous in biology; they occur in all kingdoms of life. Despite their prevalence, they have often been somewhat neglected in studies of structure-dynamics-function relationships of macromolecules to which they are attached, with the exception of nucleic acids. This is largely due to the inherent difficulties of not only studying the conformational dynamics of sugars using experimental methods but indeed also resolving their static structures. Molecular dynamics (MD) simulations offer a route to the prediction of conformational ensembles and the time-dependent behavior of sugars and glycosylated macromolecules. However, at the all-atom level of detail, MD simulations are often too computationally demanding to allow a systematic investigation of molecular interactions in systems of interest. To overcome this, large scale simulations of complex biological systems have profited from advances in coarse-grained (CG) simulations. Perhaps the most widely used CG force field for biomolecular simulations is Martini. Here, we present a parameter set for glucose- and mannose-based disaccharides for Martini 3. The generation of the CG parameters from atomistic trajectories is automated as fully as possible, and where not possible, we provide details of the protocol used for manual intervention.

Free Energy, Rates, and Mechanism of Transmembrane Dimerization in Lipid Bilayers from Dynamically Unbiased Molecular Dynamics Simulations

The assembly of proteins in membranes plays a key role in many crucial cellular pathways. Despite their importance, characterizing transmembrane assembly remains challenging for experiments and simulations. Equilibrium molecular dynamics simulations do not cover the time scales required to sample the typical transmembrane assembly. Hence, most studies rely on enhanced sampling schemes that steer the dynamics of transmembrane proteins along a collective variable that should encode all slow degrees of freedom. However, given the complexity of the condensed-phase lipid environment, this is far from trivial, with the consequence that free energy profiles of dimerization can be poorly converged. Here, we introduce an alternative approach, which relies only on simulating short, dynamically unbiased paths, avoiding using collective variables or biasing forces. By merging all paths, we obtain free energy profiles, rates, and mechanisms of transmembrane dimerization with the same set of simulations. We showcase our algorithm by sampling the spontaneous association and dissociation of a transmembrane protein in a lipid bilayer, the popular coarse-grained Martini force field. Our algorithm represents a promising way to investigate assembly processes in biologically relevant membranes, overcoming some of the challenges of conventional methods.

Effects of Lipid Headgroups on the Mechanical Properties and In Vitro Cellular Internalization of Liposomes

We performed all-atom and coarse-grained simulations of lipid bilayer mixtures of the unsaturated lipid DOPC, with saturated lipids having the same 18-carbon acyl tails and different headgroups, to understand their mechanical properties. The secondary lipids were DSPG, DSPA, DSPS, DSPC, and DSPE. The DOPC:DSPG system with 65:35 molar ratio was the softest, with area compressibility modulus KA ∼ 22% smaller than the pure DOPC value. Raising the mole % of DOPC leads to increases in KA, yet at any given composition the KA trend is DSPG < DSPA < DSPS < DSPC < DSPE. Lipid-lipid interactions are weaker in DOPC:DSPG mixtures and stronger in DSPE systems. The head and phosphate groups of the secondary lipids DSPG, DSPA, and DSPS interact strongly with salt ions. Adding secondary lipids leads to DOPC having more ordered acyl tails relative to pure DOPC systems. No evidence of phase separation or inhomogeneities was observed in our simulations. We synthesized three liposomal formulations, L-DOPC (pure DOPC) and L-DOPC/DSPG and L-DOPC/DSPA, both with 15 mol % of secondary lipid. L-DOPC/DSPA had approximately 3- and 2-times higher in vitro internalization by normal epithelial (EpH4-Ev) and metastatic breast cancer (4T1) cells, compared with L-DOPC. The uptake of L-DOPC/DSPG by EpH4-Ev cells was almost 2-fold compared to L-DOPC, but both liposomes had similar uptakes by cancer cells. As L-DOPC/DSPG and L-DOPC/DSPA have similar KA values, we presumed that the mechanical properties, possibly in combination with the higher negative surface charges in L-DOPC/DSPA and differences in effective liposome diameters and diffusivities, contributed to these observations.

How well do empirical molecular mechanics force fields model the cholesterol condensing effect?

Membrane properties are determined in part by lipid composition, and cholesterol plays a large role in determining these properties. Cellular membranes show a diverse range of cholesterol compositions, the effects of which include alterations to cellular biomechanics, lipid raft formation, membrane fusion, signaling pathways, metabolism, pharmaceutical therapeutic efficacy, and disease onset. In addition, cholesterol plays an important role in non-cellular membranes, with its concentration in the skin lipid matrix being implicated in several skin diseases. In phospholipid membranes, cholesterol increases the tail ordering of neighboring lipids, decreasing the membrane lateral area and increasing the thickness. This reduction in the lateral area, known as the cholesterol condensing effect, results from cholesterol-lipid mixtures deviating from ideal mixing. Capturing the cholesterol condensing effect is crucial for molecular dynamics simulations as it directly affects the accuracy of predicted membrane properties, which are essential for understanding membrane function. We present a comparative analysis of cholesterol models across several popular force fields: CHARMM36, Slipids, Lipid17, GROMOS 53A6L, GROMOS-CKP, MARTINI 2, MARTINI 3, and ELBA. The simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) membranes with varying cholesterol concentrations were conducted to calculate the partial-molecular areas of cholesterol and other condensing parameters, which are compared to the experimental data for validation. While all tested force fields predict small negative deviations from ideal mixing in cholesterol-DOPC membranes, only all-atom force fields capture the larger deviations expected in DMPC membranes. United-atom and coarse-grained models under-predict this effect, condensing fewer neighboring lipids by smaller magnitudes, resulting in too small deviations from ideal mixing. These results suggest that all-atom force fields, particularly CHARMM36 or Slipids, should be used for accurate simulations of cholesterol-containing membranes.

Phospholipid Transport Across the Bacterial Periplasm Through the Envelope-spanning Bridge YhdP

The outer membrane of Gram-negative bacteria provides a formidable barrier, essential for both pathogenesis and antimicrobial resistance. Biogenesis of this complex structure necessitates the transport of phospholipids across the cell envelope. Recently, YhdP was implicated as a major protagonist in the trafficking of inner membrane phospholipids to the outer membrane; however the molecular mechanism of YhdP mediated transport remains elusive. Here, utilising AlphaFold, we observe YhdP to form an elongated assembly of 60 β-strands that curve to form a continuous hydrophobic groove. This architecture is consistent with our negative stain electron microscopy data which reveals YhdP to be approximately 250 Å in length and thus sufficient to span the bacterial cell envelope. Furthermore, molecular dynamics simulations and bacterial growth assays indicate essential helical regions at the N- and C-termini of YhdP, that may embed into the inner and outer membranes respectively, reinforcing its envelope spanning nature. Our in vivo crosslinking data reveal phosphate-containing substrates captured along the length of the YhdP groove, providing direct evidence that YhdP interacts with a phosphate-containing substrate, which we propose to be phospholipids. This finding is congruent with our molecular dynamics simulations which demonstrate the propensity for inner membrane lipids to spontaneously enter the groove of YhdP. Collectively, our results support a model in which YhdP bridges the cell envelope, providing a hydrophobic environment for the transport of phospholipids to the outer membrane.

An in silico framework to visualize how cancer-associated mutations influence structural plasticity of the chemokine receptor CCR3

G protein Coupled Receptors (GPCRs) are the largest family of cell surface receptors in humans. Somatic mutations in GPCRs are implicated in cancer progression and metastasis, but mechanisms are poorly understood. Emerging evidence implicates perturbation of intra-receptor activation pathway motifs whereby extracellular signals are transmitted intracellularly. Recently, sufficiently sensitive methodology was described to calculate structural strain as a function of missense mutations in AlphaFold-predicted model structures, which was extensively validated on experimental and predicted structural datasets. When paired with Molecular Dynamics (MD) simulations, these tools provide a facile approach to screen mutations in silico. We applied this framework to calculate the structural and dynamic effects of cancer-associated mutations in the chemokine receptor CCR3, a Class A GPCR involved in cancer and autoimmune disorders. Residue-residue contact scoring refined effective strain results, highlighting significant remodeling of inter- and intra-motif contacts along the highly conserved GPCR activation pathway network. We then integrated AlphaFold-derived predicted Local Distance Difference Test scores with per-residue Root Mean Square Fluctuations and activation pathway Contact Analysis (CONAN) from coarse grain MD simulations to identify statistically significant changes in receptor dynamics upon mutation. Finally, analysis of negative control mutants suggests false positive results in AlphaFold pipelines should be considered but can be mitigated with stricter control of statistical analysis. Our results indicate selected mutants influence structural plasticity of CCR3 related to ligand interaction, activation, and G protein coupling, using a framework that could be applicable to a wide range of biochemically relevant protein targets following further validation.

Targeting the Weak Spot: Preferential Disruption of Bacterial Poles by Janus Nanoparticles

The interaction between nanoparticles (NPs) and bacterial cell envelopes is crucial for designing effective antibacterial materials against multi-drug-resistant pathogens. However, the current understanding assumes a uniform bacterial cell wall. This study challenges that assumption by investigating how bacterial cell wall curvature impacts antibacterial NP action. Focusing on Janus NPs, which feature segregated hydrophobic and polycationic ligands and previously demonstrated high efficacy against diverse bacteria, we found that these NPs preferentially target and disrupt bacterial poles. Experimental and computational approaches reveal that curvature at E. coli poles induces conformational changes in lipopolysaccharide (LPS) polymers on the outer membrane, exposing underlying lipids for NP-mediated disruption. We establish that curvature-induced targeting by Janus NPs depends on the outer membrane composition and is most pronounced at physiologically relevant LPS densities. This work demonstrates that high-curvature regions of bacterial cell walls are "weak spots" for Janus NPs, thereby aiding the development of more effective targeted therapies.

VPS4A is the selective receptor for lipophagy in mice and humans

Lipophagy is a ubiquitous mechanism for degradation of lipid droplets (LDs) in lysosomes. Autophagy receptors selectively target organelles for lysosomal degradation. The selective receptor for lipophagy remains elusive. Using mouse liver phosphoproteomics and human liver transcriptomics, we identify vacuolar-protein-sorting-associated protein 4A (VPS4A), a member of a large family AAA+ ATPases, as a selective receptor for lipophagy. We show that phosphorylation of VPS4A on Ser95,97 and its localization to LDs in response to fasting drives lipophagy. Imaging/three-dimensional (3D) reconstruction and biochemical analyses reveal the concomitant degradation of VPS4A and LDs in lysosomes in an autophagy-gene-7-sensitive manner. Either silencing VPS4A or targeting VPS4AS95,S97 phosphorylation or VPS4A binding to LDs or LC3 blocks lipophagy without affecting other forms of selective autophagy. Finally, VPS4A levels and markers of lipophagy are markedly reduced in human steatotic livers-revealing a fundamental role of VPS4A as the lipophagy receptor in mice and humans.

Local Stress in Cylindrically Curved Lipid Membrane: Insights into Local Versus Global Lateral Fluidity Models

Lipid membranes, which are fundamental to cellular function, undergo various mechanical deformations. Accurate modeling of these processes necessitates a thorough understanding of membrane elasticity. The lateral shear modulus, a critical parameter describing membrane resistance to lateral stresses, remains elusive due to the membrane's fluid nature. Two contrasting hypotheses, local fluidity and global fluidity, have been proposed. While the former suggests a zero local lateral shear modulus anywhere within lipid monolayers, the latter posits that only the integral of this modulus over the monolayer thickness vanishes. These differing models lead to distinct estimations of other elastic moduli and affect the modeling of biological processes, such as membrane fusion/fission and membrane-mediated interactions. Notably, they predict distinct local stress distributions in cylindrically curved membranes. The local fluidity model proposes isotropic local lateral stress, whereas the global fluidity model predicts anisotropy due to anisotropic local lateral stretching of lipid monolayers. Using molecular dynamics simulations, this study directly investigates these models by analyzing local stress in a cylindrically curved membrane. The results conclusively demonstrate the existence of static local lateral shear stress and anisotropy in local lateral stress within the monolayers of the cylindrical membrane, strongly supporting the global fluidity model. These findings have significant implications for the calculation of surface elastic moduli and offer novel insights into the fundamental principles governing lipid membrane elasticity.

Origins of synergy in multilipid lubrication

Lipid bilayers, ubiquitous in living systems, form lubricious boundary layers in aqueous media, with broad relevance for biolubrication, especially in mechanically stressed environments such as articular cartilage in joints, as well as for modifying material interfacial properties. Model studies have revealed efficient lubricity by single-component lipid bilayers; synovial joints, however (e.g. hips and knees), comprise over a hundred different lipids, raising the question of whether this is natural redundancy or whether it confers any lubrication benefits. Here, we examine lubrication by progressively more complex mixtures of lipids representative of those in joints, using a surface forces balance at physiologically relevant salt concentrations and pressures. We find that different combinations of such lipids differ very significantly in the robustness of the bilayers to hemifusion under physiological loads (when lubrication breaks down), indicating a clear lubrication synergy afforded by multiple lipid types in the bilayers. Insight into the origins of this synergy is provided by detailed molecular dynamics simulations of potential profiles for the formation of stalks, the essential precursors of hemifusion, between bilayers of the different lipid mixtures used in the experiments. These reveal how bilayer hemifusion-and thus lubrication breakdown-depends on the detailed lipid bilayer composition, through the corresponding separation into domains that are better able to resist stalk formation. Our results shed light on the role of lipid-type proliferation in biolubrication synergy, point to improved treatment modalities for common joint diseases such as osteoarthritis, and indicate how lipid-based interfacial modification in a materials context may be optimized.

Lipid organization by the Caveolin-1 complex

Caveolins are lipid-binding proteins that can organize membrane remodeling and oligomerize into the 8S complex. The CAV1-8S complex comprises a disk-like structure, about 15 nm in diameter, with a central beta barrel. Further oligomerization of 8S complexes remodels the membrane into caveolae vessels, with a dependence on cholesterol concentration. However, the molecular mechanisms behind membrane remodeling and cholesterol filtering are still not understood. Performing atomistic molecular dynamics simulations in combination with advanced sampling techniques, we describe how the CAV1-8S complex bends the membrane and accumulates cholesterol. Here, our simulations show an enhancing effect by the palmitoylations of CAV1, and we predict that the CAV1-8S complex can extract cholesterol molecules from the lipid bilayer and accommodate them in its beta barrel. Through backmapping to the all-atom level, we also conclude that the Martini v.2 coarse-grained force field overestimates membrane bending, as the atomistic simulations exhibit only very localized bending.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Polymyxin B1 in the Escherichia coli inner membrane: A complex story of protein and lipopolysaccharide-mediated insertion

The rise in multi-drug resistant Gram-negative bacterial infections has led to an increased need for "last-resort" antibiotics such as polymyxins. However, the emergence of polymyxin-resistant strains threatens to bring about a post-antibiotic era. Thus, there is a need to develop new polymyxin-based antibiotics, but a lack of knowledge of the mechanism of action of polymyxins hinders such efforts. It has recently been suggested that polymyxins induce cell lysis of the Gram-negative bacterial inner membrane (IM) by targeting trace amounts of lipopolysaccharide (LPS) localized there. We use multiscale molecular dynamics (MD), including long-timescale coarse-grained (CG) and all-atom (AA) simulations, to investigate the interactions of polymyxin B1 (PMB1) with bacterial IM models containing phospholipids (PLs), small quantities of LPS, and IM proteins. LPS was observed to (transiently) phase separate from PLs at multiple LPS concentrations, and associate with proteins in the IM. PMB1 spontaneously inserted into the IM and localized at the LPS-PL interface, where it cross-linked lipid headgroups via hydrogen bonds, sampling a wide range of interfacial environments. In the presence of membrane proteins, a small number of PMB1 molecules formed interactions with them, in a manner that was modulated by local LPS molecules. Electroporation-driven translocation of PMB1 via water-filled pores was favored at the protein-PL interface, supporting the 'destabilizing' role proteins may have within the IM. Overall, this in-depth characterization of PMB1 modes of interaction reveals how small amounts of mislocalized LPS may play a role in pre-lytic targeting and provides insights that may facilitate rational improvement of polymyxin-based antibiotics.

Increased burden of rare risk variants across gene expression networks predisposes to sporadic Parkinson's disease

Alpha-synuclein (αSyn) is an intrinsically disordered protein that accumulates in the brains of patients with Parkinson's disease and forms intraneuronal inclusions called Lewy Bodies. While the mechanism underlying the dysregulation of αSyn in Parkinson's disease is unclear, it is thought that prionoid cell-to-cell propagation of αSyn has an important role. Through a high throughput screen, we recently identified 38 genes whose knock down modulates αSyn propagation. Follow up experiments were undertaken for two of those genes, TAX1BP1 and ADAMTS19, to study the mechanism with which they regulate αSyn homeostasis. We used a recently developed M17D neuroblastoma cell line expressing triple mutant (E35K+E46K+E61K) "3K" αSyn under doxycycline induction. 3K αSyn spontaneously forms inclusions that show ultrastructural similarities to Lewy Bodies. Experiments using that cell line showed that TAX1BP1 and ADAMTS19 regulate how αSyn interacts with lipids and phase separates into inclusions, respectively, adding to the growing body of evidence implicating those processes in Parkinson's disease. Through RNA sequencing, we identified several genes that are differentially expressed after knock-down of TAX1BP1 or ADAMTS19. Burden analysis revealed that those differentially expressed genes (DEGs) carry an increased frequency of rare risk variants in Parkinson's disease patients versus healthy controls, an effect that was independently replicated across two separate cohorts (GP2 and AMP-PD). Weighted gene co-expression network analysis (WGCNA) showed that the DEGs cluster within modules in regions of the brain that develop high degrees of αSyn pathology (basal ganglia, cortex). We propose a novel model for the genetic architecture of sporadic Parkinson's disease: increased burden of risk variants across genetic networks dysregulates pathways underlying αSyn homeostasis, thereby leading to pathology and neurodegeneration.

Faster but Not Sweeter: A Model of Escherichia coli Re-level Lipopolysaccharide for Martini 3 and a Martini 2 Version with Accelerated Kinetics

Lipopolysaccharide (LPS) is a complex glycolipid molecule that is the main lipidic component of the outer leaflet of the outer membrane of Gram-negative bacteria. It has very limited lateral motion compared to phospholipids, which are more ubiquitous in biological membranes, including in the inner leaflet of the outer membrane of Gram-negative bacteria. The slow-moving nature of LPS can present a hurdle for molecular dynamics simulations, given that the (pragmatically) accessible timescales to simulations are currently limited to microseconds, during which LPS displays some conformational dynamics but hardly any lateral diffusion. Thus, it is not feasible to observe phenomena such as insertion of molecules, including antibiotics/antimicrobials, directly into the outer membrane from the extracellular side nor to observe LPS dissociating from proteins via molecular dynamics using currently available models at the atomistic and more coarse-grained levels of granularity. Here, we present a model of deep rough LPS compatible with the Martini 2 coarse-grained force field with scaled down nonbonded interactions to enable faster diffusion. We show that the faster-diffusing LPS model is able to reproduce the salient biophysical properties of the standard models, but due to its faster lateral motion, molecules are able to penetrate deeper into membranes containing the faster model. We show that the fast ReLPS model is able to reproduce experimentally determined patterns of interaction with outer membrane proteins while also allowing for LPS to associate and dissociate with proteins within microsecond timescales. We also complete the Martini 3 LPS toolkit for Escherichia coli by presenting a (standard) model of deep rough LPS for this force field.

Chain-Length Dependence of Peptide-Lipid Bilayer Interaction Strength and Binding Kinetics: A Combined Theoretical and Experimental Approach

Physical interactions between polypeptide chains and lipid membranes underlie critical cellular processes. Yet, despite fundamental importance, key mechanistic aspects of these interactions remain elusive. Bulk experiments have revealed a linear relationship between free energy and peptide chain length in a model system, but does this linearity extend to the interaction strength and to the kinetics of lipid binding? To address these questions, we utilized a combination of coarse-grained molecular dynamics (CG MD) simulations, analytical modeling, and atomic force microscopy (AFM)-based single molecule force spectroscopy. Following previous bulk experiments, we focused on interactions between short hydrophobic peptides (WLn, n = 1, ..., 5) with 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayers, a simple system that probes peptide primary structure effects. Potentials of mean force extracted from CG MD recapitulated the linearity of free energy with the chain length. Simulation results were quantitatively connected to bulk biochemical experiments via a single scaling factor of order unity, corroborating the methodology. Additionally, CG MD revealed an increase in the distance to the transition state, a result that weakens the dependence of the dissociation force on the peptide chain length. AFM experiments elucidated rupture force distributions and, through modeling, intrinsic dissociation rates. Taken together, the analysis indicates a rupture force plateau in the WLn-POPC system, suggesting that the final rupture event involves the last 2 or 3 residues. In contrast, the linear dependence on chain length was preserved in the intrinsic dissociation rate. This study advances the understanding of peptide-lipid interactions and provides potentially useful insights for the design of peptides with tailored membrane-interacting properties.

Lipid Organization by the Caveolin-1 Complex

Caveolins are lipid-binding proteins that can organize membrane remodeling and oligomerize into the 8S-complex. The CAV1 8S-complex comprises a disk-like structure, about 15nm in diameter, with a central beta barrel. Further oligomerization of 8S-complexes remodels the membrane into caveolae vessels, with a dependence on cholesterol concentration. However, the molecular mechanisms behind membrane remodeling and cholesterol filtering are still not understood. Performing atomistic Molecular Dynamics simulations in combination with advanced sampling techniques, we describe how the CAV1-8S complex bends the membrane and accumulates cholesterol. Here, our simulations show an enhancing effect by the palmitoylations of CAV1, and we predict that the CAV1-8S complex can extract cholesterol molecules from the lipid bilayer and accommodate them in its beta barrel. Through backmapping to the all-atom level we also conclude that the Martini v2 coarse-grained forcefield overestimates membrane bending, as the atomistic simulations exhibit only very localized bending.

A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory complex I

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations aimed at investigating new CI variants. One of the primary CI variants associated with the Leber hereditary optic neuropathy (m.14484T>C/MT-ND6) was used as a test case and was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant positioned in the E-channel region, which is fundamental for CI function, stiffens the enzyme dynamics. Moreover, a new mechanism for the transition between π- and α-conformation in the helix carrying the primary variant is proposed. This may have implications for the E-channel opening/closing mechanism. Finally, our findings show that one of the rare variants, located next to the primary one, further worsens the stiffening, while the other rare variant does not affect CI function. This approach may be extended to other variants candidate to exert a pathogenic impact on CI dynamics, or to investigate the interaction of multiple variants.

Molecular Dynamics Simulations Help Determine the Molecular Mechanisms of Lasioglossin-III and Its Variant Peptides' Membrane Interfacial Interactions

Lasioglossin-III (LL-III) is a potent broad-spectrum antimicrobial peptide used in diverse antimicrobial applications. In this work, coarse-grained and all-atom molecular dynamics simulation strategies were used in tandem to interpret the molecular mechanisms involved in the interfacial dynamics of LL-III and its recombinant variants during interactions with diverse cell membrane systems. Our results indicate that the membrane charges act as the driving force for initiating the membrane-peptide interactions, while the hydrophobic or van der Waals forces help to reinforce the membrane-peptide bindings. The optimized charge-hydrophobicity ratio of the LL-III peptides helps ensure their high specificity toward bacterial membranes compared to mammalian membrane systems, which also helps explain our experimental observations. Overall, we hope that our work gives new insight into the antimicrobial action of LL-III peptides and that the adopted simulation strategy will help other scientists and engineers extract maximal information from complex molecular simulations using minimal computational power.

There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning

We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.

Manufacturing process of liposomal Formation: A coarse-grained molecular dynamics simulation

A method of producing liposomes has been previously developed using a continuous manufacturing technology that involves a co-axial turbulent jet in co-flow. In this study, coarse-grained molecular dynamics (CG-MD) simulations were used to gain a deeper understanding of how the self-assembly process of liposomes is affected by the material attributes (such as the concentration of ethanol) and the process parameters (such as temperature), while also providing detailed information on a nano-scale molecular level. Specifically, the CG-MD simulations yield a comprehensive internal view of the structure and formation mechanisms of liposomes containing DPPC, DPPG, and cholesterol molecules. The importance of this work is that structural details on the molecular level are proposed, and such detail is not possible to obtain through experimental studies alone. The assessment of structural properties, including the area per lipid, diffusion coefficient, and order parameters, indicated that a thicker bilayer was observed at higher ethanol concentrations, while a thinner bilayer was present at higher temperatures. These conditions led to more water penetrating the interior of the bilayer and an unstable structure, as indicated by a larger contact area between lipids and water, and a higher coefficient of lipid lateral diffusion. However, stable liposomes were found through these evaluations at lower ethanol concentrations and/or lower process temperatures. Furthermore, the CG-MD model was further compared and validated with experimental and computational data including liposomal bilayer thickness and area per lipid measurements.

Inhibition and transport mechanisms of the ABC transporter hMRP5

Human multidrug resistance protein 5 (hMRP5) effluxes anticancer and antivirus drugs, driving multidrug resistance. To uncover the mechanism of hMRP5, we determine six distinct cryo-EM structures, revealing an autoinhibitory N-terminal peptide that must dissociate to permit subsequent substrate recruitment. Guided by these molecular insights, we design an inhibitory peptide that could block substrate entry into the transport pathway. We also identify a regulatory motif, comprising a positively charged cluster and hydrophobic patches, within the first nucleotide-binding domain that modulates hMRP5 localization by engaging with membranes. By integrating our structural, biochemical, computational, and cell biological findings, we propose a model for hMRP5 conformational cycling and localization. Overall, this work provides mechanistic understanding of hMRP5 function, while informing future selective hMRP5 inhibitor development. More broadly, this study advances our understanding of the structural dynamics and inhibition of ABC transporters.

Construction of Multiscale Dissipative Particle Dynamics (DPD) Models from Other Coarse-Grained Models

We present a general scheme for converting coarse-grained models into Dissipative Particle Dynamics (DPD) models. We build the corresponding DPD models by analogy with the de novo DPD coarse-graining scheme suggested by Groot and Warren (J. Chem. Phys., 1997). Electrostatic interactions between charged DPD particles are represented though the addition of a long-range Slater Coulomb potential as suggested by González-Melchor et al. (J. Chem. Phys., 2006). The construction is illustrated by converting MARTINI models for various proteins into a DPD representation, but it not restricted to the usual potential form in the MARTINI model-viz., Lennard-Jones potentials. We further extended the DPD scheme away from the typical use of homogeneous particle sizes, therefore faithfully representing the variations in the particle sizes seen in the underlying MARTINI model. The accuracy of the resulting construction of our generalized DPD models with respect to several structural observables has been benchmarked favorably against all-atom and MARTINI models for a selected set of peptides and proteins, and variations in the scales of the coarse-graining of the water solvent.

Building complex membranes with Martini 3

The Martini model is a popular force field for coarse-grained simulations. Membranes have always been at the center of its development, with the latest version, Martini 3, showing great promise in capturing more and more realistic behavior. In this chapter we provide a step-by-step tutorial on how to construct starting configurations, run initial simulations and perform dedicated analysis for membrane-based systems of increasing complexity, including leaflet asymmetry, curvature gradients and embedding of membrane proteins.

Dissecting the mechanism of atlastin-mediated homotypic membrane fusion at the single-molecule level

Homotypic membrane fusion of the endoplasmic reticulum (ER) is mediated by dynamin-like GTPase atlastin (ATL). This fundamental process relies on GTP-dependent domain rearrangements in the N-terminal region of ATL (ATLcyto), including the GTPase domain and three-helix bundle (3HB). However, its conformational dynamics during the GTPase cycle remain elusive. Here, we combine single-molecule FRET imaging and molecular dynamics simulations to address this conundrum. Different from the prevailing model, ATLcyto can form a loose crossover dimer upon GTP binding, which is tightened by GTP hydrolysis for membrane fusion. Furthermore, the α-helical motif between the 3HB and transmembrane domain, which is embedded in the surface of the lipid bilayer and self-associates in the crossover dimer, is required for ATL function. To recycle the proteins, Pi release, which disassembles the dimer, activates frequent relative movements between the GTPase domain and 3HB, and subsequent GDP dissociation alters the conformational preference of the ATLcyto monomer for entering the next reaction cycle. Finally, we found that two disease-causing mutations affect human ATL1 activity by destabilizing GTP binding-induced loose crossover dimer formation and the membrane-embedded helix, respectively. These results provide insights into ATL-mediated homotypic membrane fusion and the pathological mechanisms of related disease.

Cooperative Gating of a K+ Channel by Unmodified Biological Anionic Lipids Viewed by Solid-State NMR Spectroscopy

Lipids adhere to membrane proteins to stimulate or suppress molecular and ionic transport and signal transduction. Yet, the molecular details of lipid-protein interaction and their functional impact are poorly characterized. Here we combine NMR, coarse-grained molecular dynamics (CGMD), and functional assays to reveal classic cooperativity in the binding and subsequent activation of a bacterial inward rectifier potassium (Kir) channel by phosphatidylglycerol (PG), a common component of many membranes. Past studies of lipid activation of Kir channels focused primarily on phosphatidylinositol bisphosphate, a relatively rare signaling lipid that is tightly regulated in space and time. We use solid-state NMR to quantify the binding of unmodified 13C-PG to the K+ channel KirBac1.1 in liposomes. This specific lipid-protein interaction has a dissociation constant (Kd) of ∼7 mol percentage PG (ΧPG) with positive cooperativity (n = 3.8) and approaches saturation near 20% ΧPG. Liposomal flux assays show that K+ flux also increases with PG in a cooperative manner with an EC50 of ∼20% ΧPG, within the physiological range. Further quantitative fitting of these data reveals that PG acts as a partial (80%) agonist with fivefold K+ flux amplification. Comparisons of NMR chemical shift perturbation and CGMD simulations at different ΧPG confirm the direct interaction of PG with key residues, several of which would not be accessible to lipid headgroups in the closed state of the channel. Allosteric regulation by a common lipid is directly relevant to the activation mechanisms of several human ion channels. This study highlights the role of concentration-dependent lipid-protein interactions and tightly controlled protein allostery in the activation and regulation of ion channels.

Bilayer tension-induced clustering of the UPR sensor IRE1

The endoplasmic reticulum acts as a protein quality control center where a range of chaperones and foldases facilitates protein folding. IRE1 is a sensory transmembrane protein that transduces signals of proteotoxic stress by forming clusters and activating a cellular program called the unfolded protein response (UPR). Recently, membrane thickness variation due to membrane compositional changes have been shown to drive IRE1 cluster formation, activating the UPR even in the absence of proteotoxic stress. Here, we demonstrate a direct relationship between bilayer tension and UPR activation based on IRE1 dimer stability. The stability of the IRE1 dimer in a (50%DOPC-50%POPC) membrane at different applied bilayer tensions was analyzed via molecular dynamics simulations. The potential of mean force for IRE1 dimerization predicts a higher concentration of IRE1 dimers for both tensed and compressed ER membranes. This study shows that IRE1 may be a mechanosensitive membrane protein and establishes a direct biophysical relationship between bilayer tension and UPR activation.

Molecular Modeling and Simulation of the Mycobacterial Cell Envelope: From Individual Components to Cell Envelope Assemblies

Unlike typical Gram-positive bacteria, the cell envelope of mycobacteria is unique and composed of a mycobacterial outer membrane, also known as the mycomembrane, a peptidoglycan layer, and a mycobacterial inner membrane, which is analogous to that of Gram-negative bacteria. Despite its importance, however, our understanding of this complex cell envelope is rudimentary at best. Thus, molecular modeling and simulation of such an envelope can benefit the scientific community by proposing new hypotheses about the biophysical properties of its different layers. In this Perspective, we present recent advances in molecular modeling and simulation of the mycobacterial cell envelope from individual components to cell envelope assemblies. We also show how modeling other types of cell envelopes, such as that of Escherichia coli, may help modeling part of the mycobacterial envelopes. We hope that the studies presented here are just the beginning of the road and more and more new modeling and simulation studies help us to understand crucial questions related to mycobacteria such as antibiotic resistance or bacterial survival in the host.

Photodynamic inactivation of bacteria: Why it is not enough to excite a photosensitizer

BACKGROUND

The development of multidrug resistance (MDR) in infectious agents is one of the most serious global problems facing humanity. Antimicrobial photodynamic therapy (APDT) shows encouraging results in the fight against MDR pathogens, including those in biofilms.

METHODS

Photosensitizers (PS), monocationic methylene blue, polycationic and polyanionic derivatives of phthalocyanines, electroneutral and polycationic derivatives of bacteriochlorin were used to study photodynamic inactivation of Gram-positive and Gram-negative planktonic bacteria and biofilms under LED irradiation. Zeta potential measurements, confocal fluorescence imaging, and coarse-grained modeling were used to evaluate the interactions of PS with bacteria. PS aggregation and photobleaching were studied using absorption and fluorescence spectroscopy.

RESULTS

The main approaches to ensure high efficiency of bacteria photosensitization are analyzed.

CONCLUSIONS

PS must maintain a delicate balance between binding to exocellular and external structures of bacterial cells and penetration through the cell wall so as not to get stuck on the way to photooxidation-sensitive structures of the bacterial cell.

Computational microbiology of bacteria: Advancements in molecular dynamics simulations

Microbiology is traditionally considered within the context of wet laboratory methodologies. Computational techniques have a great potential to contribute to microbiology. Here, we describe our loose definition of "computational microbiology" and provide a short survey focused on molecular dynamics simulations of bacterial systems that fall within this definition. It is our contention that increased compositional complexity and realistic levels of molecular crowding within simulated systems are key for bridging the divide between experimental and computational microbiology.

Probing the Structural Topology and Dynamic Properties of gp28 Using Continuous Wave Electron Paramagnetic Resonance Spectroscopy

Lysis of Gram-negative bacteria by dsDNA phages is accomplished through either the canonical holin-endolysin pathway or the pinholin-SAR endolysin pathway. During lysis, the outer membrane (OM) is disrupted, typically by two-component spanins or unimolecular spanins. However, in the absence of spanins, phages use alternative proteins called Disruptin to disrupt the OM. The Disruptin family includes the cationic antimicrobial peptide gp28, which is found in the virulent podophage φKT. In this study, EPR spectroscopy was used to analyze the dynamics and topology of gp28 incorporated into a lipid bilayer, revealing differences in mobility, depth parameter, and membrane interaction among different segments and residues of the protein. Our results indicate that multiple points of helix 2 and helix 3 interact with the phospholipid membrane, while others are solvent-exposed, suggesting that gp28 is a surface-bound peptide. The CW-EPR power saturation data and helical wheel analysis confirmed the amphipathic-helical structure of gp28. Additionally, course-grain molecular dynamics simulations were further used to develop the structural model of the gp28 peptide associated with the lipid bilayers. Based on the data obtained in this study, we propose a structural topology model for gp28 with respect to the membrane. This work provides important insights into the structural and dynamic properties of gp28 incorporated into a lipid bilayer environment.

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.

Structural basis of dimerization of chemokine receptors CCR5 and CXCR4

G protein-coupled receptors (GPCRs) are prominent drug targets responsible for extracellular-to-intracellular signal transduction. GPCRs can form functional dimers that have been poorly characterized so far. Here, we show the dimerization mechanism of the chemokine receptors CCR5 and CXCR4 by means of an advanced free-energy technique named coarse-grained metadynamics. Our results reproduce binding events between the GPCRs occurring in the minute timescale, revealing a symmetric and an asymmetric dimeric structure for each of the three investigated systems, CCR5/CCR5, CXCR4/CXCR4, and CCR5/CXCR4. The transmembrane helices TM4-TM5 and TM6-TM7 are the preferred binding interfaces for CCR5 and CXCR4, respectively. The identified dimeric states differ in the access to the binding sites of the ligand and G protein, indicating that dimerization may represent a fine allosteric mechanism to regulate receptor activity. Our study offers structural basis for the design of ligands able to modulate the formation of CCR5 and CXCR4 dimers and in turn their activity, with therapeutic potential against HIV, cancer, and immune-inflammatory diseases.

Structural basis for triacylglyceride extraction from mycobacterial inner membrane by MFS transporter Rv1410

Mycobacterium tuberculosis is protected from antibiotic therapy by a multi-layered hydrophobic cell envelope. Major facilitator superfamily (MFS) transporter Rv1410 and the periplasmic lipoprotein LprG are involved in transport of triacylglycerides (TAGs) that seal the mycomembrane. Here, we report a 2.7 Å structure of a mycobacterial Rv1410 homologue, which adopts an outward-facing conformation and exhibits unusual transmembrane helix 11 and 12 extensions that protrude ~20 Å into the periplasm. A small, very hydrophobic cavity suitable for lipid transport is constricted by a functionally important ion-lock likely involved in proton coupling. Combining mutational analyses and MD simulations, we propose that TAGs are extracted from the core of the inner membrane into the central cavity via lateral clefts present in the inward-facing conformation. The functional role of the periplasmic helix extensions is to channel the extracted TAG into the lipid binding pocket of LprG.

Computational investigation on lipid bilayer disruption induced by amphiphilic Janus nanoparticles: combined effect of Janus balance and charged lipid concentration

Janus nanoparticles (NPs) with charged/hydrophobic compartments have garnered attention for their potential antimicrobial activity. These NPs have been shown to disrupt lipid bilayers in experimental studies, yet the underlying mechanisms of this disruption at the particle-membrane interface remain unclear. To address this knowledge gap, the present study conducts a computational investigation to systematically examine the disruption of lipid bilayers induced by amphiphilic Janus NPs. The focus of this study is on the combined effects of the hydrophobicity of the Janus NP, referred to as the Janus balance, defined as the ratio of hydrophilic to hydrophobic surface coverage, and the concentration of charged phospholipids on the interactions between Janus NPs and lipid bilayers. Computational simulations were conducted using a coarse-grained molecular dynamics (MD) approach. The results of these MD simulations reveal that while the area change of the bilayer increases monotonically with the Janus balance, the effect of charged lipid concentration in the membrane is not easy to be predicted. Specifically, it was found that the concentration of negatively charged lipids is directly proportional to the intensity of membrane disruption. Conversely, positively charged lipids have a negligible effect on membrane defects. This study provides molecular insights into the significant role of Janus balance in the disruption of lipid bilayers by Janus NPs and supports the selectivity of Janus NPs for negatively charged lipid membranes. Furthermore, the anisotropic properties of Janus NPs were found to play a crucial role in their ability to disrupt the membrane via the combination of hydrophobic and electrostatic interactions. This finding is validated by testing the current Janus NP design on a bacterial membrane-mimicking model. This computational study may serve as a foundation for further studies aimed at optimizing the properties of Janus NPs for specific antimicrobial applications.

Cardiolipin occupancy profiles of YidC paralogs reveal the significance of respective TM2 helix residues in determining paralog-specific phenotypes

YidC belongs to an evolutionarily conserved family of insertases, YidC/Oxa1/Alb3, in bacteria, mitochondria, and chloroplasts, respectively. Unlike Gram-negative bacteria, Gram-positives including Streptococcus mutans harbor two paralogs of YidC. The mechanism for paralog-specific phenotypes of bacterial YidC1 versus YidC2 has been partially attributed to the differences in their cytoplasmic domains. However, we previously identified a W138R gain-of-function mutation in the YidC1 transmembrane helix 2. YidC1W138R mostly phenocopied YidC2, yet the mechanism remained unknown. Primary sequence comparison of streptococcal YidCs led us to identify and mutate the YidC1W138 analog, YidC2S152 to W/A, which resulted in a loss of YidC2- and acquisition of YidC1-like phenotype. The predicted lipid-facing side chains of YidC1W138/YidC2S152 led us to propose a role for membrane phospholipids in specific-residue dependent phenotypes of S. mutans YidC paralogs. Cardiolipin (CL), a prevalent phospholipid in the S. mutans cytoplasmic membrane during acid stress, is encoded by a single gene, cls. We show a concerted mechanism for cardiolipin and YidC2 under acid stress based on similarly increased promoter activities and similar elimination phenotypes. Using coarse grain molecular dynamics simulations with the Martini2.2 Forcefield, YidC1 and YidC2 wild-type and mutant interactions with CL were assessed in silico. We observed substantially increased CL interaction in dimeric versus monomeric proteins, and variable CL occupancy in YidC1 and YidC2 mutant constructs that mimicked characteristics of the other wild-type paralog. Hence, paralog-specific amino acid- CL interactions contribute to YidC1 and YidC2-associated phenotypes that can be exchanged by point mutation at positions 138 or 152, respectively.

Phospholipid transport to the bacterial outer membrane through an envelope-spanning bridge

The outer membrane of Gram-negative bacteria provides a formidable barrier, essential for both pathogenesis and antimicrobial resistance. Biogenesis of the outer membrane requires the transport of phospholipids across the cell envelope. Recently, YhdP was implicated as a major protagonist in the transport of phospholipids from the inner membrane to the outer membrane however the molecular mechanism of YhdP mediated transport remains elusive. Here, utilising AlphaFold, we observe YhdP to form an elongated assembly of 60 β strands that curve to form a continuous hydrophobic groove. This architecture is consistent with our negative stain electron microscopy data which reveals YhdP to be approximately 250 Å in length and thus sufficient to span the bacterial cell envelope. Furthermore, molecular dynamics simulations and in vivo bacterial growth assays indicate essential helical regions at the N- and C-termini of YhdP, that may embed into the inner and outer membranes respectively, reinforcing its envelope spanning nature. Our in vivo crosslinking data reveal phosphate-containing substrates captured along the length of the YhdP groove, providing direct evidence that YhdP transports phospholipids. This finding is congruent with our molecular dynamics simulations which demonstrate the propensity for inner membrane lipids to spontaneously enter the groove of YhdP. Collectively, our results support a model in which YhdP bridges the cell envelope, providing a hydrophobic environment for the transport of phospholipids to the outer membrane.

Permeability of antioxidants through a lipid bilayer model with coarse-grained simulations

Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

Exploring the Molecular Dynamics of a Lipid-A Vesicle at the Atom Level: Morphology and Permeation Mechanism

Lipid-A was previously shown to spontaneously aggregate into a vesicle via the hybrid particle field approach. We assess the validity of the proposed vesiculation mechanism by simulating the resulting lipid-A vesicle at the atom level. The spatial confinement imposed by the vesicle geometry on the conformation and packing of lipid-A induces significant heterogeneity of physical properties in the inner and outer leaflets. It also induces tighter molecular packing and lower acyl chain order compared to the lamellar arrangement. Around 5% of water molecules passively permeates the vesicle membrane inward and outward. The permeation is facilitated by interactions with water molecules that are transported across the membrane by a network of electrostatic interactions with the hydrogen bond donors/acceptors in the N-acetylglucosamine ring and upper region of the acyl chains of lipid-A. The permeation process takes place at low rates but still at higher frequencies than observed for the lamellar arrangement of lipid-A. These findings not only substantiate the proposed lipid-A vesiculation mechanism but also reveal the complex structural dynamics of an important nonlamellar arrangement of lipid-A.

pyCHARMM: Embedding CHARMM Functionality in a Python Framework

CHARMM is rich in methodology and functionality as one of the first programs addressing problems of molecular dynamics and modeling of biological macromolecules and their partners, e.g., small molecule ligands. When combined with the highly developed CHARMM parameters for proteins, nucleic acids, small molecules, lipids, sugars, and other biologically relevant building blocks, and the versatile CHARMM scripting language, CHARMM has been a trendsetting platform for modeling studies of biological macromolecules. To further enhance the utility of accessing and using CHARMM functionality in increasingly complex workflows associated with modeling biological systems, we introduce pyCHARMM, Python bindings, functions, and modules to complement and extend the extensive set of modeling tools and methods already available in CHARMM. These include access to CHARMM function-generated variables associated with the system (psf), coordinates, velocities and forces, atom selection variables, and force field related parameters. The ability to augment CHARMM forces and energies with energy terms or methods derived from machine learning or other sources, written in Python, CUDA, or OpenCL and expressed as Python callable routines is introduced together with analogous functions callable during dynamics calculations. Integration of Python-based graphical engines for visualization of simulation models and results is also accessible. Loosely coupled parallelism is available for workflows such as free energy calculations, using MBAR/TI approaches or high-throughput multisite λ-dynamics (MSλD) free energy methods, string path optimization calculations, replica exchange, and molecular docking with a new Python-based CDOCKER module. CHARMM accelerated platform kernels through the CHARMM/OpenMM API, CHARMM/DOMDEC, and CHARMM/BLaDE API are also readily integrated into this Python framework. We anticipate that pyCHARMM will be a robust platform for the development of comprehensive and complex workflows utilizing Python and its extensive functionality as well as an optimal platform for users to learn molecular modeling methods and practices within a Python-friendly environment such as Jupyter Notebooks.

Effect of cetrimonium carrier micelles on bacterial membranes and extracellular DNA, an in silico study

Microorganisms do not live as dispersed single cells but rather they form aggregates with extracellular polymeric substances at interfaces. Biofilms are considered efficient life forms because they shield bacteria from biocides and collect dilute nutrients. This is a big concern in industry since the microorganisms can colonize a wide range of surfaces, accelerating material deterioration, colonizing medical devices, contaminating ultrapure drinking water, increasing energy costs and creating focus of infection. Conventional biocides that target a specific component of the bacteria are not effective in the presence of biofilms. Efficient biofilm inhibitors are based on a multitarget approach interacting with the bacteria and the biofilm matrix. Their rationale design requires a thorough understanding of inhibitory mechanisms that are still largely lacking today. Herein we uncover via molecular modelling the inhibition mechanism of cetrimonium 4-OH cinnamate (CTA-4OHcinn). Simulations show that CTA-4OH micelles can disrupt symmetric and asymmetric bilayers, representative of inner and outer bacterial membranes, following three stages: adsorption, assimilation, and defect formation. The main driving force for micellar attack is electrostatic interactions. In addition to disrupting the bilayers, the micelles work as carriers facilitating the trapping of 4OH cinnamate anions within the bilayer upper leaflet and overcoming electrostatic repulsion. The micelles also interact with extracellular DNA (e-DNA), which is one of the main components of biofilms. It is observed that CTA-4OHcinn forms spherical micelles on the DNA backbone; which hinders their ability to pack. This is demonstrated by modelling the DNA along the hbb histone-like protein, showing that in the presence of CTA-4OHcinn, DNA does not pack properly around hbb. The abilities of CTA-4OHcinn to cause cell death through membrane disruption and to disperse a mature, multi-species biofilm are also confirmed experimentally.

Nonlocal Effects of Antibiotic-Resistance-Causing Mutations Reveal an Alternative Region for Targeting on FtsW-Penicillin-Binding Protein 3 Complex of Haemophilus influenzae

Currently prescribed antibiotics target the catalytic sites of wild-type bacterial proteins; however, bacteria adopt mutations at this site, eventually leading to the emergence of resistance. Therefore, the identification of alternative drug binding sites is crucial, which requires knowledge of the dynamics of the mutant protein. Here, we set out to investigate the impact of a high-resistance-causing triple mutation (S385T + L389F + N526K) on the dynamics of a prioritized resistant pathogen, Haemophilus influenzae, using computational techniques. We studied penicillin-binding protein 3 (PBP3) and its complex with FtsW, which display resistance toward β-lactam antibiotics. We showed that mutations displayed local and nonlocal effects. In terms of the former, the orientation of the β-sheet, which surrounds the active site of PBP3, was impacted and the catalytic site was exposed to the periplasmic region. In addition, the flexibility of the β3-β4 loop, which modulates the catalysis of the enzyme, increased in the mutant FtsW-PBP3 complex. As for nonlocal effects, the dynamics of the pedestal domain (N-terminal periplasmic modulus (N-t)), i.e., the opening of the fork, was different between the wild-type and mutant enzymes. We showed the closed fork caused a greater number of residues to participate in the hypothesized allosteric communication network connecting N-t to the transpeptidase domain in the mutant enzyme. Finally, we demonstrated that the closed fork results in more favorable binding with β-lactam antibiotics, particularly cefixime, suggesting that small therapeutics that can stabilize the closed fork of mutant PBP3 may lead to the development of more effective molecules to combat resistant bacteria.

HylleraasMD: A Domain Decomposition-Based Hybrid Particle-Field Software for Multiscale Simulations of Soft Matter

We present HylleraasMD (HyMD), a comprehensive implementation of the recently proposed Hamiltonian formulation of hybrid particle-field molecular dynamics. The methodology is based on a tunable, grid-independent length-scale of coarse graining, obtained by filtering particle densities in reciprocal space. This enables systematic convergence of energies and forces by grid refinement, also eliminating nonphysical force aliasing. Separating the time integration of fast modes associated with internal molecular motion from slow modes associated with their density fields, we enable the first time-reversible, energy-conserving hybrid particle-field simulations. HyMD comprises the optional use of explicit electrostatics, which, in this formalism, corresponds to the long-range potential in particle-mesh Ewald. We demonstrate the ability of HyMD to perform simulations in the microcanonical and canonical ensembles with a series of test cases, comprising lipid bilayers and vesicles, surfactant micelles, and polypeptide chains, comparing our results to established literature. An on-the-fly increase of the characteristic coarse-grain length significantly speeds up dynamics, accelerating self-diffusion and leading to expedited aggregation. Exploiting this acceleration, we find that the time scales involved in the self-assembly of polymeric structures can lie in the tens to hundreds of picoseconds instead of the multimicrosecond regime observed with comparable coarse-grained models.

Bilayer lipid membrane formation on surface assemblies with sparsely distributed tethers

A combined computational and experimental study of small unilamellar vesicle (SUV) fusion on mixed self-assembled monolayers (SAMs) terminated with different deuterated tether moieties (-(CD₂)₇CD₃ or -(CD₂)₁₅CD₃) is reported. Tethered bilayer lipid membrane (tBLM) formation of synthetic 1-stearoyl-2-oleoyl-sn-glycero-3-phosphocholine was initially probed on SAMs with controlled tether (d-alkyl tail) surface densities and lateral molecular packing using quartz crystal microbalance with dissipation monitoring (QCM-D). Long time-scale coarse-grained molecular dynamics (MD) simulations were then employed to elucidate the mechanisms behind the interaction between the SUVs and the different phases formed by the -(CD₂)₇CD₃ and -(CD₂)₁₅CD₃ tethers. Furthermore, a series of real time kinetics was recorded under different osmotic conditions using QCM-D to determine the accumulated lipid mass and for probing the fusion process. It is shown that the key factors driving the SUV fusion and tBLM formation on this type of surfaces involve tether insertion into the SUVs along with vesicle deformation. It is also evident that surface densities of the tethers as small as a few mol% are sufficient to obtain stable tBLMs with a high reproducibility. The described "sparsely tethered" tBLM system can be advantageous in studying different biophysical phenomena, such as membrane protein insertion, effects of receptor clustering, and raft formation.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Soft Matter under Pressure: Pushing Particle-Field Molecular Dynamics to the Isobaric Ensemble

Hamiltonian hybrid particle-field molecular dynamics is a computationally efficient method to study large soft matter systems. In this work, we extend this approach to constant-pressure (NPT) simulations. We reformulate the calculation of internal pressure from the density field by taking into account the intrinsic spread of the particles in space, which naturally leads to a direct anisotropy in the pressure tensor. The anisotropic contribution is crucial for reliably describing the physics of systems under pressure, as demonstrated by a series of tests on analytical and monatomic model systems as well as realistic water/lipid biphasic systems. Using Bayesian optimization, we parametrize the field interactions of phospholipids to reproduce the structural properties of their lamellar phases, including area per lipid, and local density profiles. The resulting model excels in providing pressure profiles in qualitative agreement with all-atom modeling, and surface tension and area compressibility in quantitative agreement with experimental values, indicating the correct description of long-wavelength undulations in large membranes. Finally, we demonstrate that the model is capable of reproducing the formation of lipid droplets inside a lipid bilayer.