HCV Antiviral Drugs Have the Potential to Adversely Perturb the Fetal-Maternal Communication Axis through Inhibition of CYP3A7 DHEA-S Oxidation.

The hepatitis C virus (HCV) poses a great risk to pregnant people and their developing fetus, yet no HCV antiviral treatment guidelines have been established. While there has been a substantial increase in the development of HCV antivirals, the effect they have on the developing fetus remains poorly defined. Many of these drugs are metabolized through the cytochrome P450 CYP3A pathway, which is mediated by cytochrome P450 3A7 (CYP3A7) in the fetus and developing infant. In this study, we sought to investigate the effect HCV antivirals have on CYP3A7 metabolism, as this CYP enzyme plays a vital role in proper fetal and neonatal development. Of the 13 HCV antivirals we investigated, 8 (∼62%) inhibited CYP3A7 metabolic activity by 50% or more at a concentration of 20 µM. Furthermore, paritaprevir, asunaprevir, simeprevir, danoprevir, and glecaprevir all had observed half-maximal inhibitory concentrations between the range of 10 and 20 µM, which is physiologically relevant in comparison with the Km of dehydroepiandrosterone-sulfate (DHEA-S) oxidation (reported to be between 5 and 20 µM). We also discovered that paritaprevir is a time-dependent inhibitor of CYP3A7, which shifts the IC50 ∼twofold from 11 µM to 5 µM. Upon further characterization, paritaprevir inactivates DHEA-S metabolism by CYP3A7, with KI and Kinact values of 4.66 µM and 0.00954 minute-1, respectively. Depending on treatment plan and off-label drug use, HCV treatment could adversely affect the fetal-maternal communication axis by blocking fetal CYP3A7 metabolism of important endogenous hormones. SIGNIFICANCE STATEMENT: The prevalence of HCV in pregnant people is estimated at between 1% and 8% of the global population, yet little to no information exists about the risk antiviral treatment poses to the developing fetus. There is a potential risk of drugs adversely affecting mother-fetal communication by inhibiting fetal hepatic CYP3A7, an integral enzyme for estriol production. We discovered that five HCV antivirals inhibited DHEA-S metabolism by CYP3A7, and paritaprevir inactivated the enzyme. Our studies demonstrate the potential threat these drugs pose to proper fetal development.

Induced Fit Describes Ligand Binding to Membrane-Associated Cytochrome P450 3A4.

Cytochrome P450 3A4 (CYP3A4) is the dominant P450 involved in human xenobiotic metabolism. Competition for CYP3A4 therefore underlies several adverse drug-drug interactions. Despite its clinical significance, the mechanisms CYP3A4 uses to bind diverse ligands remain poorly understood. Highly monodisperse CYP3A4 embedded in anionic lipoprotein nanodiscs containing an equal mixture of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) were used to determine which of the limiting kinetic schemes that include protein conformational change, conformational selection (CS) or induced fit (IF), best described the binding of four known irreversible inhibitors. Azamulin, retapamulin, pleuromutilin, and mibrefadil binding to CYP3A4 nanodiscs conformed to a single-site binding model. Exponential fits of stopped-flow UV-visible absorption spectroscopy data supported multiple-step binding mechanisms. Trends in the rates of relaxation to equilibrium with increasing ligand concentrations were ambiguous as to whether IF or CS was involved; however, global fitting and consideration of the rate constants favored an IF mechanism. In the case of mibrefadil, a transient complex was observed in the stopped-flow UV-visible experiment, definitively assigning the presence of IF in ligand binding. While these studies only consider a small region of CYP3A4's vast ligand space, they provide kinetic evidence that CYP3A4 can use an IF mechanism. SIGNIFICANCE STATEMENT: CYP3A4 is capable of oxidizing numerous xenobiotics, including many drugs. Such promiscuity could not be achieved without conformational changes to accommodate diverse substrates. It is unknown whether conformational heterogeneity is present before (conformational selection) or after (induced fit) ligand binding. Stopped-flow measurements of suicide inhibitors binding to nanodisc-embedded CYP3A4 combined with sophisticated numerical analyses support that induced fit better describes ligand binding to this important enzyme.

PSI relieves the pressure of membrane fusion.

Some plant proteases contain a latent sequence known as the plant-specific insert (PSI) that, upon release from the full protease sequence, initiates membrane fusion to defend from pathogens. However, the mechanism by which it exerts its effects has been unclear. Zhao et al. report an elegant integration of biophysical experiments and molecular dynamics simulations to reveal events leading up to PSI-mediated membrane fusion. Their results demonstrate a pH-dependent monomer-to-dimer transition, clear evidence of membrane association, and probable structures of prefusion intermediates. These data expand our understanding of the elusive PSIs and may provide new directions for antimicrobial development.

Dynamics and Location of the Allosteric Midazolam Site in Cytochrome P4503A4 in Lipid Nanodiscs.

Promiscuous and allosteric drug interactions with cytochrome P450 3A4 (CYP3A4) are ubiquitous but incompletely understood at the molecular level. A classic allosteric CYP3A4 drug interaction includes the benzodiazepine midazolam (MDZ). MDZ exhibits homotropic and heterotropic allostery when metabolized to 1'-hydroxy and 4-hydroxy metabolites in varying ratios. The combination of hydrogen-deuterium exchange mass spectrometry (HDX-MS) and Gaussian accelerated molecular dynamics (GaMD) simulations of CYP3A4 in lipid nanodiscs and in a lipid bilayer, respectively, reveals MDZ-dependent changes in dynamics in a membrane environment. The F-, G-, and intervening helices, as well as the loop preceding the ß1-sheets, display the largest observed changes in HDX. The GaMD suggests a potential allosteric binding site for MDZ in the F'- and G'-regions, which undergo significant increases in HDX at near-saturating MDZ concentrations. The HDX-MS and GaMD results confirm that changes in dynamics are most significant near the developing consensus allosteric site, and these changes are distinct from those observed previously with the nonallosteric inhibitor ketoconazole. The results suggest that the allosteric MDZ remains mobile in its binding site at the Phe-cluster. The results further suggest that this binding site remains dynamic or changes the depth of insertion in the membrane.

Dynamics and Mechanism of Binding of Androstenedione to Membrane-Associated Aromatase.

Aromatase (CYP19A1) catalyzes the synthesis of estrogens from androgens and is an invaluable target of pharmacotherapy for estrogen-dependent cancers. CYP19A1 is also one of the most primordial human CYPs and, to the extent that its fundamental dynamics are conserved, is highly relevant to understanding those of the more recently evolved and promiscuous enzymes. A complementary approach employing molecular dynamics simulations and hydrogen-deuterium exchange mass spectrometry (HDX-MS) was employed to interrogate the changes in CYP19A1 dynamics coupled to binding androstenedione (ASD). Gaussian-accelerated molecular dynamics and HDX-MS agree that ASD globally suppresses CYP19A1 dynamics. Bimodal HDX patterns of the B'-C loop potentially arising from at least two conformations are present in free 19A1 only, supporting the possibility that conformational selection is operative. Random-acceleration molecular dynamics and adaptive biasing force simulations illuminate ASD's binding pathway, predicting ASD capture in the lipid headgroups and a pathway to the active site shielded from solvent. Intriguingly, the predicted access channel in 19A1 aligns well with the steroid binding sites of other human sterol-oxidizing CYPs.

Membrane-embedded substrate recognition by cytochrome P450 3A4.

Cytochrome P450 3A4 (CYP3A4) is the dominant xenobiotic-metabolizing enzyme in the liver and intestine and is involved in the disposition of more than 50% of drugs. Because of its ability to bind multiple substrates, its reaction kinetics are complex, and its association with the microsomal membrane confounds our understanding of how this enzyme recognizes and recruits diverse substrates. Testosterone (TST) hydroxylation is the prototypical CYP3A4 reaction, displaying positive homotropic cooperativity with three binding sites. Here, exploiting the capability of accelerated molecular dynamics (aMD) to sample events in the millisecond regime, I performed >25-µs aMD simulations in the presence of three TST molecules. These simulations identified high-occupancy surface-binding sites as well as a pathway for TST ingress into the CYP3A4 active site originating in the membrane. Adaptive biasing force analysis of the latter pathway revealed a metastable intermediate that could constitute a third binding site at high TST concentrations. Prompted by the observation that interactions between TST and the G'-helix mobilize the ligand into the active site, a free-energy analysis of TST distribution in the membrane was conducted and revealed that the depth of the G'-helix is optimal for extracting TST. In summary, these simulations confirm separate, but adjacent substrate-binding sites within the enzyme and the existence of an auxiliary TST-binding site. The broader impact of these simulations is that they support a mechanism in which cytochromes P450 directly recruit membrane-solubilized substrates.

Noncovalent interactions dominate dynamic heme distortion in cytochrome P450 4B1.

Cytochrome P450 4B1 (4B1) functions in both xenobiotic and endobiotic metabolism. An ester linkage between Glu-310 in 4B1 and the 5-methyl group of heme facilitates preferential hydroxylation of terminal (ω) methyl groups of hydrocarbons (HCs) and fatty acids compared with ω-1 sites bearing weaker C-H bonds. This preference is retained albeit diminished 4-fold for the E310A mutant, but the reason for this is unclear. Here, a crystal structure of the E310A-octane complex disclosed that noncovalent interactions maintain heme deformation in the absence of the ester linkage. Consistent with the lower symmetry of the heme, resonance Raman (RR) spectroscopy revealed large enhancements of RR peaks for high-spin HC complexes of 4B1 and the E310A mutant relative to P450 3A4. Whereas these enhancements were diminished in RR spectra of a low-spin 4B1-N-hydroxy-N'-(4-butyl-2-methylphenyl)formamidine complex, a crystal structure indicated that this inhibitor does not alter heme ruffling. RR spectra of Fe2+-CO HC complexes revealed larger effects of HC length in E310A than in 4B1, suggesting that reduced rigidity probably underlies increased E310A-catalyzed (ω-1)-hydroxylation. Diminished effects of the HC on the position of the Fe-CO stretching mode in 4B1 suggested that the ester linkage limits substrate access to the CO. Heme ruffling probably facilitates autocatalytic ester formation by reducing inhibitory coordination of Glu-310 with the heme iron. This also positions the 5-methyl for a reaction with the proposed glutamyl radical intermediate and potentially enhances oxo-ferryl intermediate reactivity for generation of the glutamyl radical to initiate ester bond formation and ω-hydroxylation.

N-Heterocyclic Carbene Capture by Cytochrome P450 3A4.

Cytochrome P450 3A4 (CYP3A4) is the dominant P450 enzyme involved in human drug metabolism, and its inhibition may result in adverse interactions or, conversely, favorably reduce the systemic elimination rates of poorly bioavailable drugs. Herein we describe a spectroscopic investigation of the interaction of CYP3A4 with N-methylritonavir, an analog of ritonavir, widely used as a pharmacoenhancer. In contrast to ritonavir, the binding affinity of N-methylritonavir for CYP3A4 is pH-dependent. At pH <7.4, the spectra are definitively type I, whereas at pH ≥7.4 the spectra have split Soret bands, including a red-shifted component characteristic of a P450-carbene complex. Variable-pH UV-visible spectroscopy binding studies with molecular fragments narrows the source of this pH dependence to its N-methylthiazolium fragment. The C2 proton of this group is acidic, and variable-pH resonance Raman spectroscopy tentatively assigns it a pKa of 7.4. Hence, this fragment of N-methylritonavir is expected to be readily deprotonated under physiologic conditions to yield a thiazol-2-ylidene, which is an N-heterocyclic carbene that has high-affinity for and is presumed to be subsequently captured by the heme iron. This mechanism is supported by time-dependent density functional theory with an active site model that accurately reproduces distinguishing features of the experimental UV-visible spectra of N-methylritonavir bound to CYP3A4. Finally, density functional theory calculations support that this novel interaction is as strong as the tightest-binding azaheterocycles found in P450 inhibitors and could offer new avenues for inhibitor development.

Evidence for an elevated aspartate pK(a) in the active site of human aromatase.

Aromatase (CYP19A1), the enzyme that converts androgens to estrogens, is of significant mechanistic and therapeutic interest. Crystal structures and computational studies of this enzyme shed light on the critical role of Asp(309) in substrate binding and catalysis. These studies predicted an elevated pK(a) for Asp(309) and proposed that protonation of this residue was required for function. In this study, UV-visible absorption, circular dichroism, resonance Raman spectroscopy, and enzyme kinetics were used to study the impact of pH on aromatase structure and androstenedione binding. Spectroscopic studies demonstrate that androstenedione binding is pH-dependent, whereas, in contrast, the D309N mutant retains its ability to bind to androstenedione across the entire pH range studied. Neither pH nor mutation perturbed the secondary structure or heme environment. The origin of the observed pH dependence was further narrowed to the protonation equilibria of Asp(309) with a parallel set of spectroscopic studies using exemestane and anastrozole. Because exemestane interacts with Asp(309) based on its co-crystal structure with the enzyme, its binding is pH-dependent. Aromatase binding to anastrozole is pH-independent, consistent with the hypothesis that this ligand exploits a distinct set of interactions in the active site. In summary, we assign the apparent pK(a) of 8.2 observed for androstenedione binding to the side chain of Asp(309). To our knowledge, this work represents the first experimental assignment of a pK(a) value to a residue in a cytochrome P450. This value is in agreement with theoretical calculations (7.7-8.1) despite the reliance of the computational methods on the conformational snapshots provided by crystal structures.

Photoswitchable anticancer activity via trans-cis isomerization of a combretastatin A-4 analog.

Combretastatin A-4 (CA4) is highly potent anticancer drug that acts as an inhibitor of tubulin polymerization. The core of the CA4 structure contains a cis-stilbene, and it is known that the trans isomer is significantly less potent. We prepared an azobenzene analog of CA4 (Azo-CA4) that shows 13-35 fold enhancement in potency upon illumination. EC50 values in the light were in the mid nM range. Due to its ability to thermally revert to less toxic trans form, Azo-CA4 also has the ability to automatically turn its activity off with time. Azo-CA4 is less potent than CA-4 because it degrades in the presence of glutathione as evidenced by UV-Vis spectroscopy and ESI-MS. Nevertheless, Azo-CA4 represents a promising strategy for switchable potency for treatment of cancer.

Ab initio dynamics of the cytochrome P450 hydroxylation reaction.

The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.

Small angle scattering reveals the orientation of cytochrome P450 19A1 in lipoprotein nanodiscs.

Human aromatase (CYP19A1), the cytochrome P450 enzyme responsible for conversion of androgens to estrogens, was incorporated into lipoprotein nanodiscs (NDs) and interrogated by small angle X-ray and neutron scattering (SAXS/SANS). CYP19A1 was associated with the surface and centered at the edge of the long axis of the ND membrane. In the absence of the N-terminal anchor, the amphipathic A'- and G'-helices were predominately buried in the lipid head groups, with the possibly that their hydrophobic side chains protrude into the hydrophobic, aliphatic tails. The prediction is like that for CYP3A4 based on SAXS employing a similar modeling approach. The orientation of CYP19A1 in a ND is consistent with our previous predictions based on molecular dynamics simulations and lends additional credibility to the notion that CYP19A1 captures substrates from the membrane.

Cytochrome P450 compound I in the plane wave pseudopotential framework: GGA electronic and geometric structure of thiolate-ligated iron(IV)-oxo porphyrin.

The cytochromes P450 constitute a ubiquitous family of metalloenzymes, catalyzing manifold reactions of biological and synthetic importance via a thiolate-ligated iron-oxo (IV) porphyrin radical species denoted compound I (Cpd I). Experimental investigations have implicated this intermediate in a broad spectrum of biophysically interesting phenomena, further augmenting the importance of a Cpd I model system. Ab initio molecular dynamics, including Car-Parrinello and path integral methods, conjoin electronic structure theory with finite temperature simulation, affording tools most valuable to approach such enzymes. These methods are typically driven by density functional theory (DFT) in a plane-wave pseudopotential framework; however, existing studies of Cpd I have been restricted to localized Gaussian basis sets. The appropriate choice of density functional and pseudopotential for such simulations is accordingly not obvious. To remedy this situation, a systematic benchmarking of thiolate-ligated Cpd I is performed using several generalized-gradient approximation (GGA) functionals in the Martins-Troullier and Vanderbilt ultrasoft pseudopotential schemes. The resultant electronic and structural parameters are compared to localized-basis DFT calculations using GGA and hybrid density functionals. The merits and demerits of each scheme are presented in the context of reproducing existing experimental and theoretical results for Cpd I.

Solvatochromism and the solvation structure of benzophenone.

Many complex molecular phenomena, including macromolecular association, protein folding, and chemical reactivity, are determined by the nuances of their electrostatic landscapes. The measurement of such electrostatic effects is nonetheless difficult, and is typically accomplished by exploiting a spectroscopic probe within the system of interest, such as through the vibrational Stark effect. Raman spectroscopy and solvatochromism afford an alternative to this method, circumventing the limitations of infrared spectroscopy, providing a lower detection limit, and permitting measurement in a native chemical environment. To explore this possibility, the solvatochromism of the C=O and aromatic C-H stretching modes of benzophenone are investigated using Raman spectroscopy. In conjunction with density functional theory calculations, these observations are sufficient to determine the probe electrostatic environment as well as contributions from halogen and hydrogen bonding. Further analysis using a detailed Kubo-Anderson lineshape model permits the detailed assignment of distinct hydrogen bonding configurations for water in the benzophenone solvation shell. These observations reinforce the use of benzophenone as an effective electrostatic probe for complex chemical systems.

Nanodisc-embedded cytochrome P450 P3A4 binds diverse ligands by distributing conformational dynamics to its flexible elements.

Cytochrome P450 3A4 (CYP3A4) metabolizes a wide range of drugs and toxins. Interactions of CYP3A4 with ligands are difficult to predict due to promiscuity and conformational flexibility. To better understand CYP3A4 conformational responses to ligands we use hydrogen deuterium exchange mass spectrometry (HDX-MS) to investigate the effect of ligands on nanodisc-embedded CYP3A4. For a subset of CYP3A4-ligand complexes, differences in the low-frequency modes derived by principal component analyses of molecular dynamics trajectories mirrored the HDX-MS results. The effects of ligands are distributed to flexible elements of CYP3A4 between stretches of secondary structure. The largest effects occur in the F- and G-helices, where most ligands increase the flexibility of the F-helix and connecting loops and decrease the flexibility of the C-term of the G-helix. Most ligands affect the E-F-G, CD and HI regions of the protein. Ligand-dependent differences are observed in the A"-A' loop, BC region, E-helix, K-ß1 region, proximal loop, and C-term loop. Correlated HDX responses were observed in the CD region and the C-term of the G-helix that were most pronounced for Type II ligands. Collectively, the HDX and molecular dynamics results suggest that CYP3A4 accommodates diverse binding partners by propagating local backbone fluctuations from the binding site onto the flexible regions of the enzyme via long-range interactions that are differentially modulated by ligands. In contrast to the paradigm wherein ligands decrease protein dynamics at their binding site, a wide range of ligands modestly increase CYP3A4 dynamics throughout the protein including effects remote from the active site.

Coupled electron transfer and proton hopping in the final step of CYP19-catalyzed androgen aromatization.

Aromatase (CYP19) catalyzes the terminal step in estrogen biosynthesis, which requires three separate oxidation reactions, culminating in an enigmatic aromatization that converts an androgen to an estrogen. A stable ferric peroxo (Fe(3+)O(2)(2-)) intermediate is seen by electron paramagnetic resonance, but its role in this complex reaction remains controversial. Combining molecular dynamics simulation and hybrid quantum mechanics/molecular mechanics, we show that ferric peroxo addition to the 19-aldehyde initiates the reaction. Stepwise cleavage of the C10-C19 and O-O bonds of the peroxohemiacetal extrudes formate and yields Compound II, which in turn desaturates the steroid through successive abstraction of the 1ß-hydrogen atom and deprotonation of the 2ß-position. Throughout the transformation, a proton is cyclically relayed between D309 and the substrate to stabilize reaction intermediates. This mechanism invokes novel oxygen intermediates and provides a unifying interpretation of past experimental mechanistic studies.

Tetrahydrofolate recognition by the mitochondrial folate transporter.

A mitochondrial carrier family (MCF) of transport proteins facilitates the transfer of charged small molecules across the inner mitochondrial membrane. The human genome has ∼50 genes corresponding to members of this family. All MCF proteins contain three repeats of a characteristic and conserved PX(D/E)XX(K/R) motif thought to be central to the mechanism of these transporters. The mammalian mitochondrial folate transporter (MFT) is one of a few MCF members, known as the P(I/L)W subfamily, that have evolved a tryptophan residue in place of the (D/E) in the second conserved motif; the function of this substitution (Trp-142) is unclear. Molecular dynamics simulations of the MFT in its explicit membrane environment identified this tryptophan, as well as several other residues lining the transport cavity, to be involved in a series of sequential interactions that coordinated the movement of the tetrahydrofolate substrate within the transport cavity. We probed the function of these residues by mutagenesis. The mutation of every residue identified by molecular dynamics to interact with tetrahydrofolate during simulated transit into the aqueous channel severely impaired folate transport. Mutation of the subfamily-defining tryptophan residue in the MFT to match the MCF consensus at this position (W142D) was incompatible with cell survival. These studies indicate that MFT Trp-142, as well as other residues lining the transporter interior, coordinate tetrahydrofolate descent and positioning of the substrate in the transporter basin. Overall, we identified residues in the walls and at the base of the transport cavity that are involved in substrate recognition by the MFT.

A GGA+U approach to effective electronic correlations in thiolate-ligated iron-oxo (IV) porphyrin.

High-valent oxo-metal complexes exhibit correlated electronic behavior on dense, low-lying electronic state manifolds, presenting challenging systems for electronic structure methods. Among these species, the iron-oxo (IV) porphyrin denoted Compound I occupies a privileged position, serving a broad spectrum of catalytic roles. The most reactive members of this family bear a thiolate axial ligand, exhibiting high activity toward molecular oxygen activation and substrate oxidation. The default approach to such systems has entailed the use of hybrid density functionals or multi-configurational/multireference methods to treat electronic correlation. An alternative approach is presented based on the GGA+U approximation to density functional theory, in which a generalized gradient approximation (GGA) functional is supplemented with a localization correction to treat on-site correlation as inspired by the Hubbard model. The electronic structure of thiolate-ligated iron-oxo (IV) porphyrin and corresponding Coulomb repulsion U are determined both empirically and self-consistently, yielding spin-distributions, state level splittings, and electronic densities of states consistent with prior hybrid functional calculations. Comparison of this detailed electronic structure with model Hamiltonian calculations suggests that the localized 3d iron moments induce correlation in the surrounding electron gas, strengthening local moment formation. This behavior is analogous to strongly correlated electronic systems such as Mott insulators, in which the GGA+U scheme serves as an effective single-particle representation for the full, correlated many-body problem.

Structures and Dynamics of Anionic Lipoprotein Nanodiscs.

Nanolipoprotein particles known as nanodiscs (NDs) have emerged as versatile and powerful tools for the stabilization of membrane proteins permitting a plethora of structural and biophysical studies. Part of their allure is their flexibility to accommodate many types of lipids and precise control of the composition. However, little is known about how variations in lipid composition impact their structures and dynamics. Herein, we investigate how the introduction of the anionic lipid POPG into POPC NDs impacts these features. Small-angle X-ray and neutron scattering (SAXS and SANS) of variable-composition NDs are complemented with molecular dynamics simulations to interrogate how increasing the concern of POPG impacts the ND shape, structure of the lipid core, and the dynamics of the popular membrane scaffold protein, MSP1D1(-). A convenient benefit of including POPG is that it eliminates D2O-induced aggregation observed in pure POPC NDs, permitting studies by SANS at multiple contrasts. SAXS and SANS data could be globally fit to a stacked elliptical cylinder model as well as an extension of the model that accounts for membrane curvature. Fitting to both models supports that the introduction of POPG results in strongly elliptical NDs; however, MD simulations predict the curvature of the membrane, thereby supporting the use of the latter model. Trends in the model-independent parameters suggest that increases in POPG reduce the conformational heterogeneity of the MSP1D1(-), which is in agreement with MD simulations that show that the incorporation of sufficient POPG suppresses disengagement of the N-terminal helix from the lipid core. These studies highlight novel structural changes in NDs in response to an anionic lipid and will inform the interpretation of future structural studies of membrane proteins embedded in NDs of mixed lipid composition.