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

Structural insights into the semiquinone form of human Cytochrome P450 reductase by DEER distance measurements between a native flavin and a spin labelled non-canonical amino acid

The flavoprotein Cytochrome P450 reductase (CPR) is the unique electron pathway from NADPH to Cytochrome P450 (CYPs). The conformational dynamics of human CPR in solution, which involves transitions from a "locked/closed" to an "unlocked/open" state, is crucial for electron transfer. To date, however, the factors guiding these changes remain unknown. By Site-Directed Spin Labelling coupled to Electron Paramagnetic Resonance spectroscopy, we have incorporated a non-canonical amino acid onto the flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) domains of soluble human CPR, and labelled it with a specific nitroxide spin probe. Taking advantage of the endogenous FMN cofactor, we successfully measured for the first time, the distance distribution by DEER between the semiquinone state FMNH• and the nitroxide. The DEER data revealed a salt concentration-dependent distance distribution, evidence of an "open" CPR conformation at high salt concentrations exceeding previous reports. We also conducted molecular dynamics simulations which unveiled a diverse ensemble of conformations for the "open" semiquinone state of the CPR at high salt concentration. This study unravels the conformational landscape of the one electron reduced state of CPR, which had never been studied before.

Single Mutations in Cytochrome P450 Oxidoreductase Can Alter the Specificity of Human Cytochrome P450 1A2-Mediated Caffeine Metabolism

A unique cytochrome P450 (CYP) oxidoreductase (CPR) sustains activities of human microsomal CYPs. Its function requires toggling between a closed conformation enabling electron transfers from NADPH to FAD and then FMN cofactors and open conformations forming complexes and transferring electrons to CYPs. We previously demonstrated that distinct features of the hinge region linking the FAD and FMN domain (FD) modulate conformer poses and their interactions with CYPs. Specific FD residues contribute in a CYP isoform-dependent manner to the recognition and electron transfer mechanisms that are additionally modulated by the structure of CYP-bound substrate. To obtain insights into the underlying mechanisms, we analyzed how hinge region and FD mutations influence CYP1A2-mediated caffeine metabolism. Activities, metabolite profiles, regiospecificity and coupling efficiencies were evaluated in regard to the structural features and molecular dynamics of complexes bearing alternate substrate poses at the CYP active site. Studies reveal that FD variants not only modulate CYP activities but surprisingly the regiospecificity of reactions. Computational approaches evidenced that the considered mutations are generally in close contact with residues at the FD-CYP interface, exhibiting induced fits during complexation and modified dynamics depending on caffeine presence and orientation. It was concluded that dynamic coupling between FD mutations, the complex interface and CYP active site exist consistently with the observed regiospecific alterations.

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.

Polymer-Nanodiscs as a Novel Alignment Medium for High-Resolution NMR-Based Structural Studies of Nucleic Acids

Residual dipolar couplings (RDCs) are increasingly used for high-throughput NMR-based structural studies and to provide long-range angular constraints to validate and refine structures of various molecules determined by X-ray crystallography and NMR spectroscopy. RDCs of a given molecule can be measured in an anisotropic environment that aligns in an external magnetic field. Here, we demonstrate the first application of polymer-based nanodiscs for the measurement of RDCs from nucleic acids. Polymer-based nanodiscs prepared using negatively charged SMA-EA polymer and zwitterionic DMPC lipids were characterized by size-exclusion chromatography, 1H NMR, dynamic light-scattering, and 2H NMR. The magnetically aligned polymer-nanodiscs were used as an alignment medium to measure RDCs from a 13C/15N-labeled fluoride riboswitch aptamer using 2D ARTSY-HSQC NMR experiments. The results showed that the alignment of nanodiscs is stable for nucleic acids and nanodisc-induced RDCs fit well with the previously determined solution structure of the riboswitch. These results demonstrate that SMA-EA-based lipid-nanodiscs can be used as a stable alignment medium for high-resolution structural and dynamical studies of nucleic acids, and they can also be applicable to study various other biomolecules and small molecules in general.

Angiotensin converting enzyme inhibitors from medicinal plants: a molecular docking and dynamic simulation approach

Angiotensin converting enzyme (ACE) is a key enzyme and mediator in the aetiology of high blood pressure (HBP) and hypertension. As one of the leading cause of untimely death worldwide, there is a lot of research and studies on the management and treatment of hypertension. The usage of medicinal plants in the management of hypertension as alternative to synthetic allopathic drugs is a common practice in folkloric and traditional medicine. Therefore, this study was aimed to investigate the ACE inhibitory activity of some medicinal plants which are commonly used in the treatment of HBP in southwestern part of Nigeria using extensive in-silico approach. Compounds identified in the plants through GC-MS technique, together with Lisinopril were docked against ACE protein. It was observed that only 40 of the compounds had binding affinity ≥ - 6.8 kcal/mol which was demonstrated by the standard drug (lisinopril). Interaction between the compounds and ACE was via conventional hydrogen, carbon hydrogen, alkyl, pi-alkyl, pi-carbon, and Van Der Wall bonds among others. Most of these compounds exhibited drug like properties, without violating majority of the physicochemical descriptors and Lipinski rule of 5. The ADMET evaluation revealed that only 2 compounds (cyclopentadecanone and oxacycloheptadecan-2-one) which were identified in Bacopa florinbunda plant were predicted non-toxic and thus were subjected to molecular dynamics and simulation with ACE. From the molecular dynamics and mechanics analysis, both cyclopentadecanone and oxacycloheptadecan-2-one showed high stability and inhibitory potentials when bound to ACE. Oxacycloheptadecan-2-one was more stable than lisinopril and cyclopentadecanone in the ligand-ACE complex; we therefore suggested its experimental and clinical validation as drug candidates for the treatment of hypertension.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-022-00135-z.

Detergent-Free Isolation of Membrane Proteins and Strategies to Study Them in a Near-Native Membrane Environment

Atomic-resolution structural studies of membrane-associated proteins and peptides in a membrane environment are important to fully understand their biological function and the roles played by them in the pathology of many diseases. However, the complexity of the cell membrane has severely limited the application of commonly used biophysical and biochemical techniques. Recent advancements in NMR spectroscopy and cryoEM approaches and the development of novel membrane mimetics have overcome some of the major challenges in this area. For example, the development of a variety of lipid-nanodiscs has enabled stable reconstitution and structural and functional studies of membrane proteins. In particular, the ability of synthetic amphipathic polymers to isolate membrane proteins directly from the cell membrane, along with the associated membrane components such as lipids, without the use of a detergent, has opened new avenues to study the structure and function of membrane proteins using a variety of biophysical and biological approaches. This review article is focused on covering the various polymers and approaches developed and their applications for the functional reconstitution and structural investigation of membrane proteins. The unique advantages and limitations of the use of synthetic polymers are also discussed.

QSAR models reveal new EPAC-selective allosteric modulators

Exchange proteins directly activated by cAMP (EPAC) are guanine nucleotide exchange factors for the small GTPases, Rap1 and Rap2. They regulate several physiological functions and mitigation of their activity has been suggested as a possible treatment for multiple diseases such as cardiomyopathy, diabetes, chronic pain, and cancer. Several EPAC-specific modulators have been developed, however studies that quantify their structure–activity relationships are still lacking. Here we propose a quantitative structure–activity relationship (QSAR) model for a series of EPAC-specific compounds. The model demonstrated high reproducibility and predictivity and the predictive ability of the model was tested against a series of compounds that were unknown to the model. The compound with the highest predicted affinity was validated experimentally through fluorescence-based competition assays and NMR experiments revealed its mode of binding and mechanism of action as a partial agonist. The proposed QSAR model can, therefore, serve as an effective screening tool to identify promising EPAC-selective drug leads with enhanced potency.

QSAR models of EPAC-specific allosteric ligands predict the affinity of a promising analogue.

A Cytochrome B5-Like Heme/Steroid Binding Domain Protein, PlCB5L1, Regulates Mycelial Growth, Pathogenicity and Oxidative Stress Tolerance in Peronophythora litchii

As an electron transport component, cytochrome b₅ is an essential component of the Class II cytochrome P450 monooxygenation system and widely present in animals, plants, and fungi. However, the roles of Cyt-b₅ domain proteins in pathogenic oomycetes remain unknown. Peronophythora litchii is an oomycete pathogen that causes litchi downy blight, the most destructive disease of litchi. In this study, we identified a gene, designated PlCB5L1, that encodes a Cyt-b₅ domain protein in P. litchii, and characterized its function. PlCB5L1 is highly expressed in the zoospores, cysts, germinated cysts, and during early stages of infection. PlCB5L1 knockout mutants showed reduced growth rate and β-sitosterol utilization. Importantly, we also found that PlCB5L1 is required for the full pathogenicity of P. litchii. Compared with the wild-type strain, the PlCB5L1 mutants exhibited significantly higher tolerance to SDS and sorbitol, but impaired tolerance to cell wall stress, osmotic stress, and oxidative stress. Further, the expression of genes involved in oxidative stress tolerance, including peroxidase, cytochrome P450, and laccase genes, were down-regulated in PlCB5L1 mutants under oxidative stress. This is the first report that a Cyt-b₅ domain protein contributes to the development, stress response, and pathogenicity in plant pathogenic oomycetes.

Lipid Nanodiscs for High-Resolution NMR Studies of Membrane Proteins

Membrane proteins (MPs) play essential roles in numerous cellular processes. Because around 70% of the currently marketed drugs target MPs, a detailed understanding of their structure, binding properties, and functional dynamics in a physiologically relevant environment is crucial for a more detailed understanding of this important protein class. We here summarize the benefits of using lipid nanodiscs for NMR structural investigations and provide a detailed overview of the currently used lipid nanodisc systems as well as their applications in solution-state NMR. Despite the increasing use of other structural methods for the structure determination of MPs in lipid nanodiscs, solution NMR turns out to be a versatile tool to probe a wide range of MP features, ranging from the structure determination of small to medium-sized MPs to probing ligand and partner protein binding as well as functionally relevant dynamical signatures in a lipid nanodisc setting. We will expand on these topics by discussing recent NMR studies with lipid nanodiscs and work out a key workflow for optimizing the nanodisc incorporation of an MP for subsequent NMR investigations. With this, we hope to provide a comprehensive background to enable an informed assessment of the applicability of lipid nanodiscs for NMR studies of a particular MP of interest.

The Central Role of Cytochrome P450 in Xenobiotic Metabolism-A Brief Review on a Fascinating Enzyme Family

Human Cytochrome P450 (CYP) enzymes constitute a superfamily of membrane-bound hemoproteins that are responsible for the metabolism of a wide variety of clinically, physiologically, and toxicologically important compounds. These heme-thiolate monooxygenases play a pivotal role in the detoxification of xenobiotics, participating in the metabolism of many structurally diverge compounds. This short-review is intended to provide a summary on the major roles of CYPs in Phase I xenobiotic metabolism. The manuscript is focused on eight main topics that include the most relevant aspects of past and current CYP research. Initially, (I) a general overview of the main aspects of absorption, distribution, metabolism, and excretion (ADME) of xenobiotics are presented. This is followed by (II) a background overview on major achievements in the past of the CYP research field. (III) Classification and nomenclature of CYPs is briefly reviewed, followed by (IV) a summary description on CYP’s location and function in mammals. Subsequently, (V) the physiological relevance of CYP as the cornerstone of Phase I xenobiotic metabolism is highlighted, followed by (VI) reviewing both genetic determinants and (VI) nongenetic factors in CYP function and activity. The last topic of the review (VIII) is focused on the current challenges of the CYP research field.

Differential effects on human cytochromes P450 by CRISPR/Cas9-induced genetic knockout of cytochrome P450 reductase and cytochrome b5 in HepaRG cells

HepaRG cells are increasingly accepted as model for human drug metabolism and other hepatic functions. We used lentiviral transduction of undifferentiated HepaRG cells to deliver Cas9 and two alternative sgRNAs targeted at NADPH:cytochrome P450 oxidoreductase (POR), the obligate electron donor for microsomal cytochromes P450 (CYP). Cas9-expressing HepaRG^VC (vector control) cells were phenotypically similar to wild type HepaRG cells and could be differentiated into hepatocyte-like cells by DMSO. Genetic POR-knockout resulted in phenotypic POR knockdown of up to 90% at mRNA, protein, and activity levels. LC–MS/MS measurement of seven CYP-activities showed differential effects of POR-knockdown with CYP2C8 being least and CYP2C9 being most affected. Further studies on cytochrome b5 (CYB5), an alternative NADH-dependent electron donor indicated particularly strong support of CYP2C8-dependent amodiaquine N-deethylation by CYB5 and this was confirmed by genetic CYB5 single- and POR/CYB5 double-knockout. POR-knockdown also affected CYP expression on mRNA and protein level, with CYP1A2 being induced severalfold, while CYP2C9 was strongly downregulated. In summary our results show that POR/NADPH- and CYB5/NADH-electron transport systems influence human drug metabolizing CYPs differentially and differently than mouse Cyps. Our Cas9-expressing HepaRG^VC cells should be suitable to study the influence of diverse genes on drug metabolism and other hepatic functions.

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.

Interaction Modes of Microsomal Cytochrome P450s with Its Reductase and the Role of Substrate Binding

The activity of microsomal cytochromes P450 (CYP) is strictly dependent on the supply of electrons provided by NADPH cytochrome P450 oxidoreductase (CPR). The variant nature of the isoform-specific proximal interface of microsomal CYPs implies that the interacting interface between the two proteins is degenerated. Recently, we demonstrated that specific CPR mutations in the FMN-domain (FD) may induce a gain in activity for a specific CYP isoform. In the current report, we confirm the CYP isoform dependence of CPR's degenerated binding by demonstrating that the effect of four of the formerly studied FD mutants are indeed exclusive of a specific CYP isoform, as verified by cytochrome c inhibition studies. Moreover, the nature of CYP's substrate seems to have a modulating role in the CPR:CYP interaction. In silico molecular dynamics simulations of the FD evidence that mutations induces very subtle structural alterations, influencing the characteristics of residues formerly implicated in the CPR:CYP interaction or in positioning of the FMN moiety. CPR seems therefore to be able to form effective interaction complexes with its structural diverse partners via a combination of specific structural features of the FD, which are functional in a CYP isoform dependent manner, and dependent on the substrate bound.

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.

An NMR based phosphodiesterase assay

We propose a phosphodiesterase assay based on 1D ¹H NMR to monitor the hydrolysis of cyclic nucleotides directly, without requiring tags or the addition of exogenous reagents. The method is suitable to measure phosphodiesterase K_M and k_cat parameters and to identify phosphodiesterase inhibitors.

The Role of the FMN-Domain of Human Cytochrome P450 Oxidoreductase in Its Promiscuous Interactions With Structurally Diverse Redox Partners

NADPH cytochrome P450 oxidoreductase (CPR) is the obligatory electron supplier that sustains the activity of microsomal cytochrome P450 (CYP) enzymes. The variant nature of the isoform-specific proximal interface of microsomal CYPs indicates that CPR is capable of multiple degenerated interactions with CYPs for electron transfer, through different binding mechanisms, and which are still not well-understood. Recently, we showed that CPR dynamics allows formation of open conformations that can be sampled by its structurally diverse redox partners in a CYP-isoform dependent manner. To further investigate the role of the CPR FMN-domain in effective binding of CPR to its diverse acceptors and to clarify the underlying molecular mechanisms, five different CPR-FMN-domain random mutant libraries were created. These libraries were screened for mutants with increased activity when combined with specific CYP-isoforms. Seven CPR-FMN-domain mutants were identified, supporting a gain in activity for CYP1A2 (P117H, G144C, A229T), 2A6 (P117L/L125V, G175D, H183Y), or 3A4 (N151D). Effects were evaluated using extended enzyme kinetic analysis, cytochrome b₅ competition, ionic strength effect on CYP activity, and structural analysis. Mutated residues were located either in or adjacent to several acidic amino acid stretches – formerly indicated to be involved in CPR:CYP interactions – or close to two tyrosine residues suggested to be involved in FMN binding. Several of the identified positions co-localize with mutations found in naturally occurring CPR variants that were previously shown to cause CYP-isoform-dependent effects. The mutations do not seem to significantly alter the geometry of the FMN-domain but are likely to cause very subtle alterations leading to improved interaction with a specific CYP. Overall, these data suggest that CYPs interact with CPR using an isoform specific combination of several binding motifs of the FMN-domain.