Permeability is an important molecular property in drug discovery, as it co-determines pharmacokinetics whenever a drug crosses the phospholipid bilayer, e.g., into the cell, in the gastrointestinal tract, or across the blood-brain barrier. Many methods for the determination of permeability have been developed, including cell line assays (CACO-2 and MDCK), cell-free model systems like parallel artificial membrane permeability assay (PAMPA) mimicking, e.g., gastrointestinal epithelia or the skin, as well as the black lipid membrane (BLM) and submicrometer liposomes. Furthermore, many in silico approaches have been developed for permeability prediction: meta-analysis of publicly available databases for permeability data (MolMeDB and ChEMBL) was performed to establish their usability. Four experimental and two computational methods were evaluated. It was shown that repeatability of the reported permeability measurement is not great even for the same method. For the PAMPA method, two different permeabilities are reported: intrinsic and apparent. They can vary in degrees of magnitude; thus, we suggest being extra cautious using literature data on permeability. When we compared data for the same molecules using different methods, the best agreement was between cell-based methods and between BLM and computational methods. Existence of unstirred water layer (UWL) permeability limits the data agreement between cell-based methods (and apparent PAMPA) with data that are not limited by UWL permeability (computational methods, BLM, intrinsic PAMPA). Therefore, different methods have different limitations. Cell-based methods provide results only in a small range of permeabilities (-8 to -4 in cm/s), and computational methods can predict a wider range of permeabilities beyond physical limitations, but their precision is therefore limited. BLM with liposomes can be used for both fast and slow permeating molecules, but its usage is more complicated than standard transwell techniques. To sum up, when working with in-house measured or published permeability data, we recommend caution in interpreting and combining them.

• Klíčová slova
• BLM, CACO-2, COSMOperm, MDCK, MolMeDB, PAMPA, PerMM, liposome, membrane, permeability,
• MeSH
• buňky MDCK MeSH
• Caco-2 buňky MeSH
• hematoencefalická bariéra metabolismus   MeSH
• lidé MeSH
• liposomy metabolismus   chemie   MeSH
• permeabilita buněčné membrány * MeSH
• permeabilita MeSH
• počítačová simulace MeSH
• psi MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• psi MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• metaanalýza MeSH
• Názvy látek
• liposomy MeSH

Lipid-mediated delivery of active pharmaceutical ingredients (API) opened new possibilities in advanced therapies. By encapsulating an API into a lipid nanocarrier (LNC), one can safely deliver APIs not soluble in water, those with otherwise strong adverse effects, or very fragile ones such as nucleic acids. However, for the rational design of LNCs, a detailed understanding of the composition-structure-function relationships is missing. This review presents currently available computational methods for LNC investigation, screening, and design. The state-of-the-art physics-based approaches are described, with the focus on molecular dynamics simulations in all-atom and coarse-grained resolution. Their strengths and weaknesses are discussed, highlighting the aspects necessary for obtaining reliable results in the simulations. Furthermore, a machine learning, i.e., data-based learning, approach to the design of lipid-mediated API delivery is introduced. The data produced by the experimental and theoretical approaches provide valuable insights. Processing these data can help optimize the design of LNCs for better performance. In the final section of this Review, state-of-the-art of computer simulations of LNCs are reviewed, specifically addressing the compatibility of experimental and computational insights.

• Klíčová slova
• ionizable lipid, lipid nanocarrier, lipid nanoparticle, liposome, molecular simulation, vesicle,
• MeSH
• léčivé přípravky chemie   aplikace a dávkování   MeSH
• lékové transportní systémy * metody   MeSH
• lidé MeSH
• lipidy * chemie   MeSH
• nanočástice chemie   MeSH
• nosiče léků * chemie   MeSH
• počítačová simulace MeSH
• simulace molekulární dynamiky MeSH
• strojové učení MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• přehledy MeSH
• Názvy látek
• léčivé přípravky MeSH
• lipidy * MeSH
• nosiče léků * MeSH

Cytochrome P450 1A2 (P450 1A2, CYP1A2) is a membrane-bound enzyme that oxidizes a broad range of hydrophobic substrates. The structure and dynamics of both the catalytic and trans-membrane (TM) domains of this enzyme in the membrane/water environment were investigated using a multiscale computational approach, including coarse-grained and all-atom molecular dynamics. Starting from the spontaneous self-assembly of the system containing the TM or soluble domain immersed in randomized dilauroyl phosphatidylcholine (DLPC)/water mixture into their respective membrane-bound forms, we reconstituted the membrane-bound structure of the full-length P450 1A2. This structure includes a TM helix that spans the membrane, while being connected to the catalytic domain by a short flexible loop. Furthermore, in this model, the upper part of the TM helix interacts directly with a conserved and highly hydrophobic N-terminal proline-rich segment of the catalytic domain; this segment and the FG loop are immersed in the membrane, whereas the remaining portion of the catalytic domain remains exposed to aqueous solution. The shallow membrane immersion of the catalytic domain induces a depression in the opposite intact layer of the phospholipids. This structural effect may help in stabilizing the position of the TM helix directly beneath the catalytic domain. The partial immersion of the catalytic domain also allows for the enzyme substrates to enter the active site from either aqueous solution or phospholipid environment via several solvent- and membrane-facing tunnels in the full-length P450 1A2. The calculated tunnel dynamics indicated that the opening probability of the membrane-facing tunnels is significantly enhanced when a DLPC molecule spontaneously penetrates into the membrane-facing tunnel 2d. The energetics of the lipid penetration process were assessed by the linear interaction energy (LIE) approximation, and found to be thermodynamically feasible.

• MeSH
• cytochrom P-450 CYP1A2 chemie   MeSH
• fosfatidylcholiny MeSH
• fosfolipidy chemie   MeSH
• katalytická doména MeSH
• katalýza MeSH
• lidé MeSH
• simulace molekulární dynamiky MeSH
• vazba proteinů MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• 1,2-dilauroylphosphatidylcholine MeSH Prohlížeč
• cytochrom P-450 CYP1A2 MeSH
• fosfatidylcholiny MeSH
• fosfolipidy MeSH

Benzo[a]pyrene (BaP) is a human carcinogen that covalently binds to DNA after activation by cytochrome P450 (P450). Here, we investigated whether NADH:cytochrome b5 reductase (CBR) in the presence of cytochrome b5 can act as sole electron donor to human P450 1A1 during BaP oxidation and replace the canonical NADPH:cytochrome P450 reductase (POR) system. We also studied the efficiencies of the coenzymes of these reductases, NADPH as a coenzyme of POR, and NADH as a coenzyme of CBR, to mediate BaP oxidation. Two systems containing human P450 1A1 were utilized: human recombinant P450 1A1 expressed with POR, CBR, epoxide hydrolase, and cytochrome b5 in Supersomes and human recombinant P450 1A1 reconstituted with POR and/or with CBR and cytochrome b5 in liposomes. BaP-9,10-dihydrodiol, BaP-7,8-dihydrodiol, BaP-1,6-dione, BaP-3,6-dione, BaP-9-ol, BaP-3-ol, a metabolite of unknown structure, and two BaP-DNA adducts were generated by the P450 1A1-Supersomes system, both in the presence of NADPH and in the presence of NADH. The major BaP-DNA adduct detected by (32)P-postlabeling was characterized as 10-(deoxyguanosin-N(2)-yl)-7,8,9-trihydroxy-7,8,9,10-tetrahydro-BaP (assigned adduct 1), while the minor adduct is probably a guanine adduct derived from 9-hydroxy-BaP-4,5-epoxide (assigned adduct 2). BaP-3-ol as the major metabolite, BaP-9-ol, BaP-1,6-dione, BaP-3,6-dione, an unknown metabolite, and adduct 2 were observed in the system using P450 1A1 reconstituted with POR plus NADPH. When P450 1A1 was reconstituted with CBR and cytochrome b5 plus NADH, BaP-3-ol was the predominant metabolite too, and an adduct 2 was also generated. Our results demonstrate that the NADH/cytochrome b5/CBR system can act as the sole electron donor both for the first and second reduction of P450 1A1 during the oxidation of BaP in vitro. They suggest that NADH-dependent CBR can replace NADPH-dependent POR in the P450 1A1-catalyzed metabolism of BaP.

• MeSH
• adukty DNA metabolismus   MeSH
• benzopyren toxicita   MeSH
• cytochrom-B(5)-reduktasa metabolismus   MeSH
• lidé MeSH
• oxidace-redukce MeSH
• Check Tag
• lidé MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• adukty DNA MeSH
• benzopyren MeSH
• cytochrom-B(5)-reduktasa MeSH

This symposium summary, sponsored by the ASPET, was held at Experimental Biology 2015 on March 29, 2015, in Boston, Massachusetts. The symposium focused on: 1) the interactions of cytochrome P450s (P450s) with their redox partners; and 2) the role of the lipid membrane in their orientation and stabilization. Two presentations discussed the interactions of P450s with NADPH-P450 reductase (CPR) and cytochrome b5. First, solution nuclear magnetic resonance was used to compare the protein interactions that facilitated either the hydroxylase or lyase activities of CYP17A1. The lyase interaction was stimulated by the presence of b5 and 17α-hydroxypregnenolone, whereas the hydroxylase reaction was predominant in the absence of b5. The role of b5 was also shown in vivo by selective hepatic knockout of b5 from mice expressing CYP3A4 and CYP2D6; the lack of b5 caused a decrease in the clearance of several substrates. The role of the membrane on P450 orientation was examined using computational methods, showing that the proximal region of the P450 molecule faced the aqueous phase. The distal region, containing the substrate-access channel, was associated with the membrane. The interaction of NADPH-P450 reductase (CPR) with the membrane was also described, showing the ability of CPR to "helicopter" above the membrane. Finally, the endoplasmic reticulum (ER) was shown to be heterogeneous, having ordered membrane regions containing cholesterol and more disordered regions. Interestingly, two closely related P450s, CYP1A1 and CYP1A2, resided in different regions of the ER. The structural characteristics of their localization were examined. These studies emphasize the importance of P450 protein organization to their function.

• MeSH
• buněčná membrána metabolismus   MeSH
• endoplazmatické retikulum metabolismus   MeSH
• interakční proteinové domény a motivy fyziologie   MeSH
• jaterní mikrozomy metabolismus   MeSH
• lidé MeSH
• sekundární struktura proteinů MeSH
• systém (enzymů) cytochromů P-450 chemie   fyziologie   MeSH
• výzkumná zpráva * MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• systém (enzymů) cytochromů P-450 MeSH

ABSTRACT: Benzo[a]pyrene (BaP) is a human carcinogen that covalently binds to DNA after metabolic activation by cytochrome P450 (CYP) enzymes. Here we investigated the efficiencies of rat hepatic microsomes and rat recombinant CYP1A1 expressed with its reductase, NADPH:CYP oxidoreductase (POR), NADH:cytochrome b5 reductase, epoxide hydrolase and/or cytochrome b5 in Supersomes™ to metabolize this carcinogen. We also studied the effectiveness of coenzymes of two of the microsomal reductases, NADPH as a coenzyme of POR, and NADH as a coenzyme of NADH:cytochrome b5 reductase, to mediate BaP metabolism in these systems. Up to eight BaP metabolites and two DNA adducts were generated by the systems, both in the presence of NADPH and NADH. Among BaP metabolites, BaP-9,10-dihydrodiol, BaP-4,5-dihydrodiol, BaP-7,8-dihydrodiol, BaP-1,6-dione, BaP-3,6-dione, BaP-9-ol, BaP-3-ol, and a metabolite of unknown structure were formed by hepatic microsomes and rat CYP1A1. One of two DNA adducts formed by examined enzymatic systems (rat hepatic microsomes and rat CYP1A1) was characterized to be 10-(deoxyguanosin-N2-yl)-7,8,9-trihydroxy-7,8,9,10-tetrahydrobenzo[a]pyrene (dG-N2-BPDE), while another adduct has similar chromatographic properties on polyethylaneimine-cellulose thin layer chromatography to a guanine adduct derived from reaction with 9-hydroxy-BaP-4,5-oxide. In the presence of either of the reductase cofactors tested, NADPH or NADH, cytochrome b5 stimulated CYP1A1-mediated formation of both BaP-DNA adducts. The results demonstrate that NADH can act as a sole electron donor for both the first and the second reduction of CYP1A1 during its reaction cycle catalyzing oxidation of BaP, and suggest that the NADH:cytochrome b5 reductase as the NADH-dependent reductase might substitute POR in this enzymatic system.

• Klíčová slova
• Coenzymes, DNA, Enzymes, High-pressure liquid chromatography,
• Publikační typ
• časopisecké články MeSH

BACKGROUND: Channels and pores in biomacromolecules (proteins, nucleic acids and their complexes) play significant biological roles, e.g., in molecular recognition and enzyme substrate specificity. RESULTS: We present an advanced software tool entitled MOLE 2.0, which has been designed to analyze molecular channels and pores. Benchmark tests against other available software tools showed that MOLE 2.0 is by comparison quicker, more robust and more versatile. As a new feature, MOLE 2.0 estimates physicochemical properties of the identified channels, i.e., hydropathy, hydrophobicity, polarity, charge, and mutability. We also assessed the variability in physicochemical properties of eighty X-ray structures of two members of the cytochrome P450 superfamily. CONCLUSION: Estimated physicochemical properties of the identified channels in the selected biomacromolecules corresponded well with the known functions of the respective channels. Thus, the predicted physicochemical properties may provide useful information about the potential functions of identified channels. The MOLE 2.0 software is available at http://mole.chemi.muni.cz.

• Publikační typ
• časopisecké články MeSH

Biomolecular channels play important roles in many biological systems, e.g. enzymes, ribosomes and ion channels. This article introduces a web-based interactive MOLEonline 2.0 application for the analysis of access/egress paths to interior molecular voids. MOLEonline 2.0 enables platform-independent, easy-to-use and interactive analyses of (bio)macromolecular channels, tunnels and pores. Results are presented in a clear manner, making their interpretation easy. For each channel, MOLEonline displays a 3D graphical representation of the channel, its profile accompanied by a list of lining residues and also its basic physicochemical properties. The users can tune advanced parameters when performing a channel search to direct the search according to their needs. The MOLEonline 2.0 application is freely available via the Internet at http://ncbr.muni.cz/mole or http://mole.upol.cz.

• MeSH
• cytochrom P-450 CYP3A chemie   MeSH
• enzymy chemie   MeSH
• internet MeSH
• iontové kanály chemie   MeSH
• konformace proteinů MeSH
• molekulární modely MeSH
• počítačová grafika MeSH
• ribozomy chemie   MeSH
• software * MeSH
• uživatelské rozhraní počítače MeSH
• velké podjednotky ribozomu archebakteriální chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• cytochrom P-450 CYP3A MeSH
• enzymy MeSH
• iontové kanály MeSH

Atomistic molecular dynamics (MD) simulations of druglike molecules embedded in lipid bilayers are of considerable interest as models for drug penetration and positioning in biological membranes. Here we analyze partitioning of coumarin in dioleoylphosphatidylcholine (DOPC) bilayer, based on both multiple, unbiased 3 μs MD simulations (total length) and free energy profiles along the bilayer normal calculated by biased MD simulations (∼7 μs in total). The convergences in time of free energy profiles calculated by both umbrella sampling and z-constraint techniques are thoroughly analyzed. Two sets of starting structures are also considered, one from unbiased MD simulation and the other from "pulling" coumarin along the bilayer normal. The structures obtained by pulling simulation contain water defects on the lipid bilayer surface, while those acquired from unbiased simulation have no membrane defects. The free energy profiles converge more rapidly when starting frames from unbiased simulations are used. In addition, z-constraint simulation leads to more rapid convergence than umbrella sampling, due to quicker relaxation of membrane defects. Furthermore, we show that the choice of RESP, PRODRG, or Mulliken charges considerably affects the resulting free energy profile of our model drug along the bilayer normal. We recommend using z-constraint biased MD simulations based on starting geometries acquired from unbiased MD simulations for efficient calculation of convergent free energy profiles of druglike molecules along bilayer normals. The calculation of free energy profile should start with an unbiased simulation, though the polar molecules might need a slow pulling afterward. Results obtained with the recommended simulation protocol agree well with available experimental data for two coumarin derivatives.

• Publikační typ
• časopisecké články MeSH