Statistical Convergence of Equilibrium Properties in Simulations of Molecular Solutes Embedded in Lipid Bilayers

Exploring the potential of silymarin-loaded nanovesicles as an effective drug delivery system for cancer therapy: in vivo, in vitro, and in silico experiments

We aimed to perform a comprehensive study on the development and characterization of silymarin (Syl)-loaded niosomes as potential drug delivery systems. The results demonstrate significant novelty and promising outcomes in terms of morphology, size distribution, encapsulation efficiency, in vitro release behavior, free energy profiles of Syl across the niosome bilayer, hydrogen bonding interactions, antimicrobial properties, cytotoxicity, and in vivo evaluations. The physical appearance, size, and morphology assessment of free niosomes and Syl-loaded niosomes indicated stable and well-formed vesicular structures suitable for drug delivery. Transmission electron microscopy (TEM) analysis revealed spherical shapes with distinct sizes for each formulation, confirming uniform distribution. Dynamic light scattering (DLS) analysis confirmed the size distribution results with higher polydispersity index for Syl-loaded niosomes. The encapsulation efficiency of Syl in the niosomes was remarkable at approximately 91%, ensuring protection and controlled release of the drug. In vitro release studies showed a sustained release profile for Syl-loaded niosomes, enhancing therapeutic efficacy over time. Free energy profiles analysis identified energy barriers hindering Syl permeation through the niosome bilayer, emphasizing challenges in drug delivery system design. Hydrogen bonding interactions between Syl and niosome components contributed to energy barriers, impacting drug permeability. Antimicrobial assessments revealed significant differences in inhibitory effects against S. aureus and E. coli. Cytotoxicity evaluations demonstrated the superior tumor-killing potential of Syl-loaded niosomes compared to free Syl. In vivo studies indicated niosome formulations' safety profiles in terms of liver and kidney parameters compared to bulk Syl, showcasing potential for clinical applications. Overall, this research highlights the promising potential of Syl-loaded niosomes as effective drug delivery systems with enhanced stability, controlled release, and improved therapeutic outcomes.

Insertion and Anchoring of HIV-1 Fusion Peptide into Complex Membrane Mimicking Human T-cell

A fundamental understanding of how HIV-1 envelope (Env) protein facilitates fusion is still lacking. The HIV-1 fusion peptide, consisting of 15 to 22 residues, is the N-terminus of the gp41 subunit of the Env protein. Further, this peptide, a promising vaccine candidate, initiates viral entry into target cells by inserting and anchoring into human immune cells. The influence of membrane lipid reorganization and the conformational changes of the fusion peptide during the membrane insertion and anchoring processes, which can significantly affect HIV-1 cell entry, remains largely unexplored due to the limitations of experimental measurements. In this work, we investigate the insertion of the fusion peptide into an immune cell membrane mimic through multiscale molecular dynamics simulations. We mimic the native T-cell by constructing a 9-lipid asymmetric membrane, along with geometrical restraints accounting for insertion in the context of gp41. To account for the slow timescale of lipid mixing while enabling conformational changes, we implement a protocol to go back and forth between atomistic and coarse-grained simulations. Our study provides a molecular understanding of the interactions between the HIV-1 fusion peptide and the T-cell membrane, highlighting the importance of conformational flexibility of fusion peptides and local lipid reorganization in stabilizing the anchoring of gp41 into the targeted host membrane during the early events of HIV-1 cell entry. Importantly, we identify a motif within the fusion peptide critical for fusion that can be further manipulated in future immunological studies.

Quantifying Induced Dipole Effects in Small Molecule Permeation in a Model Phospholipid Bilayer

The cell membrane functions as a semipermeable barrier that governs the transport of materials into and out of cells. The bilayer features a distinct dielectric gradient due to the amphiphilic nature of its lipid components. This gradient influences various aspects of small molecule permeation and the folding and functioning of membrane proteins. Here, we employ polarizable molecular dynamics simulations to elucidate the impact of the electronic environment on the permeation process. We simulated eight distinct amino-acid side chain analogs within a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer using the Drude polarizable force field (FF). Our approach includes both unbiased and umbrella sampling simulations. By using a polarizable FF, we sought to investigate explicit dipole responses in relation to local electric fields along the membrane normal. We evaluate molecular dipole moments, which exhibit variation based on their localization within the membrane, and compare the outcomes with analogous simulations using the nonpolarizable CHARMM36 FF. This comparative analysis aims to discern characteristic differences in the free energy surfaces of permeation for the various amino-acid analogs. Our results provide the first systematic quantification of the impact of employing an explicitly polarizable FF in this context compared to the fixed-charge convention inherent to nonpolarizable FFs, which may not fully capture the influence of the membrane dielectric gradient.

Leveraging Machine Learning-Guided Molecular Simulations Coupled with Experimental Data to Decipher Membrane Binding Mechanisms of Aminosterols

Understanding the molecular mechanisms of the interactions between specific compounds and cellular membranes is essential for numerous biotechnological applications, including targeted drug delivery, elucidation of the drug mechanism of action, pathogen identification, and novel antibiotic development. However, estimation of the free energy landscape associated with solute binding to realistic biological systems is still a challenging task. In this work, we leverage the Time-lagged Independent Component Analysis (TICA) in combination with neural networks (NN) through the Deep-TICA approach for determining the free energy associated with the membrane insertion processes of two natural aminosterol compounds, trodusquemine (TRO), and squalamine (SQ). These compounds are particularly noteworthy because they interact with the outer layer of neuron membranes, protecting them from the toxic action of misfolded proteins involved in neurodegenerative disorders, in both their monomeric and oligomeric forms. We demonstrate how this strategy could be used to generate an effective collective variable for describing solute absorption in the membrane and for estimating free energy landscape of translocation via on-the-fly probability enhanced sampling (OPES) method. In this context, the computational protocol allowed an exhaustive characterization of the aminosterol entry pathway into a neuron-like lipid bilayer. Furthermore, it provided accurate prediction of membrane binding affinities, in close agreement with the experimental binding data obtained by using fluorescently labeled aminosterols and large unilamellar vesicles (LUVs). The findings contribute significantly to our understanding of aminosterol entry pathways and aminosterol-lipid membrane interactions. Finally, the computational methods deployed in this study further demonstrate considerable potential for investigating membrane binding processes.

Interaction of Uperin Peptides with Model Membranes: Molecular Dynamics Study

The interaction of antimicrobial and amyloid peptides with cell membranes is a critical step in their activities. Peptides of the uperin family obtained from the skin secretion of Australian amphibians demonstrate antimicrobial and amyloidogenic properties. All-atomic molecular dynamics and an umbrella sampling approach were used to study the interaction of uperins with model bacterial membrane. Two stable configurations of peptides were found. In the bound state, the peptides in helical form were located right under the head group region in parallel orientation with respect to the bilayer surface. Stable transmembrane configuration was observed for wild-type uperin and its alanine mutant in both alpha-helical and extended unstructured forms. The potential of mean force characterized the process of peptide binding from water to the lipid bilayer and its insertion into the membrane, and revealed that the transition of uperins from the bound state to the transmembrane position was accompanied by the rotation of peptides and passes through the energy barrier of 4-5 kcal/mol. Uperins have a weak effect on membrane properties.

Simulated Tempering-Enhanced Umbrella Sampling Improves Convergence of Free Energy Calculations of Drug Membrane Permeation

Molecular dynamics simulations have been widely used to study solute permeation across biological membranes. The potential of mean force (PMF) for solute permeation is typically computed using enhanced sampling techniques such as umbrella sampling (US). For bulky drug-like permeants, however, obtaining converged PMFs remains challenging and often requires long simulation times, resulting in an unacceptable computational cost. Here, we augmented US with simulated tempering (ST), an extended-ensemble technique that consists in varying the temperature of the system along a pre-defined temperature ladder. Simulated tempering-enhanced US (STeUS) was employed to improve the convergence of PMF calculations for the permeation of methanol and three common drug molecules. To obtain sufficient sampling of the umbrella histograms, which were computed only from the ground temperature, we modified the simulation time fraction spent at the ground temperature between 1/K and 50%, where K is the number of ST temperature states. We found that STeUS accelerates convergence, when compared to standard US, and that the benefit of STeUS is system-dependent. For bulky molecules, for which standard US poorly converged, the application of ST was highly successful, leading to a more than fivefold accelerated convergence of the PMFs. For the small methanol solute, for which conventional US converges moderately, the application of ST is only beneficial if 50% of the STeUS simulation time is spent at the ground temperature. This study establishes STeUS as an efficient and simple method for PMF calculations, thereby strongly reducing the computational cost of routine high-throughput studies of drug permeability.

Modeling the Effect of Hydrophobicity on the Passive Permeation of Solutes across a Bacterial Model Membrane

Passive diffusion across biomembranes is an important mechanism of permeation for multiple drugs, including antibiotics. However, this process is frequently neglected while studying drug uptake and, in our view, warrants further investigation. Here, we apply molecular dynamics simulations to investigate the impact of changes in molecular hydrophobicity on the permeability of a series of inhibitors of the quorum sensing of Pseudomonas aeruginosa, previously discovered by us, across a membrane model. Overall, we show that permeation across this membrane model does not correlate with the molecule's hydrophobicity. We demonstrate that using a simple model for permeation, based on the difference between the maximum and minimum of the free energy profile, outperforms the inhomogeneous solubility-diffusion model, yielding a permeability ranking that better agrees with the experimental results, especially for hydrophobic permeants. The calculated differences in permeability could not explain differences in in bacterio activity. Nevertheless, substantial differences in molecular orientation along the permeation pathway correlate with the in bacterio activity, emphasizing the importance of analyzing, at an atomistic level, the permeation pathway of these solutes.

Lipid Composition Is Critical for Accurate Membrane Permeability Prediction of Cyclic Peptides by Molecular Dynamics Simulations

Cyclic peptides have attracted attention as a promising pharmaceutical modality due to their potential to selectively inhibit previously undruggable targets, such as intracellular protein-protein interactions. Poor membrane permeability is the biggest bottleneck hindering successful drug discovery based on cyclic peptides. Therefore, the development of computational methods that can predict membrane permeability and support elucidation of the membrane permeation mechanism of drug candidate peptides is much sought after. In this study, we developed a protocol to simulate the behavior in membrane permeation steps and estimate the membrane permeability of large cyclic peptides with more than or equal to 10 residues. This protocol requires the use of a more realistic membrane model than a single-lipid phospholipid bilayer. To select a membrane model, we first analyzed the effect of cholesterol concentration in the model membrane on the potential of mean force and hydrogen bonding networks along the direction perpendicular to the membrane surface as predicted by molecular dynamics simulations using cyclosporine A. These results suggest that a membrane model with 40 or 50 mol % cholesterol was suitable for predicting the permeation process. Subsequently, two types of membrane models containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 40 and 50 mol % cholesterol were used. To validate the efficiency of our protocol, the membrane permeability of 18 ten-residue peptides was predicted. Correlation coefficients of R > 0.8 between the experimental and calculated permeability values were obtained with both model membranes. The results of this study demonstrate that the lipid membrane is not just a medium but also among the main factors determining the membrane permeability of molecules. The computational protocol proposed in this study and the findings obtained on the effect of membrane model composition will contribute to building a schematic view of the membrane permeation process. Furthermore, the results of this study will eventually aid the elucidation of design rules for peptide drugs with high membrane permeability.

Interactions between Rhodamine Dyes and Model Membrane Systems—Insights from Molecular Dynamics Simulations

Background: rhodamines are dyes widely used as fluorescent tags in cell imaging, probing of mitochondrial membrane potential, and as P-glycoprotein model substrates. In all these applications, detailed understanding of the interaction between rhodamines and biomembranes is fundamental. Methods: we combined atomistic molecular dynamics (MD) simulations and fluorescence spectroscopy to characterize the interaction between rhodamines 123 and B (Rh123 and RhB, respectively) and POPC bilayers. Results: while the xanthene moiety orients roughly parallel to the membrane plane in unrestrained MD simulations, variations on the relative position of the benzoic ring (below the xanthene for Rh123, above it for RhB) were observed, and related to the structure of the two dyes and their interactions with water and lipids. Subtle distinctions were found among different ionization forms of the probes. Experimentally, RhB displayed a lipid/water partition coefficient more than two orders of magnitude higher than Rh123, in agreement with free energy profiles obtained from umbrella sampling MD. Conclusions: this work provided detailed insights on the similarities and differences in the behavior of bilayer-inserted Rh123 and RhB, related to the structure of the probes. The much higher affinity of RhB for the membranes increases the local concentration and explains its higher apparent affinity for P-glycoprotein reconstituted in model membranes.

Molecular dynamics simulations of ethanol permeation through single and double-lipid bilayers

Permeation of small molecules through membranes is a fundamental biological process, and molecular dynamics simulations have proven to be a promising tool for studying the permeability of membranes by providing a precise characterization of the free energy and diffusivity. In this study, permeation of ethanol through three different membranes of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS), PO-phosphatidylethanolamine (POPE), and PO-phosphatidylcholine (POPC) is studied. Permeabilities are calculated and compared with two different approaches based on Fick's first law and the inhomogeneous solubility-diffusion model. Microsecond simulation of double bilayers of these membranes provided a direct measurement of permeability by a flux-based counting method. These simulations show that a membrane of POPC has the highest permeability, followed by POPE and POPS. Due to the membrane-modulating properties of ethanol, the permeability increases as functions of concentration and saturation of the inner leaflet in a double bilayer setting, as opposed to the customary definition as a proportionality constant. This concentration dependence is confirmed by single bilayer simulations at different ethanol concentrations ranging from 1% to 18%, where permeability estimates are available from transition-based counting and the inhomogeneous solubility-diffusion model. We show that the free energy and diffusion profiles for ethanol lack accuracy at higher permeant concentrations due to non-Markovian kinetics caused by collective behavior. In contrast, the counting method provides unbiased estimates. Finally, the permeabilities obtained from single bilayer simulations are combined to represent natural gradients felt by a cellular membrane, which accurately models the non-equilibrium effects on ethanol permeability from single bilayer simulations in equilibrium.

Lipophilicity of Cationic Ligands Promotes Irreversible Adsorption of Nanoparticles to Lipid Bilayers

A mechanistic understanding of the influence of the surface properties of engineered nanomaterials on their interactions with cells is essential for designing materials for applications such as bioimaging and drug delivery as well as for assessing nanomaterial safety. Ligand-coated gold nanoparticles have been widely investigated because their highly tunable surface properties enable investigations into the effect of ligand functionalization on interactions with biological systems. Lipophilic ligands have been linked to adverse biological outcomes through membrane disruption, but the relationship between ligand lipophilicity and membrane interactions is not well understood. Here, we use a library of cationic ligands coated on 2 nm gold nanoparticles to probe the impact of ligand end group lipophilicity on interactions with supported phosphatidylcholine lipid bilayers as a model for cytoplasmic membranes. Nanoparticle adsorption to and desorption from the model membranes were investigated by quartz crystal microbalance with dissipation monitoring. We find that nanoparticle adsorption to model membranes increases with ligand lipophilicity. The effects of ligand structure on gold nanoparticle attachment were further analyzed using atomistic molecular dynamics simulations, which showed that the increase in ligand lipophilicity promotes ligand intercalation into the lipid bilayer. Together, the experimental and simulation results could be described by a two-state model that accounts for the initial attachment and subsequent conversion to a quasi-irreversibly bound state. We find that only nanoparticles coated with the most lipophilic ligands in our nanoparticle library undergo conversion to the quasi-irreversible state. We propose that the initial attachment is governed by interaction between the ligands and phospholipid tail groups, whereas conversion into the quasi-irreversibly bound state reflects ligand intercalation between phospholipid tail groups and eventual lipid extraction from the bilayer. The systematic variation of ligand lipophilicity enabled us to demonstrate that the lipophilicity of cationic ligands correlates with nanoparticle-bilayer adsorption and suggested that changing the nonpolar ligand R group promotes a mechanism of ligand intercalation into the bilayer associated with irreversible adsorption.

Rational design of antimicrobial peptides targeting Gram-negative bacteria

Membrane-targeting host antimicrobial peptides (AMPs) can kill or inhibit the growth of Gram-negative bacteria. However, the evolution of resistance among microbes poses a substantial barrier to the long-term utility of the host AMPs. Combining experiment and molecular dynamics simulations, we show that terminal carboxyl capping enhances both membrane insertion and antibacterial activity of an AMP called P1. Furthermore, we show that a bacterial strain with evolved resistance to this peptide becomes susceptible to P1 variants with either backbone capping or lysine-to-arginine substitutions. Our results suggest that cocktails of closely related AMPs may be useful in overcoming evolved resistance.

Molecular dynamics simulations of bacterial outer membrane lipid extraction: Adequate sampling?

The outer membrane of Gram-negative bacteria is almost exclusively composed of lipopolysaccharide in its outer leaflet, whereas the inner leaflet contains a mixture of phospholipids. Lipopolysaccharide diffuses at least an order of magnitude slower than phospholipids, which can cause issues for molecular dynamics simulations in terms of adequate sampling. Here, we test a number of simulation protocols for their ability to achieve convergence with reasonable computational effort using the MARTINI coarse-grained force-field. This is tested in the context both of potential of mean force (PMF) calculations for lipid extraction from membranes and of lateral mixing within the membrane phase. We find that decoupling the cations that cross-link the lipopolysaccharide headgroups from the extracted lipid during PMF calculations is the best approach to achieve convergence comparable to that for phospholipid extraction. We also show that lateral lipopolysaccharide mixing/sorting is very slow and not readily addressable even with Hamiltonian replica exchange. We discuss why more sorting may be unrealistic for the short (microseconds) timescales we simulate and provide an outlook for future studies of lipopolysaccharide-containing membranes.

Passive Diffusion of Ciprofloxacin and its Metalloantibiotic: A Computational and Experimental study

Fluoroquinolones (FQ) are antibiotics widely used in clinical practise, but the development of bacterial resistance to these drugs is currently a critical public health problem. In this context, ternary copper complexes of FQ (CuFQPhen) have been studied as a potential alternative. In this study, we compared the passive diffusion across the lipid bilayer of one of the most used FQ, ciprofloxacin (Cpx), and its ternary copper complex, CuCpxPhen, that has shown previous promising results regarding antibacterial activity and membrane partition. A combination of spectroscopic studies and molecular dynamics simulations were used and two different model membranes tested: one composed of anionic phospholipids, and the other composed of zwitterionic phospholipids. The obtained results showed a significantly higher membrane permeabilization activity, larger partition, and a more favourable free energy landscape for the permeation of CuCpxPhen across the membrane, when compared to Cpx. Furthermore, the computational results indicated a more favourable translocation of CuCpxPhen across the anionic membrane, when compared to the zwitterionic one, suggesting a higher specificity towards the former. These findings are important to decipher the influx mechanism of CuFQPhen in bacterial cells, which is crucial for the ultimate use of CuFQPhen complexes as an alternative to FQ to tackle multidrug-resistant bacteria.

Selecting Collective Variables and Free-Energy Methods for Peptide Translocation across Membranes

The selective permeability of cellular membranes is a crucial property for controlled transport into and out of cells. Molecules that can bypass the cellular machinery and spontaneously translocate across membranes could be used as therapeutics or drug carriers. Peptides are a prominent class of such molecules, which include natural and man-developed antimicrobial and cell-penetrating peptides. However, the necessary peptide properties for translocation remain elusive. Computer simulations could uncover these properties once we have a good collective variable (CV) that accurately describes the translocation process. Here, we developed a new CV, which includes a description of peptide insertion, local membrane deformation, and peptide internal degrees of freedom related to its charged groups. By comparison of CVs, we demonstrated that all these components are necessary for an accurate description of peptide translocation. Moreover, the advantages and disadvantages of three common methods for free-energy calculations with our CV were evaluated using the MARTINI coarse-grained model: umbrella sampling, umbrella sampling with replica exchange, and metadynamics. The developed CV leads to the reliable and effective calculation of the free energy of peptide translocation, and thus, it could be useful in the design of spontaneously translocating peptides.

Molecular simulations of lipid membrane partitioning and translocation by bacterial quorum sensing modulators

Quorum sensing (QS) is a bacterial communication process mediated by both native and non-native small-molecule quorum sensing modulators (QSMs), many of which have been synthesized to disrupt QS pathways. While structure-activity relationships have been developed to relate QSM structure to the activation or inhibition of QS receptors, less is known about the transport mechanisms that enable QSMs to cross the lipid membrane and access intracellular receptors. In this study, we used atomistic MD simulations and an implicit solvent model, called COSMOmic, to analyze the partitioning and translocation of QSMs across lipid bilayers. We performed umbrella sampling at atomistic resolution to calculate partitioning and translocation free energies for a set of naturally occurring QSMs, then used COSMOmic to screen the water-membrane partition and translocation free energies for 50 native and non-native QSMs that target LasR, one of the LuxR family of quorum-sensing receptors. This screening procedure revealed the influence of systematic changes to head and tail group structures on membrane partitioning and translocation free energies at a significantly reduced computational cost compared to atomistic MD simulations. Comparisons with previously determined QSM activities suggest that QSMs that are least likely to partition into the bilayer are also less active. This work thus demonstrates the ability of the computational protocol to interrogate QSM-bilayer interactions which may help guide the design of new QSMs with engineered membrane interactions.

The importance of intramolecular hydrogen bonds on the translocation of the small drug piracetam through a lipid bilayer

The number of hydrogen bond donors and acceptors is a fundamental molecular descriptor to predict the oral bioavailability of small drug candidates. In fact, the most widely used oral bioavailability rules (such as the Lipinsky's rule-of-five and the Veber rules) make use of this molecular descriptor. It is generally assumed that hydrogen bond donors and acceptors impact on passive diffusion across cell membranes, a fundamental event during drug absorption and distribution. Although the relationship between the number of these motifs and the probability of having good oral bioavailability has been studied and described for more than 20 years, little attention has been given to their spatial distribution in the molecule. In this paper, we used molecular dynamics to describe the effect of intramolecular hydrogen bonding on the passive diffusion of a small drug (piracetam) through a lipid membrane. The results indicated that the formation of an intramolecular hydrogen bond decreases the barrier for translocation by ca. 4 kcal mol-1 and increases the permeability of the tested molecule, partially compensating the desolvation penalty arising from the penetration of the drug into the biological membrane core. This effect was apparent in simulations where the formation of this interaction was prevented with the help of modified potentials, and in simulations with a similar compound to piracetam that was not able to form this intramolecular hydrogen bond due to a larger distance between the hydrogen bond donor and acceptor groups. These results were also supported by coarse-grained methods, which are becoming an important resource for sampling a larger chemical space of molecules, with reduced computational effort. Furthermore, entropy and enthalpy derived profiles were also obtained as the compounds translocated across the membrane, suggesting that, even though the process of formation of internal hydrogen bonds is entropically unfavorable, the enthalpic gain is such that the formation of these interactions is beneficial for the passive diffusion across cell membranes.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

Unassisted N-acetyl-phenylalanine-amide transport across membrane with varying lipid size and composition: kinetic measurements and atomistic molecular dynamics simulation

Biological membranes are essential to preserve structural integrity and regulate functional properties through the permeability of nutrients, pharmaceutical drugs, and neurotransmitters of a living cell. The movement of acetylated and amidated phenylalanine (NAFA) across 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane bilayers is investigated to probe physical transport. The rate of transport is measured experimentally applying parallel artificial membrane permeation assay (PAMPA). At the physiological temperature, 310 K, the measured time constants in the neutral pH were ∼6 h in DOPC and ∼3 h in POPC, while in a more acidic condition, at a pH 4.8, the time constants were ∼8 h in both lipids. Computationally, we have expanded our transport study of three aromatic dipeptides across a bilayer composed of DOPC¹⁸. In this study, we have examined the effects of lipid composition and bilayer size on the passive transport of NAFA by simulating the dipeptide in three different bilayers, with 50 DOPC lipids, 50 POPC lipids, and 40 POPC molecules. Specifically, atomistic molecular dynamics simulations with umbrella sampling were used to calculate the potential of mean force for the passive permeation of NAFA across the bilayers. Diffusion constants were then calculated by numerically solving the Smoluchowski equation. Permeability coefficients and mean first passage times were then calculated. Structural properties - Ramachandran plots, sidechain torsions, peptide insertion angles, radial distribution functions, and proximal peptide water molecules - were also examined to determine the effect of system size and lipid type. In terms of systems size, we observed a small decrease in the highest barrier of the potential of mean force and fewer sampled sidechain dihedral angle conformations with 40 versus 50 POPC lipids due to weaker membrane deformations within a smaller lipid bilayer. In terms of lipid type, DOPC contains two monounsaturated acyl chains compared to only one such acyl chain in POPC; therefore, DOPC bilayers are less ordered and more easily deformed, as seen by a much broader potential of mean force profile. The NAFA in DOPC lipid also transitioned to an internally hydrogen-bonded backbone conformation at lower membrane depths than in POPC. Similarly, as for other aromatic dipeptides, NAFA tends to insert into the membrane sidechain-first, remains mostly desolvated in the membrane center, and exhibits slow reorientations within the bilayer in both DOPC and POPC. With a joint experimental and computational study we have gained a new insight into the rate of transport and the underlying microscopic mechanism in different lipid bilayer conditions of the simplest hydrophobic aromatic dipeptide.Communicated by Ramaswamy H. Sarma.

The cut-off effect of n-alcohols in lipid rafts: A lipid-dependent phenomenon☆

n-Aliphatic alcohols act as anesthetics only up to a certain chain length, beyond which its biological activity disappears. This is known as the 'cut-off' phenomenon. Although the most accepted explanation is based on action sites in membrane proteins, it is not well understood why alcohols alter their functions. The structural dependence of these protein receptors to lipid domains known as 'lipid rafts', suggests a new approach to tackle the puzzling phenomenon. In this work, by performing molecular dynamic simulations (MDS) to explore the lipid role, we provide relevant molecular details about the membrane-alcohol interaction at the cut-off point regime. Since the high variability of the cut-off points found on protein receptors in neurons may be a consequence of differences in the lipid composition surrounding such proteins, our results could have a clear-cut importance.

NMR Structure and Dynamics Studies of Yeast Respiratory Supercomplex Factor 2

The Saccharomyces cerevisiae respiratory supercomplex factor 2 (Rcf2) is a 224-residue protein located in the mitochondrial inner membrane where it is involved in the formation of supercomplexes composed of cytochrome bc₁ and cytochrome c oxidase. We previously demonstrated that Rcf2 forms a dimer in dodecylphosphocholine micelles, and here we report the solution NMR structure of this Rcf2 dimer. Each Rcf2 monomer has two soluble α helices and five putative transmembrane (TM) α helices, including an unexpectedly charged TM helix at the C terminus, which mediates dimer formation. The NOE contacts indicate the presence of inter-monomer salt bridges and hydrogen bonds at the dimer interface, which stabilize the Rcf2 dimer structure. Moreover, NMR chemical shift change mapping upon lipid titrations as well as molecular dynamics analysis reveal possible structural changes upon embedding Rcf2 into a native lipid environment. Our results contribute to the understanding of respiratory supercomplex formation and regulation.

Effect of Protein Flexibility from Coarse-Grained Elastic Network Parameterizations on the Calculation of Free Energy Profiles of Ligand Binding

The characterization of the affinity and binding mechanism of specific molecules to a protein active site is scientifically and industrially relevant for many applications. In principle, this information can be obtained using molecular dynamics (MD) simulations by calculating the free energy profile of the process. However, this is a computationally demanding calculation. Currently, coarse-grained (CG) force fields are very well implemented for MD simulations of biomolecular systems. These computationally efficient force fields are a major advantage to the study of large model systems and/or those requiring long simulation times. The Martini model is currently one of the most popular CG force fields for these systems. For the specific case of protein simulations, to correctly maintain the macromolecular three-dimensional structure, the Martini model needs to include an elastic network (EN). In this work, the effect of protein flexibility, as induced by three EN models compatible with the Martini force field, was tested on the calculation of free energy profiles for protein-ligand binding. The EN models used were ElNeDyn, GoMartini, and GEN. The binding of triolein (TOG) and triacetin (TAG) to a lipase protein (thermomyces lanuginosa lipase-TLL) was used as a case study. The results show that inclusion of greater flexibility in the CG parameterization of proteins is of high importance in the calculation of the free energy profiles of protein-ligand systems. However, care must be taken in order to avoid unjustified large protein deformations. In addition, due to molecular flexibility there may be no absolute need for the center of the ligand to reach the center of the protein-binding site. The calculation of the energy profile to a distance of about 0.5 nm from the active site center can be sufficient to differentiate the affinity of different ligands to a protein.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

A Minimal Membrane Metal Transport System: Dynamics and Energetics of mer Proteins

The mer operon in bacteria encodes a set of proteins and enzymes that impart resistance to environmental mercury toxicity by importing Hg2+ and reducing it to volatile Hg(0). Because the reduction occurs in the cytoplasm, mercuric ions must first be transported across the cytoplasmic membrane by one of a few known transporters. MerF is the smallest of these, containing only two transmembrane helices and two pairs of vicinal cysteines that coordinate mercuric ions. In this work, we use molecular dynamics simulations to characterize the dynamics of MerF in its apo and Hg2+ -bound states. We find that the apo state positions one of the cysteine pairs closer to the periplasmic side of the membrane, while in the bound state the same pair approaches the cytoplasmic side. This finding is consistent with the functional requirement of accepting Hg2+ from the periplasmic space, sequestering it on acceptance, and transferring it to the cytoplasm. Conformational changes in the TM helices facilitate the functional interaction of the two cysteine pairs. Free-energy calculations provide a barrier of 16 kcal/mol for the association of the periplasmic Hg2+ -bound protein MerP with MerF and 7 kcal/mol for the subsequent association of MerF's two cysteine pairs. Despite the significant conformational changes required to move the binding site across the membrane, coarse-grained simulations of multiple copies of MerF support the expectation that it functions as a monomer. Our results demonstrate how conformational changes and binding thermodynamics could lead to such a small membrane protein acting as an ion transporter. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Effects of lipid composition on membrane distribution and permeability of natural quinones

Natural quinones are amphiphilic molecules that function as mobile charge carriers in biological energy transduction. Their distribution and permeation across membranes are important for binding to enzymatic complexes and for proton translocation. Here, we employ molecular dynamics simulations and free energy calculations with a carefully calibrated classical force-field to probe quinone distribution and permeation in a multi-component bilayer trying to mimic the composition of membranes involved in bioenergetic processes. Ubiquinone, ubiquinol, plastoquinone and menaquinone molecules with short and long isoprenoid tails are simulated. We find that penetration of water molecules bound to the polar quinone head increases considerably in the less ordered and porous bilayer formed by di-linoleoyl (18:2) phospholipids, resulting in a lower free energy barrier for quinone permeation and faster transversal diffusion. In equilibrium, quinone and quinol heads localize preferentially near lipid glycerol groups, but do not perform specific contacts with lipid polar heads. Quinone distribution is not altered significantly by the quinone head, tail and lipid composition in comparison to a single-component bilayer. This study highlights the role of lipid acyl chain unsaturation for permeation and transversal diffusion of polar molecules across biological membranes.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

The asymmetry of plasma membranes and their cholesterol content influence the uptake of cisplatin

The composition of the plasma membrane of malignant cells is thought to influence the cellular uptake of cisplatin and to take part in developing resistance to this widespread anti-cancer drug. In this work we study the permeation of cisplatin through the model membranes of normal and cancer cells using molecular dynamics simulations. A special attention is paid to lipid asymmetry and cholesterol content of the membranes. The loss of lipid asymmetry, which is common for cancer cells, leads to a decrease in their permeability to cisplatin by one order of magnitude in comparison to the membranes of normal cells. The change in the cholesterol molar ratio from 0% to 33% also decreases the permeability of the membrane by approximately one order of magnitude. The permeability of pure DOPC membrane is 5-6 orders of magnitude higher than one of the membranes with realistic lipid composition, which makes it as an inadequate model for the studies of drug permeability.

Iron and sulfur oxidation pathways of Acidithiobacillus ferrooxidans

Acidithiobacillus ferrooxidans is a gram-negative, autotrophic and rod-shaped bacterium. It can gain energy through the oxidation of Fe(II) and reduced inorganic sulfur compounds for bacterial growth when oxygen is sufficient. It can be used for bio-leaching and bio-oxidation and contributes to the geobiochemical circulation of metal elements and nutrients in acid mine drainage environments. The iron and sulfur oxidation pathways of A. ferrooxidans play key roles in bacterial growth and survival under extreme circumstances. Here, the electrons transported through the thermodynamically favourable pathway for the reduction to H₂O (downhill pathway) and against the redox potential gradient reduce to NAD(P)(H) (uphill pathway) during the oxidation of Fe(II) were reviewed, mainly including the electron transport carrier, relevant operon and regulation of its expression. Similar to the electron transfer pathway, the sulfur oxidation pathway of A. ferrooxidans, related genes and operons, sulfur oxidation mechanism and sulfur oxidase system are systematically discussed.

Different effects of cholesterol on membrane permeation of arginine and tryptophan revealed by bias-exchange metadynamics simulations

Experiments have shown that cholesterol influences the membrane permeability of small molecules, amino acids, and cell-penetrating peptides. However, their exact translocation mechanisms under the influence of cholesterol remain poorly understood. Given the practical importance of cell-penetrating peptides and the existence of varied cholesterol contents in different cell types, it is necessary to examine the permeation of amino acids in cholesterol-containing membranes at atomic level of details. Here, bias-exchange metadynamics simulations were employed to investigate the molecular mechanism of the membrane permeation of two amino acids Arg and Trp important for cell-penetrating peptides in the presence of different concentrations of cholesterol. We found that the free energy barrier of Arg+ (the protonated form) permeation increased linearly as the cholesterol concentration increased, whereas the barrier of Trp permeation had a rapid increase from 0 mol. % to 20 mol. % cholesterol-containing membranes and nearly unchanged from 20 mol. % to 40 mol. % cholesterol-containing membranes. Arg0 becomes slightly more stable than Arg+ at the center of the dipalmitoylphosphatidylcholine (DPPC) membrane with 40 mol. % cholesterol concentrations. As a result, Arg+ has a similar permeability as Trp at 0 mol. % and 20 mol. % cholesterol, but a significantly lower permeability than Trp at 40 mol. % cholesterol. This difference is caused by the gradual reduction of water defects for Arg+ as the cholesterol concentration increases but lack of water defects for Trp in cholesterol-containing membranes. Strong but different orientation dependence between Arg+ and Trp permeations is observed. These results provide an improved microscopic understanding of amino-acid permeation through cholesterol-containing DPPC membrane systems.

Coarse-grained model of titrating peptides interacting with lipid bilayers

Molecular-level computer simulations of peptide aggregation, translocation, and protonation at and in biomembranes are impeded by the large time and length scales involved. We present a computationally efficient, coarse-grained, and solvent-free model for the interaction between lipid bilayers and peptides. The model combines an accurate description of mechanical membrane properties with a new granular representation of the dielectric mismatch between lipids and the aqueous phase. All-atom force fields can be easily mapped onto the coarse-grained model, and parameters for coarse-grained monopeptides accurately extrapolate to membrane permeation free energies for the corresponding dipeptides and tripeptides. Acid-base equilibria of titratable amino acid residues are further studied using a constant-pH ensemble, capturing protonation state changes upon membrane translocation. Important differences between histidine, lysine, and arginine are observed, which are in good agreement with experimental observations.

The importance of the membrane interface as the reference state for membrane protein stability

The insertion of nascent polypeptide chains into lipid bilayer membranes and the stability of membrane proteins crucially depend on the equilibrium partitioning of polypeptides. For this, the transfer of full sequences of amino-acid residues into the bilayer, rather than individual amino acids, must be understood. Earlier studies have revealed that the most likely reference state for partitioning very hydrophobic sequences is the membrane interface. We have used μs-scale simulations to calculate the interface-to-transmembrane partitioning free energies ΔG_S→TM for two hydrophobic carrier sequences in order to estimate the insertion free energy for all 20 amino acid residues when bonded to the center of a partitioning hydrophobic peptide. Our results show that prior single-residue scales likely overestimate the partitioning free energies of polypeptides. The correlation of ΔG_S→TM with experimental full-peptide translocon insertion data is high, suggesting an important role for the membrane interface in translocon-based insertion. The choice of carrier sequence greatly modulates the contribution of each single-residue mutation to the overall partitioning free energy. Our results demonstrate the importance of quantifying the observed full-peptide partitioning equilibrium, which is between membrane interface and transmembrane inserted, rather than combining individual water-to-membrane amino acid transfer free energies.

Molecular transport through membranes: Accurate permeability coefficients from multidimensional potentials of mean force and local diffusion constants

Estimating the permeability coefficient of small molecules through lipid bilayer membranes plays an important role in the development of effective drug candidates. In silico simulations can produce acceptable relative permeability coefficients for a series of small molecules; however, the absolute permeability coefficients from simulations are usually off by orders of magnitude. In addition to differences between the lipid bilayers used in vitro and in silico, the poor convergence of permeation free energy profiles and over-simplified diffusion models have contributed to these discrepancies. In this paper, we present a multidimensional inhomogeneous solubility-diffusion model to study the permeability of a small molecule drug (trimethoprim) passing through a POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) lipid bilayer. Our approach improves the permeation model in three ways: First, the free energy profile (potential of mean force, PMF) is two-dimensional in two key coordinates rather than simply one-dimensional along the direction normal to the bilayer. Second, the 2-D PMF calculation has improved convergence due to application of the recently developed transition-tempered metadynamics with randomly initialized replicas, while third, the local diffusivity coefficient was calculated along the direction of the minimum free energy path on the two-dimensional PMF. The permeability is then calculated as a line integral along the minimum free energy path of the PMF. With this approach, we report a considerably more accurate permeability coefficient (only 2-5 times larger than the experimental value). We also compare our approach with the common practice of computing permeability coefficients based only on the translation of the center of mass of the drug molecule. Our paper concludes with a discussion of approaches for minimizing the computational cost for the purpose of more rapidly screening a large number of drug candidate molecules.

Simulations of lipid bilayers using the CHARMM36 force field with the TIP3P-FB and TIP4P-FB water models

The CHARMM36 force field for lipids is widely used in simulations of lipid bilayers. The CHARMM family of force fields were developed for use with the mTIP3P water model. This water model has an anomalously high dielectric constant and low viscosity, which limits its accuracy in the calculation of quantities like permeability coefficients. The TIP3P-FB and TIP4P-FB water models are more accurate in terms of the dielectric constant and transport properties, which could allow more accurate simulations of systems containing water and lipids. To test whether the CHARMM36 lipid force field is compatible with the TIP3P-FB and TIP4P-FB water models, we have performed simulations of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers. The calculated headgroup area, compressibility, order parameters, and X-ray form factors are in good agreement with the experimental values, indicating that these improved water models can be used with the CHARMM36 lipid force field without modification when calculating membrane physical properties. The water permeability predicted by these models is significantly different; the mTIP3P-model diffusion in solution and at the lipid-water interface is anomalously fast due to the spuriously low viscosity of mTIP3P-model water, but the potential of mean force of permeation is higher for the TIP3P-FB and TIP4P-FB models due to their high excess chemical potentials. As a result, the rates of water permeation calculated the FB water models are slower than the experimental value by a factor of 15-17, while simulations with the mTIP3P model only underestimate the water permeability by a factor of 3.

Reactive Oxygen and Nitrogen Species at Phospholipid Bilayers: Peroxynitrous Acid and Its Homolysis Products

Peroxynitrite is a powerful and long-lived oxidant generated in vivo. Peroxynitrous acid (ONOOH), its protonated form, may penetrate into phospholipid bilayers and undergo homolytic cleavage to nitrogen dioxide (·NO₂) and hydroxyl radicals (·OH), causing severe nitro-oxidative damage. The membrane environment is thought to influence ONOOH reactions, but the mechanisms remain speculative. Most experimental techniques lack the level of resolution required to keep track of the motion of very reactive species and their interactions with the membrane. Here, we performed molecular dynamics simulations of the permeation, interactions, and dynamics of ONOOH and its homolysis products in the phospholipid membrane environment. We started by developing an ONOOH model that successfully accounted for its conformational equilibria and solvation energies. Membrane permeation of ONOOH was accompanied by conformational changes. ONOOH exhibited a strong tendency to bind to and accumulate at the membrane headgroup region. There, ONOOH homolysis led to ·NO₂ radicals, which in turn partitioned to the membrane interior. About one-third of the ·OH radicals readily escaped to the aqueous phase within 1 ns. However, a significant number of ·OH radicals became trapped at the lipid headgroup region for a longer period. The possible implications for membrane-based nitration and oxidation processes were discussed.

In silico screening of drug-membrane thermodynamics reveals linear relations between bulk partitioning and the potential of mean force

The partitioning of small molecules in cell membranes-a key parameter for pharmaceutical applications-typically relies on experimentally available bulk partitioning coefficients. Computer simulations provide a structural resolution of the insertion thermodynamics via the potential of mean force but require significant sampling at the atomistic level. Here, we introduce high-throughput coarse-grained molecular dynamics simulations to screen thermodynamic properties. This application of physics-based models in a large-scale study of small molecules establishes linear relationships between partitioning coefficients and key features of the potential of mean force. This allows us to predict the structure of the insertion from bulk experimental measurements for more than 400 000 compounds. The potential of mean force hereby becomes an easily accessible quantity-already recognized for its high predictability of certain properties, e.g., passive permeation. Further, we demonstrate how coarse graining helps reduce the size of chemical space, enabling a hierarchical approach to screening small molecules.

Computed Free Energies of Peptide Insertion into Bilayers are Independent of Computational Method

We show that the free energy of inserting hydrophobic peptides into lipid bilayer membranes from surface-aligned to transmembrane inserted states can be reliably calculated using atomistic models. We use two entirely different computational methods: high temperature spontaneous peptide insertion calculations as well as umbrella sampling potential-of-mean-force (PMF) calculations, both yielding the same energetic profiles. The insertion free energies were calculated using two different protein and lipid force fields (OPLS protein/united-atom lipids and CHARMM36 protein/all-atom lipids) and found to be independent of the simulation parameters. In addition, the free energy of insertion is found to be independent of temperature for both force fields. However, we find major difference in the partitioning kinetics between OPLS and CHARMM36, likely due to the difference in roughness of the underlying free energy surfaces. Our results demonstrate not only a reliable method to calculate insertion free energies for peptides, but also represent a rare case where equilibrium simulations and PMF calculations can be directly compared.

Metronidazole within phosphatidylcholine lipid membranes: New insights to improve the design of imidazole derivatives

Metronidazole is a imidazole derivative with antibacterial and antiprotozoal activity. Despite its therapeutic efficacy, several studies have been developing new imidazole derivatives with lower toxicity. Considering that drug-membrane interactions are key factors for drugs pharmacokinetic and pharmacodynamic properties, the aim of this work is to provide new insights into the structure-toxicity relationship of metronidazole within phosphatidylcholine membranes. For that purpose, lipid membrane models (liposomes and monolayers) composed of dipalmitoylphosphatidylcholine were used. Experimental techniques (determination of partition coefficients and Langmuir isotherm measurements) were combined with molecular dynamics simulations. Different pHs and lipid phases were evaluated to enable a better extrapolation for in vivo conditions. The partition of metronidazole depends on the pH and on the biphasic system (octanol/water or DPPC/water system). At pH 1.2, metronidazole is hydrophilic. At pH 7.4, metronidazole disturbs the order and the packing of phospholipids. For this toxic effect, the hydroxyl group of the side chain of metronidazole is crucial by interacting with the water embedded in the membrane and with the phosphate group and the apolar chains of phospholipids.

Dimer Interface of the Human Serotonin Transporter and Effect of the Membrane Composition

The oligomeric state of membrane proteins has recently emerged in many cases as having an effect on their function. However, the intrinsic dynamics of their spatial organization in cells and model systems makes it challenging to characterize. Here we use molecular dynamics (MD) simulations at multiple resolutions to determine the dimer conformation of the human serotonin transporter (hSERT). From self-assembly simulations we predict dimer candidates and subsequently quantify their relative strength. We use umbrella sampling (US) replica exchange MD simulations for which we present extensive analysis of their efficiency and improved sampling compared to regular US MD simulations. The data shows that the most stable hSERT dimer interface is symmetrical and involves transmembrane helix 12 (TM12), similar to the crystal structure of the bacterial homologue LeuT, but with a slightly different orientation. We also describe the supramolecular organization of hSERT from a 250 μs self-assembly simulation. Finally, the effects of the presence of phosphatidylinositol bisphosphate or cholesterol in the membrane model has been quantified for the TM12-TM12 predicted interface. Collectively, the presented data bring new insight to the area of protein and lipid interplay in biological membranes.