Alpha7 helix plays an important role in the conformational stability of PTP1B

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix α7), which forms an allosteric binding site 20 ? away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

Alpha7 Helix Plays an Important Role in the Conformational Stability of PTP1B

Abstract

The C-terminus of Protein Tyrosine Phosphatase 1B (PTP1B) includes an α-helix (α7), which forms an allosteric binding site 20 Å away from the active site. This helix is specific to PTP1B and its truncation decreases the catalytic activity significantly. Here, molecular dynamics (MD) simulations in the presence and absence of α7 were performed to investigate the role played by α7. The highly mobile α7 was found to maintain its contacts with loop 11 (L11)- α3 helix throughout the simulations. The interactions of Tyr152 on L11, Tyr176, Thr177 on the catalytically important WPD loop and Ser190 on α3 are important for the conformational stability and the concerted motions of the regions surrounding the WPD loop. In the absence of α7, L11 and WPD loop move away from their crystal structure conformations, resulting in the loss of the interactions in this region, and a decrease in the residue displacement correlations in the vicinity of WPD loop. Therefore, we suggest that one of the functionally important roles of α7 may be to limit the L11 and α3 motions, and, facilitate the WPD loop motions. Truncation of α7 in PTP1B is found to affect distant regions as well, such as the substrate recognition site and the phosphate binding-loop (P-loop), changing the conformations of these regions significantly. Our results show that the PTP1B specific α7 is important for the conformation and dynamics of the WPD loop, and also may play a role in ligand binding.     

Effects of glycosylation on the conformation and dynamics of O-linked glycoproteins: carbon-13 NMR studies of ovine submaxillary mucin

Carbon-13 NMR spectroscopic studies of native and sequentially deglycosylated ovine submaxillary mucin (OSM) have been performed to examine the effects of glycosylation on the conformation and dynamics of the peptide core of O-linked glycoproteins. OSM is a large nonglobular glycoprotein in which nearly one-third of the amino acid residues are Ser and Thr which are glycosylated by the alpha-Neu-NAc(2-6)alpha-GalNAc- disaccharide. The beta-carbon resonances of glycosylated Ser and Thr residues in intact and asialo mucin display considerable chemical shift heterogeneity which, upon the complete removal of carbohydrate, coalesces to single sharp resonances. This chemical shift heterogeneity is due to peptide sequence variability and is proposed to reflect the presence of sequence-dependent conformations of the peptide core. These different conformations are thought to be determined by steric interactions of the GalNAc residue with adjacent peptide residues. The absence of chemical shift heterogeneity in apo mucin is taken to indicate a loss in the peptide-carbohydrate steric interactions, consistent with a more relaxed random coiled structure. On the basis of the 13C relaxation behavior (T1 and NOE) the dynamics of the alpha-carbons appear to be unique to each amino acid type and glycosylation state, with alpha-carbon mobilities decreasing in the order Gly greater than Ala = Ser greater than Thr much greater than monoglycosylated Ser/Thr approximately greater than disaccharide linked Ser/Thr.(ABSTRACT TRUNCATED AT 250 WORDS)     

Partial glycosylation of Asn2181 in human factor V as a cause of molecular and functional heterogeneity. Modulation of glycosylation efficiency by mutagenesis of the consensus sequence for N-linked glycosylation.

Coagulation factor V (FV) circulates in two forms, FV1 and FV2, having slightly different molecular masses and phospholipid-binding properties. The aim was to determine whether this heterogeneity is due to the degree of glycosylation of Asn(2181). FVa1 and FVa2 were isolated and digested with endoglycosidase PNGase F. As judged by Western blotting, the FVa2 light chain contained two N-linked carbohydrates, whereas FVa1 contained three. Wild-type FV and three mutants, Asn(2181)Gln, Ser(2183)Thr, and Ser(2183)Ala, were expressed in COS1 cells, activated by thrombin, and analyzed by Western blotting. Wild-type FVa contained the 71 kDa-74 kDa doublet, whereas the Asn(2181)Gln and Ser(2183)Ala mutants contained only the 71 kDa light chain. In contrast, the Ser(2183)Thr mutant gave a 74 kDa light chain. This demonstrated that the third position in the Asn-X-Ser/Thr consensus affects glycosylation efficiency, Thr being associated with a higher degree of glycosylation than Ser. The Ser(2183)Thr mutant FVa was functionally indistinguishable from plasma-purified FVa1, whereas Asn(2181)Gln and Ser(2183)Ala mutants behaved like FVa2. Thus, the carbohydrate at Asn(2181) impaired the interaction between FVa and the phospholipid membrane, an interpretation consistent with a structural analysis of a three-dimensional model of the C2 domain and the position of a proposed phospholipid-binding site. In conclusion, we show that the FV1-FV2 heterogeneity is caused by differential glycosylation of Asn(2181) related to the presence of a Ser rather than a Thr at the third position in the consensus sequence of glycosylation.     

Role of glycosylation on the ensemble of conformations in the MUC1 immunodominant epitope

MUC1 is a membrane glycoprotein, which in adenocarninomas is overexpressed and exhibits truncated O‐glycosylation. Overexpression and altered glycosylation make MUC1 into a candidate for immunotherapy. Monoclonal antibodies directed against MUC1 frequently bind an immunodominant epitope that contains a single site for O‐glycosylation. Glycosylation with tumor carbohydrate antigens such as the Tn‐antigen (GalNAc‐O‐Ser/Thr) results in antibodies binding with higher affinity. One proposed model to explain the enhanced affinity of antibodies for the glycosylated antigen is that the addition of a carbohydrate alters the conformational properties, favoring a binding‐competent state. The conformational effects associated with Tn glycosylation of the MUC1 antigen was investigated using solution‐state NMR and molecular dynamics. NMR experiments revealed distinct substructures of the glycosylated MUC1 peptides compared with the unglycosylated peptide. Molecular dynamics simulations of the MUC1 glycopeptide and peptide revealed distinguishing differences in their conformational preferences. Furthermore, the glycopeptide displayed a smaller conformational sampling compared with the peptide, suggesting that the glycopeptide sampled a narrower conformational space and is less dynamic. A comparison of the computed ensemble of conformations assuming random distribution, NMR models, and molecular dynamics simulations indicated that the MUC1 glycopeptide and aglycosylated peptide sampled structurally distinctly ensembles and that these ensembles were different from that of the random coil. Together, these data support the hypothesis that that conformational pre‐selection could be an essential feature of these peptides that dictates the binding affinities to MUC1 specific antibodies.     

Ser/Thr Motifs in Transmembrane Proteins: Conservation Patterns and Effects on Local Protein Structure and Dynamics

We combined systematic bioinformatics analyses and molecular dynamics simulations to assess the conservation patterns of Ser and Thr motifs in membrane proteins, and the effect of such motifs on the structure and dynamics of α-helical transmembrane (TM) segments. We find that Ser/Thr motifs are often present in β-barrel TM proteins. At least one Ser/Thr motif is present in almost half of the sequences of α-helical proteins analyzed here. The extensive bioinformatics analyses and inspection of protein structures led to the identification of molecular transporters with noticeable numbers of Ser/Thr motifs within the TM region. Given the energetic penalty for burying multiple Ser/Thr groups in the membrane hydrophobic core, the observation of transporters with multiple membrane-embedded Ser/Thr is intriguing and raises the question of how the presence of multiple Ser/Thr affects protein local structure and dynamics. Molecular dynamics simulations of four different Ser-containing model TM peptides indicate that backbone hydrogen bonding of membrane-buried Ser/Thr hydroxyl groups can significantly change the local structure and dynamics of the helix. Ser groups located close to the membrane interface can hydrogen bond to solvent water instead of protein backbone, leading to an enhanced local solvation of the peptide.     

Molecular basis for the differential sensitivity of rat and human α9α10 nAChRs to α-conotoxin RgIA

J. Neurochem. (2012) 122, 1137-1144. ABSTRACT: The α9α10 nicotinic acetylcholine receptor (nAChR) may be a potential target in pathophysiology of the auditory system, chronic pain, and breast and lung cancers. Alpha-conotoxins, from the predatory marine snail Conus, are potent nicotinic antagonists, some of which are selective for the α9α10 nAChR. Here, we report a two order of magnitude species difference in the potency of α-conotoxin RgIA for the rat versus human α9α10 nAChR. We investigated the molecular mechanism of this difference. Heterologous expression of the rat α9 with the human α10 subunit in Xenopus oocytes resulted in a receptor that was blocked by RgIA with potency similar to that of the rat α9α10 nAChR. Conversely, expression of the human α9 with that of the rat α10 subunit resulted in a receptor that was blocked by RgIA with potency approaching that of the human α9α10 receptor. Systematic substitution of residues found in the human α9 subunit into the homologous position in the rat α9 subunit revealed that a single point mutation, Thr56 to Ile56, primarily accounts for this species difference. Remarkably, although the α9 nAChR subunit has previously been reported to provide the principal (+) binding face for binding of RgIA, Thr56 is located in the (-) complementary binding face.     

Structural effects of divalent calcium cations on the α7 nicotinic acetylcholine receptor: A molecular dynamics simulation study

The α7 subtype of neuronal nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein that is vital to various neurological functions, including modulation of neurotransmitter release. A relatively high concentration of extracellular Ca²⁺ in the neuronal environment is likely to exert substantial structural and functional influence on nAChRs, which may affect their interactions with agonists and antagonists. In this work, we employed atomistic molecular dynamics (MD) simulations to examine the effects of elevated Ca²⁺ on the structure and dynamics of α7 nAChR embedded in a model phospholipid bilayer. Our results suggest that the presence of Ca²⁺ in the α7 nAChR environment results in closure of loop C-in the extracellular ligand-binding domain, a motion normally associated with agonist binding and receptor activation. Elevated Ca²⁺ also alters the conformation of key regions of the receptor, including the inter-helical loops, pore-lining helices and the “gate” residues, and causes partial channel opening in the absence of an agonist, leading to an attendant reduction in the free energy of Ca²⁺ permeation through the pore as elucidated by umbrella sampling simulations. Overall, the structural and permeability changes in α7 nAChR suggest that elevated Ca²⁺ induces a partially activated receptor state that is distinct from both the resting and the agonist-activated states. These results are consistent with the notion that divalent ions can serve as a potentiator of nAChRs, resulting in a higher rate of receptor activation (and subsequent desensitization) in the presence of agonists, with possible implications for diseases involving calcium dysregulation.     

Molecular dynamics simulations on pars intercerebralis major peptide-C (PMP-C) reveal the role of glycosylation and disulfide bonds in its enhanced structural stability and function

Fucosylation of Thr 9 in pars intercerebralis major peptide-C (PMP-C) enhances its structural stability and functional ability as a serine protease inhibitor. In order to understand the role of disulfide bonds and glycosylation on the structure and function of PMP-C, we have carried out multiple explicit solvent molecular dynamics (MD) simulations on fucosylated and non-fucosylated forms of PMP-C, both in the presence and absence of the disulfide bonds. Our simulations revealed that there were no significant structural changes in the native disulfide bonded forms of PMP-C due to fucosylation. On the other hand, the non-fucosylated form of PMP-C without disulfide bonds showed larger deviations from the starting structure than the fucosylated form. However, the structural deviations were restricted to the terminal regions while core β-sheet retained its hydrogen bonded structure even in absence of disulfide bonds as well as fucosylation. Interestingly, fucosylation of disulfide bonded native PMP-C led to a decreased thermal flexibility in the residue stretch 29-32 which is known to interact with the active site of the target proteases. Our analysis revealed that disulfide bonds covalently connect the residue stretch 29-32 to the central β-sheet of PMP-C and using a novel network of side chain interactions and disulfide bonds fucosylation at Thr 9 is altering the flexibility of the stretch 29-32 located at a distal site. Thus, our simulations explain for the first time, how presence of disulfide bonds between conserved cysteines and fucosylation enhance the function of PMP-C as a protease inhibitor.     

Molecular Dynamics Simulations on Pars Intercerebralis Major Peptide-C (PMP-C) Reveal the Role of Glycosylation and Disulfide Bonds in its Enhanced Structural Stability and Function

Abstract

Fucosylation of Thr 9 in pars intercerebralis major peptide-C (PMP-C) enhances its structural stability and functional ability as a serine protease inhibitor. In order to understand the role of disulfide bonds and glycosylation on the structure and function of PMP-C, we have carried out multiple explicit solvent molecular dynamics (MD) simulations on fucosylated and non-fucosylated forms of PMP-C, both in the presence and absence of the disulfide bonds. Our simulations revealed that there were no significant structural changes in the native disulfide bonded forms of PMP-C due to fucosylation. On the other hand, the non-fucosylated form of PMP-C without disulfide bonds showed larger deviations from the starting structure than the fucosylated form. However, the structural deviations were restricted to the terminal regions while core β-sheet retained its hydrogen bonded structure even in absence of disulfide bonds as well as fucosylation. Interestingly, fucosylation of disulfide bonded native PMP-C led to a decreased thermal flexibility in the residue stretch 29–32 which is known to interact with the active site of the target proteases. Our analysis revealed that disulfide bonds covalently connect the residue stretch 29–32 to the central β-sheet of PMP-C and using a novel network of side chain interactions and disulfide bonds fucosylation at Thr 9 is altering the flexibility of the stretch 29–32 located at a distal site. Thus, our simulations explain for the first time, how presence of disulfide bonds between conserved cysteines and fucosylation enhance the function of PMP-C as a protease inhibitor.     

Exploration of the conformational landscape in pregnane X receptor reveals a new binding pocket

Ligand‐regulated pregnane X receptor (PXR), a member of the nuclear receptor superfamily, plays a central role in xenobiotic metabolism. Despite its critical role in drug metabolism, PXR activation can lead to adverse drug‐drug interactions and early stage metabolism of drugs. Activated PXR can induce cancer drug resistance and enhance the onset of malignancy. Since promiscuity in ligand binding makes it difficult to develop competitive inhibitors targeting PXR ligand binding pocket (LBP), it is essential to identify allosteric sites for effective PXR antagonism. Here, molecular dynamics (MD) simulation studies unravelled the existence of two different conformational states, namely “expanded” and “contracted”, in apo PXR ligand binding domain (LBD). Ligand binding events shifted this conformational equilibrium and locked the LBD in a single “ligand‐adaptable” conformational state. Ensemble‐based computational solvent mapping identified a transiently open potential small molecule binding pocket between α5 and α8 helices, named “α8 pocket”, whose opening‐closing mechanism directly correlated with the conformational shift in LBD. A virtual hit identified through structure‐based virtual screening against α8 pocket locks the pocket in its open conformation. MD simulations further revealed that the presence of small molecule at allosteric site disrupts the LBD dynamics and locks the LBD in a “tightly‐contracted” conformation. The molecular details provided here could guide new structural studies to understand PXR activation and antagonism.     

Control of cation permeation through the nicotinic receptor channel

We used molecular dynamics (MD) simulations to explore the transport of single cations through the channel of the muscle nicotinic acetylcholine receptor (nAChR). Four MD simulations of 16 ns were performed at physiological and hyperpolarized membrane potentials, with and without restraints of the structure, but all without bound agonist. With the structure unrestrained and a potential of −100 mV, one cation traversed the channel during a transient period of channel hydration; at −200 mV, the channel was continuously hydrated and two cations traversed the channel. With the structure restrained, however, cations did not traverse the channel at either membrane potential, even though the channel was continuously hydrated. The overall results show that cation selective transport through the nAChR channel is governed by electrostatic interactions to achieve charge selectivity, but ion translocation relies on channel hydration, facilitated by a trans-membrane field, coupled with dynamic fluctuations of the channel structure.     

New insight into the architecture of oxy‐anion pocket in unliganded conformation of GAT domains: A MD‐simulation study

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N‐terminal “glutaminase” (GAT) and C‐terminal “synthetase” domain. The enzyme is identified as a potential target for anti‐cancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H‐bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)‐substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non‐catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site. Proteins 2016; 84:360–373. © 2016 Wiley Periodicals, Inc.     

Multiple binding modes for palmitate to barley lipid transfer protein facilitated by the presence of proline 12

Molecular dynamics simulations have been used to characterise the binding of the fatty acid ligand palmitate in the barley lipid transfer protein 1 (LTP) internal cavity. Two different palmitate binding modes (1 and 2), with similar protein–ligand interaction energies, have been identified using a variety of simulation strategies. These strategies include applying experimental protein–ligand atom–atom distance restraints during the simulation, or protonating the palmitate ligand, or using the vacuum GROMOS 54B7 force‐field parameter set for the ligand during the initial stages of the simulations. In both the binding modes identified the palmitate carboxylate head group hydrogen bonds with main chain amide groups in helix A, residues 4 to 19, of the protein. In binding mode 1 the hydrogen bonds are to Lys 11, Cys 13, and Leu 14 and in binding mode 2 to Thr 15, Tyr 16, Val 17, Ser 24 and also to the OH of Thr 15. In both cases palmitate binding exploits irregularity of the intrahelical hydrogen‐bonding pattern in helix A of barley LTP due to the presence of Pro 12. Simulations of two variants of barley LTP, namely the single mutant Pro12Val and the double mutant Pro12Val Pro70Val, show that Pro 12 is required for persistent palmitate binding in the LTP cavity. Overall, the work identifies key MD simulation approaches for characterizing the details of protein–ligand interactions in complexes where NMR data provide insufficient restraints.     

Enhanced self-association of mucins possessing the T and Tn carbohydrate cancer antigens at the single-molecule level

Mucins are linear O-glycosylated glycoproteins involved in inflammation, cell adhesion, and tumorigenesis. Cancer-associated mucins often possess increased expression of the T (Galβ1,3GalNAcαThr/Ser) and Tn (GalNAcαThr/Ser) cancer antigens, which are diagnostic markers for several cancers, including colon cancer. We have used AFM based single-molecule forced unbinding under near physiological conditions to investigate the self-interactions between porcine submaxillary mucin (PSM) as well as between PSM analogs possessing various carbohydrates including the T- and Tn-antigen. Distributions of unbinding forces and corresponding force loading rates were determined for force loading rates from 0.18 nN/s to 39 nN/s, and processed to yield most probable unbinding forces f* and lifetimes of the interactions. Parameter f* varied in the range 27 to 50 pN at force loading rates of about 2 nN/s among the various mucins. All mucin samples investigated showed self-interaction, but the tendency was greatest for PSM displaying only the Tn-antigen (Tn-PSM) or a mixture of Tn-, T-antigen, and the trisaccharide Fucα1,2Galβ1,3GalNAc (Tri-PSM). Weaker self-interactions were observed for native PSM (Fd-PSM), which consists of a nearly equal mixture of the longer core 1 blood group A tetrasaccharide (GalNAcα1,3(Fucα1,2)Galβ1,3GalNAcαSer/Thr) and Tn-antigen. The data are consistent with the truncated Tn and T glycans enhancing self-interaction of the mucins. These carbohydrate cancer antigens may, thus, play an active role in the disease by constitutively activating mucin and mucin-type receptors by self-association on cells.     

Structural requirements for additional N-linked carbohydrate on recombinant human erythropoietin

N-Linked glycosylation is a post-translational event whereby carbohydrates are added to secreted proteins at the consensus sequence Asn-Xaa-Ser/Thr, where Xaa is any amino acid except proline. Some consensus sequences in secreted proteins are not glycosylated, indicating that consensus sequences are necessary but not sufficient for glycosylation. In order to understand the structural rules for N-linked glycosylation, we introduced N-linked consensus sequences by site-directed mutagenesis into the polypeptide chain of the recombinant human erythropoietin molecule. Some regions of the polypeptide chain supported N-linked glycosylation more effectively than others. N-Linked glycosylation was inhibited by an adjacent proline suggesting that sequence context of a consensus sequence could affect glycosylation. One N-linked consensus sequence (Asn123-Thr125) introduced into a position close to the existing O-glycosylation site (Ser126) had an additional O-linked carbohydrate chain and not an additional N-linked carbohydrate chain suggesting that structural requirements in this region favored O-glycosylation over N-glycosylation. The presence of a consensus sequence on the protein surface of the folded molecule did not appear to be a prerequisite for oligosaccharide addition. However, it was noted that recombinant human erythropoietin analogs that were hyperglycosylated at sites that were normally buried had altered protein structures. This suggests that carbohydrate addition precedes polypeptide folding.     

Molecular dynamics simulations of the Bcl-2 protein to predict the structure of its unordered flexible loop domain

B-cell lymphoma (Bcl-2) protein is an anti-apoptotic member of the Bcl-2 family. It is functionally demarcated into four Bcl-2 homology (BH) domains: BH1, BH2, BH3, BH4, one flexible loop domain (FLD), a transmembrane domain (TM), and an X domain. Bcl-2’s BH domains have clearly been elucidated from a structural perspective, whereas the conformation of FLD has not yet been predicted, despite its important role in regulating apoptosis through its interactions with JNK-1, PKC, PP2A phosphatase, caspase 3, MAP kinase, ubiquitin, PS1, and FKBP38. Many important residues that regulate Bcl-2 anti-apoptotic activity are present in this domain, for example Asp34, Thr56, Thr69, Ser70, Thr74, and Ser87. The structural elucidation of the FLD would likely help in attempts to accurately predict the effect of mutating these residues on the overall structure of the protein and the interactions of other proteins in this domain. Therefore, we have generated an increased quality model of the Bcl-2 protein including the FLD through modeling. Further, molecular dynamics (MD) simulations were used for FLD optimization, to predict the flexibility, and to determine the stability of the folded FLD. In addition, essential dynamics (ED) was used to predict the collective motions and the essential subspace relevant to Bcl-2 protein function. The predicted average structure and ensemble of MD-simulated structures were submitted to the Protein Model Database (PMDB), and the Bcl-2 structures obtained exhibited enhanced quality. This study should help to elucidate the structural basis for Bcl-2 anti-apoptotic activity regulation through its binding to other proteins via the FLD.     

Molecular Dynamics Simulations of L-Lactate Dehydrogenase: Conformation of a Mobile Loop and Influence of the Tetrameric Protein Environment

L-lactate dehydrogenase (LDH) catalyses the interconversion of pyruvate and L-lactate in the presence of the coenzyme NADH. Molecular dynamics (MD) simulations have been performed for LDH complexed with NADH and the pyruvate-analogue inhibitor oxamate with the aim of characterising important influences on maintaining the geometry and hydrogen bond network of the active site. Two features in particular were found to dominate. First, the tetrameric protein environment is found to play a significant role in maintaining the active-site geometry. Simulations of the monomer alone reproduce the crystallographic structure poorly, and at least part of the neighboring subunit is necessary to prevent water penetration into the active site and to provide rigidity to the α1G-α2G helix immediately adjacent to the active site. These results offer one explanation for the observation that the monomer is not biologically active. Second, the conformation of Arg109 (part of the mobile loop which closes over the active site) is shown to play a key role in maintaining the active-site geometry. In some simulations, a torsional rotation in the side chain of Arg109 results in the breaking of crystallographic hydrogen bonds which are important for polarising the carbonyl bond of the substrate. This conformational change appears to be a trigger for the opening of the mobile loop. Long-range nonbonded interactions are found to be influential in maintaining the proper crystallographic conformation of Arg109. Thus, we conclude that to adequately model LDH, at least part of the neighboring subunit must be included in the MD simulations and nonbonded interactions must be properly represented to ensure that Arg109 remains in the crystallographic conformation. Out of a set of simulation protocols tested here, one meets both these criteria and will be used for the generation of starting structures for future reaction-mechanism calculations.     

Structure prediction and R115866 binding study of human CYP26A1: homology modelling,fold recognition,molecular docking and MD simulations

Two three-dimensional (3D) models of human cytochrome P450 26A1 (CYP26A1) were constructed using the programs Modeller and Sybyl-GeneFold, respectively. After refinement by molecular mechanics and molecular dynamics (MD) simulations, the two models were validated by structure analysis-validation online server. Subsequently, a flexible docking study was performed on the model constructed by GeneFold with the potent and specific inhibitor R115866 to examine the enzyme–inhibitor interactions. From the docking results, we can see R115866 interacts with amino acid residues at the active site by multiple hydrophobic interactions including the side chains of His111, Trp112, Ser115, Val116, Leu125, Ser126, Leu221, Phe222, Glu296, Phe299, Gly300, Glu303, Thr304, Pro371 and the cofactor heme. Trp112 and Thr304 form hydrogen bonds with R115866 and play important roles in stabilising the complex. This constructed CYP26A1 model may provide an opportunity to understand the action mode of the enzyme and could be useful in designing novel retinoic acid metabolism blocking agents (RAMBAs).