A computer perspective of membranes: molecular dynamics studies of lipid bilayer systems

Molecular Dynamics Simulation on the Conformational Change of a pH-Switchable Lipid

pH-sensitive lipids are important components of lipid nanoparticles, which enable the targeted delivery and controlled release of drugs. Understanding the mechanism of pH-triggered drug release at the molecular level is important for the rational design of ionizable lipids. Based on a recently reported pH-switchable lipid, named SL2, molecular dynamics (MD) simulations were employed to explore the microscopic mechanism behind the membrane destabilization induced by the conformational change of pH-switchable lipids. The simulated results showed that, at neutral pH, the neutral SL2 lipids assembled with other components (helper lipids and cholesterol) to form a structurally ordered bilayer structure. At this moment, the two hydrocarbon chains of SL2 were closely aligned and inserted in an orderly manner inside of the membrane. With a decrease in pH, the protonation of the pyridinium ring caused a large degree of molecular structural change. The pyridinium ring preferred to form intramolecular H-bonds with the methoxy groups and intermolecular H-bonds with water, resulting in the flip of the pyridinium ring. Meanwhile, due to the structural flip, the two alkane chains showed a more open state, which perturbed the arrangement of molecules within the membrane. The perturbations caused local collapse of the membrane and the formation of water molecule channels, which contributed to the pH-induced drug release. Our results verified the experimentally proposed mechanism at the molecular level and provided more complementary information, which are expected to have deeper insights into the pH-triggered drug release.

Membrane fusion by dengue virus: The first step

Flaviviruses include important human pathogens such as Dengue, Zika, West Nile, Yellow fever, Japanese encephalitis, and Tick-borne encephalitis viruses as well as some emerging viruses that affect millions of people worldwide. They fuse their membrane with the late endosomal one in a pH-dependent way and therefore the merging of the membranes is one of the main goals for obtaining new antivirals. The envelope E protein, a membrane fusion protein, is accountable for fusion and encompasses different domains involved in the fusion mechanism, including the fusion peptide segment. In this work we have used molecular dynamics to study the interaction of the distal end of domain II of the DENV envelope E protein with a membrane like the late endosomal membrane in order to observe the initiation of membrane fusion carried out by a number of trimers of the DENV envelope E protein interacting with a complex biomembrane and demonstrate its feasibility. Our results demonstrate the likelihood of membrane disorganization and pore formation by trimer complex organization, the amino acids responsible for such condition and the secondary structure arrangements needed for such fundamental process. At the same time, we define new targets of the envelope E protein sequence which could permit designing potent antiviral bioactive molecules.

Localization, aggregation, and interaction of glycyrrhizic acid with the plasma membrane

Glycyrrhizic acid (GLA) is the most important bioactive constituent of licorize root and exhibits antiviral, antimicrobial, anti-oxidant, anti-inflammatory, anti-allergic, and antitumor activities. GLA has an amphiphilic nature consisting of two hydrophilic and one hydrophobic part, and its mechanism of action could be mediated by its incorporation into the membrane. Furthermore, GLA presents two different forms, protonated (GLA) and deprotonated (GLAD), and has been suggested that their location inside the membrane could be different. Since GLA could be a source against many types of diseases, we have localized the GLA molecule in the presence of a complex membrane and established the detailed interactions of GLA with lipids using all-atom molecular dynamics. Our outcomes sustain that GLA/GLAD tend to locate amid the CHOL oxygen atom and the phospholipid phosphates, preferably perpendicular to the membrane surface, increasing membrane fluidity. Interestingly, GLA and GLAD tend to be surrounded by specific phospholipids, different for each type of molecule. Outstandingly, both GLA and GLAD tend to spontaneously associate in solution forming aggregates, precluding them from inserting into the membrane and, therefore, interacting with it. Consequently, some of the biological properties of GLA/GLAD could be credited to the alteration of the membrane biophysical properties by interacting with specific lipids. However, the formation of an aggregate in solution could hinder its bioactive properties and should be considered a suited vehicle when prepared to be used in biological or clinical assays.

Enhancing the Assembly Properties of Bottom-Up Coarse-Grained Phospholipids

A plethora of key biological events occur at the cellular membrane where the large spatiotemporal scales necessitate dimensionality reduction or coarse-graining approaches over conventional all-atom molecular dynamics simulation. Constructing coarse-grained descriptions of membranes systematically from statistical mechanical principles has largely remained challenging due to the necessity of capturing amphipathic self-assembling behavior in coarse-grained models. We show that bottom-up coarse-grained lipid models can possess metastable morphological behavior and that this potential metastability has ramifications for accurate development and training. We in turn develop a training algorithm which evades metastability issues by linking model training to self-assembling behavior, and demonstrate its robustness via construction of solvent-free coarse-grained models of various phospholipid membranes, including lipid species such as phosphatidylcholines, phosphatidylserines, sphingolipids, and cholesterol. The resulting coarse-grained lipid models are orders of magnitude faster than their atomistic counterparts while retaining structural fidelity and constitute a promising direction for the development of coarse-grained models of realistic cell membranes.

Structure-dependent destructive adsorption of organophosphate flame retardants on lipid membranes

The widespread use of organophosphate flame retardants (OPFRs), a serious type of pervasive environmental contaminants, has led to a global concern regarding their diverse toxicities to living beings. Using a combination of experimental and theoretical approaches, we systematically studied the adsorption, accumulation, and influence of a series of OPFRs on the lipid membranes of bacteria and cells. Our results revealed that OPFRs can aggregate in lipid membranes, leading to the destruction of membrane integrity. During this process, the molecular structure of the OPFRs is a dominant factor that significantly influences the strength of their interaction with the lipid membrane, resulting in varying degrees of biotoxicity. Triphenyl phosphate (TPHP), owing to its large molecular size and strong hydrophobicity, causes severe membrane disruption through the formation of nanoclusters. The corresponding severe toxicity originates from the phase transitions of the lipid membranes. In contrast, smaller OPFRs such as triethyl phosphate (TEP) and tris(2-chloroethyl) phosphate (TCEP) have weaker hydrophobicity and induce minimal membrane disturbance and ineffective damage. In vivo, gavage of TPHP induced more severe barrier damage and inflammatory infiltration in mice than TEP or TCEP, confirming the higher toxicity of TPHP. Overall, our study elucidates the structure-dependent adsorption of OPFRs onto lipid membranes, highlighting their destructive interactions with membranes as the origin of OPFR toxicity.

The Membrane Dipole Potential and the Roles of Interfacial Water and Lipid Hydrocarbon Chains

Understanding membrane charge transport processes, including the actions of ion channels, pumps, carriers, and membrane-active peptides, requires a description of the electrostatics of the lipid bilayer. We have simulated a library of different lipid chemistries to reveal the impact of the headgroup, glycerol backbone, and hydrocarbon chains on the membrane dipole potential. We found a strong dependence of the potential on lipid packing, but this was not caused by the packing of lipid polar components, due to cancellation of their electric fields by electrolyte. In contrast, lipid tail contributions were determined by area per lipid, arising from two countering effects. Increased area per lipid leads to chain tilting that increases methylene dipole projections to strengthen the electric field within the bilayer, while at the same time decreasing the electric field from terminal methyl groups. Moreover, electric fields from some nonterminal groups and the terminal methyl group can extend beyond the bilayer center and be canceled by the opposing leaflet. This interleaflet field annulment explains the experimental reduction in dipole potential for unsaturated and branched lipid bilayers, by as much as ∼200 mV, as well as experiments that substitute chain carbons with sulfur. Replacing ester with ether groups (eliminating two carbonyl groups) causes a significant reduction in potential, also by ∼200 mV, in agreement with experiment. We show that the effect can be largely attributed to the loss of aligned water molecules in the glycerol backbone region, lowering the potential inside the bilayer core. When only one of the two carbonyls is removed (using a hybrid ester-ether lipid or a single-chain lipid), most of this reduction in potential was lost, with the single carbonyl group able to maintain full hydration in the interfacial region. While headgroup chemistry can have a major effect (by as much as ±100 mV relative to phosphatidylcholine), anionic headgroups either decrease or increase the dipole potential, with the variation involving perturbation in hydrogen-bonded water molecules and changes in packing of lipid tails. Overall, these results suggest that membrane electrostatics are dominated by aligned water molecules at the polar-hydrocarbon interface and, surprisingly, by the charge distribution of the nonpolar lipid tails, and not the packing of headgroup and glycerol carbonyl dipoles.

Bisphenol F and Bisphenol S in a Complex Biomembrane: Comparison with Bisphenol A

Bisphenols are a group of endocrine-disrupting chemicals used worldwide for the production of plastics and resins. Bisphenol A (BPA), the main bisphenol, exhibits many unwanted effects. BPA has, currently, been replaced with bisphenol F (BPF) and bisphenol S (BPS) in many applications in the hope that these molecules have a lesser effect on metabolism than BPA. Since bisphenols tend to partition into the lipid phase, their place of choice would be the cellular membrane. In this paper, I carried out molecular dynamics simulations to compare the localization and interactions of BPA, BPF, and BPS in a complex membrane. This study suggests that bisphenols tend to be placed at the membrane interface, they have no preferred orientation inside the membrane, they can be in the monomer or aggregated state, and they affect the biophysical properties of the membrane lipids. The properties of bisphenols can be attributed, at least in part, to their membranotropic effects and to the modulation of the biophysical membrane properties. The data support that both BPF and BPS, behaving in the same way in the membrane as BPA and with the same capacity to accumulate in the biological membrane, are not safe alternatives to BPA.

Location and interaction of idebenone and mitoquinone in a membrane similar to the inner mitochondrial membrane. Comparison with ubiquinone 10

Oxygen is essential for aerobic life on earth but it is also the origin of harmful reactive oxygen species (ROS). Ubiquinone is par excellence the endogenous cellular antioxidant, but a very hydrophobic one. Because of that, other molecules have been envisaged, such as idebenone (IDE) and mitoquinone (MTQ), molecules having the same redox active benzoquinone moiety but higher solubility. We have used molecular dynamics to determine the location and interaction of these molecules, both in their oxidized and reduced forms, with membrane lipids in a membrane similar to that of the mitochondria. Both IDE and reduced IDE (IDOL) are situated near the membrane interface, whereas both MTQ and reduced MTQ (MTQOL) locate in a position adjacent to the phospholipid hydrocarbon chains. The quinone moieties of both ubiquinone 10 (UQ10) and reduced UQ10 (UQOL10) in contraposition to the same moieties of IDE, IDOL, MTQ and MTQOL, located near the membrane interphase, whereas the isoprenoid chains remained at the middle of the hydrocarbon chains. These molecules do not aggregate and their functional quinone moieties are located in the membrane at different depths but near the hydrophobic phospholipid chains whereby protecting them from ROS harmful effects.

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

Validating Fractalkine receptor as a target and identifying candidates for drug discovery against type 2 diabetes

Type 2 diabetes mellitus (T2DM) is one of the most common chronic diseases employing abnormal levels of insulin. Enhancing the insulin production is greatly aided by the regulatory mechanisms of the Fractalkine receptor (CX3CR1) system in islet β-cell function. However, elements including a high-fat diet, obesity, and ageing negatively impact the expression of CX3CR1 in islets. CX3CL1/CX3CR1 receptor-ligand complex is now recognized as a novel therapeutic target. It suggests that T2DM-related β-cell dysfunction may result from lower amount of these proteins. We analyzed the differential expression of CX3CR1 gene samples taken from persons with T2DM using data obtained from the Gene Expression Omnibus database. Homology modeling enabled us to generate the three-dimensional structure of CX3CR1 and a possible binding pocket. The optimized CX3CR1 structure was subjected to rigorous screening against a massive library of 693 million drug-like molecules from the ZINC15 database. This screening process led to the identification of three compounds with strong binding affinity at the identified binding pocket of CX3CR1. To further evaluate the potential of these compounds, molecular dynamics simulations were conducted over a 50 ns time scale to assess the stability of the protein-ligand complexes. These simulations revealed that ZINC000032506419 emerged as the most promising drug-like compound among the three potent molecules. The discovery of ZINC000032506419 holds exciting promise as a potential therapeutic agent for T2D and other related metabolic disorders. These findings pave the way for the development of effective medications to address the complexities of T2DM and its associated metabolic diseases.

Designing Reaction Coordinate for Ion-Induced Pore-Assisted Mechanism of Halide Ions Permeation through Lipid Bilayer by Umbrella Sampling

Ion permeation mechanism through lipid membranes helps to understand cellular processes. We propose new reaction coordinates that allow ions to permeate according to their water affinity and interaction with the hydrophilic layer. Simulations were done for three different halides (F-, Cl-, and I-) in two different lipid bilayers, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dinervonoyl-sn-glycero-3-phosphocholine (DNPC). It is found that the involvement of the water molecules decreases the free energy barrier. The ions were found to follow different pathways for permeation. Formation of proper pores required a collaboration effort of the hydration shell water molecules and the hydrophilic lipid layer, which was favored in the case of Cl- ions. The optimum charge density and good water affinity of Cl- with respect to F- and I- ions helped to form the pore. The effect was prominently seen in the case of DNPC membrane because of its higher hydrophobic thickness. The umbrella sampling results were compared with other methods such as the Markov state model (MSM) and well-tempered metadynamics (WT-metaD).

Effect of quaternary ammonium surfactants on biomembranes using molecular dynamics simulation

Research conducted both prior to and after the emergence of the COVID-19 pandemic reveals a notable rise in human exposure to cleaning products, hand sanitizers, and personal care items. Moreover, there has been a corresponding increase in the environmental release of these chemicals. Cleaning and disinfecting products often contain quaternary ammonium compounds (QACs) with alkyl chains as long as 8-12 carbon atoms. The attachment of quaternary ammonium surfactants to the membrane resulted in the deformation of the bilayer and membrane disruption. Before interactions with cell membranes, these surfactant molecules may form different aggregates depending on their architecture. Interaction of surfactant monomers or clusters with the cell membrane changes the physiochemical properties of the biomembranes. To investigate this interaction and its influence on membrane properties, we conducted molecular dynamics simulations of cationic quaternary ammonium surfactants interacting with dipalmitoylphosphatidylcholine (DPPC) membranes. Our simulations revealed significant interactions between the surfactants and the phospholipids, leading to substantial alterations in the structure of the bilayer. The results are compared with the simulated anionic (SDS) and nonionic surfactants/bilayer systems. Various aspects were considered, including the aggregation process, migration behavior, and eventual equilibrium of these molecules at the interface between the membrane and water. This analysis used various techniques such as density profiles, distribution functions, cluster analysis, order parameters, hydrogen bonding (H-bonding), and mean-square displacements. The results indicate that while surfactants with shorter alkyl tails (N-(2-hydroxyethyl)-N,N-dimethyloctan-1-aminium chloride (HEDMOAC)) make strong hydrogen bonds with the phosphate group and ester oxygen of the phosphatidylcholine bilayer and enter toward the bilayer in the monomer form, surfactants of longer alkyl tails aggregated on the membrane head-water interface and interact minimally with the head groups of the DPPC bilayer. For DDEDMEAC, a quaternary ammonium surfactant with a hydrophobic alkyl chain consisting of two decanoate groups, alteration of the structural and dynamical properties of the bilayer is expected to be governed by two different factors. First, the structural order of DPPC increases as surfactant aggregates interact with the membrane head group. Second, the decrease in the order of the bilayer occurs due to the insertion of surfactant monomers within the hydrophobic region of the bilayer. Strong interactions between constituents of tetraoctylammonium bromide (TOABr) and lipid head groups lead to a reduction in interlipid interactions and order, which further results in increased porosity of cellular membranes. Understanding the extent of these interactions plays a pivotal role in the toxicological assessment of these surfactants.

The impact of physiologically relevant temperatures on physical properties of thylakoid membranes: a molecular dynamics study

Thylakoid membranes are energy-converting membranes with a unique lipid composition. Though the membranes are primarily composed of proteins, their photosynthetic function is strongly influenced by the lipid constituents. Here we characterize, with molecular dynamics (MD) simulations, lipid bilayers with compositions representative of plant thylakoid membranes. We determine, in a wide range of temperatures, the physical parameters of the model membranes which are relevant for the photosynthetic function. We found a marked impact of temperature on membrane permeability due to a combination of increased compressibility and curvature of the membrane at elevated temperatures. With increasing temperatures, we observed increasingly smeared transmembrane density profiles of the membrane forming lipid headgroups predicting increased membrane flexibility. The diffusion coefficient of the lipids increased with temperature without apparent specificity for lipid species. Instead of a comprehensive experimental dataset in the relevant temperature range, we quantitatively compared and validated our MD results with MD simulations on a dipalmitoylphosphatidylcholine model system.

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.

Dynamics of fluid bilayer vesicles: Soft meshes and robust curvature energy discretization

Continuum models like the Helfrich Hamiltonian are widely used to describe fluid bilayer vesicles. Here we study the molecular dynamics compatible dynamics of the vertices of two-dimensional meshes representing the bilayer, whose in-plane motion is only weakly constrained. We show (i) that Jülicher's discretization of the curvature energy offers vastly superior robustness for soft meshes compared to the commonly employed expression by Gommper and Kroll and (ii) that for sufficiently soft meshes, the typical behavior of fluid bilayer vesicles can emerge even if the mesh connectivity remains fixed throughout the simulations. In particular, soft meshes can accommodate large shape transformations, and the model can generate the typical ℓ^{-4} signal for the amplitude of surface undulation modes of nearly spherical vesicles all the way up to the longest wavelength modes. Furthermore, we compare results for Newtonian, Langevin, and Brownian dynamics simulations of the mesh vertices to demonstrate that the internal friction of the membrane model is negligible, making it suitable for studying the internal dynamics of vesicles via coupling to hydrodynamic solvers or particle-based solvent models.

SARS-CoV-2 Protein S Fusion Peptide Is Capable of Wrapping Negatively-Charged Phospholipids

COVID-19, caused by SARS-CoV-2, which is a positive-sense, single-stranded RNA enveloped virus, emerged in late 2019 and was declared a worldwide pandemic in early 2020 causing more than 600 million infections so far and more than 6 million deaths in the world. Although new vaccines have been implemented, the pandemic continues to impact world health dramatically. Membrane fusion, critical for the viral entry into the host cell, is one of the main targets for the development of novel antiviral therapies to combat COVID-19. The S2 subunit of the viral S protein, a class I membrane fusion protein, contains the fusion domain which is directly implicated in the fusion mechanism. The knowledge of the membrane fusion mechanism at the molecular level will undoubtedly result in the development of effective antiviral strategies. We have used all-atom molecular dynamics to analyse the binding of the SARS-CoV-2 fusion peptide to specific phospholipids in model membranes composed of only one phospholipid plus cholesterol in the presence of either Na⁺ or Ca²⁺. Our results show that the fusion peptide is capable of binding to the membrane, that its secondary structure does not change significantly upon binding, that it tends to preferentially bind electronegatively charged phospholipids, and that it does not bind cholesterol at all. Understanding the intricacies of the membrane fusion mechanism and the molecular interactions involved will lead us to the development of antiviral molecules that will allow a more efficient battle against these viruses.

Structural analysis of human G-protein-coupled receptor 17 ligand binding sites

The human G protein coupled membrane receptor (GPR17), the sensor of brain damage, is identified as a biomarker for many neurological diseases. In human brain tissue, GPR17 exist in two isoforms, long and short. While cryo-electron microscopy technology has provided the structure of the long isoform of GPR17 with Gi complex, the structure of the short isoform and its activation mechanism remains unclear. Recently, we theoretically modeled the structure of the short isoform of GPR17 with Gi signaling protein and identified novel ligands. In the present work, we demonstrated the presence of two distinct ligand binding sites in the short isoform of GPR17. The molecular docking of GPR17 with endogenous (UDP) and synthetic ligands (T0510.3657, MDL29950) found the presence of two distinct binding pockets. Our observations revealed that endogenous ligand UDP can bind stronger in two different binding pockets as evidenced by glide and autodock vina scores, whereas the other two ligand's binding with GPR17 has less docking score. The analysis of receptor-UDP interactions shows complexes' stability in the lipid environment by 100 ns atomic molecular dynamics simulations. The amino acid residues VAL83, ARG87, and PHE111 constitute ligand binding site 1, whereas site 2 constitutes ASN67, ARG129, and LYS232. Root mean square fluctuation analysis showed the residues 83, 87, and 232 with higher fluctuations during molecular dynamics simulation in both binding pockets. Our findings imply that the residues of GPR17's two binding sites are crucial, and their interaction with UDP reveals the protein's hidden signaling and communication properties. Furthermore, this finding may assist in the development of targeted therapies for the treatment of neurological diseases.

Bio-membranes: Picosecond to second dynamics and plasticity as deciphered by solid state NMR

Since the first membrane models in the 1970s, the concept of biological membranes has evolved considerably. The membrane is now seen as a very complex mixture whose dynamic behavior is even more complex. Solid-state NMR is well suited for such studies as it can probe the movements of the membrane from picoseconds to seconds. Two NMR observables can be used: motionally averaged spectra and relaxation times. They bring information on order parameters, phase transitions, correlation times, activation energies and membrane elasticity. Spectra are used to determine the nature of the membrane phase. The order parameters can be measured directly from spectra that are dominated by quadrupolar, dipolar and chemical shielding magnetic interactions and allow describing the lipid membrane as being very rigid at the glycerol and chain level and very fluid at its center and surface. Correlation times and activation energies can be measured for intramolecular motions (pico to nanoseconds), molecular motions (nano to 100 ns) and collective modes of membrane deformation (microseconds). Sterols modulate membrane phases, order parameters, correlation times and membrane elasticity. In general terms, sterols tend to act to reduce the impact of environmental changes on molecular order and dynamics. They can be described as regulators of membrane dynamics by keeping them in a state of dynamics that changes very little when the temperature or other factors change. The presence of such large-scale membrane dynamics is proposed as a means of adapting to evolutionary constraints.

Bergamottin: location, aggregation and interaction with the plasma membrane

Bioactive furanocoumarins, a group of natural secondary metabolites common in higher plants, are recognized for their benefits to human health and have been shown to have numerous biological properties. However, the knowledge of its biomolecular mechanism is not known. One of the main furanocoumarins is bergamottin (BGM), which is characterized by a planar three-ringed structure and a hydrocarbon chain, which give BGM its high lipid/water partition coefficient. Because of that, and although the biological mechanism of BGM is not known, BGM bioactive properties could be ascribed to its potential to interact with the biological membrane, modulating its structure, changing its dynamics and at the same time that it might interact with lipids. For our goal, we have applied molecular dynamics to determine the position of BGM in a complex membrane and discern the possibility of certain interactions with membrane lipids. Our findings establish that BGM tends to locate in the middle of the hydrocarbon layer of the membrane, inserts in between the hydrocarbon chains of the phospholipids in an oblique position with respect to the membrane plane, increasing the fluidity of the membrane. Significantly, BGM tends to be surrounded by POPC molecules but exclude the molecule of CHOL. Outstandingly, BGM molecules associate spontaneously creating aggregates, which does not preclude them from interacting with and inserting into the membrane. The bioactive properties of BGM could be ascribed to its membranotropic effects and support the improvement of these molecules as therapeutic molecules, giving place to new opportunities for potential medical improvements.Communicated by Ramaswamy H. Sarma.

Computing Individual Area per Head Group Reveals Lipid Bilayer Dynamics

Lipid bilayers express a range of phases from solid-like to gel-like to liquid-like as a function of temperature and lipid surface concentration. The area occupied per lipid head group serves as one useful indicator of the bilayer phase, in conjunction with the two-dimensional radial distribution function (i.e., structure factor) within the bilayer. Typically, the area per head group is determined by dividing the bilayer area equally among all head groups. Such an approach is less satisfactory for a multicomponent set of diverse lipids. In this work, area determination is performed on a lipid-by-lipid basis by attributing to a lipid the volume that surrounds each atom. Voronoi tessellation provides this division of the interfacial region on a per-atom basis. The method is applied to a multicomponent system of water, NaCl, and 19 phospholipid types that was devised recently [Langmuir2022, 38, 9481-9499] as a computational representation of the Gram-positive Staphylococcus aureus phospholipid bilayer. Results demonstrate that lipids and water molecules occupy similar extents of area within the interfacial region; ascribing area only to head groups implicitly incorporates assumptions about head group hydration. Results further show that lipid tails provide non-negligible contributions to area on the membrane side of the bilayer-water interface. Results for minimum and maximum area of individual lipids reveal that spontaneous fluctuations displace head groups more than 10 Å from the interfacial region during an NPT simulation at 310 K, leading to a zero contribution to total area at some times. Total area fluctuations and fluctuations per individual lipid relax with a correlation time of ∼10 ns. The method complements density profile as an approach to quantify the structure and dynamics of computational lipid bilayers.

Protein drift-diffusion dynamics and phase separation in curved cell membranes and dendritic spines: Hybrid discrete-continuum methods

We develop methods for investigating protein drift-diffusion dynamics in heterogeneous cell membranes and the roles played by geometry, diffusion, chemical kinetics, and phase separation. Our hybrid stochastic numerical methods combine discrete particle descriptions with continuum-level models for tracking the individual protein drift-diffusion dynamics when coupled to continuum fields. We show how our approaches can be used to investigate phenomena motivated by protein kinetics within dendritic spines. The spine geometry is hypothesized to play an important biological role regulating synaptic strength, protein kinetics, and self-assembly of clusters. We perform simulation studies for model spine geometries varying the neck size to investigate how phase-separation and protein organization is influenced by different shapes. We also show how our methods can be used to study the roles of geometry in reaction-diffusion systems including Turing instabilities. Our methods provide general approaches for investigating protein kinetics and drift-diffusion dynamics within curved membrane structures.

Molecular oxygen as a probe molecule in EPR spin-labeling studies of membrane structure and dynamics

Molecular oxygen (O2) is the perfect probe molecule for membrane studies carried out using the saturation recovery EPR technique. O2 is a small, paramagnetic, hydrophobic enough molecule that easily partitions into a membrane's different phases and domains. In membrane studies, the saturation recovery EPR method requires two paramagnetic probes: a lipid-analog nitroxide spin label and an oxygen molecule. The experimentally derived parameters of this method are the spin-lattice relaxation times (T 1s) of spin labels and rates of bimolecular collisions between O2 and the nitroxide fragment. Thanks to the long T 1 of lipid spin labels (from 1 to 10 μs), the approach is very sensitive to changes of the local (around the nitroxide fragment) O2 diffusion-concentration product. Small variations in the lipid packing affect O2 solubility and O2 diffusion, which can be detected by the shortening of T 1 of spin labels. Using O2 as a probe molecule and a different lipid spin label inserted into specific phases of the membrane and membrane domains allows data about the lateral arrangement of lipid membranes to be obtained. Moreover, using a lipid spin label with the nitroxide fragment attached to its head group or a hydrocarbon chain at different positions also enables data about molecular dynamics and structure at different membrane depths to be obtained. Thus, the method can be used to investigate not only the lateral organization of the membrane (i.e., the presence of membrane domains and phases), but also the depth-dependent membrane structure and dynamics, and, hence, the membrane properties in three dimensions.

Carbon Nanodots from an In Silico Perspective

Carbon nanodots (CNDs) are the latest and most shining rising stars among photoluminescent (PL) nanomaterials. These carbon-based surface-passivated nanostructures compete with other related PL materials, including traditional semiconductor quantum dots and organic dyes, with a long list of benefits and emerging applications. Advantages of CNDs include tunable inherent optical properties and high photostability, rich possibilities for surface functionalization and doping, dispersibility, low toxicity, and viable synthesis (top-down and bottom-up) from organic materials. CNDs can be applied to biomedicine including imaging and sensing, drug-delivery, photodynamic therapy, photocatalysis but also to energy harvesting in solar cells and as LEDs. More applications are reported continuously, making this already a research field of its own. Understanding of the properties of CNDs requires one to go to the levels of electrons, atoms, molecules, and nanostructures at different scales using modern molecular modeling and to correlate it tightly with experiments. This review highlights different in silico techniques and studies, from quantum chemistry to the mesoscale, with particular reference to carbon nanodots, carbonaceous nanoparticles whose structural and photophysical properties are not fully elucidated. The role of experimental investigation is also presented. Hereby, we hope to encourage the reader to investigate CNDs and to apply virtual chemistry to obtain further insights needed to customize these amazing systems for novel prospective applications.

Generation and Computational Characterization of a Complex Staphylococcus aureus Lipid Bilayer

Studies indicate a crucial cell membrane role in the antibiotic resistance of Staphylococcus aureus. To simulate its membrane structure and dynamics, a complex molecular-scale computational representation of the S. aureus lipid bilayer was developed. Phospholipid types and their amounts were optimized by reverse Monte Carlo to represent characterization data from the literature, leading to 19 different phospholipid types that combine three headgroups [phosphatidylglycerol, lysyl-phosphatidylglycerol (LPG), and cardiolipin] and 10 tails, including iso- and anteiso-branched saturated chains. The averaged lipid bilayer thickness was 36.7 Å, and area per headgroup was 67.8 Ų. Phosphorus and nitrogen density profiles showed that LPG headgroups tended to be bent and oriented more parallel to the bilayer plane. The water density profile showed that small amounts reached the membrane center. Carbon density profiles indicated hydrophobic interactions for all lipids in the middle of the bilayer. Bond vector order parameters along each tail demonstrated different C-H ordering even within distinct lipids of the same type; however, all tails followed similar trends in average order parameter. These complex simulations further revealed bilayer insights beyond those attainable with monodisperse, unbranched lipids. Longer tails often extended into the opposite leaflet. Carbon at and beyond a branch showed significantly decreased ordering compared to carbon in unbranched tails; this feature arose in every branched lipid. Diverse tail lengths distributed these disordered methyl groups throughout the middle third of the bilayer. Distributions in mobility and ordering reveal diverse properties that cannot be obtained with monodisperse lipids.

Fast Prediction of Binding Affinities of SARS-CoV-2 Spike Protein and Its Mutants with Antibodies through Intermolecular Interaction Modeling-Based Machine Learning

Since the introduction of the novel SARS-CoV-2 virus (COVID-19) in late 2019, various new variants have appeared with mutations that confer resistance to the vaccines and monoclonal antibodies that were developed in response to the wild-type virus. As we continue through the pandemic, an accurate and efficient methodology is needed to help predict the effects certain mutations will have on both our currently produced therapeutics and those that are in development. Using published cryo-electron microscopy and X-ray crystallography structures of the spike receptor binding domain region with currently known antibodies, in the present study, we created and cross-validated an intermolecular interaction modeling-based multi-layer perceptron machine learning approach that can accurately predict the mutation-caused shifts in the binding affinity between the spike protein (wild-type or mutant) and various antibodies. This validated artificial intelligence (AI) model was used to predict the binding affinity (Kd) of reported SARS-CoV-2 antibodies with various variants of concern, including the most recently identified "Deltamicron" (or "Deltacron") variant. This AI model may be employed in the future to predict the Kd of developed novel antibody therapeutics to overcome the challenging antibody resistance issue and develop structural bases for the effects of both current and new mutants of the spike protein. In addition, the similar AI strategy and approach based on modeling of the intermolecular interactions may be useful in development of machine learning models predicting binding affinities for other protein-protein binding systems, including other antibodies binding with their antigens.

Procyanidin C1 Location, Interaction, and Aggregation in Two Complex Biomembranes

Procyanidins are known for their many benefits to human health and show a plethora of biological effects. One of the most important procyanidin is the procyanidin trimer C1 (PC1). Due to its relatively high lipid–water partition coefficient, the properties of PC1 could be attributed to its capability to interact with the biomembrane, to modulate its structure and dynamics, and to interact with lipids and proteins, however, its biological mechanism is not known. We have used all-atom molecular dynamics in order to determine the position of PC1 in complex membranes and the presence of its specific interactions with membrane lipids, having simulated a membrane mimicking the plasma membrane and another mimicking the mitochondrial membrane. PC1 has a tendency to be located at the membrane interphase, with part of the molecule exposed to the water solvent and part of it reaching the first carbons of the hydrocarbon chains. It has no preferred orientation, and it completely excludes the CHOL molecule. Remarkably, PC1 has a tendency to spontaneously aggregate, forming high-order oligomers. These data suggest that its bioactive properties could be attributed to its membranotropic effects, which therefore supports the development of these molecules as therapeutic molecules, which would open new opportunities for future medical advances.

Phospholipid-porphyrin conjugates: deciphering the driving forces behind their supramolecular assemblies

Phospholipid-porphyrin conjugates (PL-Por) are nowadays considered as a unique class of building blocks that can self-assemble into supramolecular structures that possess multifunctional properties and enhanced optoelectronics characteristics compared to their disassembled counterparts. However, despite their versatile properties, little is known about the impact of the packing parameter of PL-Por conjugates on their assembling mechanism and their molecular organization inside these assemblies. To gain a better understanding on their assembling properties, we synthesized two new series of PL-Por conjugates with different alkyl sn2-chain lengths linked via an amide bond to either pheophorbide-a (Ph_xLPC) or pyropheophorbide-a (Pyr_xLPC). By combining a variety of experimental techniques with molecular dynamics (MD) simulations, we investigated both the assembling and optical properties of the PL-Por either self-assembled or when incorporated into lipid bilayers. We demonstrated that independently of the linker length, Ph_xLPC assembled into closed ovoid structures, whereas Pyr_xLPC formed rigid open sheets. Interestingly, Pyr_xLPC assemblies displayed a significant red shift and narrowing of the Q-band indicating the formation of ordered J-aggregates. The MD simulations highlighted the central role of the interaction between porphyrin cores rather than the length difference between the two phospholipid chains in controlling the structure of the lipid bilayer membranes and thus their optical properties. Indeed, while Ph_xLPC have the tendency to form inter-leaflet π-stacked dimers, Pyr_xLPC conjugates formed dimers within the same leaflet. Altogether, this work could be used as guidelines for the design of new PL-Por conjugates that self-assemble into bilayer-like supramolecular structures with tunable morphology and optical properties.

Active site geometry stabilization of a presenilin homolog by the lipid bilayer promotes intramembrane proteolysis

Cleavage of membrane proteins in the lipid bilayer by intramembrane proteases is crucial for health and disease. Although different lipid environments can potently modulate their activity, how this is linked to their structural dynamics is unclear. Here, we show that the carboxy-peptidase-like activity of the archaeal intramembrane protease PSH, a homolog of the Alzheimer’s disease-associated presenilin/γ-secretase is impaired in micelles and promoted in a lipid bilayer. Comparative molecular dynamics simulations revealed that important elements for substrate binding such as transmembrane domain 6a of PSH are more labile in micelles and stabilized in the lipid bilayer. Moreover, consistent with an enhanced interaction of PSH with a transition-state analog inhibitor, the bilayer promoted the formation of the enzyme’s catalytic active site geometry. Our data indicate that the lipid environment of an intramembrane protease plays a critical role in structural stabilization and active site arrangement of the enzyme-substrate complex thereby promoting intramembrane proteolysis.

Cutting proteins into pieces is a crucial process in the cell, allowing several important processes to take place, including cell differentiation (which allows cells to develop into specific types), cell death, protein quality control, or even where in the cell a protein will end up. However, the specialized proteins that carry out this task, known as proteases, can also be involved in the development of disease. For example, in the brain, a protease called γ-secretase cuts up the amyloid-β protein precursor, producing toxic forms of amyloid-β peptides that are widely believed to cause Alzheimer’s disease.

Proteases like γ-secretase carry out their role in the membrane, the layer of fats (also known as lipids) that forms the outer boundary of the cell. The environment in this area of the cell can influence the activity of proteases, but it is poorly understood how this happens.

One way to address this question would be to compare the activity of γ-secretase in the lipid environment of the membrane to its activity when it is entirely surrounded by different molecules, such as detergent molecules. Unfortunately, γ-secretase is not active when it is removed from its lipid environment by a detergent, making it difficult to perform this comparison. To overcome this issue, Feilen et al. chose to study PSH, a protease similar to γ-secretase that produces the same amyloid-β peptides but remains active in detergent.

When Feilen et al. mixed PSH with lipid molecules like those found in the membrane and amyloid-β precursor protein, PSH produced amyloid-β peptides including those that are thought to cause Alzheimer’s. However, when a detergent was substituted for the lipid molecules this led to longer amyloid-β peptides than usual, indicating that PSH was not able to cut proteins as effectively. The change in environment appeared to reduce PSH’s ability to progressively trim small segments from the peptides.

Computer modelling of the protease’s structure in lipids versus detergent supported the experimental findings: the model predicted that the areas of PSH important for recognizing and cutting other proteins would be more stable in the membrane compared to the detergent.

These results indicate that the cell membrane plays a vital role in the stability of the active regions of proteases that are cleaving in this environment. In the future, this could help to better understand how changes to the lipid molecules in the membrane may contribute to the activity of γ-secretase and its role in Alzheimer’s disease.

Generalized Methodology for the Quick Prediction of Variant SARS-CoV-2 Spike Protein Binding Affinities with Human Angiotensin-Converting Enzyme II

Variants of the SARS-CoV-2 virus continue to remain a threat 2 years from the beginning of the pandemic. As more variants arise, and the B.1.1.529 (Omicron) variant threatens to create another wave of infections, a method is needed to predict the binding affinity of the spike protein quickly and accurately with human angiotensin-converting enzyme II (ACE2). We present an accurate and convenient energy minimization/molecular mechanics Poisson-Boltzmann surface area methodology previously used with engineered ACE2 therapeutics to predict the binding affinity of the Omicron variant. Without any additional data from the variants discovered after the publication of our first model, the methodology can accurately predict the binding of the spike/ACE2 variant complexes. From this methodology, we predicted that the Omicron variant spike has a Kd of ∼22.69 nM (which is very close to the experimental Kd of 20.63 nM published during the review process of the current report) and that spike protein of the new "Stealth" Omicron variant (BA.2) will display a Kd of ∼12.9 nM with the wild-type ACE2 protein. This methodology can be used with as-yet discovered variants, allowing for quick determinations regarding the variant's infectivity versus either the wild-type virus or its variants.

Effect of Ionic Strength on Ibuprofenate Adsorption on a Lipid Bilayer of Dipalmitoylphosphatidylcholine from Molecular Dynamics Simulations

In this work, the free energy change in the process of transferring ibuprofenate from the bulk solution to the center of a model of the dipalmitoylphosphatidylcholine bilayer at different NaCl concentrations was calculated. Two minima were found in the free energy profile: a local minimum, located in the vicinity of the membrane, and the global free energy minimum, found near the headgroup region. The downward shift of free energy minima with increasing NaCl concentration is consistent with the results of previous works. Conversely, the upward shift of the free energy maximum with increasing ionic strength is due to the competition of sodium ions and lipids molecules to coordinate with ibuprofenate and neutralize its charge. In addition, normal molecular dynamics simulations were performed to study the effects of the ibuprofenate on the lipid bilayer and in the presence of a high ibuprofenate concentration. The effect of ionic strength on the properties of the lipid bilayer and on lipid-drug interactions was analyzed. The area per lipid shrinking with increasing ionic strength, volume of lipids, and thickness of the bilayer is consistent with the experimental results. At a very high ibuprofenate concentration, the lipid bilayer dehydrates, and it consequently transforms into the gel phase, thus blocking the permeation.

Self-assembly of ultra-small-sized carbon nanoparticles in lipid membrane disrupts its integrity

Although nanomaterials are widely studied in biomedical applications, the major concern of nanotoxicity still exists. Therefore, numerous studies have been conducted on the interactions of various biomolecules with various types of nanomaterials, including carbon nanotubes, graphene, fullerene etc. However, the size effect of nanomaterials is poorly documented, especially ultra-small particles. Here, the interactions of the smallest carbon nanoparticle (NP), C₂₈, with the cell membrane were studied using molecular dynamics (MD) simulations. The results show that similar to fullerene C₆₀, the C₂₈ NPs can self-assemble into stable clusters in water. Inside the membrane, the C₂₈ NPs are more prone to aggregate to form clusters than C₆₀ NPs. The reason for C₂₈ aggregation is characterized by the potential of mean force (PMF) and can be explained by the polarized nature of C₂₈ NPs while the acyl chains of lipids are nonpolar. At the C₂₈ cluster regions, the thickness of the membrane is significantly reduced by the C₂₈ aggregation. Accordingly, the membrane loses its structural integrity, and translocation of water molecules through it was observed. Thus, these results predict a stronger cytotoxicity to cells than C₆₀ NPs. The present findings might shed light on the understanding of the cytotoxicity of NPs with different sizes and would be helpful for the potential biomedical applications of carbon NPs, especially as antibacterial agents.

Molecular insights into the differential efflux mechanism of Rv1634 protein, a multidrug transporter of major facilitator superfamily in Mycobacterium tuberculosis

Currently, multidrug-resistant tuberculosis (MDR-TB) is a public health crisis and a major health security threat globally. In Mycobacterium tuberculosis (Mtb), major facilitator superfamily (MFS) is the largest group of secondary active transporters. Along with the transport of their natural substrates, MFS proteins were involved in a drug efflux mechanism that ultimately lead to resistance against available anti-TB drugs in Mtb. In the present study, the three-dimensional structure model of an MFS protein, Rv1634, a probable multidrug transporter from Mtb, was generated using homology modeling. The protein structure model was found in inward-open conformation having 14 transmembrane helices. In addition, a central transport channel was deduced across the protein, and a single binding pocket was identified halfway through the central cavity by structural alignment with the homologous protein structures. Further, Rv1634 protein was studied based on the differential structural behavior of apo and ligand-bound forms. All the protein systems were inserted into a phospholipid bilayer to characterize the conformational dynamics of the protein using molecular dynamics (MD) simulations. Detailed analysis of the MD trajectories showed the diverse substrate specificity of the binding pocket for the antibiotics that caused differential movement in the ciprofloxacin and norfloxacin, to which Mtb strains have now become resistant. The expulsion of the drugs outside the bacterial cell occurs through the alternating-access mechanism of N and C-terminal domains, which is intriguing and essential to the understanding the drug resistance mechanism in pathogenic bacteria.

Auxin transport mechanism of membrane transporter encoded by AEC gene of Bacillus licheniformis isolated from metagenome of Tapta Kund Hotspring of Uttrakhand, India

Members of group Bacillus are most widely occurring microbes in agricultural soil and they affect crop health in various ways. They directly stimulate plant growth either by augmenting nutrients availability, invigorating plants' defence mechanisms; repressing soil-borne phytopathogens or by producing growth-regulating hormones like auxins and cytokinins. It is a well known fact that indole-3- acetic acid (a type of auxin) is a vital biologically active phytohormone excreted by certain Bacillus species, but its molecular mechanism has not yet been described. In this study, the auxin efflux carrier gene is isolated from the metagenome of the Tapta Kund hot spring, Uttrakhand, India. In addition, auxin efflux carrier (AEC) transporter protein of Bacillus licheniformis is modeled and the 318 amino acid residues long protein was found homologous to the apical sodium-dependent bile acid transporter (ASBT) of Yersinia frederiksnii, with 10 transmembrane segments (TM1-10) split into different domains: a panel domain defined by TM1, 2, 6 and 7; and a core domain defined by TM3-5 and 8-10. Finally, the predicted Bacillus licheniformis AEC protein has also been phylogenetically evaluated and its detailed molecular transport mechanism was worked out using molecular dynamics simulation analysis. Conclusively, this study demonstrates the efflux mechanism of the substrate, Indole 3- acetic acid by AEC transporter protein.

Fast Prediction of Binding Affinities of the SARS-CoV-2 Spike Protein Mutant N501Y (UK Variant) with ACE2 and Miniprotein Drug Candidates

A recently identified variant of SARS-CoV-2 virus, known as the United Kingdom (UK) variant (lineage B.1.1.7), has an N501Y mutation on its spike protein. SARS-CoV-2 spike protein binds with angiotensin-converting enzyme 2 (ACE2), a key protein for the viral entry into the host cells. Here, we report an efficient computational approach, including the simple energy minimizations and binding free energy calculations, starting from an experimental structure of the binding complex along with experimental calibration of the calculated binding free energies, to rapidly and reliably predict the binding affinities of the N501Y mutant with human ACE2 (hACE2) and recently reported miniprotein and hACE2 decoy (CTC-445.2) drug candidates. It has been demonstrated that the N501Y mutation markedly increases the ACE2-spike protein binding affinity (Kd) from 22 to 0.44 nM, which could partially explain why the UK variant is more infectious. The miniproteins are predicted to have ∼10,000- to 100,000-fold diminished binding affinities with the N501Y mutant, creating a need for design of novel therapeutic candidates to overcome the N501Y mutation-induced drug resistance. The N501Y mutation is also predicted to decrease the binding affinity of a hACE2 decoy (CTC-445.2) binding with the spike protein by ∼200-fold. This convenient computational approach along with experimental calibration may be similarly used in the future to predict the binding affinities of potential new variants of the spike protein.

Functional Group Distributions, Partition Coefficients, and Resistance Factors in Lipid Bilayers Using Site Identification by Ligand Competitive Saturation

Small molecules such as metabolites and drugs must pass through the membrane of the cell, a barrier primarily comprising phospholipid bilayers and embedded proteins. To better understand the process of passive diffusion, knowledge of the ability of various functional groups to partition across bilayers and the associated energetics would be of utility. In the present study, the site identification by ligand competitive saturation (SILCS) methodology has been applied to sample the distributions of a diverse set of chemical solutes representing the functional groups of small molecules across phospholipid bilayers composed of 0.9:0.1 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/cholesterol and a mixture of 0.52:0.18:0.3 1,2-dioleoyl-sn-glycero-3-phospho-l-serine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol used in parallel artificial membrane permeability assay experiments. A combination of oscillating chemical potential grand canonical Monte Carlo and molecular dynamics in the SILCS simulations was applied to achieve solute sampling through the bilayers and surrounding aqueous environment from which the distribution of solutes and the functional groups they represent were obtained. Results show differential distribution of aliphatic versus aromatic groups with the former having increased sampling in the center of the bilayers versus in the region of the glycerol linker for the latter. Variations in the distribution of different polar groups are evident, with large differences between negative acetate and positive methylammonium with accumulation of the polar-neutral and acetate solutes above the bilayer head groups. Conversion of the distributions to absolute free energies allows for a detailed understanding of energetics of functional groups in different regions of the bilayers and for calculation of absolute free-energy profiles of multifunctional drug-like molecules across the bilayers from which partition coefficients and resistance factors suitable for insertion into the homogenous solubility-diffusion equation for calculation of permeability were obtained. Comparisons of the calculated bilayer/solution partition coefficients with 1-octanol/water experimental data for both drug-like molecules and the solutes show overall good agreement, validating the calculated distributions and associated absolute free-energy profiles.

Effect of pH and temperature on the aggregation behaviour of dirhamnolipid biosurfactant. An experimental and molecular dynamics study

HYPOTHESIS

Pseudomonas aeruginosa dirhamnolipid (diRL) has been shown to form aggregates of different size and structure, under various conditions. Due to the presence of a carboxyl group in the molecule, it is expected that pH would strongly affect this aggregation behaviour. In addition, preliminary observations of temperature-induced changes in the states of aggregation of diRL supported the need of further investigation.

EXPERIMENTS

A systematic experimental study, using differential scanning calorimetry (DSC), small-angle Xray diffraction (SAXD), and Fourier-transform infrared spectroscopy (FTIR), has been carried out to characterize pH and temperature driven changes in the aggregation behavior of diRL biosurfactant. Molecular dynamics (MD) simulations, supported by the experimental results, allowed depicting molecular details on formation of diRL membranes and other aggregated structures under various physicochemical conditions.

FINDINGS

DiRL could adopt fairly organized multilayered structures (membranes) at low pH and temperature, which became highly disordered upon increasing either of these parameters. The effect of pH on the gauche/all-trans conformer ratio of the diRL acyl chains was not of significance, whereas temperature-induced effects were observed. For the first time it is described that diRL underwent an endothermic thermotropic transition with T_c = 34 °C as observed by DSC, at pH 4.5 (protonated diRL), but not at pH 7.4 (unprotonated diRL). FTIR confirmed these findings, showing a significant additional disordering of the all-trans acyl chains upon increasing temperature around that same value in the protonated form, an effect not observed for the dissociated form of the biosurfactant. In addition, at pH 7.4, changing temperature did not modify the hydration state of the polar moiety of diRL, whereas at pH 4.5 a significant decrease in the hydration state around 34 °C took place. SAXD data showed that protonated diRL formed multilayered structures at 20 °C, which converted into poorly correlated layers at 50 °C. MD simulations supported these findings, showing that the membrane-like structures formed by protonated diRL at 20 °C became unstable at higher temperatures, tending to form other structures, which could be micelles or other type of layered structures, whereas the negatively charged form of diRL organized in micelle-type aggregates in the whole range of temperature under study.

Influence of the gallate moiety on the interactions between green tea polyphenols and lipid membranes elucidated by molecular dynamics simulations

Previous studies suggested that naturally occurring EGCG primarily acted on the bacterial cell membrane then damaged the membrane and the gallate moiety in EGCG was very important to its anti-bacterial activity. However, the detailed mechanisms were still poorly understood. In this paper, EGCG and EGC were selected to study the great contribution of gallate moiety on the anti-bacterial activities of polyphenols. The results indicated that EGCG could penetrate deeper into the POPG lipid bilayer and possess more potent structure-perturbing potency on the POPG lipid bilayer than EGC. We also found that EGCG had the ability to form hydrogen bonds with the deeper inside oxygen atoms in the POPG lipid bilayer and the gallate moiety was the key functional group for EGCG forming hydrogen bonds with the POPG lipid bilayer. Moreover, results from the binding free energy analysis demonstrated that the gallate moiety made great contribution to the high affinity between EGCG and the POPG lipid bilayer. We believed that these findings could yield useful insights into the influence mechanisms of gallate moiety on the anti-bacterial activities of polyphenols.

Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion

The development of the CHARMM lipid force field (FF) can be traced back to the early 1990s with its current version denoted CHARMM36 (C36). The parametrization of C36 utilized high-level quantum mechanical data and free energy calculations of model compounds before parameters were manually adjusted to yield agreement with experimental properties of lipid bilayers. While such manual fine-tuning of FF parameters is based on intuition and trial-and-error, automated methods can identify beneficial modifications of the parameters via their sensitivities and thereby guide the optimization process. This work introduces a semi-automated approach to reparametrize the CHARMM lipid FF with consistent inclusion of long-range dispersion through the Lennard-Jones particle-mesh Ewald (LJ-PME) approach. The optimization method is based on thermodynamic reweighting with regularization with respect to the C36 set. Two independent optimizations with different topology restrictions are presented. Targets of the optimizations are primarily liquid crystalline phase properties of lipid bilayers and the compression isotherm of monolayers. Pair correlation functions between water and lipid functional groups in aqueous solution are also included to address headgroup hydration. While the physics of the reweighting strategy itself is well-understood, applying it to heterogeneous, complex anisotropic systems poses additional challenges. These were overcome through careful selection of target properties and reweighting settings allowing for the successful incorporation of the explicit treatment of long-range dispersion, and we denote the newly optimized lipid force field as C36/LJ-PME. The current implementation of the optimization protocol will facilitate the future development of the CHARMM and related lipid force fields.

Tuning the transdermal transport by application of external continuous electric field: a coarse-grained molecular dynamics study

The control of skin permeability to specific substances (e.g., medications, vitamins, and nutrients) through stratum corneum is a challenge. Iontophoresis is an option in spite of the lack of a detailed understanding of the underlying molecular mechanism. In the present work, the simulations concerning application of an external continuous electric field to stratum corneum, in a range of low intensity (0-24 mV nm^-1), were carried out using the coarse-grained molecular dynamics approach. Using a set of random seed replicas of the starting configuration, we observed that in the range of electric field intensity of 22-23 mV nm^-1, water-rich lipid vesicles were formed in 20% of cases. Pores appeared in the remaining 80%. We argue that lipids undergo fast re-orientations under electric field inducing mechanical instability, which originates the pores. We presented a simple electrostatic model to interpret the results where the mismatch between electrical permittivities of the membrane and external media and the gradient of the local electric field in the membrane surface, govern the time scales and electric fields for vesicle formation. Our results indicate that just 10% difference between electrical permittivities of the membrane and external media decreases 1/6 the minimal time required for vesicle formation. The minimal electric field required decreases 10 times. The control and tunning of formation of biologically compatible vesicles, capable of transporting substances under low-intensity electric fields, has a promising application in fields such as drug therapy and dermo-cosmetics allowing the use of hydrophilic substances in dermal applications.

Influence of Lipid Composition on the Insertion Process of Glyphosate into Membranes: A Thermodynamic Study

In this work, molecular dynamics simulations were applied to investigate the influence of lipid composition of the model membrane on the insertion of glyphosate (in its charged state, GLYP2-). The profiles of free energy, entropy and enthalpy were obtained through umbrella sampling calculations, for lipid bilayers composed by only 1,2-dipalmitoyl-sn-glycerol-3-phosphocholine (DPPC), only 1,2-dipalmitoyl-sn-glycerol-3-phosphoserine (DPPS) or a symmetric binary mixture of DPPC and DPPS. In general, the location, the values of minima and maxima of the free energy, and the trend of free energy profiles are influenced by the lipid composition of the lipid bilayer. The driving force in the glyphosate insertion process depends on the lipid composition of the membrane model. If the lipid bilayer is composed solely of DPPS or DPPC, GLYP2- insertion is driven by a favorable enthalpic change. However, if the membrane is composed of a mixture of both lipids, this process is driven by a favorable entropic change. In the lipid bilayer containing DPPS, the glyphosate was found to penetrate hydrated and coordinated with Na+ ions, in contrast to the pure zwitterionic lipid bilayer which penetrated only hydrated. This effect is independent of the concentration of sodium ions present in the bulk solution.

Aggregation of 25-hydroxycholesterol in a complex biomembrane. Differences with cholesterol

25-Hydroxycholesterol (25HC), one of the most important oxysterol molecules, can be used by cells to fight bacterial and viral infections but the mechanism that defines its biological effects are unknown. Using molecular dynamics, we have aimed to describe the orientation and location of 25HC in the membrane as well as the interactions it might have with lipids. We have studied two complex model membrane systems, one similar to the late endosome membrane and the other one to the plasma membrane. Our results reinforce that 25HC is inserted in the membrane in a relative stable location similar to but not identical to cholesterol. 25HC fluctuates in the membrane to a much greater degree than cholesterol, but the effect of 25HC on the phospholipid order parameters is not significantly different. One of the most notable facts about 25HC is that, unlike cholesterol, this molecule tends to aggregate, forming dimers, trimers and higher-order aggregates. These aggregates are formed spontaneously through the formation of hydrogen bonds between the two 25HC atoms, the formation of hydrogen bonds being independent of the studied system. Remarkably, no contacts or hydrogen bonds are observed between 25HC and cholesterol molecules, as well as between cholesterol molecules themselves at any time. It would be conceivable that 25HC, by forming high order aggregates without significantly altering the membrane properties, would modify the way proteins interact with the membrane and henceforth form a true innate antiviral molecule.

Energetics and structures of the tilted phases of fatty acid Langmuir monolayers

Langmuir monolayers are monomolecular deep films composed of amphiphilic molecules which are typically confined to a water/air interface in a bi-dimensional structure. Due to the important applications in many research areas, they have been studied for many years. Their phase diagrams present several condensed phases, showing untilted or tilted structures at low values of surface pressure. In this paper, we present a novel density functional study on tilted phases of different fatty acid Langmuir monolayers. By means of this study, a further understanding of the physical chemistry properties and the nature of the formation of tilted monolayers can be achieved. Our calculations reveal that, regardless of the number of carbon atoms which form the apolar chain, the transversal (or conventional in the case of untilted phases) unit cell shows similar dimensions, ca. 4.9 × 6.8 Å, which is in fair agreement with the range of the observed data. The energy variation of the unit cell as a function of the inclination of the molecules, reveals an abrupt increase in values larger than 45° and 36° for NN- and NNN-tilt, respectively, in fair agreement with the experimental observation of L2h (NN) and L2' (NNN) phases of fatty acids. All of the fatty acids explored (from 10 to 19 carbon atoms) yield similar results. Finally, the energetics and structural changes of the monolayer along the variation of the area per molecule, obtained by enlarging in a-, b- or both axes of the untilted unit cell, have been explored. This study reveals that the untilted phases are energetically more stable at low values of area per molecule (high surface concentration), as it is expected. When the area per molecule values are increased, tilted phases (along NN or NNN-direction) with b/a ratio typical of herringbone (HB) or pseudo-herringbone (PHB) structures are found in the lowest energy configurations, which depend on how the distortion of the untilted unit cell is performed. For example, HB structures are the most stable when the molecules tilt along the enlarged axis of the untilted unit cell (a or b), meanwhile unit cell structures characteristic of PHB configurations occur in the opposite cases and at larger values of the area per molecule (low surface concentrations). All these predictions are in good agreement with the GIXD observations of the different phases of the phase diagram of fatty acid Langmuir monolayers.

Saddle-curvature instability of lipid bilayer induced by amphipathic peptides: a molecular model

Amphipathic peptides that partition into lipid bilayers affect the curvature elastic properties of their host. Some of these peptides are able to shift the Gaussian modulus to positive values, thus triggering an instability with respect to the formation of saddle curvatures. To characterize the generic aspects of the underlying mechanism, we employ a molecular lipid model that accounts for the interfacial tension between the polar and apolar regions of the membrane, for interactions between the lipid headgroups, and for the energy to stretch or compress the hydrocarbon chains. Peptides are modeled as cylinders that partition into the host membrane in a parallel orientation where they diminish the space available to the lipid headgroups and chains. The penetration depth into the membrane is determined by the angular size of the peptide's hydrophilic region. We demonstrate that only peptides with a small angular size of their hydrophilic region have an intrinsic tendency to render the Gaussian modulus more positive, and we identify conditions at which the Gaussian modulus adopts a positive sign upon increasing the peptide concentration. Our model allows us to also incorporate electrostatic interactions between cationic peptides and anionic lipids on the level of the linear Debye-Hückel model. We show that electrostatic interactions tend to shift the Gaussian modulus toward more positive values. Steric and electrostatic lipid-peptide interactions jointly decrease the effective interaction strength in the headgroup region of the host membrane thus suggesting a generic mechanisms of how certain amphipathic peptides are able to induce the formation of saddle curvatures.

Bicontinuous cubic phases in biological and artificial self-assembled systems

Nature has created innumerable life forms with miraculous hierarchical structures and morphologies that are optimized for different life events through evolution over billions of years. Bicontinuous cubic structures, which are often described by triply periodic minimal surfaces (TPMSs) and their constant mean curvature (CMC)/parallel surface companions, are of special interest to various research fields because of their complex form with unique physical functionalities. This has prompted the scientific community to fully understand the formation, structure, and properties of these materials. In this review, we summarize and discuss the formation mechanism and relationships of the relevant biological structures and the artificial self-assembly systems. These structures can be formed through biological processes with amazing regulation across a great length scales; nevertheless, artificial construction normally produces the structure corresponding to the molecular size and shape. Notably, the block copolymeric system is considered to be an applicable and attractive model system for the study of biological systems due to their versatile design and rich phase behavior. Some of the phenomena found in these two systems are compared and discussed, and this information may provide new ideas for a comprehensive understanding of the relationship between molecular shape and resulting interface curvature and the self-assembly process in living organisms. We argue that the co-polymeric system may serve as a model to understand these biological systems and could encourage additional studies of artificial self-assembly and the creation of new functional materials.

Location, Orientation and Aggregation of Bardoxolone-ME, CDDO-ME, in a Complex Phospholipid Bilayer Membrane

Bardoxolone methyl (CDDO-Me), a synthetic derivative of the naturally occurring triterpenoid oleanolic acid, displays strong antioxidant, anticancer and anti-inflammatory activities, according to different bibliographical sources. However, the understanding of its molecular mechanism is missing. Furthermore, CDDO-Me has displayed a significant cytotoxicity against various types of cancer cells. CDDO-Me has a noticeable hydrophobic character and several of its effects could be attributed to its ability to be incorporated inside the biological membrane and therefore modify its structure and specifically interact with its components. In this study, we have used full-atom molecular dynamics to determine the location, orientation and interactions of CDDO-Me in phospholipid model membranes. Our results support the location of CDDO-Me in the middle of the membrane, it specifically orients so that the cyano group lean towards the phospholipid interface and it specifically interacts with particular phospholipids. Significantly, in the membrane the CDDO-Me molecules specifically interact with POPE and POPS. Moreover, CDDO-Me does not aggregates in the membrane but it forms a complex conglomerate in solution. The formation of a complex aggregate in solution might hamper its biological activity and therefore it should be taken into account when intended to be used in clinical assays. This work should aid in the development of these molecules opening new avenues for future therapeutic developments.