
==== Front
J Phys Chem B
J Phys Chem B
jp
jpcbfk
The Journal of Physical Chemistry. B
1520-6106
1520-5207
American Chemical Society

39225398
10.1021/acs.jpcb.4c03411
Article
Employing Metadynamics to Predict the Membrane Partitioning of Carboxy-2H-Azirine Natural Products
https://orcid.org/0000-0003-2051-0561
Daly Clyde A. Jr. *
https://orcid.org/0000-0001-8890-8985
Seebald Leah M.
Wolk Emma
Department of Chemistry, Haverford College, 370 Lancaster Ave., Haverford, Pennsylvania 19041, United States
* Email: cdaly2@haverford.edu. Tel.: 610-896-1000.
03 09 2024
12 09 2024
128 36 87718781
24 05 2024
27 08 2024
21 08 2024
© 2024 The Authors. Published by American Chemical Society
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

Natural products containing the carboxy-2H-azirine moiety are an exciting target for investigation due to their broad-spectrum antimicrobial activity and new chemical space they afford for novel therapeutic development. The carboxy-2H-azirine moiety, including those appended to well-characterized chemical scaffolds, is understudied, which creates a challenge for understanding potential modes of inhibition. In particular, some known natural product carboxy-2H-azirines have long hydrophobic tails, which could implicate them in membrane-associated processes. In this study, we examined a small set of carboxy-2H-azirine natural products with varied structural features that could alter membrane partitioning. We compared the predicted membrane partitioning and alignment of these compounds to those of established membrane embedders with similar chemical scaffolds. To accomplish this, we developed parameters within the framework of the CHARMM36 force field for the 2H-azirine functional group and performed metadynamics simulations of the partitioning into a model bacterial membrane from aqueous solution. We determined that the carboxy-2H-azirine functional group is strongly hydrophilic, imbuing the long-chain natural products with amphipathicity similar to the known membrane-embedding molecules to which they were compared. For the long-chain analogs, the carboxy-2H-azirine head group stays within 1 nm of the phosphate layer, while the hydrophobic tail sits within the membrane. The carboxy-2H-azirine lacking the long alkyl chain instead partitions completely into aqueous solution.

Division of Chemistry 10.13039/100000165 CHE-1229354 Haverford College 10.13039/100019684 NA Division of Chemistry 10.13039/100000165 CHE-2018427 Division of Chemistry 10.13039/100000165 CHE-1662030 document-id-old-9jp4c03411
document-id-new-14jp4c03411
ccc-price
Special Issue

Published as part of The Journal of Physical Chemistry Bspecial issue “The Dynamic Structure of the Lipid Bilayer and Its Modulation by Small Molecules”.
==== Body
pmcIntroduction

The alarming global increase in antimicrobial resistance has created a pressing need to expand the chemical space of antimicrobials to mitigate infectious diseases.1 Natural products containing a carboxy-2H-azirine moiety have garnered attention for demonstrating moderate inhibition against resistant microbes,2−8 which provides a strong starting scaffold that can be optimized through structure–activity relationships (SAR).3,4,6 Despite these intriguing inhibition results, the biological target of inhibition is currently unknown. The carboxy-2H-azirine is a unique and understudied moiety, and its structure does not readily inspire confident assumptions about physiological behavior. To clarify the properties of this moiety and begin to build an understanding of its antimicrobial activity, we set out to predict the membrane partitioning behavior of molecules with a carboxy-2H-azirine functional group in a bacterial lipid bilayer. This includes investigating the orientation of these molecules when embedded in a lipid bilayer, as well as how these molecules protrude in relation to the polar phospholipid head groups of the lipid bilayer.

Structural features of natural products can often supply hints to infer biological targets. For example, many natural products known to inhibit membrane protein activity or disrupt membrane integrated processes have a distinctive amphipathic scaffold.9 This correlation between amphipathic structure and inhibition of a membrane target may be hypothesized to exist for the carboxy-2H-azirine natural products. Many of these carboxy-2H-azirine natural products comprise a long hydrophobic alkyl chain ranging from 16 to 18 carbons, suggestive of integration into a membrane lipid bilayer. On the other hand, the physiochemical properties of just the carboxy-2H-azirine moiety are not nearly as well-established as products. While the long alkyl chains of some of these compounds are suggestive of lipid embedding, the influence of the carboxy-2H-azirine on these membrane interactions is unclear because the hydrophilicity of the carboxy-2H-azirine has not been characterized.

To provide predictive clarity, we first examine structurally similar molecules to help guide predictions of the physicochemical properties of a carboxy-2H-azirine interacting with a membrane lipid bilayer. Computational tools such as molecular dynamics (MD) have been used to make predictions of membrane partitioning, permeability, and small-molecule embedding.11−15 While extensive force fields for biomolecules and many small molecules exist, there has been little effort to parametrize carboxy-2H-azirines. Most computational modeling of this moiety has been performed using electronic structure methods where force fields are not needed.16−21 However the processes at question in this work cannot be addressed by electronic structure methods. Some studies have used automatic parameter generation without adjustment to obtain force fields for small azirine-containing molecules.21−23 We hope to address the limitations of these studies by avoiding parameter assignment by pure analogy since there are few well-characterized and analogous moieties to carboxy-2H-azirines.25,26

Typical MD simulations can fully characterize the behavior of molecular processes with timescales between 1 ps and 1 μs.11,12,27 Because of the slow dynamics of membranes, straightforward atomistic simulations tend not to be able to reliably capture the equilibrium behavior of arbitrary small molecules relative to the bilayer.11,12,28 However, enhanced sampling methods can overcome these limitations.29,30 In these methods, the probability that configurations are sampled by the simulation is manipulated such that the full space of some collective variable (CV) is observed. Metadynamics can be used to sample a particular CV and reconstruct the potential of mean force (PMF).31,32 Here, biases are added to the potential energy surface to encourage sampling of the collective variable; the accumulated biases can be used at the end of the simulation to reconstruct the potential of mean force.31,32

We use metadynamics as a tool to quantify the propensity of these carboxy-2H-azirine natural products to integrate in representative bacterial lipid bilayers. To accomplish this, we will obtain the PMF for moving the molecules into and out of the membrane. The PMF will allow us to calculate the free energy for membrane partitioning, thus determining if the molecules are thermodynamically stable in the membrane or in the aqueous environment. We also investigate the orientation and solvent exposure of the molecules of interest by unbiasing the statistics from our biased simulations.

These simulations are performed for molecules 1-Az-H, 4-Z-Dy-Me, and 5-E-Dy-Me and two structurally similar molecules known to exhibit membrane embedding: stearic acid (2-St-H) and sphingosine (3-Sp) (Figure 1). Compounds 2-St-H, 3-Sp, 4-Z-Dy-Me, and 5-E-Dy-Me all share a linear carbon chain comprising 18 carbons. More specifically, 4-Z-Dy-Me, 5-E-Dy-Me, and 3-Sp are unsaturated between C4 and C5. Regarding this unsaturation, 4-Z-Dy-Me has the Z configuration, providing contrast to 5-E-Dy-ME and 3-Sp, which have the E configuration. Compound 3-Sp is saturated between C4 and C5, so this compound was used as a comparative model for membrane partitioning of different linear C18 alkyl chains. If the carboxy-2H-azirine group is overall hydrophilic, then we should expect to see similar membrane partitioning and alignment behavior for all four long-chain molecules in our investigated set. The greatest uncertainty in this prediction arises from the difference in oxidation at C1 of 4-Z-Dy-Me and 5-E-Dy-Me compared to 3-Sp, in combination with the understudied properties of the 2H-azirine heterocycle, especially related to hydrophilicity. These molecular comparisons will allow us to ensure that our metadynamics methods have properly ascertained the probability of embedded binding and allow us to determine if the physiochemical properties of these molecules are typical of membrane-embedding molecules.

Figure 1 Structures of the molecules investigated in this study. Carbons are numbered in blue according to IUPAC rules. Notable similarities include the alkyl chain length, the C4–C5 unsaturation, and the nitrogen and oxygen heteroatoms at C1 and C2.

From the set of compounds used for this work, compounds 1-Az-H and 2-St-H have carboxyl groups, whereas 4-Z-Dy-Me and 5-E-Dy-Me have a methyl ester. Carboxyl groups deprotonate in aqueous biological conditions (pH ∼7.4) to yield carboxylates, which are especially hydrophilic. However, inside a membrane, the less hydrophilic protonated states would likely predominate. Naïvely, this implies that both protonation states should be considered for these molecules. We chose to focus on the protonated states because they would less strongly bias the head groups to be hydrophilic, allowing the properties of the 2H-azirine moiety to more strongly determine the outcomes of our study.

Computational Methods

Force Fields

The CHARMM36 force field was used to model all molecular interactions.34−37 However, force fields for molecules containing the carboxy-2H-azirine moiety are unavailable in CHARMM. CGenFF can be used to automatically generate all of the parameters needed to run an MD simulation of a new small organic molecule.25,38 The parameters are determined by analogy to molecules already parametrized in the CHARMM force field and are assigned “penalties” rating the appropriateness of the analogy. A penalty larger than 50 indicates an extremely poor analogy and thus an unreliable parameter, which should be manually determined.25

In order to obtain initial atomic point charges for the central azirine ring of the carboxy-2H-azirine moiety, we first focused on the simple 2H-azirine molecule (i.e., C2H3N). The structure of this molecule is shown in Figure 2. First, we used CGenFF to obtain initial force field parameters. The charge for the C2 atom had a low penalty, but the N and C1 charges required manual parametrization. To parametrize charges which are compatible with the CHARMM force field, one must fit interactions with TIP3P water to quantum chemical potential energy surfaces for interactions with water.39 To perform these calculations, we used the program FFParam.39 This program was used to optimize the N and C1 atomic charges by fitting them to quantum mechanical water interaction energies (MP2/6-31G*) computed using Psi4 1.5.36,40−42 As a check, the fitted charges were compared to a CHELPG charge analysis of 2H-azirine in Q-Chem 5.4 at MP2/6-31G*, and we found that both sets of charges were similar. This includes the C1 atom, which was not changed from the CGenFF guess.43,44

Figure 2 Ball and stick structures of the two central molecules used to develop our azirine force field. Carbon atoms are gray, nitrogen atoms are blue, hydrogen atoms are white, and oxygen atoms are red.

Next, we parametrized 1-Az-H as our minimal example of a carboxy-2H-azirine-containing molecule. This structure is also shown in Figure 2. After using CGenFF to obtain the initial force field, we observed that many of the high penalty interactions included C1. Using constrained optimizations in Psi4 1.4 (for bonds and angles) and Q-Chem 5.4 (for dihedrals), we obtained potential energy surfaces for all high penalty (i.e., >50) interactions except for dihedrals involving hydrogen, which FFParam does not recommend parametrizing.39,42,44 These potential energy surfaces were obtained at the MP2/6-31G* level of electronic structure theory, and all degrees of freedom besides the one being scanned were optimized at each step. All degrees of freedom were scanned starting from the optimized geometry. Bonds were scanned in 17 steps from 15 pm less than the optimized length to 15 pm more than the optimized length, and the optimized length was the central ninth step. In a similar way, angles were scanned in 11 steps from 15° below to 15° above the optimum angle. For the dihedrals, full scans from −160° to + 160° were attempted. In several cases, only a portion of the potential energy surface could be collected due to high ring strain or interatomic clashes during optimization. However, the portion of the potential energy surface collected was always sufficient to fit the interaction of interest; the probability of the molecule exploring the unscanned coordinate values in an ordinary 300 K MD simulation was negligible. Dihedrals were fit to the appropriate CHARMM sum-of-cosines functions, carefully monitoring how few could be used to achieve a satisfactory fit of the potential energy surface. One out-of-plane improper dihedral angle was scanned from −5° to 5° and fit to a harmonic interaction.

The charges for 1-Az-H were obtained in the following way. First, CGenFF was used to obtain initial charges. A CHELPG analysis (MP2/6-31G*) was used to help constrain the charge optimization space. Comparing the CGenFF penalties, the CHELPG analyses, and the charges obtained for 2H-azirine, we believed that the charges for N and C1 would be similar between the two molecules. If the charges from 2H-azirine for the analogous atoms are used for N and C1 in 1-Az-H and the CGenFF values are used for all other charges, then the overall charge is −0.209. According to CHELPG, C2 in 1-Az-H would be expected to become more positive than C2 in 2H-azirine, and CGenFF agrees with a penalty less than 50. The CGenFF penalty for C4, however, was greater than 50, and the CHELPG analysis implied that it should be more positive than CGenFF suggested. Thus, +0.209 atomic units of charge were added to this atom. This leaves the overall molecule neutral. The final obtained point charges correlate well with synthetic observations of the electrophilicity and nucleophilicity of the atoms in 2H-azirine-bearing molecules.45,46 The N and C1 are predicted to be nucleophilic, and C2 is predicted to be electrophilic.

The high penalty interactions for 4-Z-Dy-Me and 5-E-Dy-Me (which are geometric isomers and so have the same force field parameters) were the same as for 1-Az-H, so the required force field parameters were simply borrowed. Outside of the carboxy-2H-azirine group, all CGenFF penalties were small for 4-Z-Dy-Me and 5-E-Dy-Me. The initial (CGenFF) and final (FFParam) values of all the force field parameters and the CHELPG analyses are given in the SI. Some examples of how the coordinates of interest were scanned and fit are provided on github at https://github.com/Daly-Lab-at-Haverford/code_examples. Similar strategies to those described in this section have been used to obtain CHARMM-compatible parameters for a wide variety of molecules.47−52

Molecular Dynamics Systems and Equilibration

Initial structures of our small molecules embedded in a membrane were constructed using CHARMM-GUI’s membrane builder.53−56 The membrane bilayer was constructed from 75% palmitoyloleoyl phosphatidylethanolamine (POPE) and 25% palmitoyloleoyl phosphatidylglycerol (POPG) and was oriented so that its surface plane was aligned with the xy plane. This mixture has previously been used as an effective model of the inner bacterial membrane.57 The x- and y-axis lengths of the membrane and the periodic box were initially 5 nm. Above and below the membrane, 4 nm of water was added; NaCl was added to neutralize the membrane and then to a total ionic strength of 0.15 M. The structure of each solute was initially obtained using Avogadro.58 A distinct system was created for each solute where the solute was placed in the center of the membrane. CHARMM-GUI was allowed to perform some minimization of the energy of the structure. An example initial membrane structure with a solute inserted is shown in Figure 3.

Figure 3 Initial structure of 5-E-Dy-Me in the membrane prior to equilibration. Water molecules (red and white) are shown in a licorice representation, and NaCl (purple and green, respectively) ions are shown as small spheres. Licorice POPE and POPG molecules (gray, white, red, and blue depending on the atom type) are slightly transparent to more easily show 2-E-Dy-Me, which is shown in a space filling representation. Phosphorus atoms (orange) are also shown in a space filling representation to clearly demark the boundary between the hydrophilic and hydrophobic portions of the simulation. On the bottom left, the directions of the x̂ and ẑ unit vectors are shown. The ŷ unit vector is perpendicular to both of these and so points toward the reader by the right-hand rule. The correspondence between colors and elements is: blue for nitrogen, red for oxygen, white for hydrogen, green for chlorine, purple for sodium, orange for phosphorus, and gray for carbon.

The membrane systems taken from CHARMM-GUI were subjected to the equilibration procedure suggested by CHARMM-GUI, which was implemented in OpenMM 7.6.59 The steps of this equilibration procedure have been outlined by Jo et al., though in the most recent version, the NPzγT (called NPAT in their work) simulation steps are each 500 ps.54,60 To produce NPzγT ensemble simulations, we held the surface tension (γ) at 0 kJ/mol/nm2, the pressure in the z-axis direction (Pz) at 1 bar, and the temperature (T) at 300 K.61 A Monte Carlo barostat with an update frequency of 150 steps was used to maintain the pressure and surface tension.62 Temperatures were controlled using Langevin dynamics with a friction coefficient of 1 ps–1.63 All simulations were performed in this ensemble. Nonbonded interactions were cut off after 1.2 nm with a smooth switching function used from 1.0 nm. Long-range electrostatics were corrected using particle mesh Ewald summation with an error tolerance of 0.0005. Configurations were sampled every 10 ps. After the equilibration recommended by CHARMM-GUI, velocities were randomized according to the Maxwell–Boltzmann distribution at 300 K, and the system was equilibrated with no constraints for 2 ns.

Well-Tempered Metadynamics

PMFs were obtained using an OpenMM implementation of well-tempered metadynamics.32 In this method, biases are added to the system along a CV of interest, encouraging the system to explore the full range of values of the CV. As specific values of the CV are sampled, the strength of the biases added at those CV values is reduced. At long times (i.e., t → ∞), the added biases will converge such that the bias along the CV stops changing. The total added bias is related to the PMF, f(z),1

where z is the CV, T is the temperature of the simulation (300 K in this case), the biases are added in such a way that the CV is sampled as if the temperature was T + ΔT, V(z, t) is the total added bias as a function of z and time, t, and C(t) is a constant, which can be ignored once convergence is reached.32 In our simulations, T + ΔT = 6000 K (i.e., the bias factor was 20). The CV was chosen to be the distance between the z-axis location of the center of mass of the molecule of interest and the z-axis location of the center of mass of all POPE and POPG molecules (i.e., the central plane of the membrane). The biases were shaped as Gaussian functions with a standard deviation of 0.3 nm. The initial bias height was 25 kJ/mol, and they were added every 100 fs. The CV was allowed to vary from −6 nm to +6 nm. All metadynamics simulations were run under NPzγT conditions (1 bar, 0 kJ/mol/nm2, and 300 K). Similar CVs and metadynamics settings have been used in previous studies of membrane embedding and partitioning.64−66 Because of the symmetry of the system, the converged PMF should be even such that f(z) = f(−z). Metadynamics simulations were run for 2 μs and then continued in 500 ns increments until three criteria were met: (1) the RMSD between f(z) and f(−z) was at most 4 kJ/mol for any specific value of z, (2) the RMS standard error over all values of z between f(z) and f(−z) was below 2.5 kJ/mol, and (3) the RMSD between the potential of mean force obtained at the end of the simulation and that obtained 50 ns prior to the end was less than 2.5 kJ/mol across all values of z. This took different amounts of simulation time for each molecule (Table 1). Due to these criteria, the ones place is treated as the final significant place value in the free energy values obtained in this work.

Table 1 Time to Convergence for Each Metadynamics Simulation

molecule	time (μs)	
1-Az-H	2.0	
2-St-H	2.0	
3-Sp	4.0	
4-Z-Dy-Me	3.0	
5-E-Dy-Me	2.0	

Membrane Partitioning

A major question in this study is the extent to which the molecules of interest embed in the membrane compared to dissolving into the water layer. This can be captured using the partition coefficient, defined as2

where CW and CB are the concentrations of solute molecules in water (W) or in the bilayer (B), respectively, , and ΔGW→B is the partitioning free energy describing moving a molecule from the water to the bilayer and is given by3

where VW and VB are the volumes of the water and bilayer regions, respectively, and the integrals are taken over these regions.67−69 In this work, we defined the bilayer region to range from −3 < z < 3 and the water region to include all other space. While the bilayer region as defined includes some volume that may intuitively be considered part of the water region, we include this volume for two reasons. First, this definition makes VW and VB equal, simplifying our calculations and interpretations. Second, the bilayer has some influence on the space outside the phosphate head groups, and we believe that this should be included in our investigation of its interactions on the molecules of interest.

Results

Potentials of Mean Force

The PMFs (Figure 4) clearly show that 1-Az-H is distinct from the other molecules. The small ΔGW→B value and the κ value near unity (Table 2) show that 1-Az-H is essentially indifferent to the choice between the bilayer and water regions. The long-chain molecules and 2-St-H, 3-Sp, 4-Z-Dy-Me, and 5-E-Dy-Me all strongly prefer to embed within the membrane over all other possibilities. There is a substantial difference between the binding strengths of 4-Z-Dy-Me and 5-E-Dy-Me; the E isomer is significantly more stable in the membrane. The partitioning values for 2-St-H are fairly similar to those for 4-Z-Dy-Me. Meanwhile, 3-Sp is similar to 5-E-Dy-Me.

Figure 4 Potentials of mean force for the z-axis distance between the membrane center of mass and the center of mass of the molecule of interest. Values are averages between f(z) and f(−z), and error bars are 95% confidence intervals. For context, a simulation snapshot is shown behind the plots. The zero of energy in this plot is the average value in the water region (|z| ≥ 3). The boundaries between the two regions are shown as black dotted lines. The correspondence between colors and elements is: blue for nitrogen, red for oxygen, white for hydrogen, green for chlorine, purple for sodium, orange for phosphorus, and gray for carbon.

Table 2 Free Energy Changes Obtained for the Partitioning of the Molecules to the Membrane in kJ/mol and Unitless Partition Coefficientsa

molecule	ΔGW→B	log κ	
1-Az-H	2	0	
2-St-H	–33	+6	
3-Sp	–46	+8	
4-Z-Dy-Me	–31	+5	
5-E-Dy-Me	–42	+7	
a Values are rounded to the nearest ones digit.

All four molecules have similar hydrophobic regions, but 3-Sp has a highly hydrophilic head group, while the head group of 2-St-H is more modestly hydrophilic. These trends can be explained if we suppose that the hydrophilicity of the carboxy-2H-azirine group is only a little weaker than that of the head group of 3-Sp. For 3-Sp and 5-E-Dy-Me, moving from water to the bilayer region allows the molecules to begin to hide their hydrophobic tails in the membrane. Moving too far into the membrane, however, is unfavorable because of the hydrophilicity of their head groups. For 2-St-H, the interfacial area at the start of the bilayer region does not provide sufficient opportunity to hide its hydrophobic tail, which is longer than for any of the other molecules. When embedded, it is able to move more deeply into the membrane because of the relatively weak hydrophilicity of its head group. 4-Z-Dy-Me has less exposed hydrophilic and hydrophobic surface area than 5-E-Dy-Me because of its shape. This weakens the penalty for remaining exposed to the liquid, making ΔGW→B more positive for 4-Z-Dy-Me.

These partitioning free energies can be compared to partitioning free energies obtained in other membrane partitioning investigations. We also include comparisons to binding free energies or differences in the PMF between a water region and a membrane region since their values are often fairly similar to partitioning free energies in this context. The PMFs of membrane components partitioning into a wide variety of membranes were computed by Rogers et al. (for phospholipids) and Ermilova and Lyubartsev (for cholesterol). Free energy changes upon removal from the bilayer ranging from −32.8 to −72.4 kJ/mol were obtained, which broadly agreed with experiment.64,70 Other studies focused on the binding of hydrophobic molecules. MacCallum and Tieleman found a free energy change of −21 kJ/mol for incorporating hexane into a bilayer, in agreement with experiment (estimated at −24.8 kJ/mol), and Bochicchio et al. obtained a free energy change near −30 kJ/mol for small hydrocarbons.66,69 Johansson and Lindahl (for amino acid side chains) and Jämbeck and Lyubartsev (for small drug molecules) examined the membrane binding of broadly amphiphilic molecules. The partitioning depends sensitively on the specific characteristics of the molecule of interest. Observed free energy changes are as low as −22.2 kJ/mol, but binding is unfavorable for many of the polar or charged amino acids (the binding free energy was not calculated).65,71 Most of the free energies obtained in these studies agree well with experiment.

The PMF for 1-Az-H obtained in this work is similar to those obtained for polar or charged amino acid side chains, having a large peak near z = 0. Interestingly, the peak in the PMF for 1-Az-H is about 3 times as large as the largest observed by Johansson and Lindahl for their side chain analogue of aspartic acid, about 12 kJ/mol; this confirms that 1-Az-H is strongly hydrophilic.71 Our ΔGW→B values for 2-St-H, 3-Sp, 4-Z-Dy-Me, and 5-E-Dy-Me fit into the high (unfavorable) end of the range for membrane components but are lower (more favorable) than values for hydrophobic molecules or small drug molecules. Our long-chain molecules of interest are similar to the investigated membrane components studied by Rogers et al. but have smaller hydrophilic and hydrophobic regions, leading to generally weaker binding.70 The hydrophilic portions of our molecules of interest lead to stronger binding than completely hydrophobic molecules, and the long 18 carbon chains of our molecules of interest lead to stronger binding than small amphiphilic molecules. Overall, the partitioning free energies obtained in this work fall into a reasonable range compared to similar free energy calculations in the literature.

Binding Probability Distributions and Orientation

The structural relationship between our molecules and the membrane is vitally important. It is not sufficient to know that a carboxy-2H-azirine-bearing molecule embeds within the membrane. We must also ask if the carboxy-2H-azirine group itself is available for further biological interactions. This question was addressed using the unbiased probability densities for all sampled configurations. After our simulations converged, we used the reweighting method of Branduardi et al. to weight the statistics of the simulation configurations observed under metadynamics.72 In this method, the unbiased probability, Pu(r), of a given simulation configuration r (representing all atomic positions) is given by4

where z(r) is the value of the CV in the configuration r of interest, Pb(r) is the biased probability with which the configuration was observed under metadynamics, V(z(r), t→∞) is the converged bias surface, and .72,73 The quantity eβV(z(r), t→∞) can be thought of as a weighting factor, which removes the effect of the bias potential, leaving only the “true” free energy surface to affect configurational probabilities. However, using metadynamics has allowed for the sampling of configurations across z even if their unbiased probability was marginal, allowing for more complete statistics.

For each molecule, orientational information was obtained by tracking the locations of their centers of mass, head groups, and tail groups. For 1-Az-H, 5-E-Dy-Me, and 4-Z-Dy-Me, the head group was represented by the nitrogen atom. For 3-Sp and 2-St-H, an oxygen atom was selected. For all our molecules of interest, the carbon of the methyl at the end of the carbon chain was selected to represent the tail end of the molecule. A depiction of these locations is shown in Figure 5. The locations of water oxygens, NaCl ions, hydrophobic membrane carbons, and phosphorus atoms were also tracked. The unbiased probability distributions of these atoms are shown as a function of the distance from the membrane central plane, z, in Figure 6A–F. In this figure, all distributions are normalized such that the integrated total probability is 1.0.

Figure 5 Vectors used to define the unbiased probability distributions shown in Figure 6. The pictured molecule of interest is 5-E-Dy-Me. Water and ions have been removed for clarity. Per the results shown below in Figure 6, the pictured configuration is highly unlikely in the real system. The correspondence between colors and elements is: blue for nitrogen, red for oxygen, white for hydrogen, orange for phosphorus, and gray for carbon.

Figure 6 (A–E) Unbiased probability density for the z-axis distance for several of the constituents of each simulation: the water oxygens, the hydrophobic carbons in the membrane, the NaCl ions, and the head groups, tail groups, and centers of mass of the molecules of interest. Part A includes color coded labels for the nonsolute portions of the simulations; the same colors are used in all parts of the figure. Each plot contains data from simulations of a different solute: (A) 1-Az-H, (B) 2-St-H, (C) 3-Sp, (D) 4-Z-Dy-Me, and (E) 5-E-Dy-Me. (F) Parallel elongation of the solutes compared to the membrane surface normal, PE.

The behavior of the membrane, water, and ions is essentially identical across simulations of each solute. Hydrophobic carbons are strongly populated near the center of the membrane. The phosphorus atoms peak at the interface between the hydrophobic and hydrophilic portions of the simulation. Because the phosphate groups are strongly negative and the POPG molecules have a net negative charge, there is a peak in the distribution of Na+ ions overlapping with the phosphorus atom distribution. For similar reasons, the Cl– ion distribution is peaked far from the membrane. These can be compared to the water molecule distribution, which is flat essentially everywhere outside the hydrophobic region of the membrane. If the membrane was not present, then we would expect the distribution of water molecules to be completely flat and identical to the distribution of NaCl ions. The water and NaCl ion distributions show decreased probability near |z| = 6 nm. This likely arises from the constant pressure control in the z-dimension; there were many configurations where the box was slightly smaller than 12 nm, so distances near ±6 nm were unphysically less likely to be observed. This has no effect on the major results of this work.

As expected from the PMFs shown in Figure 4, the center of mass of 1-Az-H is preferentially located outside the membrane (Figure 6A). There are two peaks in its probability distribution: The largest occurs just outside the phosphorus peak, indicating favorable interfacial interactions. The smaller peak occurs just as the Cl– ions start to become more probable than the baseline set by water molecules. Tracking the nitrogen and methyl carbon of 1-Az-H allows us to extract orientational information about the molecule. The distributions for these atoms are nearly identical to that for the center of mass. This can only be true if there is no orientational preference for 1-Az-H. Focusing on the largest peak near the phosphate groups, the peak for the nitrogen at |z| = 2.56 nm is very slightly closer to the phosphate groups than the peak for the methyl carbon at |z| = 2.62 nm. There are three atoms in 1-Az-H with substantial positive charges, which might be expected to interact favorably with phosphates: H5, C4, and C2 (Figure 2). Considering the geometry of 1-Az-H, it might be that the carboxy group is engaging in hydrogen bonds with phosphates. However, even this orientational preference is slight. Without a long hydrocarbon tail, 1-Az-H is hydrophilic.

The distributions for 4-Z-Dy-Me and 5-E-Dy-Me are similar (Figure 6D,E). For both, the most probable location for the center of mass is near |z| = 1 nm. The head group is located near the phosphate groups. In fact, the phosphorus and nitrogen distributions overlap substantially. However, the nitrogen atoms are unlikely to be located beyond the phosphorus atoms. Moving out of the membrane, the head group probability goes to zero before the phosphate probability does. The methyl tails are consistently located at the center of the membrane. All the probability distributions for 5-E-Dy-Me are somewhat sharper than the equivalents for 4-Z-Dy-Me, as should be expected from the relative ΔGW→B values for the two molecules. Together, we see that 4-Z-Dy-Me and 5-E-Dy-Me are broadly aligned with the membrane lipid molecules with their head groups located at the start of the hydrophilic region. This strongly implies that their head groups are hydrophilic. Together with the results from 1-Az-H, we can conclude that the carboxy-2H-azirine moiety is hydrophilic.

The two remaining molecules, 3-Sp and 2-St-H, have similar distributions but with important differences (Figure 6B,C). Both are aligned with the membrane and have their head groups pointed toward the hydrophilic region. However, the head group of 3-Sp overlaps more than any other molecule of interest with the membrane phosphorus atoms. For this molecule, we also see that the tail group probability distribution is relatively broad. The center of mass peak is relatively sharp for 3-Sp, which corresponds well with the especially negative ΔGW→B observed for that solute. This is all consistent with the idea that the head group of 3-Sp is highly hydrophilic and is more hydrophilic than any other membrane-embedding molecule of interest. The head group distribution for 2-St-H overlaps the least with the phosphorus atom distribution. In fact, there is a small peak in the head group distribution at the center of the membrane. The tail group distribution is more strongly peaked at the center of the membrane than any other molecule of interest; at the same time, there is a small probability that the tail group is located near the phosphate groups. The center of mass of 2-St-H is consistently located between the head and tail groups. This means that in most cases, the molecule is aligned with the membrane phospholipids; however, 2-St-H is also occasionally antialigned with the membrane. This is consistent with the idea that the head group of 2-St-H is more hydrophobic than the head group of any other molecule of interest. Together, 3-Sp and 2-St-H can be seen as behavioral end points that 4-Z-Dy-Me and 5-E-Dy-Me fall between.

The above results already strongly imply that the molecules 2-St-H, 3-Sp, 4-Z-Dy-Me, and 5-E-Dy-Me are most often embedded in the membrane and oriented with their head group pointed toward the solvent and with their methyl tails pointed into the center of the membrane, a strong orientational preference. To confirm this, we computed a direct measure of alignment with the membrane bilayer. We assumed that the membrane surface normal unit vector was always identical to the z-axis unit vector, ẑ . Then, we defined the parallel elongation of the solutes with the membrane, PE, as5

where is the vector pointing from the methyl tail to the head group atom in a particular molecule. The denominator of eq 5 is the maximum molecule length ever observed in our simulations. Our metadynamics simulations sample a wide variety of molecular configurations, so this is a slight overestimate for the length of each molecule at “full” elongation. The numerator is the length of the molecule in a specific configuration, but only in the z-axis direction. When PE = 1, the molecule is fully elongated and is perfectly parallel to the z-axis. In other words, PE = 1 describes perfect alignment with a single leaflet of the membrane. The PE = 0 end point can describe several possibilities. For instance, the molecule could be completely elongated but could have a 90° angle with respect to the membrane surface normal. Alternatively, PE will be close to zero if the molecule has a small radius of gyration. In any case, values of PE near zero describe a distinct lack of alignment with the membrane. The unbiased probability distribution of alignment values observed in our simulations is given in Figure 6F. The results clearly show that all molecules of interest are most often parallel to the normal vector and elongated except 1-Az-H. All PE values for 1-Az-H are equally probable, as would be expected from the results shown in Figure 6A. The probability of PE values especially close to 1 is lower than might be expected because the maximum molecular lengths used to define perfect elongation require longer than equilibrium bond lengths, which are generally energetically unfavorable.

Discussion

This work was performed with compounds 1-Az-H, 4-Z-Dy-Me, and 5-E-Dy-Me as representatives of the natural product carboxy-2H-azirine set. We selected 4-Z-Dy-Me and 5-E-Dy-Me because these compounds are geometric isomers, and such structural differences are known to impact membrane-dependent cell functions by altering membrane fluidity.74,75 Thus, membrane embedding may be expected to differ between these two isomers, possibly leading to distinct mechanisms of inhibition. Comparing the linear alkyl compounds composed of an 18-carbon chain, 2-St-H, 3-Sp, 4-Z-Dy-Me, and 5-E-Dy-Me, against a short-alkyl-chain compound composed of a 4-carbon chain, 1-Az-H, allowed us to investigate the effect of carbon chain length on membrane embedding. Long alkyl chains in small-molecule inhibitors can be distinct drivers of membrane integration,33,76 but the alkyl chain length and chain saturation can also influence binding affinity with a membrane-bound molecular target.77,78 We also sought to identify any perturbations that the carboxy-2H-azirine could have on the ability of long-chain molecules to embed in a membrane. Many of these properties depend sensitively on the hydrophilicity or hydrophobicity of the carboxy-2H-azirine group.

Our results strongly predict that, at equilibrium, molecules 2-St-H, 3-Sp, 4-Z-Dy-Me, and 5-E-Dy-Me will display embedded binding to a lipid membrane, while 1-Az-H will not. Even so, 1-Az-H can display interfacial binding (Figure 6A). For the other molecules in our set, moving from the water to the membrane is always favorable. For molecules 2-St-H, 3-Sp, 4-Z-Dy-Me, and 5-E-Dy-Me, the reverse process of moving from the membrane to water is highly unfavorable and essentially will never be observed in experiment. The large free energy changes observed in this study give us the confidence to predict that partitioning to the membrane will remain favorable for molecules 2-St-H, 3-Sp, 4-Z-Dy-Me, and 5-E-Dy-Me even if the composition of the membrane is changed substantially from our particular mixture of 75% POPE and 25% POPG.

We also predict that the carboxy-2H-azirine moiety is sufficiently hydrophilic to act in much the same way as the amino and carboxy groups present in 3-Sp and 2-St-H, respectively. The charge distribution in our force field, supported by quantum chemical CHELPG population analysis, is consistent with other hydrophilic compounds. For instance, the charges adopted for the N (−0.357) and the C2 (+0.427) in the 2H-azirine ring are similar in magnitude to those of the H atoms in TIP3P water (+0.417). We also find that 1-Az-H, the smallest carboxy-2H-azirine natural product, strongly prefers aqueous solution to the hydrophobic space within the membrane. The carboxy-2H-azirine groups in 5-E-Dy-Me and 4-Z-Dy-Me also prefer aqueous solution, causing these molecules to orient with those moieties facing the water layer and to overlap substantially with the phosphate head groups in the membrane. This behavior is similar to that observed for 3-Sp and 2-St-H. Like these nonazirine-containing compounds, the hydrophobic tails of 5-E-Dy-Me and 4-Z-Dy-Me point toward the center of the membrane, and the molecules align with the membrane.

Conclusions

In this work, we have made predictions of the membrane partitioning of carboxy-2H-azirine-containing molecules and determined that the moiety itself is likely hydrophilic, similar to an amino or hydroxyl group. Substantial membrane complexities encompassing phospholipids, cholesterol, proteins, and carbohydrates that give the membrane a fluid character could hypothetically be added to our model and may exist in the real systems.79 However, even in our simplified representation, we see the same behavior between molecules that are known to embed in complex membranes (3-Sp and 2-St-H) and long-chain carboxy-2H-azirines (5-E-Dy-Me and 4-Z-Dy-Me). This and the large free energy changes upon embedding that we observe allow us to speculate that similar partitioning would be observed if long-chain carboxy-2H-azirines were introduced to membranes containing substantial complexity, and we welcome future simulations and experiments testing this speculation.

Like other long-chain molecules with a hydrophilic head group, the mixture of polar and nonpolar interactions leads to long-chain molecules embedding in lipid bilayers with the hydrophilic carboxy-2H-azirine group colocated with the polar phospholipid head groups of the membrane bilayer. Gaining this understanding of the physiochemical properties of carboxy-2H-azirines utilizing metadynamics has been a crucial first step for guiding the elucidation of the biological mechanisms of inhibition from these molecules. Currently, biochemical experimentation using these molecules is hampered by the synthetic feasibility of generating enough material for comprehensive microbial assays with the complete panel of carboxy-2H-azirine natural products. This computational work allows us to make tentative predictions about potential biological targets for carboxy-2H-azirine inhibition based on membrane-embedding behaviors, hopefully guiding additional experiments on these molecules. Moreover, such an understanding opens up this unexplored chemical space, offering a distinct scaffold with promising potential for therapeutic development.

It is clear from our metadynamics simulations that 1-Az-H is an expected outlier and must be separately investigated. Lacking a long alkyl chain, 1-Az-H demonstrates complete aqueous partitioning due to the isolated hydrophilicity with the minimal carboxy-2H-azirine scaffold. This mode of binding suggests a distinct mechanism of action compared to those with longer alkyl chains; for instance, 1-Az-H may target a soluble or a membrane embedded process. In contrast, the membrane partitioning for 4-Z-Dy-Me and 5-E-Dy-Me suggests that the whole set of long-alkyl-chain carboxy-2H-azirine analogs may target biological processes at, or near, the lipid membrane.80 This work has directed our first steps toward a membrane-centric investigation for uncovering the biological targets of the carboxy-2H-azirine natural products.

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcb.4c03411.Additional details confirming wide configurational sampling in our metadynamics simulations, the initial CGenFF force field files for all of the small molecules, Q-Chem output files containing CHELPG analyses of 2H-azirine and 6-Az-H, and the final force field files obtained using FFParam (ZIP)

(Figures S1–S4) Distribution of X and Y axis values sampled for the center of mass of each solute, orientational distributions of solute molecules outside the bilayer, and asymmetry of unbiased orientations (PDF)

Supplementary Material

jp4c03411_si_001.zip

jp4c03411_si_002.pdf

The authors declare no competing financial interest.

Acknowledgments

Computational resources were provided in part by the MERCURY consortium (https://mercuryconsortium.org/) under NSF grants CHE-1229354, CHE-1662030, and CHE-2018427. C.A.D. and L.M.S. also thank the Haverford College for startup funding and computational resources.
==== Refs
References

Butler M. S. ; Blaskovich M. A. ; Cooper M. A. Antibiotics in the Clinical Pipeline in 2013. J. Antibiot. 2013, 66 (10 ), 571–591. 10.1038/ja.2013.86.
Miller T. W. ; Tristram E. W. ; Wolf F. J. Azirinomycin. II. Isolation and Chemical Characterization as 3-Methyl-2(2H) Azirinecarboxylic Acid. J Antibiot (Tokyo) 1971, 24 (1 ), 48–50. 10.7164/antibiotics.24.48.5541332
Skepper C. K. ; Dalisay D. S. ; Molinski T. F. Synthesis and Antifungal Activity of (−)-(Z)-Dysidazirine. Org. Lett. 2008, 10 (22 ), 5269–5271. 10.1021/ol802065d.18937483
Keffer J. L. ; Plaza A. ; Bewley C. A. Motualevic Acids A–F, Antimicrobial Acids from the Sponge Siliquariaspongia Sp. Org. Lett. 2009, 11 (5 ), 1087–1090. 10.1021/ol802890b.19191563
Gunasekera S. P. ; Kokkaliari S. ; Ratnayake R. ; Sauvage T. ; Dos Santos L. A. H. ; Luesch H. ; Paul V. J. Anti-Inflammatory Dysidazirine Carboxylic Acid from the Marine Cyanobacterium Caldora Sp. Collected from the Reefs of Fort Lauderdale, Florida. Molecules 2022, 27 (5 ), 1717 10.3390/molecules27051717.35268819
Stapley E. O. ; Hendlin D. ; Jackson M. ; Miller A. K. ; Hernandez S. ; Mata J. M. Azirinomycin. I. Microbial Production and Biological Characteristics. J Antibiot (Tokyo) 1971, 24 (1 ), 42–47. 10.7164/antibiotics.24.42.5541331
Skepper C. K. ; Molinski T. F. Long-Chain 2 H -Azirines with Heterogeneous Terminal Halogenation from the Marine Sponge Dysidea Fragilis. J. Org. Chem. 2008, 73 (7 ), 2592–2597. 10.1021/jo702435s.18321120
Molinski T. F. ; Ireland C. M. Dysidazirine, a Cytotoxic Azacyclopropene from the Marine Sponge Dysidea Fragilis. J. Org. Chem. 1988, 53 (9 ), 2103–2105. 10.1021/jo00244a049.
Yin H. ; Flynn A. D. Drugging Membrane Protein Interactions. Annu. Rev. Biomed. Eng. 2016, 18 (1 ), 51–76. 10.1146/annurev-bioeng-092115-025322.26863923
Marrink S. J. ; Corradi V. ; Souza P. C. T. ; Ingólfsson H. I. ; Tieleman D. P. ; Sansom M. S. P. Computational Modeling of Realistic Cell Membranes. Chem. Rev. 2019, 119 (9 ), 6184–6226. 10.1021/acs.chemrev.8b00460.30623647
Enkavi G. ; Javanainen M. ; Kulig W. ; Róg T. ; Vattulainen I. Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance. Chem. Rev. 2019, 119 (9 ), 5607–5774. 10.1021/acs.chemrev.8b00538.30859819
Shinoda W. Permeability across Lipid Membranes. Biochimica et Biophysica Acta (BBA) - Biomembranes 2016, 1858 (10 ), 2254–2265. 10.1016/j.bbamem.2016.03.032.27085977
Venable R. M. ; Krämer A. ; Pastor R. W. Molecular Dynamics Simulations of Membrane Permeability. Chem. Rev. 2019, 119 (9 ), 5954–5997. 10.1021/acs.chemrev.8b00486.30747524
Hannesschlaeger C. ; Horner A. ; Pohl P. Intrinsic Membrane Permeability to Small Molecules. Chem. Rev. 2019, 119 (9 ), 5922–5953. 10.1021/acs.chemrev.8b00560.30951292
Liu Y. ; Guan P. ; Wang Y. ; Liu L. ; Cao J. Mechanistic Insight into Decomposition of 2H-Azirines: Electronic Structure Calculations and Dynamics Simulations. J. Phys. Chem. A 2015, 119 (1 ), 67–78. 10.1021/jp511208p.25489983
Cao J. Photoinduced Reactions of Both 2-Formyl-2H-Azirine and Isoxazole: A Theoretical Study Based on Electronic Structure Calculations and Nonadiabatic Dynamics Simulations. The Journal of Chemical Physics 2015, 142 (24 ), 244302 10.1063/1.4922742.26133423
Dickerson C. E. ; Bera P. P. ; Lee T. J. Characterization of Azirine and Its Structural Isomers. J. Phys. Chem. A 2018, 122 (45 ), 8898–8904. 10.1021/acs.jpca.8b07788.30350988
Cho H.-G. Matrix Infrared Spectra and DFT Computations of 2H-Azirine Produced from Acetonitrile by Laser-Ablation Plume Radiation. Bulletin of the Korean Chemical Society 2014, 35 (7 ), 2093–2096. 10.5012/bkcs.2014.35.7.2093.
Zhu B. ; Zeng X. 3-Fluoro-2H-Azirine: Generation, Characterization, and Photochemistry. J. Phys. Chem. A 2023, 127 (50 ), 10591–10599. 10.1021/acs.jpca.3c06076.38063135
Wakchaure P. D. ; Ganguly B. Tuning the Electronic Effects in Designing Ligands for the Inhibition of Rotamase Activity of FK506 Binding Protein. Theor. Chem. Acc. 2021, 140 (1 ), 5 10.1007/s00214-020-02717-6.
Duarte V. C. M. ; Faustino H. ; Alves M. J. ; Gil Fortes A. ; Micaelo N. Asymmetric Diels–Alder Cycloadditions of d-Erythrose 1,3-Butadienes to Achiral t-Butyl 2H-Azirine 3-Carboxylate. Tetrahedron: Asymmetry 2013, 24 (18 ), 1063–1068. 10.1016/j.tetasy.2013.05.015.
Serratos I. N. ; Olayo R. ; Millán-Pacheco C. ; Morales-Corona J. ; Vicente-Escobar J. O. ; Soto-Estrada A. M. ; Córdoba-Herrera J. G. ; Uribe O. ; Gómez-Quintero T. ; Arroyo-Ornelas M. Á. ; et al. Modeling Integrin and Plasma-Polymerized Pyrrole Interactions: Chemical Diversity Relevance for Cell Regeneration. Sci Rep 2019, 9 (1 ), 7009 10.1038/s41598-019-43286-4.31065077
Vanommeslaeghe K. ; MacKerell A. D. Jr. Automation of the CHARMM General Force Field (CGenFF) I: Bond Perception and Atom Typing. J. Chem. Inf. Model. 2012, 52 (12 ), 3144–3154. 10.1021/ci300363c.23146088
Dodda L. S. ; Cabeza de Vaca I. ; Tirado-Rives J. ; Jorgensen W. L. LigParGen Web Server: An Automatic OPLS-AA Parameter Generator for Organic Ligands. Nucleic Acids Res. 2017, 45 (W1 ), W331–W336. 10.1093/nar/gkx312.28444340
Leach A. R. Molecular Modelling: Principles and Applications; 2nd ed.; Prentice Hall: Harlow, England; New York, 2001.
Zwier M. C. ; Chong L. T. Reaching Biological Timescales with All-Atom Molecular Dynamics Simulations. Current Opinion in Pharmacology 2010, 10 (6 ), 745–752. 10.1016/j.coph.2010.09.008.20934381
Spiwok V. ; Sucur Z. ; Hosek P. Enhanced Sampling Techniques in Biomolecular Simulations. Biotechnol. Adv. 2015, 33 (6, Part 2 ), 1130–1140. 10.1016/j.biotechadv.2014.11.011.25482668
Allison J. R. Computational Methods for Exploring Protein Conformations. Biochem. Soc. Trans. 2020, 48 (4 ), 1707–1724. 10.1042/BST20200193.32756904
Bussi G. ; Laio A. Using Metadynamics to Explore Complex Free-Energy Landscapes. Nat Rev Phys 2020, 2 (4 ), 200–212. 10.1038/s42254-020-0153-0.
Barducci A. ; Bussi G. ; Parrinello M. Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method. Phys. Rev. Lett. 2008, 100 (2 ), 020603 10.1103/PhysRevLett.100.020603.18232845
Skepper C. K. ; Dalisay D. S. ; Molinski T. F. Synthesis and Chain-Dependent Antifungal Activity of Long-Chain 2H-Azirine-Carboxylate Esters Related to Dysidazirine. Bioorg. Med. Chem. Lett. 2010, 20 (6 ), 2029–2032. 10.1016/j.bmcl.2010.01.068.20171099
Klauda J. B. ; Venable R. M. ; Freites J. A. ; O’Connor J. W. ; Tobias D. J. ; Mondragon-Ramirez C. ; Vorobyov I. ; MacKerell A. D. Jr ; Pastor R. W. Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. J. Phys. Chem. B 2010, 114 (23 ), 7830–7843. 10.1021/jp101759q.20496934
Venable R. M. ; Sodt A. J. ; Rogaski B. ; Rui H. ; Hatcher E. ; MacKerell A. D. ; Pastor R. W. ; Klauda J. B. CHARMM All-Atom Additive Force Field for Sphingomyelin: Elucidation of Hydrogen Bonding and of Positive Curvature. Biophys. J. 2014, 107 (1 ), 134–145. 10.1016/j.bpj.2014.05.034.24988348
Jorgensen W. L. ; Chandrasekhar J. ; Madura J. D. ; Impey R. W. ; Klein M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. The Journal of Chemical Physics 1983, 79 (2 ), 926–935. 10.1063/1.445869.
Beglov D. ; Roux B. Finite Representation of an Infinite Bulk System: Solvent Boundary Potential for Computer Simulations. The Journal of Chemical Physics 1994, 100 (12 ), 9050–9063. 10.1063/1.466711.
Vanommeslaeghe K. ; Hatcher E. ; Acharya C. ; Kundu S. ; Zhong S. ; Shim J. ; Darian E. ; Guvench O. ; Lopes P. ; Vorobyov I. ; et al. CHARMM General Force Field: A Force Field for Drug-like Molecules Compatible with the CHARMM All-Atom Additive Biological Force Fields. J. Comput. Chem. 2010, 31 (4 ), 671–690. 10.1002/jcc.21367.19575467
Kumar A. ; Yoluk O. ; MacKerell A. D. Jr. FFParam: Standalone Package for CHARMM Additive and Drude Polarizable Force Field Parametrization of Small Molecules. J. Comput. Chem. 2020, 41 (9 ), 958–970. 10.1002/jcc.26138.31886576
Møller C. ; Plesset M. S. Note on an Approximation Treatment for Many-Electron Systems. Phys. Rev. 1934, 46 (7 ), 618–622. 10.1103/PhysRev.46.618.
Binkley J. S. ; Pople J. A. ; Hehre W. J. Self-Consistent Molecular Orbital Methods. 21. Small Split-Valence Basis Sets for First-Row Elements. J. Am. Chem. Soc. 1980, 102 (3 ), 939–947. 10.1021/ja00523a008.
Smith D. G. A. ; Burns L. A. ; Simmonett A. C. ; Parrish R. M. ; Schieber M. C. ; Galvelis R. ; Kraus P. ; Kruse H. ; Di Remigio R. ; Alenaizan A. ; et al. PSI4 1.4: Open-Source Software for High-Throughput Quantum Chemistry. The Journal of Chemical Physics 2020, 152 (18 ), 184108 10.1063/5.0006002.32414239
Breneman C. M. ; Wiberg K. B. Determining Atom-Centered Monopoles from Molecular Electrostatic Potentials. The Need for High Sampling Density in Formamide Conformational Analysis. J. Comput. Chem. 1990, 11 (3 ), 361–373. 10.1002/jcc.540110311.
Epifanovsky E. ; Gilbert A. T. B. ; Feng X. ; Lee J. ; Mao Y. ; Mardirossian N. ; Pokhilko P. ; White A. F. ; Coons M. P. ; Dempwolff A. L. ; et al. Software for the Frontiers of Quantum Chemistry: An Overview of Developments in the Q-Chem 5 Package. The Journal of Chemical Physics 2021, 155 (8 ), 084801 10.1063/5.0055522.34470363
Xu F. ; Zeng F.-W. ; Luo W.-J. ; Zhang S.-Y. ; Huo J.-Q. ; Li Y.-P. 2H-Azirines: Recent Progress in Synthesis and Applications. Eur. J. Org. Chem. 2024, 27 , e202301292 10.1002/ejoc.202301292.
Do Carmo Carreiras M. ; Marco-Contelles J. Heterocyclic Targets in Advanced Organic Synthesis; Research Signpost: Kerala, India, 2011.
Pavlova A. ; Parks J. M. ; Gumbart J. C. Development of CHARMM-Compatible Force-Field Parameters for Cobalamin and Related Cofactors from Quantum Mechanical Calculations. J. Chem. Theory Comput. 2018, 14 (2 ), 784–798. 10.1021/acs.jctc.7b01236.29334459
Bouchouireb Z. ; Thany S. H. ; Le Questel J.-Y. Development of CHARMM Compatible Force Field Parameters and Molecular Dynamics Simulations for the Pesticide Flupyradifurone. J. Comput. Chem. 2024, 45 (7 ), 377–391. 10.1002/jcc.27245.37966816
Yu W. ; He X. ; Vanommeslaeghe K. ; MacKerell A. D. Jr. Extension of the CHARMM General Force Field to Sulfonyl-Containing Compounds and Its Utility in Biomolecular Simulations. J. Comput. Chem. 2012, 33 (31 ), 2451–2468. 10.1002/jcc.23067.22821581
Xu Y. ; Vanommeslaeghe K. ; Aleksandrov A. ; MacKerell A. D. Jr. ; Nilsson L. Additive CHARMM Force Field for Naturally Occurring Modified Ribonucleotides. J. Comput. Chem. 2016, 37 (10 ), 896–912. 10.1002/jcc.24307.26841080
Hatcher E. ; Guvench O. ; MacKerell A. D. Jr CHARMM Additive All-Atom Force Field for Aldopentofuranoses, Methyl-Aldopentofuranosides, and Fructofuranose. J. Phys. Chem. B 2009, 113 (37 ), 12466–12476. 10.1021/jp905496e.19694450
Mayne C. G. ; Saam J. ; Schulten K. ; Tajkhorshid E. ; Gumbart J. C. Rapid Parameterization of Small Molecules Using the Force Field Toolkit. J. Comput. Chem. 2013, 34 (32 ), 2757–2770. 10.1002/jcc.23422.24000174
Jo S. ; Lim J. B. ; Klauda J. B. ; Im W. CHARMM-GUI Membrane Builder for Mixed Bilayers and Its Application to Yeast Membranes. Biophys. J. 2009, 97 (1 ), 50–58. 10.1016/j.bpj.2009.04.013.19580743
Jo S. ; Kim T. ; Im W. Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations. PLOS ONE 2007, 2 (9 ), e880 10.1371/journal.pone.0000880.17849009
Lee J. ; Patel D. S. ; Ståhle J. ; Park S.-J. ; Kern N. R. ; Kim S. ; Lee J. ; Cheng X. ; Valvano M. A. ; Holst O. ; et al. CHARMM-GUI Membrane Builder for Complex Biological Membrane Simulations with Glycolipids and Lipoglycans. J. Chem. Theory Comput. 2019, 15 (1 ), 775–786. 10.1021/acs.jctc.8b01066.30525595
Wu E. L. ; Cheng X. ; Jo S. ; Rui H. ; Song K. C. ; Dávila-Contreras E. M. ; Qi Y. ; Lee J. ; Monje-Galvan V. ; Venable ; et al. CHARMM-GUI Membrane Builder toward Realistic Biological Membrane Simulations. J. Comput. Chem. 2014, 35 (27 ), 1997–2004. 10.1002/jcc.23702.25130509
Murzyn K. ; Róg T. ; Pasenkiewicz-Gierula M. Phosphatidylethanolamine-Phosphatidylglycerol Bilayer as a Model of the Inner Bacterial Membrane. Biophys. J. 2005, 88 (2 ), 1091–1103. 10.1529/biophysj.104.048835.15556990
Hanwell M. D. ; Curtis D. E. ; Lonie D. C. ; Vandermeersch T. ; Zurek E. ; Hutchison G. R. Avogadro: An Advanced Semantic Chemical Editor, Visualization, and Analysis Platform. J Cheminform 2012, 4 (1 ), 17 10.1186/1758-2946-4-17.22889332
Eastman P. ; Swails J. ; Chodera J. D. ; McGibbon R. T. ; Zhao Y. ; Beauchamp K. A. ; Wang L.-P. ; Simmonett A. C. ; Harrigan M. P. ; Stern C. D. ; et al. OpenMM 7: Rapid Development of High Performance Algorithms for Molecular Dynamics. PLOS Computational Biology 2017, 13 (7 ), e1005659 10.1371/journal.pcbi.1005659.28746339
Park S. ; Choi Y. K. ; Kim S. ; Lee J. ; Im W. CHARMM-GUI Membrane Builder for Lipid Nanoparticles with Ionizable Cationic Lipids and PEGylated Lipids. J. Chem. Inf. Model. 2021, 61 (10 ), 5192–5202. 10.1021/acs.jcim.1c00770.34546048
Smith D. J. ; Klauda J. B. ; Sodt A. J. Simulation Best Practices for Lipid Membranes [Article v1.0]. Living Journal of Computational Molecular Science 2019, 1 (1 ), 5966–5966. 10.33011/livecoms.1.1.5966.36204133
Chow K.-H. ; Ferguson D. M. Isothermal-Isobaric Molecular Dynamics Simulations with Monte Carlo Volume Sampling. Comput. Phys. Commun. 1995, 91 (1 ), 283–289. 10.1016/0010-4655(95)00059-O.
Izaguirre J. A. ; Sweet C. R. ; Pande V. S. Multiscale Dynamics of Macromolecules Using Normal Mode Langevin. In Biocomputing 2010; World Scientific: 2009; pp 240–251. 10.1142/9789814295291_0026.
Ermilova I. ; Lyubartsev A. P. Cholesterol in Phospholipid Bilayers: Positions and Orientations inside Membranes with Different Unsaturation Degrees. Soft Matter 2019, 15 (1 ), 78–93. 10.1039/C8SM01937A.
Jämbeck J. P. M. ; Lyubartsev A. P. Exploring the Free Energy Landscape of Solutes Embedded in Lipid Bilayers. J. Phys. Chem. Lett. 2013, 4 (11 ), 1781–1787. 10.1021/jz4007993.26283109
Bochicchio D. ; Panizon E. ; Ferrando R. ; Monticelli L. ; Rossi G. Calculating the Free Energy of Transfer of Small Solutes into a Model Lipid Membrane: Comparison between Metadynamics and Umbrella Sampling. The Journal of Chemical Physics 2015, 143 (14 ), 144108 10.1063/1.4932159.26472364
Boggara M. B. ; Krishnamoorti R. Partitioning of Nonsteroidal Antiinflammatory Drugs in Lipid Membranes: A Molecular Dynamics Simulation Study. Biophys. J. 2010, 98 (4 ), 586–595. 10.1016/j.bpj.2009.10.046.20159155
Ben-Tal N. ; Honig B. ; Bagdassarian C. K. ; Ben-Shaul A. Association Entropy in Adsorption Processes. Biophys. J. 2000, 79 (3 ), 1180–1187. 10.1016/S0006-3495(00)76372-7.10968982
MacCallum J. L. ; Tieleman D. P. Computer Simulation of the Distribution of Hexane in a Lipid Bilayer: Spatially Resolved Free Energy, Entropy, and Enthalpy Profiles. J. Am. Chem. Soc. 2006, 128 (1 ), 125–130. 10.1021/ja0535099.16390139
Rogers J. R. ; Espinoza Garcia G. ; Geissler P. L. Membrane Hydrophobicity Determines the Activation Free Energy of Passive Lipid Transport. Biophys. J. 2021, 120 (17 ), 3718–3731. 10.1016/j.bpj.2021.07.016.34302793
Johansson A. C. V. ; Lindahl E. Position-Resolved Free Energy of Solvation for Amino Acids in Lipid Membranes from Molecular Dynamics Simulations. Proteins: Structure, Function, and Bioinformatics 2008, 70 (4 ), 1332–1344. 10.1002/prot.21629.
Branduardi D. ; Bussi G. ; Parrinello M. Metadynamics with Adaptive Gaussians. J. Chem. Theory Comput. 2012, 8 (7 ), 2247–2254. 10.1021/ct3002464.26588957
Schäfer T. M. ; Settanni G. Data Reweighting in Metadynamics Simulations. J. Chem. Theory Comput. 2020, 16 (4 ), 2042–2052. 10.1021/acs.jctc.9b00867.32192340
Junker F. ; Ramos J. L. Involvement of the Cis/Trans Isomerase Cti in Solvent Resistance of Pseudomonas Putida DOT-T1E. J. Bacteriol. 1999, 181 (18 ), 5693–5700. 10.1128/JB.181.18.5693-5700.1999.10482510
Kabelitz N. ; Santos P. M. ; Heipieper H. J. Effect of Aliphatic Alcohols on Growth and Degree of Saturation of Membrane Lipids in Acinetobacter Calcoaceticus. FEMS Microbiol Lett 2003, 220 (2 ), 223–227. 10.1016/S0378-1097(03)00103-4.12670684
Alberts B. ; Johnson A. ; Lewis J. ; Raff M. ; Roberts K. ; Walter P. The Lipid Bilayer. In Molecular Biology of the Cell; 4 th edition; Garland Science: 2002.
Entova S. ; Guan Z. ; Imperiali B. Investigation of the Conserved Reentrant Membrane Helix in the Monotopic Phosphoglycosyl Transferase Superfamily Supports Key Molecular Interactions with Polyprenol Phosphate Substrates. Arch. Biochem. Biophys. 2019, 675 , 108111 10.1016/j.abb.2019.108111.31563509
Price N. P. J. ; Jackson M. A. ; Hartman T. M. ; Brändén G. ; Ek M. ; Koch A. A. ; Kennedy P. D. Branched Chain Lipid Metabolism As a Determinant of the N-Acyl Variation of Streptomyces Natural Products. ACS Chem Biol 2021, 16 (1 ), 116–124. 10.1021/acschembio.0c00799.33411499
Nicolson G. L. ; Ferreira de Mattos G. A Brief Introduction to Some Aspects of the Fluid–Mosaic Model of Cell Membrane Structure and Its Importance in Membrane Lipid Replacement. Membranes 2021, 11 (12 ), 947 10.3390/membranes11120947.34940448
Choo M. Z. Y. ; Chai C. L. L. The Polypharmacology of Natural Products in Drug Discovery and Development. In Natural Products; Annual Reports in Medicinal Chemistry; Elsevier: 2023; Vol. 61 , pp 55–100. 10.1016/bs.armc.2023.10.002.
