Insights from molecular dynamics simulations into pH-dependent enantioselective hydrolysis of ibuprofen esters by Candida rugosa lipase

An interesting observation was found during our continued studies on the hydrolysis of ibuprofen esters by Candida rugosa lipase (CRL). An important role is played by pH in the stereospecific hydrolysis of these esters. The flap region of CRL plays a significant role in the access of the substrate to the active site of the enzyme. At pH 5.6, 48% of the methyl ester and 5% of the butyl ester of ibuprofen were hydrolysed in 5.5 h, whereas at pH 7.2, 9% of methyl ester and 45% of the butyl ester of ibuprofen was hydrolysed in a identical reaction time using CRL. This lead us to assume that CRL prefers the methyl ester of ibuprofen as a substrate at an acidic pH and the butyl ester of ibuprofen at a neutral pH. Therefore, in order to understand the role of pH in the substrate selection by CRL for the esters of ibuprofen we used the crystallographic coordinates of the open form of the CRL (1CRL) for molecular dynamics (MD) simulations under acidic and neutral conditions for 2 ns using GROMACS. The final structures obtained after simulation in acidic and neutral conditions were compared with the energy-minimized structure, and the root-mean-square deviations (r.m.s.ds) were calculated. The r.m.s.d. of the CRL flap at neutral pH was found to be greater than that of the CRL flap at acidic pH. The extent to which the flap opens at neutral pH allowed the bulkier substrate, the butyl ester of ibuprofen, to diffuse into the active site and provides the best enzyme-substrate fit for this specific substrate. At acidic pH there is a decreased opening of the flap thereby accommodating a more compact substrate, namely the methyl ester of ibuprofen. Thus, simulation experiments using MD provide reasonable insight for the pH-dependent substrate selectivity of CRL in aqueous environments.     

The role of tyrosine 71 in modulating the flap conformations of BACE1

β-Site amyloid precursor protein cleaving enzyme 1 (BACE1) is a potential target for treating Alzheimer's disease. BACE1's binding site is partially covered by a flexible loop on its N-terminal domain, known as the "flap," which has been found in several conformations in crystal structures of BACE1 and other aspartyl proteases. The side chain of the invariant residue Tyr71 on the flap adopts several rotameric orientations, leading to our hypothesis that the orientation of this residue dictates the movement and conformations available to the flap. We investigated this hypothesis by performing 220 ns of molecular dynamics simulations of bound and unbound wild-type BACE1 as well as the unbound Y71A mutant. Our findings indicate that the flap exhibits various degrees of mobility and adopts different conformations depending on the Tyr71 orientation. Surprisingly, the "self-inhibited" form is stable in our simulations, making it a reasonable target for drug design. The alanine mutant, lacking a large side chain at position 71, displays significant differences in flap dynamics from wild type, freely sampling very open and closed conformations. Our simulations show that Tyr71, in addition to its previously determined functions in catalysis and substrate binding, has the important role of modulating flap conformations in BACE1.     

Structure and dynamics of Candida rugosa lipase: the role of organic solvent

The effect of organic solvent on the structure and dynamics of proteins was investigated by multiple molecular dynamics simulations (1 ns each) of Candida rugosa lipase in water and in carbon tetrachloride. The choice of solvent had only a minor structural effect. For both solvents the open and the closed conformation of the lipase were near to their experimental X-ray structures (C rms deviation 1–1.3 Å). However, the solvents had a highly specific effect on the flexibility of solvent-exposed side chains: polar side chains were more flexible in water, but less flexible in organic solvent. In contrast, hydrophobic residues were more flexible in organic solvent, but less flexible in water. As a major effect solvent changed the dynamics of the lid, a mobile element involved in activation of the lipase, which fluctuated as a rigid body about its average position. While in water the deviations were about 1.6 Å, organic solvent reduced flexibility to 0.9 Å. This increase rigidity was caused by two salt bridges (Lys85–Asp284, Lys75–Asp79) and a stable hydrogen bond (Lys75–Asn 292) in organic solvent. Thus, organic solvents stabilize the lid but render the side chains in the hydrophobic substrate-binding site more mobile. Figure Superimposition of open (black, PDB entry 1CRL) and closed (gray, PDB entry 1TRH) conformers of C. rugosa lipase. The mobile lid is indicatedThis revised version was published online in October 2004 with corrections to the Graphical Abstract.     

Investigation of flap flexibility of β-secretase using molecular dynamic simulations

Flap motif and its dynamics were extensively reported in aspartate proteases, e.g. HIV proteases and plasmepsins. Herein, we report the first account of flap dynamics amongst different conformations of β-secretase using molecular dynamics simulation. Various parameters were proposed and a selected few were picked which could appropriately describe the flap motion. Three systems were studied, namely Free (BACE_Free) and two ligand-bound conformations, which belonged to space groups P6₁22 (BACE_Bound1) and C222₁ (BACE_Bound2), respectively and four parameters (distance between the flaps tip residue, Thr72 and Ser325, d₁; dihedral angle, ? (Thr72-Asp32-Asp228-Ser325); TriCα angles, θ₁ (Thr72-Asp32-Ser325), and θ₂ (Thr72-Asp228-Ser325)) were proposed to understand the change in dynamics of flap domain and the extent of flap opening and closing. Analysis of, θ₂, d₁, θ₁ and ? confirmed that the BACE_Free adopted semi-open, open and closed conformations with slight twisting during flap opening. However, BACE_Bound1 (P6122) showed an adaptation to open conformation due to lack of hydrogen bond interaction between the ligand and flap tip residue. A slight flap twisting, ? (lateral twisting) was observed for BACE_Bound1 during flap opening which correlates with the opening of BACE_Free. Contradictory to the BACE_Bound1, the BACE_Bound2 locked the flap in a closed conformation throughout the simulation due to formation of a stable hydrogen bond interaction between the flap tip residue and ligand. Analyses of all three systems highlight that d₁, θ₂ and ? can be precisely used to describe the extent of flap opening and closing concurrently with snapshots along the molecular dynamics trajectory across several conformations of β-secretase.     

A molecular dynamics study comparing a wild-type with a multiple drug resistant HIV protease: differences in flap and aspartate 25 cavity dimensions

HIV proteases can develop resistance to therapeutic drugs by mutating specific residues, but still maintain activity with their natural substrates. To gain insight into why mutations confer such resistance, long ( approximately 70 ns) Molecular Dynamics simulations in explicit solvent were performed on a multiple drug resistant (MDR) mutant (with Asn25 in the crystal structure mutated in silico back to the catalytically active Asp25) and a wild type (WT) protease. HIV proteases are homodimers, with characteristic flap tips whose conformations and dynamics are known to be important influences of ligand binding to the aspartates that form the catalytic center. The WT protease undergoes a transition between 25 and 35 ns that is absent in the MDR protease. The origin of this distinction is investigated using principal component analysis, and is related to differences in motion mainly in the flap region of each monomer. Trajectory analysis suggests that the WT transition arises from a concerted motion of the flap tip distances to their catalytic aspartate residues, and the distance between the two flap tips. These distances form a triangle that in the WT expands the active site from an initial (semi-open) form to an open form, in a correlated manner. In contrast, the MDR protease remains in a more closed configuration, with uncorrelated fluctuations in the distances defining the triangle. This contrasting behavior suggests that the MDR mutant achieves its resistance to drugs by making its active site less accessible to inhibitors. The migration of water to the active site aspartates is monitored. Water molecules move in and out of the active site and individual waters hydrogen bond to both aspartate carboxylate oxygens, with residence times in the ns time regime.     

Interchain hydrophobic clustering promotes rigidity in HIV-1 protease flap dynamics: new insights from Molecular Dynamics

The dynamics of HIV-1 protease (HIV-pr), a drug target for HIV infection, has been studied extensively by both computational and experimental methods. The flap dynamics of HIV-pr is considered to be more important for better ligand binding and enzymatic actions. Moreover, it has been demonstrated that the drug-induced mutations can change the flap dynamics of HIV-pr affecting the binding affinity of the ligands. Therefore, detailed understanding of flap dynamics is essential for designing better inhibitors. Previous computational investigations observed significant variation in the flap opening in nanosecond time scale indicating that the dynamics is highly sensitive to the simulation protocols. To understand the sensitivity of the flap dynamics on the force field and simulation protocol, molecular dynamics simulations of HIV-pr have been performed with two different AMBER force fields, ff99 and ff02. Two different trajectories (20?ns each) were obtained using the ff99 and ff02 force field. The results showed polarizable force field (ff02) make the flap tighter than the nonpolarizable force field (ff99). Some polar interactions and hydrogen bonds involving flap residues were found to be stronger with ff02 force field. The formation of interchain hydrophobic cluster (between flap tip of one chain and active site wall of another chain) was found to be dominant in the semi-open structures obtained from the simulations irrespective of the force field. It is proposed that an inhibitor, which will promote this interchain hydrophobic clustering, may make the flaps more rigid, and presumably the effect of mutation would be small on ligand binding.     

A poke in the eye: inhibiting HIV-1 protease through its flap-recognition pocket

A novel mechanism of inhibiting HIV-1 protease (HIVp) is presented. Using computational solvent mapping to identify complementary interactions and the Multiple Protein Structure method to incorporate protein flexibility, we generated a receptor-based pharmacophore model of the flexible flap region of the semiopen, apo state of HIVp. Complementary interactions were consistently observed at the base of the flap, only within a cleft with a specific structural role. In the closed, bound state of HIVp, each flap tip docks against the opposite monomer, occupying this cleft. This flap-recognition site is filled by the protein and cannot be identified using traditional approaches based on bound, closed structures. Virtual screening and dynamics simulations show how small molecules can be identified to complement this cleft. Subsequent experimental testing confirms inhibitory activity of this new class of inhibitor. This may be the first new inhibitor class for HIVp since dimerization inhibitors were introduced 17 years ago.     

Dynamic flaps in HIV-1 protease adopt unique ordering at different stages in the catalytic cycle

The flexibility of HIV-1 protease flaps is known to be essential for the enzymatic activity. Here we attempt to capture a multitude of conformations of the free and substrate-bound HIV-1 protease that differ drastically in their flap arrangements. The substrate binding process suggests the opening of active site gate in conjunction with a reversal of flap tip ordering, from the native semiopen state. The reversed-flap, open-gated enzyme readily transforms to a closed conformation after proper placement of the substrate into the binding cleft. After substrate processing, the closed state protease which possessed opposite flap ordering relative to the semiopen state, encounters another flap reversal via a second open conformation that facilitates the evolution of native semiopen state of correct flap ordering. The complicated transitional pathway, comprising of many high and low energy states, is explored by combining standard and activated molecular dynamics (MD) simulation techniques. The study not only complements the existing findings from X-ray, NMR, EPR, and MD studies but also provides a wealth of detailed information that could help the structure-based drug design process.     

Molecular docking based virtual screening of natural compounds as potential BACE1 inhibitors: 3D QSAR pharmacophore mapping and molecular dynamics analysis

Beta-site APP cleaving enzyme1 (BACE1) catalyzes the rate determining step in the generation of Aβ peptide and is widely considered as a potential therapeutic drug target for Alzheimer’s disease (AD). Active site of BACE1 contains catalytic aspartic (Asp) dyad and flap. Asp dyad cleaves the substrate amyloid precursor protein with the help of flap. Currently, there are no marketed drugs available against BACE1 and existing inhibitors are mostly pseudopeptide or synthetic derivatives. There is a need to search for a potent inhibitor with natural scaffold interacting with flap and Asp dyad. This study screens the natural database InterBioScreen, followed by three-dimensional (3D) QSAR pharmacophore modeling, mapping, in silico ADME/T predictions to find the potential BACE1 inhibitors. Further, molecular dynamics of selected inhibitors were performed to observe the dynamic structure of protein after ligand binding. All conformations and the residues of binding region were stable but the flap adopted a closed conformation after binding with the ligand. Bond oligosaccharide interacted with the flap as well as catalytic dyad via hydrogen bond throughout the simulation. This led to stabilize the flap in closed conformation and restricted the entry of substrate. Carbohydrates have been earlier used in the treatment of AD because of their low toxicity, high efficiency, good biocompatibility, and easy permeability through the blood–brain barrier. Our finding will be helpful in identify the potential leads to design novel BACE1 inhibitors for AD therapy.     

Solvent-induced lid opening in lipases: A molecular dynamics study

In most lipases, a mobile lid covers the substrate binding site. In this closed structure, the lipase is assumed to be inactive. Upon activation of the lipase by contact with a hydrophobic solvent or at a hydrophobic interface, the lid opens. In its open structure, the substrate binding site is accessible and the lipase is active. The molecular mechanism of this interfacial activation was studied for three lipases (from Candida rugosa, Rhizomucor miehei, and Thermomyces lanuginosa) by multiple molecular dynamics simulations for 25 ns without applying restraints or external forces. As initial structures of the simulations, the closed and open structures of the lipases were used. Both the closed and the open structure were simulated in water and in an organic solvent, toluene. In simulations of the closed lipases in water, no conformational transition was observed. However, in three independent simulations of the closed lipases in toluene the lid gradually opened. Thus, pathways of the conformational transitions were investigated and possible kinetic bottlenecks were suggested. The open structures in toluene were stable, but in water the lid of all three lipases moved towards the closed structure and partially unfolded. Thus, in all three lipases opening and closing was driven by the solvent and independent of a bound substrate molecule.     

Rapid structural fluctuations of the free HIV protease flaps in solution: Relationship to crystal structures and comparison with predictions of dynamics calculations

Crystal structures have shown that the HIV-1 protease flaps, domains that control access to the active site, are closed when the active site is occupied by a ligand. Although flap structures ranging from closed to semi-open are observed in the free protease, crystal structures reveal that even the semi-open flaps block access to the active site, indicating that the flaps are mobile in solution. The goals of this paper are to characterize the secondary structure and fast (sub-ns) dynamics of the flaps of the free protease in solution, to relate these results to X-ray structures and to compare them with predictions of dynamics calculations. To this end we have obtained nearly complete backbone and many sidechain signal assignments of a fully active free-protease construct that is stabilized against autoproteolysis by three point mutations. The secondary structure of this protein was characterized using the chemical shift index, measurements of (3h)J(NC') couplings across hydrogen bonds, and NOESY connectivities. Analysis of these measurements indicates that the protease secondary structure becomes irregular near the flap tips, residues 49-53. Model-free analysis of (15)N relaxation parameters, T(1), T(2) (T(1rho)) and (15)N-[(1)H] NOE, shows that residues in the flap tips are flexible on the sub-ns time scale, in contrast with previous observations on the inhibitor-bound protease. These results are compared with theoretical predictions of flap dynamics and the possible biological significance of the sub-ns time scale dynamics of the flap tips is discussed.     

The molecular basis of CRL4DDB2/CSA ubiquitin ligase architecture, targeting, and activation

The DDB1-CUL4-RBX1 (CRL4) ubiquitin ligase family regulates a diverse set of cellular pathways through dedicated substrate receptors (DCAFs). The DCAF DDB2 detects UV-induced pyrimidine dimers in the genome and facilitates nucleotide excision repair. We provide the molecular basis for DDB2 receptor-mediated cyclobutane pyrimidine dimer recognition in chromatin. The structures of the fully assembled DDB1-DDB2-CUL4A/B-RBX1 (CRL4(DDB2)) ligases reveal that the mobility of the ligase arm creates a defined ubiquitination zone around the damage, which precludes direct ligase activation by DNA lesions. Instead, the COP9 signalosome (CSN) mediates the CRL4(DDB2) inhibition in a CSN5 independent, nonenzymatic, fashion. In turn, CSN inhibition is relieved upon DNA damage binding to the DDB2 module within CSN-CRL4(DDB2). The Cockayne syndrome A DCAF complex crystal structure shows that CRL4(DCAF(WD40)) ligases share common architectural features. Our data support a general mechanism of ligase activation, which is induced by CSN displacement from CRL4(DCAF) on substrate binding to the DCAF.     

Closing of the flaps of HIV-1 protease induced by substrate binding: a model of a flap closing mechanism in retroviral aspartic proteases

The active site of aspartic proteases is covered by one or more flaps, which control access to the active site and play a significant role in the binding of the substrate. An extensive conformational change of the flaps takes place upon binding of substrate to the active site. A long molecular dynamics simulation was performed on the complex consisting of a peptide (CA-p2) from a natural substrate cleavage site of the gag/pol polyprotein placed in the active site of HIV-1 protease (PR) with an open flap conformation. During the simulation, the substrate induced the closing of the flaps into the closed conformation in an asymmetrical way through a hydrophobic intermediate state cluster. The nature of the residues of HIV-1 PR identified to be important in the flap closing mechanism is conserved across known structures of retroviral aspartic proteases family. The flap closing mechanism described in HIV-1 PR is proposed to be a general model for flap closing in retroviral aspartic proteases.     

Computational studies of tryptophanyl-tRNA synthetase: activation of ATP by induced-fit

Catalysis of amino acid activation by Bacillus stearothermophilus tryptophanyl-tRNA synthetase involves three allosteric states: (1) Open; (2) closed pre-transition state (PreTS); and (3) closed products (Product). The interconversions of these states entail significant domain motions driven by ligand binding. We explore the application of molecular dynamics simulations to investigate ligand-linked conformational stability changes associated with this catalytic cycle. Multiple molecular dynamics trajectories (5 ns) for 11 distinct liganded and unliganded monomer configurations show that the PreTS conformation is unstable in the absence of ATP, reverting within approximately 600 ps nearly to the Open conformation. In contrast, Open and Product state trajectories were stable, even without ligands, confirming the previous suggestion that catalysis entails destabilization of the protein conformation, driven by ATP-binding energies developed as the PreTS state assembles during induced-fit. The simulations suggest novel mechanistic details associated with both induced-fit (Open-PreTS) and catalysis (PreTS-Product). Notably, Mg2+ -ATP interactions are coupled to interactions between ATP and active-site lysine side-chains via mechanisms that cannot be captured under the molecular mechanics approximations, and which therefore require restraining potentials for stable simulation. Simulations of Mg2+. ATP-bound PreTS complexes with restraining potentials and with a virtual K111A mutant confirm that these coupling interactions are necessary to sustain the PreTS conformation and, in turn, provide a new model for how the PreTS conformation activates ATP for catalysis. These results emphasize the central role of the PreTS state as a high-energy intermediate structure along the catalytic pathway and suggest that Mg2+ and the KMSKS loop function cooperatively during catalysis.     

Molecular dynamics simulations of the conformational changes of the glutamate receptor ligand-binding core in the presence of glutamate and kainate

Excitatory synaptic transmission is mediated by ionotropic glutamate receptors (iGluRs) through the induced transient opening of transmembrane ion channels. The three-dimensional structure of the extracellular ligand-binding core of iGluRs shares the overall features of bacterial periplasmic binding proteins (PBPs). In both families of proteins, the ligand-binding site is arranged in two domains separated by a cleft and connected by two peptide stretches. PBPs undergo a typical hinge motion of the two domains associated with ligand binding that leads to a conformational change from an open to a closed form. The common architecture suggests a similar closing mechanism in the ligand-binding core of iGluRs induced by the binding of specific agonists. Starting from the experimentally determined kainate-bound closed form of the S1S2 GluR2 construct, we have studied by means of molecular dynamics simulations the opening motion of the ligand-binding core in the presence and in the absence of both glutamate and kainate. Our results suggest that the opening/closing interdomain hinge motions are coupled to conformational changes in the insertion region of the transmembrane segments. These changes are triggered by the interaction of the agonists with the essential Glu 209 residue. A plausible mechanism for the coupling of agonist binding to channel gating is discussed.     

A Mesorhizobium lipopolysaccharide (LPS) specific lectin (CRL) from the roots of nodulating host plant, Cicer arietinum

A 30 kDa rabbit erythrocyte agglutinating glycoprotein isolated and characterized from the roots of Cicer arietinum and designated as cicer root lectin (CRL). Hemagglutination activity of CRL is strongly inhibited by cell surface LPS of nodulating cicer specific Rhizobium. CRL agglutinates mesorhizobial cells and not Escherichia coli or yeast cells. It binds to immobilized LPS of cicer specific Rhizobium only. The primary structure of CRL as predicted by peptide mass fingerprinting by MALDI-TOF (matrix-assisted laser desorption/ionization-time of flight) indicated ∼54% amino acid sequence homology with C. arietinum seedling lectin (Accession no. gi/3204123) and ∼26% with C. arietinum (Accession no. gi/110611256), and Pisum sativum (Accession nos. gi/230612, gi/6729956, gi/126148) lectins. These suggested CRL to be a member of vegetative tissue lectin. Circular dichroism analysis indicated that the secondary structure of CRL consists of 48% β-sheets, 26% random coils, and 11% α-helix. CRL has six isoforms of closely associated molecular mass with differential acidic pI of 5.30, 5.20, 5.15, 5.05, 5.00, 4.80. Identity of these isoforms was confirmed from their binding with cicer specific Rhizobium LPS. All the isoforms of CRL are differentially glycosylated as identified by deglycosylation and monosaccharide analysis using high-performance anion exchange chromatography with pulsed amperometric detection (HPAEC-PAD). All these results suggest that unlike other plant lectins CRL is a LPS-binding lectin.     

CRL4在增殖的肝细胞系中上调乙型肝炎病毒抗原的分泌表达

乙型肝炎病毒（hepatitis B virus,HBV）编码的X蛋白（hepatitis B virus X protein,HBx）对HBV感染的起始和维持至关重要。HBx可能作为病毒来源的接头分子,介导Cullin-RING E3泛素连接酶4（cullin-RING ubiquitin E3 ligase 4,CRL4）复合物对染色体外DNA限制因子SMC5/6的降解。最近研究发现CRL4接头分子DNA损伤结合蛋白1（DNA damage-binding protein 1,DDB1）可不依赖与HBx的相互作用而直接上调病毒的表达和复制。本研究基于HBx基因删除（X-null）的HBV重组cccDNA（recombinant covalently closed circular DNA,rcccDNA）模型系统,在多种体外培养肝细胞系中证实上述发现。有意思的是,CRL4刺激rcccDNA^X-null转染细胞抗原分泌表达的效应能被血清饥饿实验抵消。应用尾静脉高压注射小鼠模型,同样发现CRL4并不上调rcccDNA^X-null在非增殖小鼠肝脏细胞中的表达。以上结果提示,细胞增殖特征与CRL4不依赖HBx上调病毒抗原表达的效应密切相关,有助于HBx生物学意义的准确分析和理解。     

Characterization of the conformational behavior of peptide Contryphan Vn: a theoretical study

In this work we report the study of a peptide, the Contryphan Vn produced by Conus ventricosus, a vermivorous cone snail living in the temperate Mediterranean sea. This cyclic peptide of nine residues is a ring closed by a Cys-Cys (Cys: cysteine) disulfide bond containing two proline (Pro) residues and two tryptophans (Trp), one of them being a D-Trp. We present a statistical mechanical characterization of the peptide, simulated in water for about 200 ns with classical molecular dynamics (MD). In recent years there has been a growing interest in the study of the mechanics and dynamics of biological molecules, and in particular for proteins and peptides, about the relationship between collective motions and the active conformations which exert the biological function. To this aim we used the essential dynamics analysis on the MD trajectory and extracted, from the total fluctuations of the molecule, the dominant dynamical modes responsible of the principal conformational transitions. The Contryphan Vn small size allowed us to investigate in details the all-atoms dynamics and the corresponding thermodynamics in conformational space defined by the most significant intramolecular motions.     

Induced fit in liver X receptor beta: a molecular dynamics-based investigation

Ligand induced fit phenomenon occurring at the ligand binding domain of the liver X receptor beta (LXRbeta) was investigated by means of molecular dynamics. Reliability of a 4-ns trajectory was tested from two distinct LXRbeta crystal complexes 1PQ6B/GW and 1PQ9B/T09 characterized by an open and a closed state of the pocket, respectively. Crossed complexes 1PQ6B/T09 and 1PQ9B/GW were then submitted to the same molecular dynamic conditions, which were able to recover LXRbeta conformations similar to the original crystallography data. Analysis of "open to closed" and "closed to open" conformational transitions pointed out the dynamic role of critical residues lining the ligand binding pocket involved in the local remodeling upon ligand binding (e.g., Phe271, Phe329, Phe340, Arg319, Glu281). Altogether, the present study indicates that the molecular dynamic protocol is a consistent approach for managing LXRbeta-related induced fit process. This protocol could therefore be used for refining ligand docking solutions of a structure-based design strategy.