Search

Phthalate acid ester hydrolases are of promise for the microbial degradation of phthalate acid ester plasticizers, however, only a small number of PAE hydrolases have been characterized. Here, the authors report the structural and functional characterization of MehpH, a monoalkyl phthalate (MBP) hydrolase that degrades MBP to phthalic acid and butanol.

Flexible thermoelectric generators (FTEGs), which can overcome the energy supply limitations of wearable devices, have received considerable attention. However, the use of toxic Te-based materials and fracture-prone electrodes constrains the application of FTEGs. In this study, a novel Ag2Se and Poly (3,4-ethylene dioxythiophene): poly (styrene sulfonate) (PEDOT:PSS)/multi-walled carbon nanotube (MWCNT) FTEG with a high output performance and good flexibility is developed. The thermoelectric columns formulated in the work are environmentally friendly and reliable. The key enabler of this work is the use of embedded EGaIn electrodes, which increase the temperature difference collected by the thermoelectric column, thereby improving the FTEG output performance. Additionally, the embedded EGaIn electrodes could be directly printed on polydimethylsiloxane (PDMS) molds without wax paper, which simplifies the preparation process of FTEGs and enhances the fabrication efficiency. The FTEG with embedded electrodes exhibits the highest output power density of 25.83 μW/cm2 and the highest output power of 10.95 μW at ΔT = 15 K. The latter is 31.6% higher than that of silver-based FTEGs and 2.5% higher than that of covered EGaIn-based FTEGs. Moreover, the prepared FTEG has an excellent flexibility (>1500 bends) and output power stability (>30 days). At high humidity and high temperature, the prepared FTEG maintains good performance. These results demonstrate that the prepared FTEGs can be used as a stable and environmentally friendly energy supply for wearable devices.

The application feedback on existing silicon nitride ceramic bearings and RCF experimental research all indicate that the primary failure mode of silicon nitride ceramic bearings is material spalling on the contact surface. Spalling failure occurs due to the initiation and propagation of cracks under rolling contact. However, silicon nitride ceramic bearings, owing to their unique manufacturing method, inevitably exhibit defects and cracks. Therefore, as silicon nitride ceramic bearings are increasingly prevalent, reducing the probability of spalling failure is crucial for extending their service life. This can only be achieved by gaining a clear understanding of the crack initiation and expansion mechanisms in silicon nitride ceramic bearings. This paper is based on silicon nitride rolling friction experiments. It involves the joint simulation of Franc3D-V8.4 and ABAQUS2020, wherein the crack front SIFs are calculated for each load contact position of the surface crack on the silicon nitride ceramic bearing ring during cyclic movement. The study also delves into the determination of the maximum effective stress intensity factors and explores the influence of the initial crack depth on the cycle life and direction of crack propagation. The research yields several valuable conclusions. The findings of this research offer theoretical guidance for formulating grinding technologies for silicon nitride rings and adjusting and controlling working parameters of silicon nitride ceramic ball bearings. These insights are crucial for enhancing the reliability and longevity of silicon nitride ceramic bearings in practical applications.

Aromatic aldehydes are important aromatic compounds for the flavour and fragrance industry. In this study, a parallel cascade combining aryl alcohol oxidase from Pleurotus eryngii (PeAAOx) and unspecific peroxygenase from the basidiomycete Agrocybe aegerita (AaeUPO) to convert aromatic primary alcohols into high-value aromatic aldehydes is proposed. Key influencing factors in the process of enzyme cascade catalysis, such as enzyme dosage, pH and temperature, were investigated. The universality of PeAAOx coupled with AaeUPO cascade catalysis for the synthesis of aromatic aldehyde flavour compounds from aromatic primary alcohols was evaluated. In a partially optimised system (comprising 30 μM PeAAOx, 2 μM AaeUPO at pH 7 and 40 °C) up to 84% conversion of 50 mM veratryl alcohol into veratryl aldehyde was achieved in a self-sufficient aerobic reaction. Promising turnover numbers of 2800 and 21,000 for PeAAOx and AaeUPO, respectively, point towards practical applicability.

A new phospholipase D from marine Moritella sp. JT01 (MsPLD) was recombinantly expressed and biochemically characterized. The optimal reaction temperature and pH of MsPLD were determined to be 35 °C and 8.0. MsPLD was stable at a temperature lower than 35 °C, and the t1/2 at 4 °C was 41 days. The crystal structure of apo-MsPLD was resolved and the functions of a unique extra loop segment on the enzyme activity were characterized. The results indicated that a direct deletion or fastening of the extra loop segment by introducing disulfide bonds both resulted in a complete loss of its activity. The results of the maximum insertion pressure indicated that the deletion of the extra loop segment significantly decreased MsPLD’s interfacial binding properties to phospholipid monolayers. Finally, MsPLD was applied to the synthesis of phosphatidic acid by using a biphasic reaction system. Under optimal reaction conditions, the conversion rate of phosphatidic acid reached 86%. The present research provides a foundation for revealing the structural–functional relationship of this enzyme.

Mining of Phospholipase D (PLD) with high activity and stability has attracted strong interest for investigation. A novel PLD from marine Moritella sp. JT01 (MsPLD) was biochemically and structurally characterized in our previous study; however, the short half-life time (t1/2) under its optimum reaction temperature seriously hampered its further applications. Herein, the disulfide bond engineering strategy was applied to improve its thermostability. Compared with wild-type MsPLD, mutant S148C-T206C/D225C-A328C with the addition of two disulfide bonds exhibited a 3.1-fold t1/2 at 35 °C and a 5.7 °C increase in melting temperature (Tm). Unexpectedly, its specific activity and catalytic efficiency (kcat/Km) also increased by 22.7% and 36.5%, respectively. The enhanced activity might be attributed to an increase in the activation entropy by displacing more water molecules by the transition state. The results of molecular dynamics simulations (MD) revealed that the introduction of double disulfide bonds rigidified the global structure of the mutant, which might cause the enhanced thermostability. Finally, the synthesis capacity of the mutant to synthesize phosphatidic acid (PA) was evaluated. The conversion rate of PA reached about 80% after 6 h reaction with wild-type MsPLD but reached 78% after 2 h with mutant S148C-T206C/D225C-A328C, which significantly reduced the time needed for the reaction to reach equilibrium. The present results pave the way for further application of MsPLD in the food and pharmaceutical industries.

Single mutations that can substantially enhance stability are highly desirable for protein engineering. However, it is generally rare for this kind of mutant to emerge from directed evolution experiments. This study used computational approaches to identify hotspots in a diacylglycerol-specific lipase for mutagenesis with functional hotspot and sequence consensus strategies, followed by ∆∆G calculations for all possible mutations using the Rosetta ddg_monomer protocol. Single mutants with significant ∆∆G changes (≤−2.5 kcal/mol) were selected for expression and characterization. Three out of seven tested mutants showed a significantly enhanced thermostability, with Q282W and A292G in the catalytic pocket and D245G located on the opposite surface of the protein. Remarkably, A292G increased the T5015 (the temperature at which 50% of the enzyme activity was lost after a 15 min of incubation) by ~7 °C, concomitant with a twofold increase in enzymatic activity at the optimal reaction temperature. Structural analysis showed that both A292 and D245 adopted unfavored dihedral angles in the wild-type (WT) enzyme. Substitution of them by glycine might release a steric strain to increase the stability. In sum, substitution by glycine might be a promising strategy to improve protein thermostability.

The efficiency and accuracy of the synthesis of structural lipids are closely related to the regiospecificity of lipases. Understanding the structural mechanism of their regiospecificity contributes to the regiospecific redesign of lipases for meeting the technological innovation needs. Here, we used a thermostable lipase from Streptomyces sp. W007 (MAS1), which has been recently reported to show great potential in industry, to gain an insight into the structural basis of its regiospecificity by molecular modelling and mutagenesis experiments. The results indicated that increasing the steric hindrance of the site for binding a non-reactive carbonyl group of TAGs could transform the non-specific MAS1 to a α-specific lipase, such as the mutants G40E, G40F, G40Q, G40R, G40W, G40Y, N45Y, H108W and T237Y (PSI > 80). In addition, altering the local polarity of the site as well as the conformational stability of its composing residues could also impact the regiospecificity. Our present study could not only aid the rational design of the regiospecificity of lipases, but open avenues of exploration for further industrial applications of lipases.

Mining of phospholipase D (PLD) with altered acyl group recognition except its head group specificity is also useful in terms of specific acyl size phospholipid production and as diagnostic reagents for quantifying specific phospholipid species. Microbial PLDs from Actinomycetes, especially Streptomyces, best fit this process requirements. In the present studies, a new PLD from marine Streptomyces klenkii (SkPLD) was purified and biochemically characterized. The optimal reaction temperature and pH of SkPLD were determined to be 60 °C and 8.0, respectively. Kinetic analysis showed that SkPLD had the relatively high catalytic efficiency toward phosphatidylcholines (PCs) with medium acyl chain length, especially 12:0/12:0-PC (67.13 S−1 mM−1), but lower catalytic efficiency toward PCs with long acyl chain (>16 fatty acids). Molecular docking results indicated that the different catalytic efficiency was related to the increased steric hindrance of long acyl-chains in the substrate-binding pockets and differences in hydrogen-bond interactions between the acyl chains and substrate-binding pockets. The enzyme displayed suitable transphosphatidylation activity and the reaction process showed 26.18% yield with L-serine and soybean PC as substrates. Present study not only enriched the PLD enzyme library but also provide guidance for the further mining of PLDs with special phospholipids recognition properties.

Cement and ionic curing agent were blended with iron tailings (ITS) and natural soil to prepare a new road base material. 7-d unconfined compressive strength (UCS) was selected as the evaluation index for determining the optimal mix ratio. The results reveal that the 7-d UCS always increases with the increase of cement content, and firstly increases with the increase of curing agent content and then tends to remain stable, exhibiting an increasing trend first and then dropping with the increase of ITS content. In lieu of this, the optimum ratio of iron tailings blended soil (ITBS) is determined as follows: ITS and soil are mixed in ratio of 50:50, and cement and soil curing agent account for 5% and 2% of the total mass, respectively. The results of splitting strength and flexural-tensile strength tests show that these strengths of the ITBS always increase with the curing age, and firstly increases and then decreases with the increase in ITS content. XRD and SEM analyses show that the ratio of ITS to soil significantly influences the micro-pore structure of the blended soil. At an ITS to soil ratio of 50%:50%, the combination of soil and ITS particles is the closest which results in the highest strength attainment.

Phospholipases D (PLDs) play important roles in different organisms and in vitro phospholipid modifications, which attract strong interests for investigation. However, the lack of PLD structural information has seriously hampered both the understanding of their structure–function relationships and the structure-based bioengineering of this enzyme. Herein, we presented the crystal structure of a PLD from the plant-associated bacteria Serratia plymuthica strain AS9 (SpPLD) at a resolution of 1.79 Å. Two classical HxKxxxxD (HKD) motifs were found in SpPLD and have shown high structural consistence with several PLDs in the same family. While comparing the structure of SpPLD with the previous resolved PLDs from the same family, several unique conformations on the C-terminus of the HKD motif were demonstrated to participate in the arrangement of the catalytic pocket of SpPLD. In SpPLD, an extented loop conformation between β9 and α9 (aa228–246) was found. Moreover, electrostatic surface potential showed that this loop region in SpPLD was positively charged while the corresponding loops in the two Streptomyces originated PLDs (PDB ID: 1F0I, 2ZE4/2ZE9) were neutral. The shortened loop between α10 and α11 (aa272–275) made the SpPLD unable to form the gate-like structure which existed specically in the two Streptomyces originated PLDs (PDB ID: 1F0I, 2ZE4/2ZE9) and functioned to stabilize the substrates. In contrast, the shortened loop conformation at this corresponding segment was more alike to several nucleases (Nuc, Zuc, mZuc, NucT) within the same family. Moreover, the loop composition between β11 and β12 was also different from the two Streptomyces originated PLDs (PDB ID: 1F0I, 2ZE4/2ZE9), which formed the entrance of the catalytic pocket and were closely related to substrate recognition. So far, SpPLD was the only structurally characterized PLD enzyme from Serratia. The structural information derived here not only helps for the understanding of the biological function of this enzyme in plant protection, but also helps for the understanding of the rational design of the mutant, with potential application in phospholipid modification.

Presently, 151 widely-diverse pyridinylimidazole-based compounds that show inhibitory activities at the TNF-α release were investigated. By using the distance comparison technique (DISCOtech), comparative molecular field analysis (CoMFA), and comparative molecular similarity index analysis (CoMSIA) methods, the pharmacophore models and the three-dimensional quantitative structure-activity relationships (3D-QSAR) of the compounds were explored. The proposed pharmacophore model, including two hydrophobic sites, two aromatic centers, two H-bond donor atoms, two H-bond acceptor atoms, and two H-bond donor sites characterizes the necessary structural features of TNF-α release inhibitors. Both the resultant CoMFA and CoMSIA models exhibited satisfactory predictability (with Q2 (cross-validated correlation coefficient) = 0.557, R2ncv (non-cross-validated correlation coefficient) = 0.740, R2pre (predicted correlation coefficient) = 0.749 and Q2 = 0.598, R2ncv = 0.767, R2pre = 0.860, respectively). Good consistency was observed between the 3D-QSAR models and the pharmacophore model that the hydrophobic interaction and hydrogen bonds play crucial roles in the mechanism of actions. The corresponding contour maps generated by these models provide more diverse information about the key intermolecular interactions of inhibitors with the surrounding environment. All these models have extended the understanding of imidazole-based compounds in the structure-activity relationship, and are useful for rational design and screening of novel 2-thioimidazole-based TNF-α release inhibitors.

Cyclin-dependent kinase 2 (CDK2) is a crucial regulator of the eukaryotic cell cycle. However it is well established that monomeric CDK2 lacks regulatory activity, which needs to be aroused by its positive regulators, cyclins E and A, or be phosphorylated on the catalytic segment. Interestingly, these activation steps bring some dynamic changes on the 3D-structure of the kinase, especially the activation segment. Until now, in the monomeric CDK2 structure, three binding sites have been reported, including the adenosine triphosphate (ATP) binding site (Site I) and two non-competitive binding sites (Site II and III). In addition, when the kinase is subjected to the cyclin binding process, the resulting structural changes give rise to a variation of the ATP binding site, thus generating an allosteric binding site (Site IV). All the four sites are demonstrated as being targeted by corresponding inhibitors, as is illustrated by the allosteric binding one which is targeted by inhibitor ANS (fluorophore 8-anilino-1-naphthalene sulfonate). In the present work, the binding mechanisms and their fluctuations during the activation process attract our attention. Therefore, we carry out corresponding studies on the structural characterization of CDK2, which are expected to facilitate the understanding of the molecular mechanisms of kinase proteins. Besides, the binding mechanisms of CDK2 with its relevant inhibitors, as well as the changes of binding mechanisms following conformational variations of CDK2, are summarized and compared. The summary of the conformational characteristics and ligand binding mechanisms of CDK2 in the present work will improve our understanding of the molecular mechanisms regulating the bioactivities of CDK2.

Thermostability and substrate specificity are important characteristics of enzymes for industrial application, which can be improved by protein engineering. SMG1 lipase from Malassezia globosa is a mono- and diacylglycerol lipase (MDL) that shows activity toward mono- and diacylglycerols, but no activity toward triacylglycerols. SMG1 lipase is considered a potential biocatalyst applied in oil/fat modification and its crystal structure revealed that an interesting residue-Asn277 may contribute to stabilize loop 273–278 and the 3104 helix which are important to enzyme characterization. In this study, to explore its role in affecting the stability and catalytic activity, mutagenesis of N277 with Asp (D), Val (V), Leu (L) and Phe (F) was conducted. Circular dichroism (CD) spectral analysis and half-life measurement showed that the N277D mutant has better thermostability. The melting temperature and half-life of the N277D mutant were 56.6 °C and 187 min, respectively, while that was 54.6 °C and 121 min for SMG1 wild type (WT). Biochemical characterization of SMG1 mutants were carried out to test whether catalytic properties were affected by mutagenesis. N277D had similar enzymatic properties as SMG1 WT, but N277F showed a different substrate selectivity profile as compared to other SMG1 mutants. Analysis of the SMG1 3D model suggested that N277D formed a salt bridge via its negative charged carboxyl group with a positively charged guanidino group of R227, which might contribute to confer N277D higher temperature stability. These findings not only provide some clues to understand the molecular basis of the lipase structure/function relationship but also lay the framework for engineering suitable MDL lipases for industrial applications.

Diacylglycerol (DAG)-like lipases are found to play an important role in the life sciences and industrial fields. A putative DAG-like lipase (MgMDL2) from Malassezia globosa was cloned and expressed in recombinant Pichia pastoris. The recombinant MgMDL2 was expressed as a glycosylated protein and purified into homogeneity by anion exchange chromatography. The activity of recombinant MgMDL2 was optimal at 15 °C and pH 6.0, and it keeps over 50% of relative activity at 5 °C, suggesting that MgMDL2 was a cold active lipase. MgMDL2 retained over 80% of initial activity after incubation at 30 and 40 °C for 2.5 h, but it was not stable at 50 °C. Incubation of methanol and ethanol at a concentration of 30% for 2 h did not affect the recombinant enzyme activity, while metal ions, including Ca2+, Mn2+ and Ni2+, sharply inhibited the MgMDL2 activity at 5 mM by 42%, 35% and 36%, respectively. MgMDL2 exhibited a preference for medium chain-length esters with highest activity toward p-nitrophenyl caprylate, while it was active on mono- and diacylglycerol but not on triacylglycerol, indicating that it was a typical DAG-like lipase. By homology modeling, Phe278 was predicted to be involved in the preference of MgMDL2 for monoacyl- and diacyl-glyceride substrates, but not triglycerides.

Inflammation is a hallmark of many diseases like diabetes, cancers, atherosclerosis and arthritis. Thus, lots of concerns have been raised toward developing novel anti-inflammatory agents. Many alternative herbal medicines possess excellent anti-inflammatory properties, yet their precise mechanisms of action are yet to be elucidated. Here, a novel systems pharmacology approach based on a large number of chemical, biological and pharmacological data was developed and exemplified by a probe herb Folium Eriobotryae, a widely used clinical anti-inflammatory botanic drug. The results show that 11 ingredients of this herb with favorable pharmacokinetic properties are predicted as active compounds for anti-inflammatory treatment. In addition, via systematic network analyses, their targets are identified to be 43 inflammation-associated proteins including especially COX2, ALOX5, PPARG, TNF and RELA that are mainly involved in the mitogen-activated protein kinase (MAPK) signaling pathway, the rheumatoid arthritis pathway and NF-κB signaling pathway. All these demonstrate that the integrated systems pharmacology method provides not only an effective tool to illustrate the anti-inflammatory mechanisms of herbs, but also a new systems-based approach for drug discovery from, but not limited to, herbs, especially when combined with further experimental validations.

Mono- and di-acylglycerol lipase has been applied to industrial usage in oil modification for its special substrate selectivity. Until now, the reported mono- and di-acylglycerol lipases from microorganism are limited, and there is no report on the mono- and di-acylglycerol lipase from bacteria. A predicted lipase (named MAJ1) from marine Janibacter sp. strain HTCC2649 was purified and biochemical characterized. MAJ1 was clustered in the family I.7 of esterase/lipase. The optimum activity of the purified MAJ1 occurred at pH 7.0 and 30 °C. The enzyme retained 50% of the optimum activity at 5 °C, indicating that MAJ1 is a cold-active lipase. The enzyme activity was stable in the presence of various metal ions, and inhibited in EDTA. MAJ1 was resistant to detergents. MAJ1 preferentially hydrolyzed mono- and di-acylglycerols, but did not show activity to triacylglycerols of camellia oil substrates. Further, MAJ1 is low homologous to that of the reported fungal diacylglycerol lipases, including Malassezia globosa lipase 1 (SMG1), Penicillium camembertii lipase U-150 (PCL), and Aspergillus oryzae lipase (AOL). Thus, we identified a novel cold-active bacterial lipase with a sn-1/3 preference towards mono- and di-acylglycerides for the first time. Moreover, it has the potential, in oil modification, for special substrate selectivity.

Recent studies revealed the role of lipase in the pathogenicity of Malassezia restricta in dandruff and seborrheic dermatitis (D/SD). The lipase from M. restricta (Mrlip1) is considered a potential target for dandruff therapy. In this work, we performed structure-based virtual screening in Zinc database to find the natural bioactive inhibitors of Mrlip1. We identified three compounds bearing superior affinity and specificity from the Traditional Chinese Medicine database (~60,000 compounds), and their binding patterns with Mrlip1 were analyzed in detail. Additionally, we performed three sets of 100 ns MD simulations of each complex in order to understand the interaction mechanism of Mrlip1 with known inhibitor RHC80267 and the newly identified compounds such as ZINC85530919, ZINC95914464 and ZINC85530320, respectively. These compounds bind to the active site cavity and cause conformational changes in Mrlip1. The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) studies suggested that the average binding energy was stronger in the case of Mrlip1-ZINC85530919 and Mrlip1-ZINC95914464. The selected natural inhibitors might act as promising lead drugs against Mrlip1. Further, the present study will contribute to various steps involved in developing and creating potent drugs for several skin diseases including dandruff.

Lipases with unique substrate specificity are highly desired in biotechnological applications. In this study, a putative marine Geobacillus sp. monoacylglycerol lipase (GMGL) encoded gene was identified by a genomic mining strategy. The gene was expressed in Escherichia coli as a His-tag fusion protein and purified by affinity chromatography with a yield of 264 mg per liter fermentation broth. The recombinant GMGL shows the highest hydrolysis activity at 60 °C and pH 8.0, and the half-life was 60 min at 70 °C. The GMGL is active on monoacylglycerol (MAG) substrate but not diacylglycerol (DAG) or triacylglycerol (TAG), and produces MAG as the single product in the esterification reaction. Modeling structure analysis showed that the catalytic triad is formed by Ser97, Asp196 and His226, and the flexible cap region is constituted by residues from Ala120 to Thr160. A mutagenesis study on Leu142, Ile145 and Ile170 located in the substrate binding tunnel revealed that these residues were related with its substrate specificity. The kcat/Km value toward the pNP-C6 substrate in mutants Leu142Ala, Ile145Ala and Ile170Phe increased to 2.3-, 1.4- and 2.2-fold as compared to that of the wild type, respectively.

Salt stress is one of the most serious abiotic factors that inhibit plant growth. Dunaliella salina has been recognized as a model organism for stress response research due to its high capacity to tolerate extreme salt stress. A proteomic approach based on isobaric tags for relative and absolute quantitation (iTRAQ) was used to analyze the proteome of D. salina during early response to salt stress and identify the differentially abundant proteins (DAPs). A total of 141 DAPs were identified in salt-treated samples, including 75 upregulated and 66 downregulated DAPs after 3 and 24 h of salt stress. DAPs were annotated and classified into gene ontology functional groups. The Kyoto Encyclopedia of Genes and Genomes pathway analysis linked DAPs to tricarboxylic acid cycle, photosynthesis and oxidative phosphorylation. Using search tool for the retrieval of interacting genes (STRING) software, regulatory protein⁻protein interaction (PPI) networks of the DAPs containing 33 and 52 nodes were built at each time point, which showed that photosynthesis and ATP synthesis were crucial for the modulation of early salinity-responsive pathways. The corresponding transcript levels of five DAPs were quantified by quantitative real-time polymerase chain reaction (qRT-PCR). These results presented an overview of the systematic molecular response to salt stress. This study revealed a complex regulatory mechanism of early salt tolerance in D. salina and potentially contributes to developing strategies to improve stress resilience.

Traditional Chinese Medicine (TCM) involves a broad range of empirical testing and refinement and plays an important role in the health maintenance for people all over the world. However, due to the complexity of Chinese herbs, a full understanding of TCM’s action mechanisms is still unavailable despite plenty of successful applications of TCM in the treatment of various diseases, including especially cardiovascular diseases (CVD), one of the leading causes of death. Thus in the present work, by incorporating the chemical predictors, target predictors and network construction approaches, an integrated system of TCM has been constructed to systematically uncover the underlying action mechanisms of TCM. From three representative Chinese herbs, i.e., Ligusticum chuanxiong Hort., Dalbergia odorifera T. Chen and Corydalis yanhusuo WT Wang which have been widely used in CVD treatment, by combinational use of drug absorption, distribution, metabolism and excretion (ADME) screening and network pharmacology techniques, we have generated 64 bioactive ingredients and identified 54 protein targets closely associated with CVD, of which 29 are common targets (52.7%) of the three herbs. The result provides new information on the efficiency of the Chinese herbs for the treatment of CVD and also explains one of the basic theories of TCM, i.e., “multiple herbal drugs can treat one disease”. The predicted potential targets were then mapped to target-disease and target-signal pathway connections, which revealed the relationships of the active ingredients with their potential targets, diseases and signal systems. This means that for the first time, the action mechanism of these three important Chinese herbs for the treatment of CVD is uncovered, by generating and identifying both their active ingredients and novel targets specifically related to CVD, which clarifies some of the common conceptions in TCM, and thus provides clues to modernize such specific herbal medicines.

Orally administered drugs must overcome several barriers before reaching their target site. Such barriers depend largely upon specific membrane transport systems and intracellular drug-metabolizing enzymes. For the first time, the P-glycoprotein (P-gp) and cytochrome P450s, the main line of defense by limiting the oral bioavailability (OB) of drugs, were brought into construction of QSAR modeling for human OB based on 805 structurally diverse drug and drug-like molecules. The linear (multiple linear regression: MLR, and partial least squares regression: PLS) and nonlinear (support-vector machine regression: SVR) methods are used to construct the models with their predictivity verified with five-fold cross-validation and independent external tests. The performance of SVR is slightly better than that of MLR and PLS, as indicated by its determination coefficient (R2) of 0.80 and standard error of estimate (SEE) of 0.31 for test sets. For the MLR and PLS, they are relatively weak, showing prediction abilities of 0.60 and 0.64 for the training set with SEE of 0.40 and 0.31, respectively. Our study indicates that the MLR, PLS and SVR-based in silico models have good potential in facilitating the prediction of oral bioavailability and can be applied in future drug design.

Chitinolytic β-N-acetyl-D-hexosaminidases, as a class of chitin hydrolysis enzyme in insects, are a potential species-specific target for developing environmentally-friendly pesticides. Until now, pesticides targeting chitinolytic β-N-acetyl-D-hexosaminidase have not been developed. This study demonstrates a combination of different theoretical methods for investigating the key structural features of this enzyme responsible for pesticide inhibition, thus allowing for the discovery of novel small molecule inhibitors. Firstly, based on the currently reported crystal structure of this protein (OfHex1.pdb), we conducted a pre-screening of a drug-like compound database with 8 × 106 compounds by using the expanded pesticide-likeness criteria, followed by docking-based screening, obtaining 5 top-ranked compounds with favorable docking conformation into OfHex1. Secondly, molecular docking and molecular dynamics simulations are performed for the five complexes and demonstrate that one main hydrophobic pocket formed by residues Trp424, Trp448 and Trp524, which is significant for stabilization of the ligand–receptor complex, and key residues Asp477 and Trp490, are respectively responsible for forming hydrogen-bonding and π–π stacking interactions with the ligands. Finally, the molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) analysis indicates that van der Waals interactions are the main driving force for the inhibitor binding that agrees with the fact that the binding pocket of OfHex1 is mainly composed of hydrophobic residues. These results suggest that screening the ZINC database can maximize the identification of potential OfHex1 inhibitors and the computational protocol will be valuable for screening potential inhibitors of the binding mode, which is useful for the future rational design of novel, potent OfHex1-specific pesticides.

Experiments were carried out to investigate the effects of various factors on the activity and conformation of recombinant leucine aminopeptidase of Bacillus kaustophilus CCRC 11223 (BkLAP) and potential utilization of BkLAP in the hydrolysis of anchovy protein. Optimal temperature and pH of BkLAP were 70 °C and 8.0 in potassium-phosphate buffer, respectively, and the activity was strongly stimulated by Ni2+, followed by Mn2+ and Co2+. Conformational studies via circular dichroism spectroscopy indicated that various factors could influence the secondary structure of BkLAP to different extents and further induce the changes in enzymatic activity. The secondary structure of BkLAP was slightly modified by Ni2+ at the concentration of 1×10−4 M, however, significant changes on the secondary structures of the enzyme were observed when Hg2+ was added to the concentration of 1×10−4 M. The potential application of BkLAP was evaluated through combination with the commercial or endogenous enzyme to hydrolysis the anchovy protein. Results showed that combining the BkLAP with other enzymes could significantly increase the degree of hydrolysis and amino acid component of hydrolysate. In this regard, BkLAP is a potential enzyme that can be used in the protein hydrolysate industry.

Lipases from microorganisms have multi-faceted properties and play an important role in ever-growing modern biotechnology and, consequently, it is of great significance to develop new ones. In the present work, a lipase gene from Candida albicans (CaLIP10) was cloned and two non-unusual CUG serine codons were mutated into universal codons, and its expression in Pichia pastoris performed optimally, as shown by response surface methodology. Optimal conditions were: initial pH of culture 6.86, temperature 25.53 °C, 3.48% of glucose and 1.32% of yeast extract. The corresponding maximal lipolytic activity of CaLIP10 was 8.06 U/mL. The purified CaLIP10 showed maximal activity at pH 8.0 and 25 °C, and a good resistance to non-ionic surfactants and polar organic solvent was noticed. CaLIP10 could effectively hydrolyze coconut oil, but exhibited no obvious preference to the fatty acids with different carbon length, and diacylglycerol was accumulated in the reaction products, suggesting that CaLIP10 is a potential lipase for the oil industry.

In order to obtain structural features of 3-arylpyrimidin-2,4-diones emerged as promising inhibitors of insect γ-aminobutyric acid (GABA) receptor, a set of ligand-/receptor-based 3D-QSAR models for 60 derivatives are generated using Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Index Analysis (CoMSIA). The statistically optimal CoMSIA model is produced with highest q2 of 0.62, r2ncv of 0.97, and r2pred of 0.95. A minor/bulky electronegative hydrophilic polar substituent at the 1-/6-postion of the uracil ring, and bulky substituents at the 3'-, 4'- and 5'-positions of the benzene ring are beneficial for the enhanced potency of the inhibitors as revealed by the obtained 3D-contour maps. Furthermore, homology modeling, molecular dynamics (MD) simulation and molecular docking are also carried out to gain a better understanding of the probable binding modes of these inhibitors, and the results show that residues Ala-183(C), Thr-187(B), Thr-187(D) and Thr-187(E) in the second transmembrane domains of GABA receptor are responsible for the H-bonding interactions with the inhibitor. The good correlation between docking observations and 3D-QSAR analyses further proves the model reasonability in probing the structural features and the binding mode of 3-arylpyrimidin-2,4-dione derivatives within the housefly GABA receptor.

An abnormal ubiquitin-proteasome is found in many human diseases, especially in cancer, and has received extensive attention as a promising therapeutic target in recent years. In this work, several in silico models have been built with two classes of proteasome inhibitors (PIs) by using 3D-QSAR, homology modeling, molecular docking and molecular dynamics (MD) simulations. The study resulted in two types of satisfactory 3D-QSAR models, i.e., the CoMFA model (Q2 = 0.462, R2pred = 0.820) for epoxyketone inhibitors (EPK) and the CoMSIA model (Q2 = 0.622, R2pred = 0.821) for tyropeptin-boronic acid derivatives (TBA). From the contour maps, some key structural factors responsible for the activity of these two series of PIs are revealed. For EPK inhibitors, the N-cap part should have higher electropositivity; a large substituent such as a benzene ring is favored at the C6-position. In terms of TBA inhibitors, hydrophobic substituents with a larger size anisole group are preferential at the C8-position; higher electropositive substituents like a naphthalene group at the C3-position can enhance the activity of the drug by providing hydrogen bond interaction with the protein target. Molecular docking disclosed that residues Thr60, Thr80, Gly106 and Ser189 play a pivotal role in maintaining the drug-target interactions, which are consistent with the contour maps. MD simulations further indicated that the binding modes of each conformation derived from docking is stable and in accord with the corresponding structure extracted from MD simulation overall. These results can offer useful theoretical references for designing more potent PIs.

Experimental pEC50s for 216 selective respiratory syncytial virus (RSV) inhibitors are used to develop classification models as a potential screening tool for a large library of target compounds. Variable selection algorithm coupled with random forests (VS-RF) is used to extract the physicochemical features most relevant to the RSV inhibition. Based on the selected small set of descriptors, four other widely used approaches, i.e., support vector machine (SVM), Gaussian process (GP), linear discriminant analysis (LDA) and k nearest neighbors (kNN) routines are also employed and compared with the VS-RF method in terms of several of rigorous evaluation criteria. The obtained results indicate that the VS-RF model is a powerful tool for classification of RSV inhibitors, producing the highest overall accuracy of 94.34% for the external prediction set, which significantly outperforms the other four methods with the average accuracy of 80.66%. The proposed model with excellent prediction capacity from internal to external quality should be important for screening and optimization of potential RSV inhibitors prior to chemical synthesis in drug development.

Hsp90 is involved in correcting, folding, maturation and activation of a diverse array of client proteins; it has also been implicated in the treatment of cancer in recent years. In this work, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular docking and molecular dynamics were performed on three different series of Hsp90 inhibitors to build 3D-QSAR models, which were based on the ligand-based or receptor-based methods. The optimum 3D-QSAR models exhibited reasonable statistical characteristics with averaging internal q2 > 0.60 and external r2pred > 0.66 for Benzamide tetrahydro-4H-carbazol-4-one analogs (BT), AT13387 derivatives (AT) and Dihydroxylphenyl amides (DA). The results revealed that steric effects contributed the most to the BT model, whereas H-bonding was more important to AT, and electrostatic, hydrophobic, H-bond donor almost contributed equally to the DA model. The docking analysis showed that Asp93, Tyr139 and Thr184 in Hsp90 are important for the three series of inhibitors. Molecular dynamics simulation (MD) further indicated that the conformation derived from docking is basically consistent with the average structure extracted from MD simulation. These results not only lead to a better understanding of interactions between these inhibitors and Hsp90 receptor but also provide useful information for the design of new inhibitors with a specific activity.

Over the years development of selective estrogen receptor (ER) ligands has been of great concern to researchers involved in the chemistry and pharmacology of anticancer drugs, resulting in numerous synthesized selective ER subtype inhibitors. In this work, a data set of 82 ER ligands with ERα and ERβ inhibitory activities was built, and quantitative structure-activity relationship (QSAR) methods based on the two linear (multiple linear regression, MLR, partial least squares regression, PLSR) and a nonlinear statistical method (Bayesian regularized neural network, BRNN) were applied to investigate the potential relationship of molecular structural features related to the activity and selectivity of these ligands. For ERα and ERβ, the performances of the MLR and PLSR models are superior to the BRNN model, giving more reasonable statistical properties (ERα: for MLR, Rtr2 = 0.72, Qte2 = 0.63; for PLSR, Rtr2 = 0.92, Qte2 = 0.84. ERβ: for MLR, Rtr2 = 0.75, Qte2 = 0.75; for PLSR, Rtr2 = 0.98, Qte2 = 0.80). The MLR method is also more powerful than other two methods for generating the subtype selectivity models, resulting in Rtr2 = 0.74 and Qte2 = 0.80. In addition, the molecular docking method was also used to explore the possible binding modes of the ligands and a relationship between the 3D-binding modes and the 2D-molecular structural features of ligands was further explored. The results show that the binding affinity strength for both ERα and ERβ is more correlated with the atom fragment type, polarity, electronegativites and hydrophobicity. The substitutent in position 8 of the naphthalene or the quinoline plane and the space orientation of these two planes contribute the most to the subtype selectivity on the basis of similar hydrogen bond interactions between binding ligands and both ER subtypes. The QSAR models built together with the docking procedure should be of great advantage for screening and designing ER ligands with improved affinity and subtype selectivity property.

This work is devoted to the prediction of a series of 208 structurally diverse PKCθ inhibitors using the Random Forest (RF) based on the Mold2 molecular descriptors. The RF model was established and identified as a robust predictor of the experimental pIC50 values, producing good external R2pred of 0.72, a standard error of prediction (SEP) of 0.45, for an external prediction set of 51 inhibitors which were not used in the development of QSAR models. By using the RF built-in measure of the relative importance of the descriptors, an important predictor—the number of group donor atoms for H-bonds (with N and O)―has been identified to play a crucial role in PKCθ inhibitory activity. We hope that the developed RF model will be helpful in the screening and prediction of novel unknown PKCθ inhibitory activity.

The effects of N-terminal (1–34 amino acids) and C-terminal (434–487 amino acids) amino acid sequences on the interfacial binding properties of Phospholipase D from Vibrio parahaemolyticus (VpPLD) were characterized by using monomolecular film technology. Online tools allowed the prediction of the secondary structure of the target N- and C-terminal VpPLD sequences. Various truncated forms of VpPLD with different N- or C-terminal deletions were designed, based on their secondary structure, and their membrane binding properties were examined. The analysis of the maximum insertion pressure (MIP) and synergy factor “a” indicated that the loop structure (1–25 amino acids) in the N-terminal segment of VpPLD had a positive effect on the binding of VpPLD to phospholipid monolayers, especially to 1,2-dimyristoyl-sn-glycero-3-phosphoserine and 1,2-dimyristoyl-sn-glycero-3-phosphocholine. The deletion affecting the N-terminus loop structure caused a significant decrease of the MIP and synergy factor a of the protein for these phospholipid monolayers. Conversely, the deletion of the helix structure (26–34 amino acids) basically had no influence on the binding of VpPLD to phospholipid monolayers. The deletion of the C-terminal amino acids 434–487 did not significantly change the binding selectivity of VpPLD for the various phospholipid monolayer tested here. However, a significant increase of the MIP value for all the phospholipid monolayers strongly indicated that the three-strand segment (434–469 amino acids) had a great negative effect on the interfacial binding to these phospholipid monolayers. The deletion of this peptide caused a significantly greater insertion of the protein into the phospholipid monolayers examined. The present study provides detailed information on the effect of the N- and C-terminal segments of VpPLD on the interfacial binding properties of the enzyme and improves our understanding of the interactions between this enzyme and cell membranes.

Nicotine and a variety of other drugs and toxins are metabolized by cytochromeP450 (CYP) 2A6. The aim of the present study was to build a quantitative structure-activityrelationship (QSAR) model to predict the activities of nicotine analogues on CYP2A6.Kernel partial least squares (K-PLS) regression was employed with the electro-topologicaldescriptors to build the computational models. Both the internal and external predictabilitiesof the models were evaluated with test sets to ensure their validity and reliability. As acomparison to K-PLS, a standard PLS algorithm was also applied on the same training andtest sets. Our results show that the K-PLS produced reasonable results that outperformed thePLS model on the datasets. The obtained K-PLS model will be helpful for the design ofnovel nicotine-like selective CYP2A6 inhibitors.

Using the classical emulsified system and the monomolecular film technique, the substrate specificity of recombinant Gibberella zeae lipase (rGZEL) that originates from Gibberella zeae was characterized in detail. Under the emulsified reaction system, both phospholipase and glycolipid hydrolytic activities were observed, except for the predominant lipase activity. The optimum conditions for different activity exhibition were also determined. Compared with its lipase activity, a little higher ratio of glycolipid hydrolytic activity (0.06) than phospholipase activity (0.02) was found. rGZEL preferred medium chain-length triglycerides, while lower activity was found for the longer-chain triglyceride. Using the monomolecular film technique, we found that the preference order of rGZEL to different phospholipids was 1,2-diacyl-sn-glycero-3-phospho-l-serine (PS) > 1,2-dioleoyl-sn-glycero-3-phospho-rac-(1-glycerol) sodium salt (PG) > 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) > l-α-phosphatidylinositol (PI) > cardiolipin (CL) > 3-sn-phosphatidic acid sodium salt (PA) > l-α-phosphatidylethanolamine (PE), while no hydrolytic activity was detected for sphingomyelin (SM). Moreover, rGZEL showed higher galactolipase activity on 1,2-distearoyimonoglactosylglyceride (MGDG). A kinetic study on the stereo- and regioselectivity of rGZEL was also performed by using three pairs of pseudodiglyceride enantiomers (DDGs). rGZEL presented higher preference for distal DDG enantiomers than adjacent ester groups, however, no hydrolytic activity to the sn-2 position of diglyceride analogs was found. Furthermore, rGZEL preferred the R configuration of DDG enantiomers. Molecular docking results were in concordance with in vitro tests.

Glycerophosphodiester phosphodiesterases (GDPD) are enzymes which degrade various glycerophosphodiesters to produce glycerol-3-phosphate and the corresponding alcohol moiety. Apart from this, a very interesting finding is that this enzyme could be used in the degradation of toxic organophosphorus esters, which has resulted in much attention on the biochemical and application research of GDPDs. In the present study, a novel GDPD from Pyrococcus furiosus DSM 3638 (pfGDPD) was successfully expressed in Escherichia coli and biochemically characterized. This enzyme hydrolyzed bis(p-nitrophenyl) phosphate, one substrate analogue of organophosphorus diester, with an optimal reaction temperature 55 °C and pH 8.5. The activity of pfGDPD was strongly dependent on existing of bivalent cations. It was strongly stimulated by Mn2+ ions, next was Co2+ and Ni2+ ions. Further investigations were conducted on its substrate selectivity towards different phospholipids. The results indicated that except of glycerophosphorylcholine (GPC), this enzyme also possessed lysophospholipase D activity toward both sn1-lysophosphatidylcholine (1-LPC) and sn2-lysophosphatidylcholine (2-LPC). Higher activity was found for 1-LPC than 2-LPC; however, no hydrolytic activity was found for phosphatidylcholine (PC). Molecular docking based on the 3D-modeled structure of pfGDPD was conducted in order to provide a structural foundation for the substrate selectivity.

Drug discovery is a risky, costly and time-consuming process depending on multidisciplinary methods to create safe and effective medicines. Although considerable progress has been made by high-throughput screening methods in drug design, the cost of developing contemporary approved drugs did not match that in the past decade. The major reason is the late-stage clinical failures in Phases II and III because of the complicated interactions between drug-specific, human body and environmental aspects affecting the safety and efficacy of a drug. There is a growing hope that systems-level consideration may provide a new perspective to overcome such current difficulties of drug discovery and development. The systems pharmacology method emerged as a holistic approach and has attracted more and more attention recently. The applications of systems pharmacology not only provide the pharmacodynamic evaluation and target identification of drug molecules, but also give a systems-level of understanding the interaction mechanism between drugs and complex disease. Therefore, the present review is an attempt to introduce how holistic systems pharmacology that integrated in silico ADME/T (i.e., absorption, distribution, metabolism, excretion and toxicity), target fishing and network pharmacology facilitates the discovery of small molecular drugs at the system level.