Probing membrane enhanced protein-protein interactions in a minimal redox complex of cytochrome-P450 and P450-reductase

On the determinants of electron transfer reorganization energy in a cytochrome P450: cytochrome b5 complex. A combined quantum mechanics and molecular dynamics simulation study

The electron transfer steps in the catalytic cycle of cytochrome P450 (CYP) enzymes, ubiquitous proteins with key roles in processes such as drug metabolism and steroidogenesis, are often rate-limiting. To predict ET rates from atomistic molecular dynamics simulations using Marcus theory, values of the reaction free energy ΔG0 and the reorganization free energy λ are required from either experiments or computations. For the reduction of cytochrome P450 17A1 (CYP17A1) by the secondary redox protein cytochrome b5 (CYb5), a critical step in the regulation of steroidogenesis, experimental measurements of λ are not available. We here describe the computation of λ for this system from a combination of molecular mechanics/molecular dynamics simulations and quantum mechanics computations. Our results show that a quantum mechanical treatment of the redox-active cofactors is necessary, even though the surrounding protein and solvent, which are modeled classically, contribute most to the reorganization energy. The values of λ computed for structural ensembles corresponding to two predicted binding modes of the proteins are 1.23 and 1.16 eV. We find that the λ values computed for the individual soluble globular domains of the two proteins sum to approximately the λ values computed for the membrane-bound CYP17A1-CYb5 complex, indicating that additivity can be invoked in a computationally efficient approach to estimating λ values for such protein-protein complexes.

Membrane Proteins in Nanodiscs: Methods and Applications

Membrane proteins, a principal class of drug targets, play indispensable roles in various biological processes and are closely associated with essential life functions. Their study, however, is complicated by their low solubility in aqueous environments and distinctive structural characteristics, necessitating a suitable native-like environment for molecular analysis. Nanodisc technology has revolutionized this field, providing biochemists with a powerful tool to stabilize membrane proteins and significantly enhance their research possibilities. This review outlines the substantial advancements in nanodisc methodologies and applications from 2018 to 2024. We cover the development of various nanodisc models, as well as structural and functional studies of membrane proteins that utilize nanodiscs, highlighting their medical applications.

Effect of Lipid Bilayer Anchoring on the Conformational Properties of the Cytochrome P450 2D6 Binding Site

Human cytochrome P450 (CYP) proteins metabolize 75% of small-molecule pharmaceuticals, which makes structure-based modeling of CYP metabolism and inhibition, bolstered by the current availability of X-ray crystal structures of CYP globular catalytic domains, an attractive prospect. Accounting for this broad metabolic capacity is a combination of the existence of multiple different CYP proteins and the capacity of a single CYP protein to metabolize multiple different small molecules. It is thought that structural plasticity and flexibility contribute to this latter property; therefore, incorporating diverse conformational states of a particular CYP is likely an important consideration in structure-based CYP metabolism and inhibition modeling. While all-atom explicit-solvent molecular dynamics simulations can be used to generate conformational ensembles under biologically relevant conditions, existing CYP crystal structures are of the globular domain only, whereas human CYPs contain N-terminal transmembrane and linker peptides that anchor the globular catalytic domain to the endoplasmic reticulum. To determine whether this can cause significant differences in the sampled binding site conformations, microsecond scale all-atom explicit-solvent molecular dynamics simulations of the CYP2D6 globular domain in an aqueous environment were compared with those of the full-length protein anchored in a POPC lipid bilayer. While bilayer-anchoring damped some structural fluctuations in the globular domain relative to the aqueous simulations, none of the affected residues included binding site pocket residues. Furthermore, clustering of molecular dynamics snapshots based on either pairwise binding site pocket RMSD or volume differences demonstrated a lack of separation of snapshots from the two simulation conditions into different clusters. These results suggest the substantially simpler and computationally cheaper aqueous simulation approach can be used to generate a relevant conformational ensemble of the CYP2D6 binding site for structure-based metabolism and inhibition modeling.

Polymer Nanodiscs and Their Bioanalytical Potential

Membrane proteins (MPs) play a pivotal role in cellular function and are therefore predominant pharmaceutical targets. Although detailed understanding of MP structure and mechanistic activity is invaluable for rational drug design, challenges are associated with the purification and study of MPs. This review delves into the historical developments that became the prelude to currently available membrane mimetic technologies before shining a spotlight on polymer nanodiscs. These are soluble nanosized particles capable of encompassing MPs embedded in a phospholipid ring. The expanding range of reported amphipathic polymer nanodisc materials is presented and discussed in terms of their tolerance to different solution conditions and their nanodisc properties. Finally, the analytical scope of polymer nanodiscs is considered in both the demonstration of basic nanodisc parameters as well as in the elucidation of structures, lipid-protein interactions, and the functional mechanisms of reconstituted membrane proteins. The final emphasis is given to the unique benefits and applications demonstrated for native nanodiscs accessed through a detergent free process.

Nanodisc reconstitution of flavin mononucleotide binding domain of cytochrome-P450-reductase enables high-resolution NMR probing

Cytochrome-P450-reductase transfers electrons to cytochrome-P450 through its flavin mononucleotide binding domain (FBD). Despite the importance of membrane-anchoring for FBD function, studies have focused on its soluble domain lacking the transmembrane-domain. Here we demonstrate that the reconstitution of FBD in nanodiscs enables high-resolution NMR measurements and renders a stable conformation.

Amphipol-facilitated elucidation of the functional tetrameric complex of full-length cytochrome P450 CYP2B4 and NADPH-cytochrome P450 oxidoreductase

Interactions of membrane-bound mammalian cytochromes P450 (CYPs) with NADPH-cytochrome P450 oxidoreductase (POR), which are required for metabolism of xenobiotics, are facilitated by membrane lipids. A variety of membrane mimetics, such as phospholipid liposomes and nanodiscs, have been used to simulate the membrane to form catalytically active CYP:POR complexes. However, the exact mechanism(s) of these interactions are unclear because of the absence of structural information of full-length mammalian CYP:POR complexes in membranes. Herein, we report the use of amphipols (APols) to form a fully functional, soluble, homogeneous preparation of full-length CYP:POR complexes amenable to biochemical and structural study. Incorporation of CYP2B4 and POR into APols resulted in a CYP2B4:POR complex with a stoichiometry of 1:1, which was fully functional in demethylating benzphetamine at a turnover rate of 37.7 ± 2.2 min^-1, with a coupling efficiency of 40%. Interestingly, the stable complex had a molecular weight (M_w) of 338 ± 22 kDa determined by multiangle light scattering, suggestive of a tetrameric complex of 2CYP2B4:2POR embedded in one APol nanoparticle. Moreover, negative stain electron microscopy (EM) validated the homogeneity of the complex and allowed us to generate a three-dimensional EM map and model consistent with the tetramer observed in solution. This first report of the full-length mammalian CYP:POR complex by transmission EM not only reveals the architecture that facilitates electron transfer but also highlights a potential use of APols in biochemical and structural studies of functional CYP complexes with redox partners.

Amphipathic Peptide-Phospholipid Nanofibers: Phospholipid Specificity and Dependence on Concentration and Temperature

The design of nanoassemblies is an important part of the development of new materials for applications in nanomedicine and biosensors. In our previous study, cysteine substitution of the apolipoprotein A-I-derived peptide 18A at residue 11, 18A[A11C], bound to 1-palmitoyl-2-oleoylphosphatidylcholine to form fibrous aggregates at a lipid-to-peptide molar ratio of ≤2 and a fiber diameter of 10-20 nm. However, the mechanisms underlying the lipid-peptide interactions that enable nanofiber formation remain unclear. Here, we evaluated the phospholipid specificity, concentration dependence, and temperature dependence of the formation of 18A[A11C]-lipid nanofibers. Nanofibers were found to form in the presence of specific phospholipids and have a constant lipid/peptide stoichiometry of 1.2 ± 0.2. Moreover, an increase in the length of the acyl chain in phosphatidylcholines was found to increase the structural stability of the nanofibers. These results indicate that specific molecular interactions between peptides and both the headgroups and acyl chains of phospholipids are involved in nanofiber formation. Furthermore, the formation and disassembly of the nanofibers were reversibly controlled by changes in temperature and concentration. The results of the present study provide an insight into the creation of nanoassembling structures.

An electron transfer competent structural ensemble of membrane-bound cytochrome P450 1A1 and cytochrome P450 oxidoreductase

Cytochrome P450 (CYP) heme monooxygenases require two electrons for their catalytic cycle. For mammalian microsomal CYPs, key enzymes for xenobiotic metabolism and steroidogenesis and important drug targets and biocatalysts, the electrons are transferred by NADPH-cytochrome P450 oxidoreductase (CPR). No structure of a mammalian CYP-CPR complex has been solved experimentally, hindering understanding of the determinants of electron transfer (ET), which is often rate-limiting for CYP reactions. Here, we investigated the interactions between membrane-bound CYP 1A1, an antitumor drug target, and CPR by a multiresolution computational approach. We find that upon binding to CPR, the CYP 1A1 catalytic domain becomes less embedded in the membrane and reorients, indicating that CPR may affect ligand passage to the CYP active site. Despite the constraints imposed by membrane binding, we identify several arrangements of CPR around CYP 1A1 that are compatible with ET. In the complexes, the interactions of the CPR FMN domain with the proximal side of CYP 1A1 are supplemented by more transient interactions of the CPR NADP domain with the distal side of CYP 1A1. Computed ET rates and pathways agree well with available experimental data and suggest why the CYP-CPR ET rates are low compared to those of soluble bacterial CYPs.

Nanodiscs: A toolkit for membrane protein science

Membrane proteins are involved in numerous vital biological processes, including transport, signal transduction and the enzymes in a variety of metabolic pathways. Integral membrane proteins account for up to 30% of the human proteome and they make up more than half of all currently marketed therapeutic targets. Unfortunately, membrane proteins are inherently recalcitrant to study using the normal toolkit available to scientists, and one is most often left with the challenge of finding inhibitors, activators and specific antibodies using a denatured or detergent solubilized aggregate. The Nanodisc platform circumvents these challenges by providing a self-assembled system that renders typically insoluble, yet biologically and pharmacologically significant, targets such as receptors, transporters, enzymes, and viral antigens soluble in aqueous media in a native-like bilayer environment that maintain a target's functional activity. By providing a bilayer surface of defined composition and structure, Nanodiscs have found great utility in the study of cellular signaling complexes that assemble on a membrane surface. Nanodiscs provide a nanometer scale vehicle for the in vivo delivery of amphipathic drugs, therapeutic lipids, tethered nucleic acids, imaging agents and active protein complexes. This means for generating nanoscale lipid bilayers has spawned the successful use of numerous other polymer and peptide amphipathic systems. This review, in celebration of the Anfinsen Award, summarizes some recent results and provides an inroad into the current and historical literature.

Dissecting the Structural Plasticity and Dynamics of Cytochrome P450 2B4 by Molecular Dynamics Simulations

The plasticity of cytochromes P450 (P450s) is known to contribute significantly to their catalytic capacity of metabolizing various substrates. Although numerous studies have been performed, factors governing the plasticity and dynamics of P450s are still not fully understood. In this study, taking CYP2B4 as an example, we dissect the protein plasticity and dynamics in different environments. CYP2B4 is featured by a high degree of plasticity, which exhibits open, closed, and intermediate states. By analyzing the CYP2B4 crystal structures, we identified the structural features for the closed, open, and intermediate states. Interestingly, formation of the dimer structure was found in the open and intermediate states. The subsequent molecular dynamics (MD) simulations of the open structure in water confirmed the importance of the dimer form in stabilizing the open conformations. MD simulations of the closed and open structures in the membrane environment and the free energies for opening the F-G cassette obtained from the umbrella sampling calculations indicate that the membrane environment is important for stabilizing the F-G cassette. The dynamical network analysis indicates that Asp105 on the B-C loop plays an important role in transiting the structure from the open to the intermediate state. Our results thus unveil the mechanisms of dimer formation and open-to-intermediate transition for CYP2B4 in the water and membrane environments.

Effector role of cytochrome P450 reductase for androstenedione binding to human aromatase

Cytochromes P450 constitute a large superfamily of monooxygenases involved in many metabolic pathways. Most of them are not self-sufficient and need a reductase protein to provide the electrons necessary for catalysis. It was shown that the redox partner plays a role in the modulation of the structure and function of some bacterial P450 enzymes. Here, the effect of NADPH-cytochrome reductase (CPR) on human aromatase (Aro) is studied for what concerns its role in substrate binding. Pre-steady-state kinetic experiments indicate that both the substrate binding rates and the percentage of spin shift detected for aromatase are increased when CPR is present. Moreover, aromatase binds the substrate through a conformational selection mechanism, suggesting a possible effector role of CPR. The thermodynamic parameters for the formation of the CPR-Aro complex were studied by isothermal titration calorimetry. The dissociation constant of the complex formation is 4.5 folds lower for substrate-free compared to the substrate-bound enzyme. The enthalpy change observed when the CPR-Aro complex forms in the absence of the substrate are higher than in its presence, indicating that more interactions are formed/broken in the former case. Taken together, our data confirm that CPR has a role in promoting aromatase conformation optimal for substrate binding.

Two Distinct Structures of Membrane-Associated Homodimers of GTP- and GDP-Bound KRAS4B Revealed by Paramagnetic Relaxation Enhancement

KRAS homo-dimerization has been implicated in the activation of RAF kinases, however, the mechanism and structural basis remain elusive. We developed a system to study KRAS dimerization on nanodiscs using paramagnetic relaxation enhancement (PRE) NMR spectroscopy, and determined distinct structures of membrane-anchored KRAS dimers in the active GTP- and inactive GDP-loaded states. Both dimerize through an α4-α5 interface, but the relative orientation of the protomers and their contacts differ substantially. Dimerization of KRAS-GTP, stabilized by electrostatic interactions between R135 and E168, favors an orientation on the membrane that promotes accessibility of the effector-binding site. Remarkably, "cross"-dimerization between GTP- and GDP-bound KRAS molecules is unfavorable. These models provide a platform to elucidate the structural basis of RAF activation by RAS and to develop inhibitors that can disrupt the KRAS dimerization. The methodology is applicable to many other farnesylated small GTPases.

Exploring the role of the membrane bilayer in the recognition of candesartan by its GPCR AT1 receptor

Cardiovascular diseases and hypertension in particular are major health risks worldwide and the improvement on their treatment will be beneficial for the human health. AT1R antagonists belong to the sartans family that targets the renin-angiotensin aldosterone system (RAAS) through blocking the hormone angiotensin II to exert its detrimental effects in pathological states. As a consequence, they are beneficial to treat hypertension, diabetes related kidney failure and hyperaemic episodes. Long unbiased Molecular Dynamics (MD) simulations are performed in order to explore candesartan's possible 2D and 3D diffusion mechanisms towards AT1R receptor. 3D diffusion mechanism is referred to the direct binding of the AT1 antagonist candesartan to the AT1R 3D structure (PDB ID: 4YAY). 2D diffusion mechanism involves first, the incorporation of candesartan in the bilayer core and then its localization on the AT1R binding cavity, through a diffusion mechanism. The obtained results indicate that membranes interact significantly with the neutral form of candesartan, which is indeed approaching the receptors' active site through diffusion via the lipids. On the other hand, the deprotonated form of the drug is interacting with AT1R's extracellular loop and fails to enter the membrane, pointing out the importance of the pH microenvironment around the receptor. To validate the calculated diffusion coefficients of the drug in the lipid bilayers 2D DOSY NMR experiments were recorded and they were in good agreement. Information on the impact that has the interaction of candesartan with the membrane is very important for the rationally design and development of potent ARBs. Thus, its conformational features as well as its localization in the membrane core have to be thoroughly explored.

Probing protein-protein and protein-substrate interactions in the dynamic membrane-associated ternary complex of cytochromes P450, b5, and reductase

Cytochrome P450 (cytP450) interacts with two redox partners, cytP450 reductase and cytochrome-b5, to metabolize substrates. Using NMR, we reveal changes in the dynamic interplay when all three proteins are incorporated into lipid nanodiscs in the absence and presence of substrates.

Conformational turn triggers regio-selectivity in the bioactivation of thiophene-contained compounds mediated by cytochrome P450

In the present work, we performed Density Functional Theory calculations to explore the bioactivation mechanism of thiophene-containing molecules mediated by P450s. For this purpose, relatively large size compounds, 2,5-diaminothiophene derivatives were selected particularly for this investigation. Here we found the successive regio-selectivity triggered by conformational turn played a significant role in the occurrence of bioactivation. 2,5-Diaminothiophene was oxidized to a 2,5-diimine thiophene-reactive intermediate by Compound I (Cpd I) through successive activations of two N-H bonds (H3-N11 and H1-N6). This reaction exhibited three special characteristics: (1) self-controlled regio-selectivity during the oxidation process. There was a large scale of conformational turn in the abstraction of the first H atom which triggers the selection of the second H for abstraction. (2) Proton-shuttle mechanism. In high spin (HS) state, proton-shuttle mechanism was observed for the abstraction of the second H atom. (3) Spin-selective manner. In protein environment, the energy barrier in HS state was much lower than that in low spin state. The novel proposed bioactivation mechanism of 2,5-diaminothiophene compounds can help us in rational design of thiophene-contained drugs avoiding the occurrence of bioactivation.