PubMed Health⌕ Search

SEARCH · PubMed Health

Results for “Molecular Docking Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78Linked to original sources

Structure and mechanism of arylamine N-acetyltransferases.

Arylamine N-acetyltransferases (NATs) are a family of phase II drug-metabolising enzymes which are important in the biotransformation of various aromatic and heterocyclic amines and hydroxylamines, arylhydrazines and arylhydrazides. NATs are present in a wide range of eukaryotes and prokaryotes. Humans have two functional NAT isoforms, both of which are highly polymorphic. The pharmacogenetics of NATs is an area which has been extensively studied. The determination of the X-ray crystal structure of NAT from Salmonella typhimurium led to the identification of the catalytically essential triad of residues: Cys-His-Asp, which is present in all functional NAT enzymes. Recent co-crystallisation data and in silico docking studies of NAT from Mycobacterium smegmatis with substrates and inhibitors have aided the identification of important contact residues within the active site. The X-ray crystal structures of four prokaryotic NAT proteins have now been determined, and these have been used to generate structural models of eukaryotic NATs, providing valuable insight into their active-site architecture. In addition to aiding crystallographic experiments, recent progress in the production of recombinant prokaryotic and eukaryotic NATs has allowed comparative studies of the kinetics and activity profiles of these enzymes. In this review we present an overview of recent structural and activity studies on NAT enzymes, and we outline how in silico methods may be used to predict NAT protein-ligand interactions based on the current knowledge.

Animals↗

In silico panning for a non-competitive peptide inhibitor.

BACKGROUND: Peptide ligands have tremendous therapeutic potential as efficacious drugs. Currently, more than 40 peptides are available in the market for a drug. However, since costly and time-consuming synthesis procedures represent a problem for high-throughput screening, novel procedures to reduce the time and labor involved in screening peptide ligands are required. We propose the novel approach of 'in silico panning' which consists of a two-stage screening, involving affinity selection by docking simulation and evolution of the peptide ligand using genetic algorithms (GAs). In silico panning was successfully applied to the selection of peptide inhibitor for water-soluble quinoprotein glucose dehydrogenase (PQQGDH). RESULTS: The evolution of peptide ligands for a target enzyme was achieved by combining a docking simulation with evolution of the peptide ligand using genetic algorithms (GAs), which mimic Darwinian evolution. Designation of the target area as next to the substrate-binding site of the enzyme in the docking simulation enabled the selection of a non-competitive inhibitor. In all, four rounds of selection were carried out on the computer; the distribution of the docking energy decreased gradually for each generation and improvements in the docking energy were observed over the four rounds of selection. One of the top three selected peptides with the lowest docking energy, 'SERG' showed an inhibitory effect with Ki value of 20 microM. PQQGDH activity, in terms of the Vmax value, was 3-fold lower than that of the wild-type enzyme in the presence of this peptide. The mechanism of the SERG blockage of the enzyme was identified as non-competitive inhibition. We confirmed the specific binding of the peptide, and its equilibrium dissociation constant (KD) value was calculated as 60 microM by surface plasmon resonance (SPR) analysis. CONCLUSION: We demonstrate an effective methodology of in silico panning for the selection of a non-competitive peptide inhibitor from small virtual peptide library. This study is the first to demonstrate the usefulness of in silico evolution using experimental data. Our study highlights the usefulness of this strategy for structure-based screening of enzyme inhibitors.

Acinetobacter calcoaceticus↗

Adenosine receptor modelling. A1/A2a selectivity.

Three-dimensional models of the A(1) and A(2a) adenosine receptors (AR) were constructed by means of a homology procedure, using bovine rhodopsin as a template. In order to validate the two models, a docking analysis of selective agonists was carried out. The study shows that A(1)/A(2a) selectivity is mainly influenced by the different ability of the two receptors to give lipophilic interactions, instead of giving different H bonds. The binding site cavity of the A(1)AR is smaller than that of the A(2a)AR, and for this reason, less bulky ligands like CPA are able to give close interactions with the A(1)AR, unlike larger ligands such as CGS-21680. The different dimensions of the binding site cavity could be due to the presence of three residues of proline, which cause a different rearrangement of the TM, thus modifying the side chain disposition inside the inter-helix channel.

Adenosine↗

Quaternary structure predictions of transmembrane proteins starting from the monomer: a docking-based approach.

BACKGROUND: We introduce a computational protocol for effective predictions of the supramolecular organization of integral transmembrane proteins, starting from the monomer. Despite the demonstrated constitutive and functional importance of supramolecular assemblies of transmembrane subunits or proteins, effective tools for structure predictions of such assemblies are still lacking. Our computational approach consists in rigid-body docking samplings, starting from the docking of two identical copies of a given monomer. Each docking run is followed by membrane topology filtering and cluster analysis. Prediction of the native oligomer is therefore accomplished by a number of progressive growing steps, each made of one docking run, filtering and cluster analysis. With this approach, knowledge about the oligomerization status of the protein is required neither for improving sampling nor for the filtering step. Furthermore, there are no size-limitations in the systems under study, which are not limited to the transmembrane domains but include also the water-soluble portions. RESULTS: Benchmarks of the approach were done on ten homo-oligomeric membrane proteins with known quaternary structure. For all these systems, predictions led to native-like quaternary structures, i.e. with Calpha-RMSDs lower than 2.5 A from the native oligomer, regardless of the resolution of the structural models. CONCLUSION: Collectively, the results of this study emphasize the effectiveness of the prediction protocol that will be extensively challenged in quaternary structure predictions of other integral membrane proteins.

Amino Acid Sequence↗

Conformational transitions induced by the binding of MgATP to the vitamin B12 ATP-binding cassette (ABC) transporter BtuCD.

ATP-binding cassette transporters use the free energy of ATP hydrolysis to transport structurally diverse molecules across prokaryotic and eukaryotic membranes. Computer simulation studies of the "real-time" dynamics of the ATP binding process in BtuCD, the vitamin B12 importer from Escherichia coli, demonstrate that the docking of ATP to the catalytic pockets progressively draws the two cytoplasmic nucleotide-binding cassettes toward each other. Movement of the cassettes into closer opposition in turn induces conformational rearrangement of alpha-helices in the transmembrane domain. The shape of the translocation pathway consequently changes in a manner that could aid the vectorial movement of vitamin B12. These results suggest that ATP binding may indeed represent the power stroke in the catalytic mechanism. Moreover, occlusion of ATP at one catalytic site is mechanically coupled to opening of the nucleotide-binding pocket at the second site. We propose that this asymmetry in nucleotide binding behavior at the two catalytic pockets may form the structural basis by which the transporter is able to alternate ATP hydrolysis from one site to the other.

ATP-Binding Cassette Transporters↗

Protein structure prediction via combinatorial assembly of sub-structural units.

Following the hierarchical nature of protein folding, we propose a three-stage scheme for the prediction of a protein structure from its sequence. First, the sequence is cut to fragments that are each assigned a structure. Second, the assigned structures are combinatorially assembled to form the overall 3D organization. Third, highly ranked predicted arrangements are completed and refined. This work focuses on the second stage of this scheme: the combinatorial assembly. We present CombDock, a combinatorial docking algorithm. CombDock gets an ordered set of protein sub-structures and predicts the inter-contacts that define their overall organization. We reduce the combinatorial assembly to a graph-theory problem, and give a heuristic polynomial solution to this computationally hard problem. We applied CombDock to various examples of structural units of two types: protein domains and building blocks, which are relatively stable sub-structures of domains. Moreover, we tested CombDock using increasingly distorted input, where the native structural units were replaced by similarly folded units extracted from homologous proteins and, in the more difficult cases, from globally unrelated proteins. The algorithm is robust, showing low sensitivity to input distortion. This suggests that CombDock is a useful tool in protein structure prediction that may be applied to large target proteins.

Algorithms↗

EF loop conformational change triggers ligand binding in beta-lactoglobulins.

Beta-lactoglobulins, belonging to the lipocalin family, are a widely studied group of proteins, characterized by the ability to solubilize and transport hydrophobic ligands, especially fatty acids. Despite many reports, the mechanism of ligand binding and the functional role of these proteins is still unclear, and many contradicting concepts are often encountered in the literature. In the present paper the comparative analysis of the binding properties of beta-lactoglobulins has been performed using sequence-derived information, structure-based electrostatic calculations, docking simulations, and NMR experiments. Our results reveal for the first time the mechanism of beta-lactoglobulin ligand binding, which is completely determined by the opening-closing of EF loop, triggered by Glu89 protonation. The alkaline shift observed for Glu89 pKa in porcine beta-lactoglobulin (pKa 9.7) with respect to the bovine species (pKa 5.5) depends upon the interplay of electrostatic effects of few nearby key residues. Porcine protein is therefore able to bind fatty acids provided that the appropriate pH solution conditions are met (pH > 8.6), where the EF loop conformational change can take place. The unusually high pH of binding detected for porcine beta-lactoglobulin seems to be functional to lipases activity. Theoretical pKa calculations extended to representative beta-lactoglobulins allowed the identification of key residues involved in structurally and functionally important electrostatic interactions. The results presented here provide a strong indication that the described conformational change is a common feature of all beta-lactoglobulins.

Amino Acid Sequence↗

Signaling through receptors and scaffolds: independent interactions reduce combinatorial complexity.

After activation, many receptors and their adaptor proteins act as scaffolds displaying numerous docking sites and engaging multiple targets. The consequent assemblage of a variety of protein complexes results in a combinatorial increase in the number of feasible molecular species presenting different states of a receptor-scaffold signaling module. Tens of thousands of such microstates emerge even for the initial signal propagation events, greatly impeding a quantitative analysis of networks. Here, we demonstrate that the assumption of independence of molecular events occurring at distinct sites enables us to approximate a mechanistic picture of all possible microstates by a macrodescription of states of separate domains, i.e., macrostates that correspond to experimentally verifiable variables. This analysis dissects a highly branched network into interacting pathways originated by protein complexes assembled on different sites of receptors and scaffolds. We specify when the temporal dynamics of any given microstate can be expressed using the product of the relative concentrations of individual sites. The methods presented here are equally applicable to deterministic and stochastic calculations of the temporal dynamics. Our domain-oriented approach drastically reduces the number of states, processes, and kinetic parameters to be considered for quantification of complex signaling networks that propagate distinct physiological responses.

Binding Sites↗

Synergistic use of compound properties and docking scores in neural network modeling of CYP2D6 binding: predicting affinity and conformational sampling.

Cytochrome P450 2D6 (CYP2D6) is used to develop an approach for predicting affinity and relevant binding conformation(s) for highly flexible binding sites. The approach combines the use of docking scores and compound properties as attributes in building a neural network (NN) model. It begins by identifying segments of CYP2D6 that are important for binding specificity, based on structural variability among diverse CYP enzymes. A family of distinct, low-energy conformations of CYP2D6 are generated using simulated annealing (SA) and a collection of 82 compounds with known CYP2D6 affinities are docked. Interestingly, docking poses are observed on the backside of the heme as well as in the known active site. Docking scores for the active site binders, along with compound-specific attributes, are used to train a neural network model to properly bin compounds as strong binders, moderate binders, or nonbinders. Attribute selection is used to preselect the most important scores and compound-specific attributes for the model. A prediction accuracy of 85+/-6% is achieved. Dominant attributes include docking scores for three of the 20 conformations in the ensemble as well as the compound's formal charge, number of aromatic rings, and AlogP. Although compound properties were highly predictive attributes (12% improvement over baseline) in the NN-based prediction of CYP2D6 binders, their combined use with docking score attributes is synergistic (net increase of 23% above baseline). Beyond prediction of affinity, attribute selection provides a way to identify the most relevant protein conformation(s), in terms of binding competence. In the case of CYP2D6, three out of the ensemble of 20 SA-generated structures are found to be the most predictive for binding.

Algorithms↗

A model for the three-dimensional structure of human plasma vitronectin from small-angle scattering measurements.

Small-angle X-ray scattering (SAXS) measurements were used to characterize vitronectin, a circulatory protein found in human plasma that functions in regulating cell adhesion and migration, as well as proteolytic cascades that affect blood coagulation, fibrinolysis, and pericellular proteolysis. SAXS measurements were taken over a 3-fold range of protein concentrations, yielding data that characterize a monodisperse system of particles with an average radius of gyration of 30.3 +/- 0.6 A and a maximum linear dimension of 110 A. Shape restoration was applied to the data to produce two models of the solution structure of the ligand-free protein. A low-resolution model of the protein was generated that indicates the protein to be roughly peanut-shaped. A better understanding of the domain structure of vitronectin resulted from low-resolution models developed from available high-resolution structures of the domains. These domains include the N-terminal domain that was determined experimentally by NMR [Mayasundari, A., Whittemore, N. A., Serpersu, E. H., and Peterson, C. B. (2004) J. Biol. Chem. 279, 29359-29366] and the docked structure of the central and C-terminal domains that were determined by computational threading [Xu, D., Baburaj, K., Peterson, C. B., and Xu, Y. (2001) Proteins: Struct., Funct., Genet. 44, 312-320]. This model provides an indication of the disposition of the central domain and C-terminal heparin-binding domains of vitronectin with respect to the N-terminal somatomedin B (SMB) domain. This model constructed from the available domain structures, which agrees with the low-resolution model produced from the SAXS data, shows the SMB domain well separated from the central and heparin-binding domains by a disordered linker (residues 54-130). Also, binding sites within the SMB domain are predicted to be well exposed to the surrounding solvent for ease of access to its various ligands.

Computational Biology↗

Insights into the selective inhibition of Candida albicans secreted aspartyl protease: a docking analysis study.

Severe fungal infections have taken precedence over other bacterial infections. Of the several fungal species, Candida albicans and others belonging to the genus Candida are responsible for several clinically important fungal infections. Emerging cases of drug resistance to the currently available drugs has limited the spectrum of currently available antifungal agents. Thus, it is imperative that new biochemical targets are identified so that better effective and selective agents can be developed. Many enzymes contribute towards the complex disease process of fungal infections; the secreted aspartyl protease (SAP), expressed both in vitro and during infection, has been implicated as one of the major virulence factors of C. albicans. Three-dimensional crystal structures of C. albicans SAP and closely related clinical isolate designated as SAP2X complexed with the same potent inhibitor A-70450 have been reported. Several analogues of A-70450 with potent C. albicans SAP2X inhibitory activity are also known. However, the structural effects of the binding of these compounds with the enzyme active site are not completely understood. Our efforts in this direction involve the docking analysis of C. albicans SAP2X inhibitors complexed with SAP2X enzyme, which is reported in this work. Docking analysis was performed on a set of molecules with differing selectivities and inhibitory potencies towards C. albicans, renin and cathepsin D. The structural effects of ligand binding were analyzed on the basis of hydrophobic and hydrogen bond interactions, binding energy analysis, interaction energies, rms deviations, etc. in the resulting energy-minimized structures of the receptor-ligand complexes. Structural analysis of the resulting models indicates that hydrophobic and hydrogen bonding interactions together with binding and interaction energies are responsible for selective inhibition of C. albicans SAP2X. Hydrophobic and hydrogen bonding interactions in the various subsites of the enzyme, contributing to both increase as well as decrease in selectivity of the molecules have been detailed. Hydrogen bonding interaction plays an important role for amino acid residues such as Gly-85, Asp-86, Asp-32, Asp-218, Tyr-225, Ala-133, and so on. Significant hydrophobic interactions with the S3, S2 and S2' subsites contribute to selectivity of the compounds. These molecular modeling analyses should, in our view, contribute for further development of selective C. albicans secreted aspartyl protease inhibitors.

Antifungal Agents↗

Data-driven docking: HADDOCK's adventures in CAPRI.

We have shown previously that given high-resolution structures of the unbound molecules, structure determination of protein complexes is possible by including biochemical and/or biophysical data as highly ambiguous distance restraints in a docking approach. We applied this method, implemented in the HADDOCK (High Ambiguity Driven DOCKing) package (Dominguez et al., J Am Chem Soc 2003;125:1731-1737), to the targets in the fourth and fifth rounds of CAPRI. Here we describe our results and analyze them in detail. Special attention is given to the role of flexibility in our docking method and the way in which this improves the docking results. We describe extensions to our approach that were developed as a direct result of our participation in CAPRI. In addition to experimental information, we also included interface residue predictions from PPISP (Protein-Protein Interaction Site Predictor; Zhou and Shan, Proteins 2001;44:336-343), a neural network method. Using HADDOCK we were able to generate acceptable structures for 6 of the 8 targets, and to submit at least 1 acceptable structure for 5 of them. Of these 5 submissions, 3 were of medium quality (Targets 10, 11, and 15) and 2 of high quality (Targets 13 and 14). In all cases, predictions were obtained containing at least 40% of the correct epitope at the interface for both ligand and receptor simultaneously.

Algorithms↗

DNA intercalative potential of marketed drugs testing positive in in vitro cytogenetics assays.

We have previously noted that the Physicians' Desk Reference (PDR) contains over 80 instances in which a drug elicited a positive genotoxic response in one or more in vitro assays, despite having no obvious structural features predictive of covalent drug/DNA interactive potential or known mechanistic basis. Furthermore, in most cases, these drugs were "missed" by computational genotoxicity-predicting models such as DEREK, MCASE and TOPKAT. We have previously reported the application of a V79 cell-based model and a 3D DNA docking model for predicting non-covalent chemical/DNA interactions. Those studies suggested that molecules that are very widely structurally diverse may be capable of intercalating into DNA. To determine whether such non-covalent drug/DNA interactions might be involved in unexpected drug genotoxicity, we evaluated, using both models where possible, 56 marketed pharmaceuticals, 40 of which were reported as being clastogenic in in vitro cytogenetics assays (chromosome aberrations/mouse lymphoma assay). As seen before, the two approaches showed good concordance (62%) and 26 of the 40 (65%) drugs exhibiting in vitro clastogenicity were predicted as intercalators by one or both methods. This finding provides support for the hypothesis that non-covalent DNA interaction may be a common mechanism of clastogenicity for many drugs having no obvious structural alerts for covalent DNA interaction.

Animals↗

Iron(II) triggered conformational changes in Escherichia coli fur upon DNA binding: a study using molecular modeling.

In order to identify the Fur dimerization domain, a three-dimensional structure of the ferric uptake regulation protein from Escherichia coli (Fur EC) was determined using homology modeling and energy minimization. The Fur monomer consists of turn- helix -turn motif on the N-terminal domain, followed by another helix-turn-helix-turn motif, and two beta-strands separated by a turn which forms the wing. The C-terminal domain, separated by a long coil from the N-terminal, and consisting of two anti parallel beta strands, and a turn-helix-turn-helix-turn motif. Residues in central domain were found to aid the dimer formation, residues 45-70 as evident in the calculated distances; this region is rich in hydrophobic residues. Most interactions occur between residues Val55, Leu53, Gln52, Glu49 and Tyr56 with closest contacts occurring at residues 49-56. These residues are part of an alpha-helix (alpha(4)) near the N-terminal. Upon raising the Fe(2+) concentration the binding of Fur dimer to DNA was enhanced, this was evident when, the Fur EC dimer was docked onto DNA "iron box" (it was found to bind the AT-rich region) and upon addition of Fe(2+) the helices near the N-terminal bound to the major groove of the DNA. Addition of high Fe(2+) concentration triggered further conformational changes in the Fur dimer as was measured by distances between the two subunits, Fe(2+) mediated the Fur binding to DNA by attaching itself to the DNA. At the same time DNA changed conformation as was evident in the distortion in the backbone and the shrinking of major groove distance from 11.4 to 9.3A. Two major Fe(2+) sites were observed on the C-terminal domain: site 1, the traditional Zn site, the cavity contains the residues Cys92, Cys95, Asp137, Asp141, Arg139, Glu 140, His 145 and His 143 at distances range from 1.3 to 2.2A. Site 2 enclave consists of His71, Ile50, Asn72, Gly97, Asp105 and Ala109 at very close proximity to Fe(2+). The closest contacts between Fur dimer and DNA at the AT-rich region were at residues Ala11, Gly12, Leu13, Pro18 and Arg19 mostly hydrophobic residues near the N-terminal domain. Close contacts repeated at His87, His88 and Arg112, and a third region near the C-terminal at Asn137, Arg 139, Glu140, Asn141, His143, Asn141 and His145. Fur dimer has three major contact regions with DNA, the first on the N-terminal domain, a second smaller region at His87, His88 and Arg112 mediated by Fe(2+) ions, and a third region on the C-terminal domain consisting mainly of hydrophobic contacts and mediated by Fe(2+) ions at high concentration.

Amino Acid Sequence↗

Free energy perturbation approach to the critical assessment of selective cyclooxygenase-2 inhibitors.

The discovery of selective cyclooxygenase-2 (COX-2) inhibitors represents a major achievement of the efforts over the past few decades to develop therapeutic treatments for inflammation. To gain insights into designing new COX-2-selective inhibitors, we address the energetic and structural basis for the selective inhibition of COX isozymes by means of a combined computational protocol involving docking experiment, force field design for the heme prothetic group, and free energy perturbation (FEP) simulation. We consider both COX-2- and COX-1-selective inhibitors taking the V523I mutant of COX-2 to be a relevant structural model for COX-1 as confirmed by a variety of experimental and theoretical evidences. For all COX-2-selective inhibitors under consideration, we find that free energies of binding become less favorable as the receptor changes from COX-2 to COX-1, due to the weakening and/or loss of hydrogen bond and hydrophobic interactions that stabilize the inhibitors in the COX-2 active site. On the other hand, COX-1-selective oxicam inhibitors gain extra stabilization energy with the change of residue 523 from valine to isoleucine because of the formations of new hydrogen bonds in the enzyme-inhibitor complexes. The utility of the combined computational approach, as a valuable tool for in silico screening of COX-2-selective inhibitors, is further exemplified by identifying the physicochemical origins of the enantiospecific selective inhibition of COX-2 by alpha-substituted indomethacin ethanolamide inhibitors.

Amino Acid Sequence↗

Soft docking an L and a D peptide to an anticholera toxin antibody using internal coordinate mechanics.

BACKGROUND: The tremendous increase in sequential and structural information is a challenge for computer-assisted modelling to predict the binding modes of interacting biomolecules. One important area is the structural understanding of protein-peptide interactions, information that is increasingly important for the design of biologically active compounds. RESULTS: We predicted the three-dimensional structure of a complex between the monoclonal antibody TE33 and its cholera-toxin-derived peptide epitope VPGSQHID. Using the internal coordinate mechanics (ICM) method of flexible docking, the bound conformation of the initially extended peptide epitope to the antibody crystal or modelled structure reproduced the known binding conformation to a root mean square deviation of between 1.9 A and 3.1 A. The predicted complexes are in good agreement with binding data obtained from substitutional analyses in which each epitope residue is replaced by all other amino acids. Furthermore, a de novo prediction of the recently discovered TE33-binding D peptide dwGsqhydp (single-letter amino acid code where D amino acids are represented by lower-case letters) explains results obtained from binding studies with 172 peptide analogues. CONCLUSIONS: Despite the difficulties arising from the huge conformational space of a peptide, this approach allowed the prediction of the correct binding orientation and the majority of essential binding features of a peptide-antibody complex.

Antibodies, Monoclonal↗

Clustering of class I HLA oligomers with CD8 and TCR: three-dimensional models based on fluorescence resonance energy transfer and crystallographic data.

Fluorescence resonance energy transfer (FRET) data, in accordance with lateral mobility measurements, suggested the existence of class I HLA dimers and oligomers at the surface of live human cells, including the B lymphoblast cell line (JY) used in the present study. Intra- and intermolecular class I HLA epitope distances were measured on JY B cells by FRET using fluorophore-conjugated Ag-binding fragments of mAbs W6/32 and L368 directed against structurally well-characterized heavy and light chain epitopes, respectively. Out-of-plane location of these epitopes relative to the membrane-bound BODIPY-PC (2-(4,4-difluoro-5-(4-phenyl-1,3-butadienyl)-4-bora-3a,4a-diaza-s-indacene-3-pentanoyl)-1-hexadecanoyl-sn-glycero-3-phosphocholine) was also determined by FRET. Computer-simulated docking of crystallographic structures of class I HLA and epitope-specific Ag-binding fragments, with experimentally determined interepitope and epitope to cell surface distances as constraints, revealed several sterically allowed and FRET-compatible class I HLA dimeric and tetrameric arrangements. Extension of the tetrameric class I HLA model with interacting TCR and CD8 resulted in a model of a supramolecular cluster that may exist physiologically and serve as a functionally significant unit for a network of CD8-HLA-I complexes providing enhanced signaling efficiency even at low MHC-peptide concentrations at the interface of effector and APCs.

Antigen-Presenting Cells↗

Structural motif-based homology modeling of CYP27A1 and site-directed mutational analyses affecting vitamin D hydroxylation.

Human CYP27A1 is a mitochondrial cytochrome P450, which is principally found in the liver and plays important roles in the biological activation of vitamin D(3) and in the biosynthesis of bile acids. We have applied a systematic analysis of hydrogen bonding patterns in 11 prokaryotic and mammalian CYP crystal structures to construct a homology-based model of CYP27A1. Docking of vitamin D(3) structures into the active site of this model identified potential substrate contact residues in the F-helix, the beta-3 sheet, and the beta-5 sheet. Site-directed mutagenesis and expression in COS-1 cells confirmed that these positions affect enzymatic activity, in some cases shifting metabolism of 1alpha-hydroxyvitamin D(3) to favor 25- or 27-hydroxylation. The results suggest that conserved hydrophobic residues in the beta-5 hairpin help define the shape of the substrate binding cavity and that this structure interacts with Phe-248 in the F-helix. Mutations directed toward the beta-3a strand suggested a possible heme-binding interaction centered on Asn-403 and a structural role for substrate contact residues Thr-402 and Ser-404.

Amino Acid Motifs↗