Exploring coumarin egress channels in human cytochrome p450 2a6 by random acceleration and steered molecular dynamics simulations

Cytochrome P450 monooxygenase systems: Diversity and plasticity for adaptive stress response

Superfamily of cytochromes P450 (CYPs) is composed of heme-thiolate-containing monooxygenase enzymes, which play crucial roles in the biosynthesis, bioactivation, and detoxification of a variety of organic compounds, both endogenic and exogenic. Majority of CYP monooxygenase systems are multi-component and contain various redox partners, cofactors and auxiliary proteins, which contribute to their diversity in both prokaryotes and eukaryotes. Recent progress in bioinformatics and computational biology approaches make it possible to undertake whole-genome and phylogenetic analyses of CYPomes of a variety of organisms. Considerable variations in sequences within and between CYP families and high similarity in secondary and tertiary structures between all CYPs along with dramatic conformational changes in secondary structure elements of a substrate binding site during catalysis have been reported. This provides structural plasticity and substrate promiscuity, which underlie functional diversity of CYPs. Gene duplication and mutation events underlie CYP evolutionary diversity and emergence of novel selectable functions, which provide the involvement of CYPs in high adaptability to changing environmental conditions and dietary restrictions. In our review, we discuss the recent advancements and challenges in the elucidating the evolutionary origin and mechanisms underlying the CYP monooxygenase system diversity and plasticity. Our review is in the view of hypothesis that diversity of CYP monooxygenase systems is translated into the broad metabolic profiles, and this has been acquired during the long evolutionary time to provide structural plasticity leading to high adaptative capabilities to environmental stress conditions.

Structure-based identification of potential substrate antagonists for isethionate sulfite-lyase enzyme of Bilophila Wadsworthia: Towards novel therapeutic intervention to curb gut-associated illness

Bilophila wadsworthia is one of the prominent sources of hydrogen sulfide (H₂S) production in appendices, excessive levels of which can result in a weaker colonic mucus barrier, inflammatory bowel disease, and colorectal cancer. Isethionate sulfite-lyase (IslA) enzyme catalyzes H2S production by cleaving CS bond in isethionate, producing acetaldehyde and sulfite. In this study, we aimed to identify potential substrate antagonists for IsIA using a structure-based drug design. Initially, pharmacophore-based computational screening of the ZINC20 database yielded 66 hits that were subjected to molecular docking targeting the isethionate binding site of IsIA. Based on striking docking scores, nine compounds showed strong interaction with critical IsIA residues (Arg189, Gln193, Glu470, Cys468, and Arg678), drug-like features, appropriate adsorption, metabolism, excretion, and excretion profile with non-toxicity. Molecular dynamics simulations uncovered the significant impact of binding the compounds on protein conformational dynamics. Finally, binding free energies revealed substantial binding affinity (ranging from -35.23 to -53.88 kcal/mol) of compounds (ZINC913876497, ZINC913856647, ZINC914263733, ZINC914137795, ZINC915757996, ZINC914357083, ZINC913934833, ZINC9143362047, and ZINC913854740) for IsIA. The compounds proposed herein through a multi-faceted computational strategy can be experimentally validated as potential substrate antagonists of B. wadsworthia's IsIA for developing new medications to curb gut-associated illness in the future.

Key regulators in the architecture of substrate access/egress channels in mammalian cytochromes P450 governing flexibility in substrate oxyfunctionalization

Cytochrome P450s (CYP) represent a superfamily of b-type hemoproteins catalyzing oxifunctionalization of a vast array of endogenous and exogenous compounds. The present review focuses on assessment of the topology of prospective determinants in substrate entry and product release channels of mammalian P450s, steering the conformational dynamics of substrate accessibility and productive ligand orientation toward the iron-oxene core. Based on a generalized, CYP3A4-related construct, the sum of critical elements from diverse target enzymes was found to cluster within the known substrate recognition sites. The majority of prevalent substrate access/egress tunnels revealed to be of fairly balanced functional importance. The hydrophobicity profile of the candidates revealed to be the most salient feature in functional interaction throughout the conduits, while bulkiness of the residues imposes steric restrictions on substrate traveling. Thus, small amino acids such as prolines and glycines serve as hinges, driving conformational flexibility in ligand passage. Similarly, bottlenecks in the tunnel architecture, being narrowest encounter points within the CYP3A4 model, have a vital function in substrate selectivity along with clusters of aromatic amino acids acting as gatekeepers. In addition, peripheral patches in conduits may house determinants modulating allosteric cooperativity between remote and central domains in the P450 structure. Remarkably, the bulk critical residues lining tunnels in the various isozymes reside in helices B'/C and F/G inclusive of their interhelical turns as well as in helix I. This suggests these regions to represent hotspots for targeted genetic engineering to tailor more sophisticated mammalian P450s exploitable in industrial, biotechnological and medicinal areas.

Determination of tyrosinase-cyanidin-3-O-glucoside and (-/+)-catechin binding modes reveal mechanistic differences in tyrosinase inhibition

Tyrosinase, exquisitely catalyzes the phenolic compounds into brown or black pigment, inhibition is used as a treatment for dermatological or neurodegenerative disorders. Natural products, such as cyanidin-3-O-glucoside and (-/+)-catechin, are considered safe and non-toxic food additives in tyrosinase inhibition but their ambiguous inhibitory mechanism against tyrosinase is still elusive. Thus, we presented the mechanistic insights into tyrosinase with cyanidin-3-O-glucoside and (-/+)-catechin using computational simulations and in vitro assessment. Initial molecular docking results predicted ideal docked poses (- 9.346 to - 5.795 kcal/mol) for tyrosinase with selected flavonoids. Furthermore, 100 ns molecular dynamics simulations and post-simulation analysis of docked poses established their stability and oxidation of flavonoids as substrate by tyrosinase. Particularly, metal chelation via catechol group linked with the free 3-OH group on the unconjugated dihydropyran heterocycle chain was elucidated to contribute to tyrosinase inhibition by (-/+)-catechin against cyanidin-3-O-glucoside. Also, predicted binding free energy using molecular mechanics/generalized Born surface area for each docked pose was consistent with in vitro enzyme inhibition for both mushroom and murine tyrosinases. Conclusively, (-/+)-catechin was observed for substantial tyrosinase inhibition and advocated for further investigation for drug development against tyrosinase-associated diseases.

Rational engineering of cofactor specificity of glutamate dehydrogenase for poly-γ-glutamic acid synthesis in Bacillus licheniformis

Poly-γ-glutamic acid (γ-PGA) is a multifunctional biopolymer mainly produced by Bacillus. The cofactor specificity of enzymes plays a critical role in regulating metabolic process and metabolite production. Here, we present a novel approach for switching cofactor specificity of glutamate dehydrogenase RocG from nicotinamide adenine dinucleotide phosphate (NADPH) to nicotinamide adenine dinucleotide (NADH) to improve γ-PGA production. Firstly, 3D structural modeling and molecular docking were performed to predict the binding modes of NADH and NADPH. Several site-specific mutants based on the conventional and Random Accelerated Molecular Dynamics simulations were obtained to alter cofactor specificity. Then, the effects of RocG variants overexpressions on γ-PGA production were evaluated. Compared to the wild-type, the mutant RocG^D276E showed highest increase in γ-PGA yield, increased by 40.50%. Meanwhile, yields of main by-products acetoin and 2,3-butandieol were decreased by 21.70% and 16.53%, respectively. Finally, the results of enzymatic properties confirmed that glutamate dehydrogenase mutant RocG^D276E exhibited the higher affinity for NADH, caused a shift in coenzyme preference from NADPH to NADH, with a catalytic efficiency comparable with NADPH-dependent RocG. Taken together, this research demonstrated that switching the cofactor preference of glutamate dehydrogenase via rational design was an effective strategy for high-level production of γ-PGA in Bacillus licheniformis.

Discovery of Novel Inhibitors From Medicinal Plants for V-Domain Ig Suppressor of T-Cell Activation

V-domain Ig suppressor of T cell activation (VISTA) is an immune checkpoint and is a type I transmembrane protein. VISTA is linked to immunotherapy resistance, and it is a potential immune therapeutic target, especially for triple-negative breast cancer. It expresses at a high concentration in regulatory T cells and myeloid-derived suppressor cells, and its functional blockade is found to delay tumor growth. A useful medicinal plant database for drug designing (MPD3), which is a collection of phytochemicals from diverse plant families, was employed in virtual screening against VISTA to prioritize natural inhibitors against VISTA. Three compounds, Paratocarpin K (PubChem ID: 14187087), 3-(1H-Indol-3-yl)-2-(trimethylazaniumyl)propanoate (PubChem ID: 3861164), and 2-[(5-Benzyl-4-ethyl-1,2,4-triazol-3-yl)sulfanylmethyl]-5-methyl-1,3,4-oxadiazole (PubChem ID: 6494266), having binding energies stronger than -6 kcal/mol were found to have two common hydrogen bond interactions with VISTA active site residues: Arg54 and Arg127. The dynamics of the compound-VISTA complexes were further explored to infer binding stability of the systems. Results revealed that the compound 14187087 and 6494266 systems are highly stable with an average RMSD of 1.31 Å. Further affirmation on the results was achieved by running MM-GBSA on the MD simulation trajectories, which re-ranked 14187087 as the top-binder with a net binding energy value of -33.33 kcal/mol. In conclusion, the present study successfully predicted natural compounds that have the potential to block the function of VISTA and therefore can be utilized further in experimental studies to validate their real anti-VISTA activity.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Phytocompounds from Himalayan Medicinal Plants as Potential Drugs to Treat Multidrug-Resistant Salmonella typhimurium: An In Silico Approach

Medicinal plants can be used as natural therapeutics to treat diseases in humans. Enteric bacteria possess efflux pumps to remove bile salts from cells to avoid potential membrane damage. Resistance to bile and antibiotics is associated with the survival of Salmonella enterica subspecies enterica serovar Typhimurium (S. typhimurium) within a host. The present study aimed to investigate the binding affinity of major phytocompounds derived from 35 medicinal plants of the North Western Himalayas with the RamR protein (PDB ID 6IE9) of S. typhimurium. Proteins and ligands were prepared using AutoDock software 1.5.6. Molecular docking was performed using AutoDock Vina and MD simulation was performed at 100 ns. Drug likeness and toxicity predictions of hit phytocompounds were evaluated using molinspiration and ProTox II online servers. Moreover, docking, drug likeness, and toxicity results revealed that among all the selected phytocompounds, beta-sitosterol exhibited the most efficacious binding affinity with RamR protein (PDB ID 6IE9) and was nontoxic in nature. MD simulation data revealed that beta-sitosterol in complex with 6IE9 can be used as an antimicrobial. Furthermore, beta-sitosterol is stable in the binding pocket of the target protein; hence, it can be further explored as a drug to inhibit resistance-nodulation-division efflux pumps.

Drug repurposing for ligand-induced rearrangement of Sirt2 active site-based inhibitors via molecular modeling and quantum mechanics calculations

Sirtuin 2 (Sirt2) nicotinamide adenine dinucleotide-dependent deacetylase enzyme has been reported to alter diverse biological functions in the cells and onset of diseases, including cancer, aging, and neurodegenerative diseases, which implicate the regulation of Sirt2 function as a potential drug target. Available Sirt2 inhibitors or modulators exhibit insufficient specificity and potency, and even partially contradictory Sirt2 effects were described for the available inhibitors. Herein, we applied computational screening and evaluation of FDA-approved drugs for highly selective modulation of Sirt2 activity via a unique inhibitory mechanism as reported earlier for SirReal2 inhibitor. Application of stringent molecular docking results in the identification of 48 FDA-approved drugs as selective putative inhibitors of Sirt2, but only top 10 drugs with docking scores > - 11 kcal/mol were considered in reference to SirReal2 inhibitor for computational analysis. The molecular dynamics simulations and post-simulation analysis of Sirt2-drug complexes revealed substantial stability for Fluphenazine and Nintedanib with Sirt2. Additionally, developed 3D-QSAR-models also support the inhibitory potential of drugs, which exclusively revealed highest activities for Nintedanib (pIC50 ≥ 5.90 µM). Conclusively, screened FDA-approved drugs were advocated as promising agents for Sirt2 inhibition and required in vitro investigation for Sirt2 targeted drug development.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Investigating the binding mechanism of piperidinyl ureas inhibitors based on the UBC12-DCN1 interaction by 3D-QSAR, molecular docking and molecular dynamics simulations

Neddylation regulates a variety of biological processes by modulating Cullin-RING E3 ubiquitin ligases (CRLs) which is considered to be an important target for human diseases. The activation of CRLs required Cullins Neddylation, which regulated by the interaction of UBC12-DCN1 complex. Here, to investigate the structure-activity relationship and binding mechanism of 41 piperidinyl ureas inhibitors based on the UBC12-DCN1 protein-protein interaction, we carried out molecular modeling studies using three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations.Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to generate 3D-QSAR models. The results show that the best CoMFA model has q²=0.736, r²_ncv=0.978, r²_pred=0.78 (CoMFA), and the best CoMSIA model has q²=0.761, r²_ncv=0.987, r²_pred=0.86. The electrostatic, hydrophobic and H-bond donor fields play important roles in the models. Molecular docking studies predict the binding mode and the interactions between the ligand and the receptor protein. Molecular dynamics simulations results reveal that the complex of the ligand and the receptor protein are stable at 300 K. The results of MM-GBSA indicated the residues of Ile1083, Ile1086, Ala1098, Val1102, Ile1105, Gln1114, Phe1164 and Leu1184 might be the key residues during the process of inhibitors bound to DCN1. This study could provide an important theoretical basis for further developing novel inhibitors design based on UBC12-DCN1 protein-protein interaction. All the results can provide us more useful information for our further drug design. Communicated by Ramaswamy H. Sarma.

Human CYP2E1-dependent mutagenicity of benzene and its hydroxylated metabolites in V79-derived cells: Suppression and enhancement by ethanol pretreatment

Benzene is a human carcinogen that requires metabolic activation. We previously observed that benzene and its hydroxylated metabolites induce micronuclei in mammalian cells expressing human CYP2E1. This study was initially aimed to study another endpoint, the induction of gene mutations by those compounds in the same cell models. A V79-derived cell line expressing human CYP2E1 and sulfotransferase (SULT) 1A1 (V79-hCYP2E1-hSULT1A1) pretreated with ethanol (a CYP2E1 stabilizer) was used in the hprt gene mutagenicity assay. Phenol, hydroquinone, catechol, and 1,2,4-trihydroxybenzene all induced gene mutations, while they were inactive, or only weakly positive (hydroquinone), in parental V79-Mz cells. Unexpectedly, benzene was non-mutagenic in both cell lines, but it became positive in V79-hCYP2E1-hSULT1A1 cells using regimes of short exposure/long recovery without ethanol pretreatment, for both gene mutations and micronuclei formation. In silico molecular simulation showed binding energies and positions favorable for each compound to be oxidized by human CYP2E1, benzene demonstrating the highest affinity. By tunnel analysis, ethanol binding did not limit benzene to pass tunnel S, which was specifically active for benzene. However, its end product, acetic acid, decreased the occurrence of tunnel S from 5.4 to 2.2% and extended the length of its bottleneck from 5.5 to 9.0 Å. With residual ethanol molecules still being present in CYP2E1 for a period of time after benzene exposure, the acetic acid formed could limit the entrance of benzene, thus inhibit its metabolic activation. In summary, ethanol may interfere with the activation of benzene to mutagenic metabolites, at least in cultured cells.

Docking studies and molecular dynamics simulations of the binding characteristics of waldiomycin and its methyl ester analog to Staphylococcus aureus histidine kinase

Bacterial histidine kinases (HKs) are considered attractive drug targets because of their ability to govern adaptive responses coupled with their ubiquity. There are several classes of HK inhibitors; however, they suffer from drug resistance, poor bioavailability, and a lack of selectivity. The 3D structure of Staphylococcus aureus HK was not isolated in high-resolution coordinates, precluding further disclosure of structure-dependent binding to the specific antibiotics. To elucidate structure-dependent binding, the 3D structure of the catalytic domain WalK of S. aureus HK was constructed using homology modeling to investigate the WalK-ligand binding mechanisms through molecular docking studies and molecular dynamics simulations. The binding free energies of the waldiomycin and its methyl ester analog were calculated using molecular mechanics/generalized born surface area scoring. The key residues for protein-ligand binding were postulated. The structural divergence responsible for the 7.4-fold higher potency of waldiomycin than that of its ester analog was clearly observed. The optimized 3D macromolecule-ligand binding modes shed light on the S. aureus HK/WalK-ligand interactions that afford a means to assess binding affinity to design new HK/WalK inhibitors.

Repurposing of Fluvastatin Against Candida albicans CYP450 Lanosterol 14 α-demethylase, a Target Enzyme for Antifungal Therapy: An In silico and In vitro Study

Background:

The incidence of fungal infections has increased significantly. Specifically the cases of candida albicans infection are increasing day by day and their resistance to clinically approved drugs is a major concern for humans. Various classes of antifungal drugs are available in the market for the treatment of these infections but unfortunately, none of them is able to treat the infection.

Objective:

Thus, in the present investigation, we have repurposed the well-known drug (Fluvastatin) in the treatment of Candida albicans infections by using in silico, in vitro and ex vivo techniques.

Results:

Firstly, we developed and validated a simple model of CYP45014α-lanosterol demethylase of Candida albicans by using crystal structure of Mycobacterium tuberculosis (1EA1). Further, fluvastatin was docked with a validated model of CYP45014α-lanosterol demethylase and revealed good binding affinity as that of fluconazole. In vitro results (Percentage growth retardation, Fungal growth kinetics, Biofilm test and Post antifungal test) have shown good antifungal activity of fluvastatin. Finally, the results of MTT assay have shown non-cytotoxic effect of fluvastatin in murine splenocytes and thymocytes.

Results:

Firstly, we developed and validated a simple model of CYP45014α-lanosterol demethylase of Candida albicans by using crystal structure of Mycobacterium tuberculosis (1EA1). Further, fluvastatin was docked with a validated model of CYP45014α-lanosterol demethylase and revealed good binding affinity as that of fluconazole. In vitro results (Percentage growth retardation, Fungal growth kinetics, Biofilm test and Post antifungal test) have shown good antifungal activity of fluvastatin. Finally, the results of MTT assay have shown non-cytotoxic effect of fluvastatin in murine splenocytes and thymocytes.

Conclusion:

However, further in vivo studies are required to confirm the complete role of fluvastatin as an antifungal agent.

Prediction of Drug-Drug Interactions Related to Inhibition or Induction of Drug-Metabolizing Enzymes

Drug-drug interaction (DDI) is the phenomenon of alteration of the pharmacological activity of a drug(s) when another drug(s) is co-administered in cases of so-called polypharmacy. There are three types of DDIs: pharmacokinetic (PK), pharmacodynamic, and pharmaceutical. PK is the most frequent type of DDI, which often appears as a result of the inhibition or induction of drug-metabolising enzymes (DME). In this review, we summarise in silico methods that may be applied for the prediction of the inhibition or induction of DMEs and describe appropriate computational methods for DDI prediction, showing the current situation and perspectives of these approaches in medicinal and pharmaceutical chemistry. We review sources of information on DDI, which can be used in pharmaceutical investigations and medicinal practice and/or for the creation of computational models. The problem of the inaccuracy and redundancy of these data are discussed. We provide information on the state-of-the-art physiologically- based pharmacokinetic modelling (PBPK) approaches and DME-based in silico methods. In the section on ligand-based methods, we describe pharmacophore models, molecular field analysis, quantitative structure-activity relationships (QSAR), and similarity analysis applied to the prediction of DDI related to the inhibition or induction of DME. In conclusion, we discuss the problems of DDI severity assessment, mention factors that influence severity, and highlight the issues, perspectives and practical using of in silico methods.

Four Major Channels Detected in the Cytochrome P450 3A4: A Step toward Understanding Its Multispecificity

We computed the network of channels of the 3A4 isoform of the cytochrome P450 (CYP) on the basis of 16 crystal structures extracted from the Protein Data Bank (PDB). The calculations were performed with version 2 of the CCCPP software that we developed for this research project. We identified the minimal cost paths (MCPs) output by CCCPP as probable ways to access to the buried active site. The algorithm of calculation of the MCPs is presented in this paper, with its original method of visualization of the channels. We found that these MCPs constitute four major channels in CYP3A4. Among the many channels proposed by Cojocaru et al. in 2007, we found that only four of them open in 3A4. We provide a refined description of these channels together with associated quantitative data.

Identification and characterization of a sterically robust phenylalanine ammonia-lyase among 481 natural isoforms through association of in silico and in vitro studies

The enzyme phenylalanine ammonia lyase (PAL) is of special importance for the treatment of phenylketonuria patients. The aim of this study was to find a stable recombinant PAL with suitable kinetic properties among all natural PAL producing species using in silico and experimental approaches. To find such a stable PAL among 481 natural isoforms, 48,000 of 3-D models were predicted using the Modeller 9.10 program and evaluated by Ramachandran plot. Correlation analysis between Ramachandran plot and the energy of different thermodynamic components indicated that this plot could be an appropriate tool to predict protein stability. Hence, PAL6 from Lotus japonicus (LjPAL6) was selected as a stable isoform. Molecular dynamic (MD) simulation for 50 ns and docking has been conducted for LjPAL6-phenylalanine complex. The best PAL-phenylalanine frame was selected by re-docking with l-phenylalanine (L-Phe) and root-mean-square deviation (RMSD) value. MD simulation showed that the complex has a good stability, depicted by the low RMSD value, binding free energy and hydrogen bindings. Docking results showed that LjPAL6 has a high affinity toward l-Phe according to the low level of binding free energy. By overexpressing Ljpal6 in E. coli BL21, a total of 33.5 mg/l of protein was obtained, which has been increased to 83.7 mg/l via the optimization of LjPAL6 production using response surface methodology. The optimal pH and temperature were 8.5 and 50 °C, respectively. LjPAL6 showed a specific activity of 42 nkat/mg protein, with K_m, K_cat and K_cat/K_m values of 0.483 mM, 7 S^-1 and 14.5 S^-1 mM^-1 for l-phe, respectively. In conclusion, finding models with the most reasonable stereo-chemical quality and lowest numbers of steric clashes would result in easier folding. Hence, in silico analyses of bulk data from natural origin will lead one to find an optimal model for in vitro studies and drug design.

Probing the binding mode and unbinding mechanism of LSD1 inhibitors by combined computational methods

Lysine specific demethylase 1 (LSD1) has emerged as a potential drug target in cancer therapy and a variety of inhibitors have been reported. We have recently reported the discovery of a series of triazole-dithiocarbamate based compounds, which were basically confirmed as cofactor flavin adenine dinucleotide (FAD)-competing inhibitors by experiments. However, the binding modes of the inhibitors to the binding site were undetermined. Here, we employed computational methods including molecular docking, classical molecular dynamics (MD) and steered molecular dynamics (SMD) simulations to investigate the potential binding modes of these inhibitors to LSD1. Based on the high correlation between the mean non-equilibrium pulling work W and experimental binding affinity, we identified the optimal binding modes of this class of compounds with LSD1. Using the optimal inhibitor binding conformation, we then performed SMD to study the ligand unbinding mechanism with a lower pulling velocity at 0.0005 nm ps-1. We found that residue Arg316 plays a crucial role in the binding/unbinding process. Furthermore, a gatekeeper residue Trp756 influences the ligand unbinding process by acting like a switch via steric hindrance but can enhance the hydrophobic interaction with the inhibitor. Hydrophobic interaction also dominated the interaction between LSD1 and the inhibitors. The pivotal residues and interactions between LSD1 and inhibitors determined from this study can be used to improve the inhibition activity of this series of inhibitors in development and to discover new scaffolds as FAD-competing inhibitors in compound screening.

The GluN2B-Glu413Gly NMDA receptor variant arising from a de novo GRIN2B mutation promotes ligand-unbinding and domain opening

N-methyl-D-aspartate (NMDA) receptors are transmembrane glutamate-binding ion channels that mediate neurotransmission in mammals. NMDA receptor subunits are tetrameric complexes of GluN1 and GluN2A-D subunits, encoded by the GRIN gene family. Of these subunits, GluN2B is suggested to be required for normal development of the central nervous system. A mutation identified in a patient with developmental delay, E413G, resides in the GluN2B ligand-binding domain and substantially reduces glutamate potency by an unknown mechanism. GluN2B Gly413, though near the agonist, is not in van der Waals contact with glutamate. Visual analysis of the GluN2B structure with the E413G mutation modeled suggests that replacement of Glu with Gly at this position increases solvent access to the ligand-binding domain. This was confirmed by molecular modeling, which showed that the ligand is more mobile in GluN2B-E413G than WT GluN2B. Evaluation of agonist occupancy using random accelerated molecular dynamics (RAMD) simulations predicts that the glutamate exits the binding-site more rapidly for GluN2B-E413G than WT receptors. This analysis was extended to other binding-site mutations, which produced qualitative agreement between experimentally determined EC50 values, deactivation time constants, and ligand motion within the binding-site. Furthermore, long sub-microsecond molecular dynamics simulations of the bi-lobed ligand-binding domain revealed that it adopted a cleft-open ligand-free state more often for GluN2B-E413G than wild-type GluN2B. This is consistent with the idea that L-glutamate binding is altered such that the ligand-binding domain occupies the open-cleft conformation associated with the closed channel.

Insights into product release dynamics through structural analyses of thymidylate kinase

Several studies on enzyme catalysis have pointed out that the product release event could be a rate limiting step. In this study, we have compared the release event of two products, Adenosine di-phosphate (ADP) and Thymidine di-phosphate (TDP) from the active-site of human and Thermus thermophilus thymidine mono-phosphate kinase (TMPK), referred to as hTMPK and ttTMPK, respectively. TMPK catalyses the conversion of Thymidine mono-phosphate (TMP) to TDP using ATP as phosphoryl donor in the presence of Mg²⁺ ion. Most of the earlier studies on this enzyme have focused on understanding substrate binding and catalysis, but the critical product release event remains elusive. Competitive binding experiments of the substrates and the products using ttTMPK apo crystals have indicated that the substrate (TMP) can replace the bound product (TDP), even in the presence of an ADP molecule. Further, the existing random accelerated molecular dynamics (RAMD) simulation program was modified to study the release of both the products simultaneously from the active site. The RAMD simulations on product-bound structures of both ttTMPK and hTMPK, revealed that while several exit patterns of the products are permissible, the sequential exit mode is the most preferred pattern for both ttTMPK and hTMPK enzymes. Additionally, the product release from the hTMPK was found to be faster and more directional as compared to ttTMPK. Structural investigation revealed that the critical changes in the residue composition in the LID-region of ttTMPK and hTMPK have an effect on the product release and can be attributed to the observed differences during product release event. Understanding of these dissimilarities is of considerable utility in designing potent inhibitors or prodrugs that can distinguish between eukaryotic and prokaryotic homologues of thymidylate kinase.

Molecular Dynamics Simulations Reveal Structural Differences among Allelic Variants of Membrane-Anchored Cytochrome P450 2D6

Cytochrome P450 2D6 (CYP2D6) is an enzyme that is involved in the metabolism of roughly 25% of all marketed drugs and therefore belongs to the most important enzymes in drug metabolism. CYP2D6 features a high degree of genetic polymorphism that can significantly affect the metabolic activity of an individual. In extreme cases, structural changes at the level of single amino acids can either increase its enzymatic activity abolishing the drug therapeutic effect or completely disable the enzyme and elevate drug plasma level potentially leading to adverse effects. In this study, starting from the crystal structure, we built a full-length membrane-anchored all-atom model of the wild-type CYP2D6 as well as five of its variants differing in the enzymatic activity. We validated our models with available experimental data and compared their structural properties with molecular dynamics simulations. The main focus of this study was to identify differences that could mechanistically explain the altered activity of the variants and improve our understanding of their functioning. We observed differences in the opening frequencies and minimal diameters of tunnels that connect the buried active site to the surrounding solvent environment. The variants CYP2D6*4 and CYP2D6*10 associated with missing or decreased activity showed less frequent opening of the tunnels compared to the wild-type. Both CYP2D6*10 and CYP2D6*17 showed a deprivation of an important ligand tunnel suggesting a feasible reason for their altered substrate specificity. Next, the altered fold at the N-terminal anchor region and the decreased active site volume caused by the amino acid mutations of the CYP2D6*4 variant offer an explanation for the absence of its metabolic activity. The mutations in CYP2D6*53 contributed to a significant enlargement of an important ligand tunnel and an extension of the active site cavity. This could explain the altered metabolic profile as well as the enhanced metabolic rates of this particular variant supporting its designation as a possible cause for the ultrarapid metabolizer phenotype. We believe these novel structural insights could advance the fields of personalized medicine and enzyme engineering. Furthermore, they could aid in guiding laboratory as well as computational experiments in the future.

Lipid molecules can induce an opening of membrane-facing tunnels in cytochrome P450 1A2

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.

ChannelsDB: database of biomacromolecular tunnels and pores

ChannelsDB (http://ncbr.muni.cz/ChannelsDB) is a database providing information about the positions, geometry and physicochemical properties of channels (pores and tunnels) found within biomacromolecular structures deposited in the Protein Data Bank. Channels were deposited from two sources; from literature using manual deposition and from a software tool automatically detecting tunnels leading to the enzymatic active sites and selected cofactors, and transmembrane pores. The database stores information about geometrical features (e.g. length and radius profile along a channel) and physicochemical properties involving polarity, hydrophobicity, hydropathy, charge and mutability. The stored data are interlinked with available UniProt annotation data mapping known mutation effects to channel-lining residues. All structures with channels are displayed in a clear interactive manner, further facilitating data manipulation and interpretation. As such, ChannelsDB provides an invaluable resource for research related to deciphering the biological function of biomacromolecular channels.

Thermodynamics of camphor migration in cytochrome P450cam by atomistic simulations

Understanding the mechanisms of ligand binding to enzymes is of paramount importance for the design of new drugs. Here, we report on the use of a novel biased molecular dynamics (MD) methodology to study the mechanism of camphor binding to cytochrome P450cam. Microsecond-long MD simulations allowed us to observe reaction coordinates characterizing ligand diffusion from the active site of cytochrome P450cam to solvent via three egress routes. These atomistic simulations were used to estimate thermodynamic quantities along the reaction coordinates and indicate diverse binding configurations. The results suggest that the diffusion of camphor along the pathway near the substrate recognition site (SRS) is thermodynamically preferred. In addition, we show that the diffusion near the SRS is triggered by a transition from a heterogeneous collection of closed ligand-bound conformers to the basin comprising the open conformations of cytochrome P450cam. The conformational change accompanying this switch is characterized by the retraction of the F and G helices and the disorder of the B' helix. These results are corroborated by experimental studies and provide detailed insight into ligand binding and conformational behavior of the cytochrome family. The presented methodology is general and can be applied to other ligand-protein systems.

Ligand diffusion in proteins via enhanced sampling in molecular dynamics

Computational simulations in biophysics describe the dynamics and functions of biological macromolecules at the atomic level. Among motions particularly important for life are the transport processes in heterogeneous media. The process of ligand diffusion inside proteins is an example of a complex rare event that can be modeled using molecular dynamics simulations. The study of physical interactions between a ligand and its biological target is of paramount importance for the design of novel drugs and enzymes. Unfortunately, the process of ligand diffusion is difficult to study experimentally. The need for identifying the ligand egress pathways and understanding how ligands migrate through protein tunnels has spurred the development of several methodological approaches to this problem. The complex topology of protein channels and the transient nature of the ligand passage pose difficulties in the modeling of the ligand entry/escape pathways by canonical molecular dynamics simulations. In this review, we report a methodology involving a reconstruction of the ligand diffusion reaction coordinates and the free-energy profiles along these reaction coordinates using enhanced sampling of conformational space. We illustrate the above methods on several ligand-protein systems, including cytochromes and G-protein-coupled receptors. The methods are general and may be adopted to other transport processes in living matter.

Molecular basis of P450 OleTJE: an investigation of substrate binding mechanism and major pathways

Cytochrome P450 OleT_JE has attracted much attention for its ability to catalyze the decarboxylation of long chain fatty acids to generate alkenes, which are not only biofuel molecule, but also can be used broadly for making lubricants, polymers and detergents. In this study, the molecular basis of the binding mechanism of P450 OleT_JE for arachidic acid, myristic acid, and caprylic acid was investigated by utilizing conventional molecular dynamics simulation and binding free energy calculations. Moreover, random acceleration molecular dynamics (RAMD) simulations were performed to uncover the most probable access/egress channels for different fatty acids. The predicted binding free energy shows an order of arachidic acid < myristic acid < caprylic acid. Key residues interacting with three substrates and residues specifically binding to one of them were identified. The RAMD results suggest the most likely channel for arachidic acid, myristic acid, and caprylic acid are 2e/2b, 2a and 2f/2a, respectively. It is suggested that the reaction is easier to carry out in myristic acid bound system than those in arachidic acid and caprylic acid bound system based on the distance of Hβ atom of substrate relative to P450 OleT_JE Compound I states. This study provided novel insight to understand the substrate preference mechanism of P450 OleT_JE and valuable information for rational enzyme design for short chain fatty acid decarboxylation.

Computational tools for enzyme improvement: why everyone can - and should - use them

This review presents computational methods that experimentalists can readily use to create smart libraries for enzyme engineering and to obtain insights into protein-substrate complexes. Computational tools have the reputation of being hard to use and inaccurate compared to experimental methods in enzyme engineering, yet they are essential to probe datasets of ever-increasing size and complexity. In recent years, bioinformatics groups have made a huge leap forward in providing user-friendly interfaces and accurate algorithms for experimentalists. These methods guide efficient experimental planning and allow the enzyme engineer to rationalize time and resources. Computational tools nevertheless face challenges in the realm of transient modern technology.

Structural-functional adaptations of porcine CYP1A1 to metabolize polychlorinated dibenzo-p-dioxins

Polychlorinated dibenzo-p-dioxins (PCDDs) are widespread by-products of human industrial activity. They accumulate in tissues of animals and humans, exerting numerous adverse effects on different systems. In living organisms, dioxins are metabolized by enzymes of the cytochrome P450 family, including CYP1A1. Particular dioxin congeners differ in their toxicity level and ability to undergo biodegradation. Since the molecular mechanisms underlying dioxin susceptibility or resistance to biodegradation are unknown, in the present study the molecular interactions between five selected dioxins and porcine CYP1A1 protein were investigated. It was found that the ability of a dioxin to undergo CYP1A1-mediated degradation is associated mainly with the number and position of chlorine atoms in the dioxin molecule. Among all examined congeners, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) demonstrated the highest affinity to CYP1A1 and, at the same time, the greatest distance to the active site of the enzyme. Interestingly, in contrast to other dioxins, the binding of the TCDD molecule to the porcine CYP1A1 active site resulted in a rapid and continuous closure of substrate channels. All the information may help to explain the extended half-life of TCDD in living organisms as well as its high toxicity.

Molecular Dynamics and Free Energy Simulations of Phenylacetate and CO2 Release from AMDase and Its G74C/C188S Mutant: A Possible Rationale for the Reduced Activity of the Latter

Arylmalonate decarboxylase (AMDase) catalyzes the decarboxylation of α-aryl-α-methyl malonates to produce optically pure α-arylpropionates of industrial and medicinal importance. Herein, atomistic molecular dynamics simulations have been carried out to delineate the mechanism of the release of product molecules phenylacetate (PAC) and carbon dioxide (CO₂), from the wild-type (WT) and its G74C/C188S mutant enzymes. Both of the product molecules follow a crystallographically characterized solvent-accessible channel to come out of the protein interior. A higher free energy barrier for the release of PAC from G74C/C188S compared to that in the WT is consistent with the experimentally observed compromised efficiency of the mutant. The release of CO₂ precedes that of PAC; free energy barriers for CO₂ and PAC release in the WT enzyme are calculated to be ∼1-2 and ∼23 kcal/mol, respectively. Postdecarboxylation, CO₂ moves toward a hydrophobic pocket formed by Pro 14, Leu 38, Leu 40, Leu 77, and the side chain of Tyr 48 which serves as its temporary "reservoir". CO₂ releases following a channel mainly decorated by apolar residues, unlike in the case of oxalate decarboxylase where polar residues mediate its transport.

Exploring ligand dissociation pathways from aminopeptidase N using random acceleration molecular dynamics simulation

Aminopeptidase N (APN) is a zinc-dependent ectopeptidase involved in cell proliferation, secretion, invasion, and angiogenesis, and is widely recognized as an important cancer target. However, the mechanisms whereby ligands leave the active site of APN remain unknown. Investigating ligand dissociation processes is quite difficult, both in classical simulation methods and in experimental approaches. In this study, random acceleration molecular dynamics (RAMD) simulation was used to investigate the potential dissociation pathways of ligand from APN. The results revealed three pathways (channels A, B and C) for ligand release. Channel A, which matches the hypothetical channel region, was the most preferred region for bestatin to dissociate from the enzyme, and is probably the major channel for the inner bound ligand. In addition, two alternative channels (channels B and C) were shown to be possible pathways for ligand egression. Meanwhile, we identified key residues controlling the dynamic features of APN channels. Identification of the dissociation routes will provide further mechanistic insights into APN, which will benefit the development of more promising APN inhibitors. Graphical Abstract The release pathways of bestatin inside active site of aminopeptidase N were simulated using RAMD simulation.

Homology model, molecular dynamics simulation and novel pyrazole analogs design of Candida albicans CYP450 lanosterol 14 α-demethylase, a target enzyme for antifungal therapy

Candida albicans infections and their resistance to clinically approved azole drugs are major concerns for human. The azole antifungal drugs inhibit the ergosterol synthesis by targeting lanosterol 14α-demethylase of cytochrome P450 family. The lack of high-resolution structural information of fungal pathogens has been a barrier for the design of modified azole drugs. Thus, a preliminary theoretical molecular dynamic study is carried out to develop and validate a simple homologous model using crystallographic structure of the lanosterol 14α-demethylase of Mycobacterium tuberculosis (PDB ID-1EA1) in which the active site residues are substituted with that of C. albicans (taxid 5476). Further, novel designed pyrazole analogs (SGS1-16) docked on chimeric 1EA1 and revealed that SGS-16 show good binding affinity through non-bonding interaction with the heme, which is different from the leading azole antifungals. The ADME-T results showed these analogs can be further explored in design of more safe and effective antifungal agents.

Memetic algorithms for ligand expulsion from protein cavities

Ligand diffusion through a protein interior is a fundamental process governing biological signaling and enzymatic catalysis. A complex topology of channels in proteins leads often to difficulties in modeling ligand escape pathways by classical molecular dynamics simulations. In this paper, two novel memetic methods for searching the exit paths and cavity space exploration are proposed: Memory Enhanced Random Acceleration (MERA) Molecular Dynamics (MD) and Immune Algorithm (IA). In MERA, a pheromone concept is introduced to optimize an expulsion force. In IA, hybrid learning protocols are exploited to predict ligand exit paths. They are tested on three protein channels with increasing complexity: M2 muscarinic G-protein-coupled receptor, enzyme nitrile hydratase, and heme-protein cytochrome P450cam. In these cases, the memetic methods outperform simulated annealing and random acceleration molecular dynamics. The proposed algorithms are general and appropriate in all problems where an accelerated transport of an object through a network of channels is studied.

Computational study on the unbinding pathways of B-RAF inhibitors and its implication for the difference of residence time: insight from random acceleration and steered molecular dynamics simulations

Random acceleration and steered molecular dynamics simulations reveal the unbinding pathway of B-RAF inhibitors and the difference in the residence time.

Mapping Mechanical Force Propagation through Biomolecular Complexes

Here we employ single-molecule force spectroscopy with an atomic force microscope (AFM) and steered molecular dynamics (SMD) simulations to reveal force propagation pathways through a mechanically ultrastable multidomain cellulosome protein complex. We demonstrate a new combination of network-based correlation analysis supported by AFM directional pulling experiments, which allowed us to visualize stiff paths through the protein complex along which force is transmitted. The results implicate specific force-propagation routes nonparallel to the pulling axis that are advantageous for achieving high dissociation forces.

In silico methods for predicting drug-drug interactions with cytochrome P-450s, transporters and beyond

Drug-drug interactions (DDIs) are associated with severe adverse effects that may lead to the patient requiring alternative therapeutics and could ultimately lead to drug withdrawal from the market if they are severe. To prevent the occurrence of DDI in the clinic, experimental systems to evaluate drug interaction have been integrated into the various stages of the drug discovery and development process. A large body of knowledge about DDI has also accumulated through these studies and pharmacovigillence systems. Much of this work to date has focused on the drug metabolizing enzymes such as cytochrome P-450s as well as drug transporters, ion channels and occasionally other proteins. This combined knowledge provides a foundation for a hypothesis-driven in silico approach, using either cheminformatics or physiologically based pharmacokinetics (PK) modeling methods to assess DDI potential. Here we review recent advances in these approaches with emphasis on hypothesis-driven mechanistic models for important protein targets involved in PK-based DDI. Recent efforts with other informatics approaches to detect DDI are highlighted. Besides DDI, we also briefly introduce drug interactions with other substances, such as Traditional Chinese Medicines to illustrate how in silico modeling can be useful in this domain. We also summarize valuable data sources and web-based tools that are available for DDI prediction. We finally explore the challenges we see faced by in silico approaches for predicting DDI and propose future directions to make these computational models more reliable, accurate, and publically accessible.

Accelerated molecular dynamics and protein conformational change: a theoretical and practical guide using a membrane embedded model neurotransmitter transporter

Molecular dynamics simulation provides a powerful and accurate method to model protein conformational change, yet timescale limitations often prevent direct assessment of the kinetic properties of interest. A large number of molecular dynamic steps are necessary for rare events to occur, which allow a system to overcome energy barriers and conformationally transition from one potential energy minimum to another. For many proteins, the energy landscape is further complicated by a multitude of potential energy wells, each separated by high free-energy barriers and each potentially representative of a functionally important protein conformation. To overcome these obstacles, accelerated molecular dynamics utilizes a robust bias potential function to simulate the transition between different potential energy minima. This straightforward approach more efficiently samples conformational space in comparison to classical molecular dynamics simulation, does not require advanced knowledge of the potential energy landscape and converges to the proper canonical distribution. Here, we review the theory behind accelerated molecular dynamics and discuss the approach in the context of modeling protein conformational change. As a practical example, we provide a detailed, step-by-step explanation of how to perform an accelerated molecular dynamics simulation using a model neurotransmitter transporter embedded in a lipid cell membrane. Changes in protein conformation of relevance to the substrate transport cycle are then examined using principle component analysis.

Computational insights into inhibitory mechanism of azole compounds against human aromatase

We investigated the inhibitory mechanism of azole aromatase inhibitors. The results showed that letrozole and imazalil prefer different unbinding pathways.

Flexible docking-based molecular dynamics/steered molecular dynamics calculations of protein-protein contacts in a complex of cytochrome P450 1A2 with cytochrome b5

Formation of transient complexes of cytochrome P450 (P450) with another protein of the endoplasmic reticulum membrane, cytochrome b5 (cyt b5), dictates the catalytic activities of several P450s. Therefore, we examined formation and binding modes of the complex of human P450 1A2 with cyt b5. Docking of soluble domains of these proteins was performed using an information-driven flexible docking approach implemented in HADDOCK. Stabilities of the five unique binding modes of the P450 1A2-cyt b5 complex yielded by HADDOCK were evaluated using explicit 10 ns molecular dynamics (MD) simulations in aqueous solution. Further, steered MD was used to compare the stability of the individual P450 1A2-cyt b5 binding modes. The best binding mode was characterized by a T-shaped mutual orientation of the porphyrin rings and a 10.7 Å distance between the two redox centers, thus satisfying the condition for a fast electron transfer. Mutagenesis studies and chemical cross-linking, which, in the absence of crystal structures, were previously used to deduce specific P450-cyt b5 interactions, indicated that the negatively charged convex surface of cyt b5 binds to the positively charged concave surface of P450. Our simulations further elaborate structural details of this interface, including nine ion pairs between R95, R100, R138, R362, K442, K455, and K465 side chains of P450 1A2 and E42, E43, E49, D65, D71, and heme propionates of cyt b5. The universal heme-centric system of internal coordinates was proposed to facilitate consistent classification of the orientation of the two porphyrins in any protein complex.

Prediction of amino acid residues participated in substrate recognition by cytochrome P450 subfamilies with broad substrate specificity

Cytochromes P450 comprise a large superfamily and several of their isoforms play a crucial role in metabolism of xenobiotics, including drugs. Although these enzymes demonstrate broad and cross-substrate specificity, different cytochrome P450 subfamilies exhibit certain selectivity for some types of substrates. Analysis of amino acid residues of the active sites of six cytochrome subfamilies (CYP1А, CYP2А, CYP2С, CYP2D, CYP2E and CYP3А) enables to define subfamily-specific patterns that consist of four residues. These residues are located on the periphery of the active sites of these cytochromes. We suggest that they can form a primary binding site at the entrance to the active site, defining cytochrome substrate recognition.

Homology modeling, molecular docking, and molecular dynamics simulations elucidated α-fetoprotein binding modes

Background

An important mechanism of endocrine activity is chemicals entering target cells via transport proteins and then interacting with hormone receptors such as the estrogen receptor (ER). α-Fetoprotein (AFP) is a major transport protein in rodent serum that can bind and sequester estrogens, thus preventing entry to the target cell and where they could otherwise induce ER-mediated endocrine activity. Recently, we reported rat AFP binding affinities for a large set of structurally diverse chemicals, including 53 binders and 72 non-binders. However, the lack of three-dimensional (3D) structures of rat AFP hinders further understanding of the structural dependence for binding. Therefore, a 3D structure of rat AFP was built using homology modeling in order to elucidate rat AFP-ligand binding modes through docking analyses and molecular dynamics (MD) simulations.

Methods

Homology modeling was first applied to build a 3D structure of rat AFP. Molecular docking and Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) scoring were then used to examine potential rat AFP ligand binding modes. MD simulations and free energy calculations were performed to refine models of binding modes.

Results

A rat AFP tertiary structure was first obtained using homology modeling and MD simulations. The rat AFP-ligand binding modes of 13 structurally diverse, representative binders were calculated using molecular docking, (MM-GBSA) ranking and MD simulations. The key residues for rat AFP-ligand binding were postulated through analyzing the binding modes.

Conclusion

The optimized 3D rat AFP structure and associated ligand binding modes shed light on rat AFP-ligand binding interactions that, in turn, provide a means to estimate binding affinity of unknown chemicals. Our results will assist in the evaluation of the endocrine disruption potential of chemicals.

Behavior of human cytochromes P450 on lipid membranes

Human cytochromes P450 (CYPs) are membrane-anchored enzymes involved in biotransformation of many marketed drugs. We constructed atomic models of six human CYPs (CYP1A2, 2A6, 2C9, 2D6, 2E1, and 3A4) anchored to a lipid bilayer to investigate the positions and orientations of CYPs on a membrane. We equilibrated the models by molecular dynamics simulations on a 100+ ns time scale. Catalytic domains of all studied CYPs were found to be partially immersed in the lipid bilayer, whereas the N-terminal part and F'/G' loop are deeply immersed. The proximal side of the enzyme faces the cytosol, whereas the distal side, where openings of substrate access and product release channels to the active site are primarily located, points toward the lipid bilayer. Access channels with openings in the vicinity of the B/C and F/G loops are typically positioned below the lipid head groups, whereas the solvent channel points toward the membrane-water interface. We found that the access channel opening positions match the preferred substrate positions, whereas the product release channel exit positions correspond closely with the positions of the products. This may indicate that membrane-anchored CYPs have evolutionarily adapted to facilitate uptake of nonpolar substrates from the membrane and uptake/release of polar substrates or products from/to the membrane-water interface.