Adaptable Small Ligand of CYP1 Enzymes for Use in Understanding the Structural Features Determining Isoform Selectivity.

Although several families of compounds have been identified as scaffolds for inhibitors of the CYP1 family, the isoform selectivity determining structural features has not been fully clarified at the molecular interaction level. We studied the CYP1 isoform selectivity for stilbenoid inhibitors using integrated induced fit docking and molecular dynamics simulations. The hydrophobic interactions with the specific phenylalanine residues in the F helix are correlated with inhibitory potency in the CYP1 family. Through this study, we found that the adaptable, small, and semirigid ligand is a promising starting point for the development of isoform-selective inhibitors and investigation of selectivity-determining features.

Energy transduction and alternating access of the mammalian ABC transporter P-glycoprotein.

ATP binding cassette (ABC) transporters of the exporter class harness the energy of ATP hydrolysis in the nucleotide-binding domains (NBDs) to power the energetically uphill efflux of substrates by a dedicated transmembrane domain (TMD). Although numerous investigations have described the mechanism of ATP hydrolysis and defined the architecture of ABC exporters, a detailed structural dynamic understanding of the transduction of ATP energy to the work of substrate translocation remains elusive. Here we used double electron-electron resonance and molecular dynamics simulations to describe the ATP- and substrate-coupled conformational cycle of the mouse ABC efflux transporter P-glycoprotein (Pgp; also known as ABCB1), which has a central role in the clearance of xenobiotics and in cancer resistance to chemotherapy. Pairs of spin labels were introduced at residues selected to track the putative inward-facing to outward-facing transition. Our findings illuminate how ATP energy is harnessed in the NBDs in a two-stroke cycle and elucidate the consequent conformational motion that reconfigures the TMD, two critical aspects of Pgp transport mechanism. Along with a fully atomistic model of the outward-facing conformation in membranes, the insight into Pgp conformational dynamics harmonizes mechanistic and structural data into a novel perspective on ATP-coupled transport and reveals mechanistic divergence within the efflux class of ABC transporters.

Mass spectrometry-based cross-linking study shows that the Psb28 protein binds to cytochrome b559 in Photosystem II.

Photosystem II (PSII), a large pigment protein complex, undergoes rapid turnover under natural conditions. During assembly of PSII, oxidative damage to vulnerable assembly intermediate complexes must be prevented. Psb28, the only cytoplasmic extrinsic protein in PSII, protects the RC47 assembly intermediate of PSII and assists its efficient conversion into functional PSII. Its role is particularly important under stress conditions when PSII damage occurs frequently. Psb28 is not found, however, in any PSII crystal structure, and its structural location has remained unknown. In this study, we used chemical cross-linking combined with mass spectrometry to capture the transient interaction of Psb28 with PSII. We detected three cross-links between Psb28 and the α- and ß-subunits of cytochrome b559, an essential component of the PSII reaction-center complex. These distance restraints enable us to position Psb28 on the cytosolic surface of PSII directly above cytochrome b559, in close proximity to the QB site. Protein-protein docking results also support Psb28 binding in this region. Determination of the Psb28 binding site and other biochemical evidence allow us to propose a mechanism by which Psb28 exerts its protective effect on the RC47 intermediate. This study also shows that isotope-encoded cross-linking with the "mass tags" selection criteria allows confident identification of more cross-linked peptides in PSII than has been previously reported. This approach thus holds promise to identify other transient protein-protein interactions in membrane protein complexes.

QSAR modeling to design selective histone deacetylase 8 (HDAC8) inhibitors.

HDAC8 inhibitors have become an attractive treatment for cancer. This study aimed to facilitate the identification of potential chemical scaffolds for the selective inhibition of histone deacetylase 8 (HDAC8) using in silico approaches. Non-linear QSAR classification and regression models of HDAC8 inhibitors were developed with support vector machine. Mean impact value-based sequential forward feature selection and grid search strategy were used for molecular descriptor selection and parameter optimization, respectively. The generated QSAR models were validated by leave-one-out cross validation and an external test set. The best QSAR classification model yielded 84 % of accuracy on the external test prediction and Matthews correlation coefficient is 0.69. The best QSAR regression model showed low root-mean-square error (0.63) and high squared correlation coefficient (0.53) for the test set. The validated QSAR models together with various drug-like properties, molecular docking and molecular dynamics simulation were sequentially used as a multi-step query in chemical database virtual screening. Finally, two hit compounds were discovered as new structural scaffolds which can be used for further in vitro and in vivo activity analyses. The strategy used in this study could be a promising computational strategy which can be utilized for other target drug design.

Binding Site Recognition and Docking Dynamics of a Single Electron Transport Protein: Cytochrome c2.

Small diffusible redox proteins facilitate electron transfer in respiration and photosynthesis by alternately binding to their redox partners and integral membrane proteins and exchanging electrons. Diffusive search, recognition, binding, and unbinding of these proteins often amount to kinetic bottlenecks in cellular energy conversion, but despite the availability of structures and intense study, the physical mechanisms controlling redox partner interactions remain largely unknown. The present molecular dynamics study provides an all-atom description of the cytochrome c2-docked bc1 complex in Rhodobacter sphaeroides in terms of an ensemble of favorable docking conformations and reveals an intricate series of conformational changes that allow cytochrome c2 to recognize the bc1 complex and bind or unbind in a redox state-dependent manner. In particular, the role of electron transfer in triggering a molecular switch and in altering water-mediated interface mobility, thereby strengthening and weakening complex formation, is described. The results resolve long-standing discrepancies between structural and functional data.

The cellular membrane as a mediator for small molecule interaction with membrane proteins.

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

New insights in the activation of human cholesterol esterase to design potent anti-cholesterol drugs.

Primary hypercholesterolemia is the root cause for major health issues like coronary heart disease and atherosclerosis. Regulating plasma cholesterol level, which is the product of biosynthesis as well as dietary intake, has become one of the major therapeutic strategies to effectively control these diseases. Human cholesterol esterase (hCEase) is an interesting target involved in the regulation of plasma cholesterol level and thus inhibition of this enzyme is highly effective in the treatment of hypercholesterolemia. This study was designed to understand the activation mechanism that enables the enzyme to accommodate long chain fatty acids and to identify the structural elements for the successful catalysis. Primarily the activation efficiencies of three different bile salts were studied and compared using molecular dynamics simulations. Based on the conformations of major surface loops, hydrogen bond interactions, and distance analyses, taurocholate was concluded as the preferred activator of the enzyme. Furthermore, the importance of two bile salt binding sites (proximal and remote) and the crucial role of 7α-OH group of the bile salts in the activation of hCEase was examined and evidenced. The results of our study explain the structural insights of the activation mechanism and show the key features of the bile salts responsible for the enzyme activation which are very useful in hypolipidemic drug designing strategies.

Multi-conformation dynamic pharmacophore modeling of the peroxisome proliferator-activated receptor γ for the discovery of novel agonists.

Activation of the peroxisome proliferator-activated receptor γ (PPARγ) is important for the treatment of type 2 diabetes and obesity through the regulation of glucose metabolism and fatty acid accumulation. Hence, the discovery of novel PPARγ agonists is necessary to overcome these diseases. In this study, a newly developed approach, multi-conformation dynamic pharmacophore modeling (MCDPM), was used for screening candidate compounds that can properly bind PPARγ. Highly populated structures obtained from molecular dynamics (MD) simulations were selected by clustering analysis. Based on these structures, pharmacophore models were generated from the ligand-binding pocket and then validated to check the rationality. Consequently, two hits were retrieved as final candidates by utilizing virtual screening and molecular docking simulations. These compounds can be used in the design of novel PPARγ agonists.

Molecular dynamics simulations of sonic hedgehog-receptor and inhibitor complexes and their applications for potential anticancer agent discovery.

The sonic hedgehog (Shh) signaling pathway is necessary for a variety of development and differentiation during embryogenesis as well as maintenance and renascence of diverse adult tissues. However, an abnormal activation of the signaling pathway is related to various cancers. In this pathway, the Shh signaling transduction is facilitated by binding of Shh to its receptor protein, Ptch. In this study, we modeled the 3D structure of functionally important key loop peptides of Ptch based on homologous proteins. Using this loop model, the molecular interactions between the structural components present in the pseudo-active site of Shh and key residues of Ptch was investigated in atomic level through molecular dynamics (MD) simulations. For the purpose of developing inhibitor candidates of the Shh signaling pathway, the Shh pseudo-active site of this interface region was selected as a target to block the direct binding between Shh and Ptch. Two different structure-based pharmacophore models were generated considering the key loop of Ptch and known inhibitor-induced conformational changes of the Shh through MD simulations. Finally two hit compounds were retrieved through a series of virtual screening combined with molecular docking simulations and we propose two hit compounds as potential inhibitory lead candidates to block the Shh signaling pathway based on their strong interactions to receptor or inhibitor induced conformations of the Shh.

Development of predictive quantitative structure-activity relationship model and its application in the discovery of human leukotriene A4 hydrolase inhibitors.

BACKGROUND: Human LTA4H catalyzes the conversion of LTA4 to LTB4 and plays a key role in innate immune responses. Inhibition of this enzyme can be a valid method in the treatment of inflammatory response exhibited through LTB4. RESULTS & DISCUSSION: The quantitative structure-activity relationship (QSAR) models were developed using genetic function approximation and validated. A training set of 26 diverse compounds and their molecular descriptors were used to develop highly correlating QSAR models. A six-descriptor model explaining the biological activity of the training and test sets with correlation values of 0.846 and 0.502, respectively, was selected as the best model and used in a database screening of drug-like Maybridge database followed by molecular docking. CONCLUSION: Based on the predicted potent inhibitory activities, expected binding mode and molecular interactions at the active site of hLTA4H final leads were selected as to be utilized in designing future hLTA4H inhibitors.

Molecular modeling study on tunnel behavior in different histone deacetylase isoforms.

Histone deacetylases (HDACs) have emerged as effective therapeutic targets in the treatment of various diseases including cancers as these enzymes directly involved in the epigenetic regulation of genes. However the development of isoform-selective HDAC inhibitors has been a challenge till date since all HDAC enzymes possess conserved tunnel-like active site. In this study, using molecular dynamics simulation we have analyzed the behavior of tunnels present in HDAC8, 10, and 11 enzymes of class I, II, and IV, respectively. We have identified the equivalent tunnel forming amino acids in these three isoforms and found that they are very much conserved with subtle differences to be utilized in selective inhibitor development. One amino acid, methionine of HDAC8, among six tunnel forming residues is different in isoforms of other classes (glutamic acid (E) in HDAC10 and leucine (L) in HDAC 11) based on which mutations were introduced in HDAC11, the less studied HDAC isoform, to observe the effects of this change. The HDAC8-like (L268M) mutation in the tunnel forming residues has almost maintained the deep and narrow tunnel as present in HDAC8 whereas HDAC10-like (L268E) mutation has changed the tunnel wider and shallow as observed in HDAC10. These results explained the importance of the single change in the tunnel formation in different isoforms. The observations from this study can be utilized in the development of isoform-selective HDAC inhibitors.

Comparative molecular modeling study of Arabidopsis NADPH-dependent thioredoxin reductase and its hybrid protein.

2-Cys peroxiredoxins (Prxs) play important roles in the protection of chloroplast proteins from oxidative damage. Arabidopsis NADPH-dependent thioredoxin reductase isotype C (AtNTRC) was identified as efficient electron donor for chloroplastic 2-Cys Prx-A. There are three isotypes (A, B, and C) of thioredoxin reductase (TrxR) in Arabidopsis. AtNTRA contains only TrxR domain, but AtNTRC consists of N-terminal TrxR and C-terminal thioredoxin (Trx) domains. AtNTRC has various oligomer structures, and Trx domain is important for chaperone activity. Our previous experimental study has reported that the hybrid protein (AtNTRA-(Trx-D)), which was a fusion of AtNTRA and Trx domain from AtNTRC, has formed variety of structures and shown strong chaperone activity. But, electron transfer mechanism was not detected at all. To find out the reason of this problem with structural basis, we performed two different molecular dynamics (MD) simulations on AtNTRC and AtNTRA-(Trx-D) proteins with same cofactors such as NADPH and flavin adenine dinucleotide (FAD) for 50 ns. Structural difference has found from superimposition of two structures that were taken relatively close to average structure. The main reason that AtNTRA-(Trx-D) cannot transfer the electron from TrxR domain to Trx domain is due to the difference of key catalytic residues in active site. The long distance between TrxR C153 and disulfide bond of Trx C387-C390 has been observed in AtNTRA-(Trx-D) because of following reasons: i) unstable and unfavorable interaction of the linker region, ii) shifted Trx domain, and iii) different or weak interface interaction of Trx domains. This study is one of the good examples for understanding the relationship between structure formation and reaction activity in hybrid protein. In addition, this study would be helpful for further study on the mechanism of electron transfer reaction in NADPH-dependent thioredoxin reductase proteins.

Structural origins for the loss of catalytic activities of bifunctional human LTA4H revealed through molecular dynamics simulations.

Human leukotriene A4 hydrolase (hLTA4H), which is the final and rate-limiting enzyme of arachidonic acid pathway, converts the unstable epoxide LTA4 to a proinflammatory lipid mediator LTB4 through its hydrolase function. The LTA4H is a bi-functional enzyme that also exhibits aminopeptidase activity with a preference over arginyl tripeptides. Various mutations including E271Q, R563A, and K565A have completely or partially abolished both the functions of this enzyme. The crystal structures with these mutations have not shown any structural changes to address the loss of functions. Molecular dynamics simulations of LTA4 and tripeptide complex structures with functional mutations were performed to investigate the structural and conformation changes that scripts the observed differences in catalytic functions. The observed protein-ligand hydrogen bonds and distances between the important catalytic components have correlated well with the experimental results. This study also confirms based on the structural observation that E271 is very important for both the functions as it holds the catalytic metal ion at its location for the catalysis and it also acts as N-terminal recognition residue during peptide binding. The comparison of binding modes of substrates revealed the structural changes explaining the importance of R563 and K565 residues and the required alignment of substrate at the active site. The results of this study provide valuable information to be utilized in designing potent hLTA4H inhibitors as anti-inflammatory agents.

Computational studies of novel chymase inhibitors against cardiovascular and allergic diseases: mechanism and inhibition.

To provide a new idea for drug design, a computational investigation is performed on chymase and its novel 1,4-diazepane-2,5-diones inhibitors that explores the crucial molecular features contributing to binding specificity. Molecular docking studies of inhibitors within the active site of chymase were carried out to rationalize the inhibitory properties of these compounds and understand their inhibition mechanism. The density functional theory method was used to optimize molecular structures with the subsequent analysis of highest occupied molecular orbital, lowest unoccupied molecular orbital, and molecular electrostatic potential maps, which revealed that negative potentials near 1,4-diazepane-2,5-diones ring are essential for effective binding of inhibitors at active site of enzyme. The Bayesian model with receiver operating curve statistic of 0.82 also identified arylsulfonyl and aminocarbonyl as the molecular features favoring and not favoring inhibition of chymase, respectively. Moreover, genetic function approximation was applied to construct 3D quantitative structure-activity relationships models. Two models (genetic function approximation model 1 r(2) = 0.812 and genetic function approximation model 2 r(2) = 0.783) performed better in terms of correlation coefficients and cross-validation analysis. In general, this study is used as example to illustrate how combinational use of 2D/3D quantitative structure-activity relationships modeling techniques, molecular docking, frontier molecular orbital density fields (highest occupied molecular orbital and lowest unoccupied molecular orbital), and molecular electrostatic potential analysis may be useful to gain an insight into the binding mechanism between enzyme and its inhibitors.

Molecular docking and dynamics simulation, receptor-based hypothesis: application to identify novel sirtuin 2 inhibitors.

Sirtuin, NAD(+)-dependent histone deacetylase enzyme, emerged as a potential therapeutic target, and modulations by small molecules could be effective drugs for various diseases. Owing to the absence of complex structure of sirtuin 2 (SIRT2), sirtinol was docked in the NAD(+) binding site and subjected to 5-nseconds molecular dynamics (MD) simulation. LigandScout was used to develop hypotheses based on 3-representative SIRT2 complex structures from MD. Three structure-based hypotheses are generated and merged to form dynamics hypothesis. The dynamics hypothesis was validated using test and decoy sets. The results showed that dynamic hypothesis represents the complementary features of SIRT2 active site. Dynamic hypothesis was used to screen ChemDiv database, and hits were filtered through ADMET, rule of five, and two different molecular docking studies. Finally, 21 molecules were selected as potent leads based on consensus score from LigandFit, Gold fitness score, binding affinity from VINA as well as based on the important interactions with critical residues in SIRT2 active site. Hence, we suggest that the dynamic hypothesis will be reliable in the identification of SIRT2 new lead as well as to reduce time and cost in the drug discovery process.

Molecular dynamics simulation study and hybrid pharmacophore model development in human LTA4H inhibitor design.

Human leukotriene A4 hydrolase (hLTA4H) is a bi-functional enzyme catalyzes the hydrolase and aminopeptidase functions upon the fatty acid and peptide substrates, respectively, utilizing the same but overlapping binding site. Particularly the hydrolase function of this enzyme catalyzes the rate-limiting step of the leukotriene (LT) cascade that converts the LTA4 to LTB4. This product is a potent pro-inflammatory activator of inflammatory responses and thus blocking this conversion provides a valuable means to design anti-inflammatory agents. Four structurally very similar chemical compounds with highly different inhibitory profile towards the hydrolase function of hLTA4H were selected from the literature. Molecular dynamics (MD) simulations of the complexes of hLTA4H with these inhibitors were performed and the results have provided valuable information explaining the reasons for the differences in their biological activities. Binding mode analysis revealed that the additional thiophene moiety of most active inhibitor helps the pyrrolidine moiety to interact the most important R563 and K565 residues. The hLTA4H complexes with the most active compound and substrate were utilized in the development of hybrid pharmacophore models. These developed pharmacophore models were used in screening chemical databases in order to identify lead candidates to design potent hLTA4H inhibitors. Final evaluation based on molecular docking and electronic parameters has identified three compounds of diverse chemical scaffolds as potential leads to be used in novel and potent hLTA4H inhibitor design.

Construction and characterization of chimeric cellulases with enhanced catalytic activity towards insoluble cellulosic substrates.

The chimeric proteins viz. CBM3-Cel9A, CBM4-Cel9A and CBM30-Cel9A, are constructed by fusion of family 3, 4, and 30 cellulose binding modules (CBMs) to N-terminus of family 9 endoglucanase (Cel9A) from Alicyclobacillus acidocaldrious. The chimeric enzymes were successfully expressed in Escherichia coli and purified to homogeneity. The chimeric enzymes showed significant increase in Avicel (8-12 folds) and filter paper (7-10 folds) degradation activities compared to Cel9A endoglucanase. Computational protein modeling and simulation on the chimeric enzymes were applied to analyze the fused CBMs effect on the increased insoluble cellulosic substrates degradation activity. Thin layer chromatography analysis of the enzymatic hydrolysis products and distribution of reducing sugars between soluble and insoluble fractions indicated processive cleavage of insoluble cellulosic substrates by the chimeras. The fused CBMs played a critical accessory role for the Cel9A catalytic domain and changed its character to facilitate the processive cleavage of insoluble cellulosic substrates.

Potential human cholesterol esterase inhibitor design: benefits from the molecular dynamics simulations and pharmacophore modeling studies.

Human pancreatic cholesterol esterase (hCEase) is one of the lipases found to involve in the digestion of large and broad spectrum of substrates including triglycerides, phospholipids, cholesteryl esters, etc. The presence of bile salts is found to be very important for the activation of hCEase. Molecular dynamic simulations were performed for the apoform and bile salt complexed form of hCEase using the co-ordinates of two bile salts from bovine CEase. The stability of the systems throughout the simulation time was checked and two representative structures from the highly populated regions were selected using cluster analysis. These two representative structures were used in pharmacophore model generation. The generated pharmacophore models were validated and used in database screening. The screened hits were refined for their drug-like properties based on Lipinski's rule of five and ADMET properties. The drug-like compounds were further refined by molecular docking simulation using GOLD program based on the GOLD fitness score, mode of binding, and molecular interactions with the active site amino acids. Finally, three hits of novel scaffolds were selected as potential leads to be used in novel and potent hCEase inhibitor design. The stability of binding modes and molecular interactions of these final hits were re-assured by molecular dynamics simulations.

Pharmacophore modeling, molecular docking, and molecular dynamics simulation approaches for identifying new lead compounds for inhibiting aldose reductase 2.

Aldose reductase 2 (ALR2), which catalyzes the reduction of glucose to sorbitol using NADP as a cofactor, has been implicated in the etiology of secondary complications of diabetes. A pharmacophore model, Hypo1, was built based on 26 compounds with known ALR2-inhibiting activity values. Hypo1 contains important chemical features required for an ALR2 inhibitor, and demonstrates good predictive ability by having a high correlation coefficient (0.95) as well as the highest cost difference (128.44) and the lowest RMS deviation (1.02) among the ten pharmacophore models examined. Hypo1 was further validated by Fisher's randomization method (95%), test set (r = 0.91), and the decoy set shows the goodness of fit (0.70). Furthermore, during virtual screening, Hypo1 was used as a 3D query to screen the NCI database, and the hit leads were sorted by applying Lipinski's rule of five and ADME properties. The best-fitting leads were subjected to docking to identify a suitable orientation at the ALR2 active site. The molecule that showed the strongest interactions with the critical amino acids was used in molecular dynamics simulations to calculate its binding affinity to the candidate molecules. Thus, Hypo1 describes the key structure-activity relationship along with the estimated activities of ALR2 inhibitors. The hit molecules were searched against PubChem to find similar molecules with new scaffolds. Finally, four molecules were found to satisfy all of the chemical features and the geometric constraints of Hypo1, as well as to show good dock scores, PLPs and PMFs. Thus, we believe that Hypo1 facilitates the selection of novel scaffolds for ALR2, allowing new classes of ALR2 inhibitors to be designed.

Molecular dynamics simulation study explaining inhibitor selectivity in different class of histone deacetylases.

Histone deacetylases (HDACs) are key regulators of gene expression and thereby compelling targets in the treatment of various cancers. Class- and isoform-selective HDAC inhibitors targeting the particular isoform to treat cancers without affecting the normal expression of other isoforms are highly desirable. Molecular dynamics simulations were performed with the set of selective inhibitors and HDAC isoforms of three different classes. The results were compared both within and across the isoforms. The hydrogen bonds between protein and inhibitors are directly correlated with the selective experimental activity. The calculated distances between important amino acids and the metal binding part of inhibitors have disclosed the optimal distance to be maintained by a selective inhibitor. In addition, the calculated non-bonded interaction energies between inhibitor and catalytic residues revealed that the subtle difference in the amino acids at the highly conserved active sites of HDAC isoforms effectively scripts the selectivity story observed experimentally. The results of this study provide valuable information in designing highly selective HDAC inhibitors.