Effects of protein flexibility on the site of metabolism prediction for CYP2A6 substrates

Structure-based prediction for the site of metabolism (SOM) of a compound metabolized by human cytochrome P450s (CYPs) is highly beneficial in drug discovery and development. However, the flexibility of the CYPs’ active site remains a huge challenge for accurate SOM prediction. Compared with other CYPs, the active site of CYP2A6 is relatively small and rigid. To address the impact of the flexibility of CYP2A6 active site residues on the SOM prediction for substrates, in this work, molecular dynamics (MD) simulations and molecular docking were used to predict the SOM of 96 CYP2A6 substrates. Substrates with known SOM were docked into the snapshot structures from MD simulations and the crystal structures of CYP2A6. Compared to the crystal structures, the protein structures obtained from MD simulations showed more accurate prediction for SOM. Our results indicated that the flexibility of the active site of CYP2A6 significantly affects the SOM prediction results. Further analysis for the 40 substrates with definite K_m values showed that the prediction accuracy for the low K_m substrates is comparable to that by ligand-based methods.     

Modeling of the three-dimensional structure of the human melanocortin 1 receptor, using an automated method and docking of a rigid cyclic melanocyte-stimulating hormone core peptide

A model is presented of the melanocortin 1 receptor (MC1R), constructed by use of an unbiased, objective method. The model is created directly from data derived from multiple sequence analysis, a low-resolution EM-projection map of rhodopsin, and the approximate membrane thickness. The model agrees well with available data concerning natural mutations of MC1Rs occurring in different species. A model is also presented of the most rigid ligand for this receptor, the cyclic pentapeptide cHFRWG, shown docked in the receptor model. The receptor–ligand complex model agrees well with available experimental data. The ligand is located between transmembrane region 1 (TM1), TM2, TM3, TM6, and TM7 of the receptor. Multiple interactions occur between ligand and receptor, including interactions with Leu-48 (TM1), Ser-52 (TM1), Glu-55 (TM1), Asn-91 (TM2), Glu-94 (TM2), Thr-95 (TM2) Ile-98 (TM2), Asp-121 (TM3), Thr-124 (TM3), Phe-257 (TM6), Phe-283 (TM7), Asn-290 (TM7), and Asp-294 (TM7) of the receptor.     

Prediction of three-dimensional structures and structural flexibilities of wild-type and mutant cytochrome P450 1A2 using molecular dynamics simulations

In this study, we investigated the effect of genetic polymorphism on the three-dimensional (3D) conformation of cytochrome P450 1A2 (CYP1A2) using molecular dynamics (MD) simulations. CYP1A2, a major drug-metabolizing enzyme among cytochrome P450 enzymes (CYPs), is known to have many variant alleles. The genetic polymorphism of CYP1A2 may cause individual differences in the pharmacokinetics of medicines. By performing 100 ns or longer MD simulations, we investigated the influence of amino acid mutation on the 3D structures and the dynamic properties of proteins. The results show that the static structures were changed by the mutations of amino acid residues, not only near the mutated residues but also in distant portions of the proteins. Moreover, the mutation of only one amino acid was shown to change the structural flexibility of proteins, which may influence the substrate recognition and enzymatic activity. Our results clearly suggest that it is necessary to investigate the dynamic property as well as the static 3D structure for understanding the change of the enzymatic activity of mutant CYP1A2.     

Glycosylation of prions and its effects on protein conformation relevant to amino acid mutations

The three-dimensional coordinates from a nuclear magnetic resonance (NMR)-averaged structure containing residues 121-226 of mouse prion were used as the starting geometry for MD of prion either with or without glycan in both mutant and wild-type forms. The following mutants were studied: Asp-178 to Asn, Thr-183 to Ala, Phe-198 to Ser, Glu-200 to Lys, and Gln-217 to Arg. NMR data vs structural models were compared to observe any major differences. Simulations of the change in protein structure with and without glycan were performed, as they cannot be tested by NMR analysis. Several mutants were expressed and analyzed for altered glycosylation and the results interpreted in terms of molecular modeling. N-linked glycosylation is likely to play an important role in prion biology as shown by visualization of glycoprotein conformation.     

Automated luminescence-based cytochrome P450 profiling using a simple, elegant robotic platform

The determination of inhibitory effects that lead compounds have on cytochrome P450 (CYP) enzymes is an important part of today's drug discovery process. Assays can be performed early in the discovery process to predict adverse drug-drug interactions caused by CYP inhibition and to minimize the costs associated with terminating candidates in late stage development or worse, removing a drug from the market after launch. For early discovery work, testing substantial numbers of compounds is desirable, thus automated "mix and read" assays are beneficial. Here, we demonstrate the automation of the CYP profiling process using a simple, yet robust robotic platform. Compound titration, as well as transfer of compounds and assay components was performed by the same automated pipetting system. IC(50)s of small molecule drugs were determined using recombinant CYP enzymes, CYP3A4, -2C9, and -2D6 and luminogenic substrates specific to each. Compounds were profiled against all three enzymes on the same 384-well assay plate.     

Structural insights into the camel milk lactoperoxidase: Homology modeling and molecular dynamics simulation studies

Lactoperoxodase (LPO) is a heme peroxidase enzyme present in mammalian milk. It is an antimicrobial protein with wide range of industrial applications. Although the three dimensional structure of LPO from various mammalian species has been reported, but its structure from camel source is still unknown. So far, the crystallization attempts have not been successful in determining camel LPO (cLPO) structure. Herein, we developed the three dimensional structure of cLPO by homology modeling approach using prime module available in Schrodinger suite. The developed model in complex with ligand hypothiocyanate (OSCN⁻) was further validated by Ramachandran plot followed by molecular dynamics (MD) simulation studies using Desmond module of Schrodinger. cLPO model exhibited overall structural similarity with template crystal structure, however, it displayed different interaction pattern of amino acid residues with ligand OSCN⁻ in comparison to template crystal structure. Moreover, the ligand binding site environment in cLPO is more polar, less hydrophobic, and harbours more number of charged residues than template crystal structure. The substrate binding pocket environment of cLPO shows a considerable difference from template crystal structure. This subsequently resulted in dissimilar behaviour of ligand during the course of MD simulation studies.     

Elucidating a chemical defense mechanism of Antarctic sponges: A computational study

In 2000, a novel secondary metabolite (erebusinone, Ereb) was isolated from the Antarctic sea sponge, Isodictya erinacea. The bioactivity of Ereb was investigated, and it was found to inhibit molting when fed to the arthropod species Orchomene plebs. Xanthurenic acid (XA) is a known endogenous molt regulator present in arthropods. Experimental studies have confirmed that XA inhibits molting by binding to either (or both) of two P450 enzymes (CYP315a1 or CYP314a1) that are responsible for the final two hydroxylations in the production of the molt-inducing hormone, 20-hydroxyecdysone (20E). The lack of crystal structures and biochemical assays for CYP315a1 or CYP314a1, has prevented further experimental exploration of XA and Ereb's molt inhibition mechanisms. Herein, a wide array of computational techniques – homology modeling, molecular dynamics simulations, binding site bioinformatics, flexible receptor–flexible ligand docking, and molecular mechanics-generalized Born surface area calculations – have been employed to elucidate the structure–function relationships between the aforementioned P450s and the two described small molecule inhibitors (Ereb and XA). Results indicate that Ereb likely targets CYP315a1 by interacting with a network of aromatic residues in the binding site, while XA may inhibit both CYP315a1 and CYP314a1 because of its aromatic, as well as charged nature.     

Structural insights into the interaction between molluscan hemocyanins and phenolic substrates: An in silico study using docking and molecular dynamics

Hemocyanin is a multimeric type-3 copper containing oxygen carrier protein that exhibits phenoloxidase-like activity and is found in selected species of arthropoda and mollusca. The phenoloxidase activity in the molluscan hemocyanins can be triggered by the proteolytic removal of the C-terminal β-rich sandwich domain of the protein or by the treatment with chemical agents like SDS, both of which enable active site access to the phenolic substrates. The mechanism by which SDS treatment enhances active site access to the substrates is however not well understood in molluscan hemocyanins. Here, using a combination of in silico molecular dynamics (MD) and docking studies on the crystal structure of Octopus dofleini hemocyanin (PDB code:1JS8), we demonstrate that the C-terminal β-domain of the protein plays a crucial role in regulating active site access to bulky phenolic substrates. Furthermore, MD simulation of hemocyanin in SDS revealed displacement of β-domain, enhanced active site access and a resulting increase in binding affinity for substrates. These observations were further validated by enzyme kinetics experiments.     

Modeling the protonation states of β-secretase binding pocket by molecular dynamics simulations and docking studies

β-secretase (BACE1) is an aspartyl protease that processes the β-amyloid peptide in the human brain in patients with Alzheimer’s disease. There are two catalytic aspartates (ASP32 and ASP228) in the active domain of BACE1. Although it is believed that the net charge of the Asp dyad is −1, the exact protonation state still remains a matter of debate. We carried out molecular dynamic (MD) simulations for the four protonation states of BACE1 proteins. We applied Glide docking studies to 21 BACE1 inhibitors against the MD extracted conformations. The dynamic results infer that the protein/ligand complex remains stable during the entire simulation course for HD32D228 model. The results show that the hydrogen bonds between the inhibitor and the Asp dyad are maintained in the 10,000th ps snapshot of HD32D228 model. Our results also reveal the significant loop residues in maintaining the active binding conformation in the HD32D228 model. Molecular docking results show that the HD32D228 model provided the best enrichment factor score, suggesting that this model was able to recognize the most active compounds. Our observations provide an evidence for the preference of the anionic state (HD32D228) in BACE1 binding site and are in accord with reported computational data. The protonation state study would provide significant information to assign the correct protonation state for structure-based drug design and docking studies targeting the BACE1 proteins as a tactic to develop potential AD inhibitors.     

In silico design: Extended molecular dynamic simulations of a new series of dually acting inhibitors against EGFR and HER2

Based on the hit structures that have been identified in our previous studies against EGFR and HER2, new potential inhibitors that share the same scaffold of the hit structures are designed and screened in silico. Insights into understanding the potential inhibitory effect of the new inhibitors against both EGFR and HER2 receptors is obtained using extended molecular dynamics (MD) simulations and different scoring techniques. The binding mechanisms and dynamics are detailed with respect to two approved inhibitors against EGFR (lapatinib) and HER2 (SYR127063). The best scoring inhibitor (T9) is chosen for additional in silico investigation against both the wild-type and T790M mutant strain of EGFR and the wild-type HER2. The results reveal that certain substitution patterns increase the stability and assure stronger binding and higher H-bond occupancy of the conserved water molecule that is commonly observed with kinase crystal structures. Furthermore, the new inhibitor (T9) forms stable interactions with the mutant strain as a direct consequence of the enhanced ability to form additional hydrogen bonding interactions with binding site residues.     

The mean hydration of carbohydrates as studied by normalized two-dimensional radial pair distributions

The hydration of carbohydrates plays a key role in many biological processes. Molecular dynamics simulations provide an effective tool for investigating the hydration of complex solutes such as carbohydrates. In this article we devise an algorithm for the calculation of two-dimensional radial pair distributions describing the probability of finding a water molecule in a site defined by two reference atoms. The normalized 2D radial pair distribution is proposed as an effective tool for investigating and comparing localized or ordered water sites around flexible molecules such as carbohydrates when analyzing molecular dynamics simulations and the utility of 2D radial pair distributions is demonstrated using sucrose as an example. In this relatively simple structure, 2D radial pair distributions were able to characterize and quantify the importance of two unique interresidue hydration sites in which a water molecule is forming a bridge between the glycopyranosyl and fructofuranosyl residues. The approach is proposed to be a valuable tool for comparing and understanding the hydration of flexible biomolecules such as carbohydrates.     

A putative bioactive conformation for the altered peptide ligand of myelin basic protein and inhibitor of experimental autoimmune encephalomyelitis [Arg91, Ala96] MBP87-99

[Arg(91), Ala(96)] MBP(87-99) is an altered peptide ligand (APL) of myelin basic protein (MBP), shown to actively inhibit experimental autoimmune encephalomyelitis (EAE), which is studied as a model of multiple sclerosis (MS). The APL has been rationally designed by substituting two of the critical residues for recognition by the T-cell receptor. A conformational analysis of the APL has been sought using a combination of 2D NOESY nuclear magnetic resonance (NMR) experiments and detailed molecular dynamics (MD) calculations, in order to comprehend the stereoelectronic requirements for antagonistic activity, and to propose a putative bioactive conformation based on spatial proximities of the native peptide in the crystal structure. The proposed structure presents backbone similarity with the native peptide especially at the N-terminus, which is important for major histocompatibility complex (MHC) binding. Primary (Val(87), Phe(90)) and secondary (Asn(92), Ile(93), Thr(95)) MHC anchors occupy the same region in space, whereas T-cell receptor (TCR) contacts (His(88), Phe(89)) have different orientation between the two structures. A possible explanation, thus, of the antagonistic activity of the APL is that it binds to MHC, preventing the binding of myelin epitopes, but it fails to activate the TCR and hence to trigger the immunologic response. NMR experiments coupled with theoretical calculations are found to be in agreement with X-ray crystallography data and open an avenue for the design and synthesis of novel peptide restricted analogues as well as peptide mimetics that rises as an ultimate goal.     

Effects of mutations on active site conformation and dynamics of RNA-dependent RNA polymerase from Coxsackievirus B3

Recent crystal structures of RNA-dependent RNA polymerase (3D^pol) from Coxsackievirus B3 (CVB3) revealed that a tyrosine mutation at Phe364 (F364Y) resulted in structures with open active site whereas a hydrophobic mutation at Phe364 (F364A) led to conformations with closed active site. Besides, the crystal structures showed that the F364W mutation had no preference between the open and closed active sites, similar to wild-type. In this paper, we present a molecular dynamics (MD) study on CVB3 3D^pol in order to address some important questions raised by experiments. First, MD simulations of F364Y and F364A were carried out to explore how these mutations at Phe364 influence active site dynamics and conformations. Second, MD simulations of wild-type and mutants were performed to discover the connection between active site dynamics and polymerase function. MD simulations reveal that the effect of mutations on active site dynamics is associated with the interaction between the structural motifs A and D in CVB3 3D^pol. Interestingly, we discover that the active site state is influenced by the formation of a hydrogen bond between backbone atoms of Ala231 (in motif A) and Ala358 (in motif D), which has never been revealed before.     

Catalytic mechanism of l,l-diaminopimelic acid with diaminopimelate epimerase by molecular docking simulations

Peptidoglycan, a key constituent of bacterial cell walls, is currently the target of broad spectrum antibiotics and a new research field involves both design and synthesis of inhibitors of its biosynthesis. Most bacteria require either lysine, or its biosynthetic precursor, diaminopimelate (meso-DAP), as a component of the peptidoglycan layer of the cell wall. In this paper, molecular modelling studies were undertaken in order to shed light on the molecular basis of interaction between (2S,6S)-diaminopimelic acid (l,l-DAP) (1) with its target enzyme DAP-epimerase, since this is a key step in the lysine biosynthetic path leading to (2R,6S)-diaminopimelic acid (meso-DAP) (2). In particular, the docking of the ligand-enzyme complex was studied by means of MD simulations and DFT computations in order to ascertain the optimal structural requirements for the epimerization reaction. Molecular dynamics simulations clearly showed that the configuration of the distal carbon C6 of l,l-DAP is critical for complex formation since both amino and carboxylate groups are involved in Hbonding interactions with the active site residues. Furthermore, the interactions occurring between the functional groups bonded to the C2 and some residues of the binding cavity immobilize the ligand in a position appropriate for the epimerization reaction, i.e., exactly in the middle of the two catalytic residues Cys73 and Cys217 as confirmed by DFT quantum mechanical computation of the Michaelis complex. All this mechanistic information could be useful for the rational design of new potential antibiotic drugs effective as inhibitors of peptidoglycan biosynthesis.     

Computational study of human phosphomannose isomerase: Insights from homology modeling and molecular dynamics simulation of enzyme bound substrate

Phosphomannose isomerase is a zinc metalloenzyme that catalyzes the reversible isomerization of mannose-6-phosphate and fructose-6-phosphate, and the three-dimensional (3D) structure of human phosphomannose isomerase has not been reported. In order to understand the catalytic mechanism, the 3D structure of the protein is built by using homology modeling based on the known crystal structure of mannose-6-phosphate isomerase from (PDB code 1PMI). The model structure is further refined by energy minimization and molecular dynamics methods. The mannose-6-phosphate-enzyme complex is developed by molecular docking and the key residues involved in the ligand binding are determined, which will facilitate the understanding of the action mode of the ligands and guide further genetic studies. Our results suggest a hydride transfer mechanism of alpha-hydrogen between the C1 and C2 positions but do not support the cis-enediol mechanism. The detailed mechanism involves, on one side, Zn2+ mediating the movement of a proton between O1 and O2, and, on the other side, the hydrophobic environment formed in part by Tyr278 promoting transfer of a hydride ion.     

Structure and dynamics of 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase

Folates are essential for life. Unlike mammals, most microorganisms must synthesize folates de novo. 6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes pyrophosphoryl transfer from ATP to 6-hydroxymethyl-7,8-dihydropterin (HP), the first reaction in folate pathway, and therefore, is an ideal target for developing novel antimicrobial agents. Because of its small size and high thermal stability, E. coli HPPK is also an excellent model enzyme for studying the mechanisms of enzymatic pyrophosphoryl transfer. We have determined the crystal structures of HPPK in the unligated form and in complex with HP, two Mg2+ ions, and AMPCPP (an ATP analog that inhibits the enzymatic reaction). Comparison of the two crystal structures reveals dramatic conformational changes of three flexible loops and many side chains and possible roles of the active site residues.     

A molecular dynamics study of the BACE1 conformational change from Apo to closed form induced by hydroxyethylamine derived compounds

BACE1 is an aspartyl protease which is a therapeutic target for Alzheimer’s disease (AD) because of its participation in the rate-limiting step in the production of Aβ-peptide, the accumulation of which produces senile plaques and, in turn, the neurodegenerative effects associated with AD. The active site of this protease is composed in part by two aspartic residues (Asp93 and Asp289). Additionally, the catalytic site has been found to be covered by an antiparallel hairpin loop called the flap. The dynamics of this flap are fundamental to the catalytic function of the enzyme. When BACE1 is inactive (Apo), the flap adopts an open conformation, allowing a substrate or inhibitor to access the active site. Subsequent interaction with the ligand induces flap closure and the stabilization of the macromolecular complex. Further, the protonation state of the aspartic dyad is affected by the chemical nature of the species entering the active site, so that appropriate selection of protonation states for the ligand and the catalytic residues will permit the elucidation of the inhibitory pathway for BACE1. In the present study, comparative analysis of different combinations of protonation states for the BACE1-hydroxyethylamine (HEA) system is reported. HEAs are potent inhibitors of BACE1 with favorable pharmacological and kinetic properties, as well as oral bioavailability. The results of Molecular Dynamics (MD) simulations and population density calculations using 8 different parameters demonstrate that the LnAsp289 configuration (HEA with a neutral amine and the Asp289 residue protonated) is the only one which permits the expected conformational change in BACE1, from apo to closed form, after flap closure. Additionally, differences in their capacities to establish and maintain interactions with residues such as Asp93, Gly95, Thr133, Asp289, Gly291, and Asn294 during this step allow differentiation among the inhibitory activities of the HEAs. The results and methodology here reported will serve to elucidate the inhibitory pathway of other families of compounds that act as BACE1 inhibitors, as well as the design of better leader compounds for the treatment of AD.     

The three-dimensional model of Dictyostelium discoideum racE based on the human rhoA-GDP crystal structure

The three-dimensional structure of racE was modeled using several homologous small G proteins, and the best model obtained using the human rhoA as modeling template is reported. The three-dimensional fold of the racE model is remarkably similar to the cellular form of human ras p21 crystal structure. Its secondary structure consists of six alpha-helices, six beta-strands and three 3(10) helices. The model retains its secondary structure after a 300 K, 300 ps molecular dynamics (MD) simulation. Important domains of the protein include its effector loop (residues 34-46), the insertion domain (residues 121-136), and the polybasic motif (between 210 and 220) not modeled in the current structure. The effector loop is inherently flexible and the structure docked with GDP exhibits the effector loop moving significantly closer to the nucleotide binding pocket, forming a tighter complex with the bound GDP. The mobility of the effector loop is conferred by a single residue 'hinge' point at residue 34Asp, also allowing the Switch I region, immediately preceding the effector loop, to be equally mobile. In comparison, the Switch II region shows average mobility. The insertion domain is highly flexible, with the insertion taking the form of a helical domain, with several charged residues forming a complex charged interface over the entire insertion region. While the GDP moiety is loosely held in the active site, the metal cation is extensively co-ordinated. The critical residue 38Thr exhibits high mobility, and is seen interacting directly with the metal ion at a distance of 2.64 A, and indirectly via an intervening water molecule. 64Gln, a key residue involved in GTP hydrolysis in ras, is seen facing the beta-phosphate group and the metal ion. Certain residues (i.e. 51Asn, 38Thr and 65Glu) exhibit unique characteristics and these residues, together with 158Val, may play important roles in the maintenance of the protein's integrity and function. There is strong consensus of secondary structural elements between models generated using various templates, such as h-rac1, h-rhoA and h-cdc42 bound to RhoGDI, all sharing only 50-55% sequence identity with racE, which suggests that this model is in all probability an accurate prediction of the true tertiary structure of racE.     

Specificity of damage recognition and catalysis of DNA repair

A common feature of DNA repair enzymes is their ability to recognize the damage independently of sequence in which they are found. The presence of a flipped out base inserted into the protein in several DNA-enzyme complexes suggests a contribution to enzyme specificity. Molecular simulations of damaged DNA indicate that the damage produces changes in DNA structure and changes the dynamics of DNA bending. The reduced bending force constant can be used by the enzyme to induce DNA bending and facilitate base flipping. We show that a thymine dimer (TD) containing DNA requires less energy to bend, lowering the barrier for base flipping. On the other hand, bending in DNA with U-G mismatch is affected only by a small amount and flipping is not enhanced significantly. T4 endonuclease V (endoV), which recognizes TD, utilizes the reduced barrier for flipping as a specific recognition element. In uracil DNA glycosylase (UDG), which recognizes U-G mismatches, base flipping is not enhanced and recognition is encoded in a highly specific binding pocket for the flipped base. Simulations of UDG and endoV in complex with damaged DNA provide insight into the essential elements of the catalytic mechanism. Calculations of pKas of active site residues in endoV and endoV-DNA complex show that the pKa, of the N-terminus is reduced from 8.01 to 6.52 while that of Glu-23 increases from 1.52 to 7.82. Thus, the key catalytic residues are in their neutral form. The simulations also show that Glu-23 is also H-bonded to O4' of the 5'-TD enhancing the nucleophilic attack on Cl and that Arg-26 enhances the hydrolysis by electrostatic stabilization but does not participate in proton transfer. In the enzyme-substrate complex of UDG, the role of electrostatic stabilization is played by His-268, whose pKa increases to 7.1 from 4.9 in the free enzyme. The pKa of Asp-145, the other important catalytic residue, remains around 4.2 in the free enzyme and in the complex. Thus, it can not act as a proton acceptor. In the complex the 3'-phosphate of uracil is stabilized next to Asp-145 by two bridging water molecules. Such a configuration activates one water molecule to act as a proton acceptor to produce a stabilizing hydronium ion and the other as a proton donor to produce the nucleophilic hydroxide. It appears that DNA glycosylases share commonalties in recognition of damage but differ in their catalytic mechanisms.