Simulated molecular evolution in a full combinatorial library

BACKGROUND: The Darwinian concept of 'survival of the fittest' has inspired the development of evolutionary optimization methods to find molecules with desired properties in iterative feedback cycles of synthesis and testing. These methods have recently been applied to the computer-guided heuristic selection of molecules that bind with high affinity to a given biological target. We describe the optimization behavior and performance of genetic algorithms (GAs) that select molecules from a combinatorial library of potential thrombin inhibitors in 'artificial molecular evolution' experiments, on the basis of biological screening results. RESULTS: A full combinatorial library of 15,360 members structurally biased towards the serine protease thrombin was synthesized, and all were tested for their ability to inhibit the protease activity of thrombin. Using the resulting large structure-activity landscape, we simulated the evolutionary selection of potent thrombin inhibitors from this library using GAs. Optimal parameter sets were found (encoding strategy, population size, mutation and cross-over rate) for this artificial molecular evolution. CONCLUSIONS: A GA-based evolutionary selection is a valuable combinatorial optimization strategy to discover compounds with desired properties without needing to synthesize and test all possible combinations (i.e. all molecules). GAs are especially powerful when dealing with very large combinatorial libraries for which synthesis and screening of all members is not possible and/or when only a small number of compounds compared with the library size can be synthesized or tested. The optimization gradient or 'learning' per individual increases when using smaller population sizes and decreases for higher mutation rates.     

Inhibitors of prolyl endopeptidase: Characterization of the pharmacophoric pattern using conformational analysis and 3D-QSAR

Summary A structure-activity study has been carried out on several compounds known as inhibitors of the serine protease prolyl endopeptidase. Conformational analysis has been done using different molecular mechanics methods such as molecular dynamics, or a randomized conformational search method. The conformers obtained were classified using geometric and energetic criteria. A pattern recognition analysis was done in order to divide conformers according to families. The resulting dominant families, for all compounds investigated, showed very similar geometric features. Based on the lowest energy conformers obtained after randomized conformational analysis, a 3D-QSAR model was established using the CoMFA approach. The validity of this model was verified by prediciting correctly the activity of other molecules not used in the construction of this model.     

Identification of the new chymotrypsin inhibitor micropeptin 996 by metabolomics-guided analysis

An untargeted metabolomics approach was used to investigate a cultured strain of Microcystis aeruginosa (UTEX LB2386) known to be a prolific producer of a diverse class of cyanopeptides. Identification of a putative new compound with a molecular weight of 996 led to the purification and structure elucidation of this new member of the micropeptin class of cyanopeptides. Micropeptin 996 displayed potent inhibition of the serine protease enzyme chymotrypisin relative to structurally related members of this class.     

Purification and characterization of Bacillus cereus protease suitable for detergent industry

An extracellular alkaline protease from an alkalophilic bacterium, Bacillus cereus, was produced in a large amount by the method of extractive fermentation. The protease is thermostable, pH tolerant, and compatible with commercial laundry detergerts. The protease purified and characterized in this study was found to be saperior to endogenous protease already present in commercial laundry detergents. The enzyme was purified to homogeneity by ammonium sulfate precipitation, concentration by ultrafiltration, anionexchange chromatography, and gel filtration. The purified enzyme had a specific activity of 3256.05 U/mg and was found to be amonomeric protein with a molecular mass of 28 and 31 kDa, as estimated by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and nondenaturing PAGE, respectively. Its maximum protease activity against casein was found to be at pH 10.5 and 50°C. Proteolytic activity of the enzyme was detected by casein and gelatin zymography, which gave a very clear protease activity zone on gel that corresponded to the band obtained on SDS-PAGE and nondenaturing PAGE with a molecular mass of nearly 31 kDa. The purified enzyme was analyzed through matrix-assisted laser desorption ionization-time-of-flight-mass spectrometry (MALDI-TOF-MS) and identified as a subtilisin class of protease. Specific serine protease inhibitors, suggesting the presence of serine residues at the active site, inhibited the enzme significantly.     

Substrate activity screening: a fragment-based method for the rapid identification of nonpeptidic protease inhibitors

A new fragment-based method for the rapid development of novel and distinct classes of nonpeptidic protease inhibitors, Substrate Activity Screening (SAS), is described. This method consists of three steps: (1) a library of N-acyl aminocoumarins with diverse, low molecular weight N-acyl groups is screened to identify protease substrates using a simple fluorescence-based assay, (2) the identified N-acyl aminocoumarin substrates are optimized by rapid analogue synthesis and evaluation, and (3) the optimized substrates are converted to inhibitors by direct replacement of the aminocoumarin with known mechanism-based pharmacophores. The SAS method was successfully applied to the cysteine protease cathepsin S, which is implicated in autoimmune diseases. Multiple distinct classes of nonpeptidic substrates were identified upon screening an N-acyl aminocoumarin library. Two of the nonpeptidic substrate classes were optimized to substrates with >8000-fold improvements in cleavage efficiency for each class. Select nonpeptidic substrates were then directly converted to low molecular weight, novel aldehyde inhibitors with nanomolar affinity to cathepsin S. This study demonstrates the unique characteristics and merits of this first substrate-based method for the rapid identification and optimization of weak fragments and provides the framework for the development of completely nonpeptidic inhibitors to many different proteases.     

Synthesis of novel bridged dinitrogen heterocycles and their evaluation as potential fragments for the design of biologically active compounds

A library of saturated bridged heterocycles based on 3,6-diazabicyclo[3.2.1]octane-2,4-dione and bispidine scaffolds (mean compound molecular weight is approximately 300 Da) with up to three stereocenters and four diversity points has been synthesized. Synthetic scaffold modifications leading to an increase in molecular complexity were studied. Well-defined stereochemical structures of both compound sets was confirmed by X-ray studies and halogenoaryl substituents were inserted appropriately for the design of novel non-basic serine protease inhibitors. Comprehensive molecular modeling has been performed for all synthesized compounds giving rationales of ligand–enzyme interactions with thrombin and trypsin. Biological testing confirmed moderate inhibitory activity of halogen-substituted saturated diazabicyclic small molecules towards thrombin.     

0.3 nm分辨率绿豆胰蛋白酶抑制剂Lys活力碎片—牛胰蛋白酶(MBILF-BTRY)复合物的晶体结构分析和分子模型

Crystal structure of MBILF-BTRY complex was successfully solved by molecular replacement method and its 3-dimensional molecular model was thereby derived. Mung bean trypsin inhibitor belongs to Bowman-Birk inhibitor group, which is by far the most complicated among the ten fundamental groups of serine protease inhibitor. Neither the 3-dimensional structure of Bowman-Birk inhibitor group nor the stereoscopic conformation of its complex with protease has ever been reported. The crystallographic data of MBILF-BTRY complex are found to be a=....     

Inclusion of conserved buried water molecules in the model structure of rat submaxillary kallikrein

A new approach to the molecular modelling of homologous serine proteases isadopted, by including a set of 21 buried waters known to be preserved inenzymes sharing the primary specificity of trypsin, in the homology modellingof rat submaxillary gland kallikrein. Buried waters – water moleculessequestered from bulk solvent within a protein matrix – appear to beintegral conserved components of all serine proteases of known structure andshould be incorporated into serine protease models built on the basis ofsequence/structural homology to this family. The absence of such waters mightinduce errors in a force field simulation, favouring the formation ofnonexistent hydrogen bonds and locally inaccurate structure. The kallikreinmodel refinement has led to the conclusion that an additional buried watershould be added to the original rigid matrix of 21 conserved water molecules.The structurally preserved protein cavities of such waters validate themodelled structure.     

'Heat-treatment aqueous two phase system' for purification of serine protease from Kesinai (Streblus asper) leaves

A 'Heat treatment aqueous two phase system' was employed for the first time to purify serine protease from kesinai (Streblus asper) leaves. In this study, introduction of heat treatment procedure in serine protease purification was investigated. In addition, the effects of different molecular weights of polyethylene glycol (PEG 4000, 6000 and 8000) at concentrations of 8, 16 and 21% (w/w) as well as salts (Na-citrate, MgSO? and K?HPO?) at concentrations of 12, 15, 18% (w/w) on serine protease partition behavior were studied. Optimum conditions for serine protease purification were achieved in the PEG-rich phase with composition of 16% PEG6000-15% MgSO?. Also, thermal treatment of kesinai leaves at 55 °C for 15 min resulted in higher purity and recovery yield compared to the non-heat treatment sample. Furthermore, this study investigated the effects of various concentrations of NaCl addition (2, 4, 6 and 8% w/w) and different pH (4, 7 and 9) on the optimization of the system to obtain high yields of the enzyme. The recovery of serine protease was significantly enhanced in the presence of 4% (w/w) of NaCl at pH 7.0. Based on this system, the purification factor was increased 14.4 fold and achieved a high yield of 96.7%.     

Binding of alkaloids into the S1 specificity pocket of α-chymotrypsin: evidence from induced circular dichroism spectra

Non-covalent binding of planar aromatic molecules into the S1 specificity pocket of the serine protease α-chymotrypsin (αCHT) can be detected by measuring induced circular dichroism (CD) spectroscopic signals. Utilizing this phenomenon, αCHT association of proflavine (PRF), the well known serine protease inhibitor has been investigated together with plant-derived compounds including isoquinoline, pyridocarbazole and indoloquinoline alkaloids, of which αCHT binding has never been reported. Non-degenerate exciton coupling between π-π* transitions of the ligand molecules and two tryptophan residues (Trp172 and Trp215) near to the binding site is proposed to be responsible for the induced CD activity. The association constants calculated from CD titration data indicated strong αCHT association of sanguninarine, ellipticine, desmethyl-isocryptolepine and isoneocryptolepine (K(a) ≈ 10(5) M(-1)) while berberine, coptisine and chelerythrine bind to the enzyme with lower, PRF-like affinity (K(a) ≈ 10(4) M(-1)). PRF-trypsin and ellipticine-trypsin binding interactions have also been demonstrated. The binding of the alkaloids into the S1 pocket of αCHT has been confirmed by CD competition experiments. Molecular docking calculations showed the inclusion of PRF as well as the alkaloid molecules in the S1 cavity where they are stabilized by hydrophobic and H-bonding interactions. These novel nonpeptidic scaffolds can be used for developing selective inhibitors of serine proteases having chymotrypsin-like folds. Furthermore, the results provide a novel, CD spectroscopic based approach for probing the ligand binding of αCHT and related proteases.     

Field-Template,QSAR, Ensemble Molecular Docking,and 3D-RISM Solvation Studies Expose Potential of FDA-Approved Marine Drugs as SARS-CoVID-2 Main Protease Inhibitors

Currently, SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) has infected people among all countries and is a pandemic as declared by the World Health Organization (WHO). SARS-CoVID-2 main protease is one of the therapeutic drug targets that has been shown to reduce virus replication, and its high-resolution 3D structures in complex with inhibitors have been solved. Previously, we had demonstrated the potential of natural compounds such as serine protease inhibitors eventually leading us to hypothesize that FDA-approved marine drugs have the potential to inhibit the biological activity of SARS-CoV-2 main protease. Initially, field-template and structure–activity atlas models were constructed to understand and explain the molecular features responsible for SARS-CoVID-2 main protease inhibitors, which revealed that Eribulin Mesylate, Plitidepsin, and Trabectedin possess similar characteristics related to SARS-CoVID-2 main protease inhibitors. Later, protein–ligand interactions are studied using ensemble molecular-docking simulations that revealed that marine drugs bind at the active site of the main protease. The three-dimensional reference interaction site model (3D-RISM) studies show that marine drugs displace water molecules at the active site, and interactions observed are favorable. These computational studies eventually paved an interest in further in vitro studies. Finally, these findings are new and indeed provide insights into the role of FDA-approved marine drugs, which are already in clinical use for cancer treatment as a potential alternative to prevent and treat infected people with SARS-CoV-2.     

Ab initio quantum mechanics-based free energy perturbation method for calculating relative solvation free energies

A free energy perturbation (FEP) method was developed that uses ab initio quantum mechanics (QM) for treating the solute molecules and molecular mechanics (MM) for treating the surroundings. Like our earlier results using AM1 semi empirical QMs, the ab initio QM/MM-based FEP method was shown to accurately calculate relative solvation free energies for a diverse set of small molecules that differ significantly in structure, aromaticity, hydrogen bonding potential, and electron density. Accuracy was similar to or better than conventional FEP methods. The QM/MM-based methods eliminate the need for time-consuming development of MM force field parameters, which are frequently required for drug-like molecules containing structural motifs not adequately described by MM. Future automation of the method and parallelization of the code for Linux 128/256/512 clusters is expected to enhance the speed and increase its use for drug design and lead optimization.     

Structure and function of complement activating enzyme complexes: C1 and MBL-MASPs

The complement system is a major effector arm of the immune defense contributing to the destruction of invading pathogens. There are three possible routes of complement cascade activation: the classical, the alternative and the lectin pathways. The activation of the classical and lectin pathways is initiated by supramolecular complexes, which resemble each other. Each complex has a recognition subunit (C1q in the classical and mannose-binding lectin (MBL) in the lectin pathway), which associates with serine protease zymogens (C1q with C1r and C1s, and MBL with MBL-associated serine proteases: MASP-1, MASP-2) to form the C1 and MBL-MASPs complexes, respectively. As the recognition subunits bind to activator structures, subsequent activation of the serine protease zymogens occurs. The precise structure of the complexes and the exact mechanism of their activation have not been solved, yet. In this review we summarize the recent advances about the structure and function of the individual subcomponents of both complexes achieved by genetic engineering, molecular modeling, physico-chemical and functional studies. Special emphasis will be laid on the serine proteases: the role of the individual domains in the assembly of the C1s-C1r-C1r-C1s tetramer and in the control of the protease activity will be discussed. We will then focus on recent functional models of the supramolecular complexes. The question of how a non-enzymatic signal (the binding of C1q or MBL to activators) can be converted into enzymatic events (activation of serine protease zymogens) will be addressed. The similarities and differences between C1 and MBL-MASPs will also be discussed.     

Large-scale motions and electrostatic properties of furin and HIV-1 protease

We present a comparative study between two members of serine and aspartic proteases complexed with a peptide substrate. The same computational setup is used to characterize the structural, electrostatic, and electronic properties for the Michaelis complex of furin, a serine protease, and of the aspartic protease from HIV-1. In both cases plane-wave density functional theory (PW-DFT) and empirical force-field-based molecular dynamics calculations are used. For furin, calculations are extended to the complex with the intermediate of the first step of the reaction. Comparisons are also made with results from recent PW-DFT investigations on both families of enzymes and with the same chemical groups in an aqueous environment. It is found that the substrate carbonyl group is more polarized in the furin complex than in the HIV-1 protease one. A further difference regards the large-scale motions of the complexes as a whole and local conformational fluctuations at the active site. The global and local fluctuations are well coupled for HIV-1 protease but not for furin. Thus, despite some chemical analogies in the first step of the reaction mechanism, furin and HIV-1 protease complexes appear to be characterized by a different interplay of electrostatics and conformational fluctuations.     

Optimization of the conditions for extraction of serine protease from kesinai plant (Streblus asper) leaves using response surface methodology

Response surface methodology (RSM) using a central composite design (CCD) was employed to optimize the conditions for extraction of serine protease from kesinai (Streblus asper) leaves. The effect of independent variables, namely temperature (42.5,47.5, X?), mixing time (2-6 min, X?), buffer content (0-80 mL, X?) and buffer pH (4.5-10.5, X?) on specific activity, storage stability, temperature and oxidizing agent stability of serine protease from kesinai leaves was investigated. The study demonstrated that use of the optimum temperature, mixing time, buffer content and buffer pH conditions protected serine protease during extraction, as demonstrated by low activity loss. It was found that the interaction effect of mixing time and buffer content improved the serine protease stability, and the buffer pH had the most significant effect on the specific activity of the enzyme. The most desirable conditions of 2.5 °C temperature, 4 min mixing time, 40 mL buffer at pH 7.5 was established for serine protease extraction from kesinai leaves.     

Conformational study on telaprevir by HPLC–DAD–MS and theoretical calculation

Telaprevir is a potent, selective, peptidomimetic inhibitor of the hepatitis C virus (HCV) NS3‐4A serine protease. it is used for the treatment of HCV infection in combination with peginterferon alfa and ribavirin. In the present work, the E–Z isomerization process of telaprevir in solution was revealed by online HPLC–DAD (diode array detector)–MS, variable‐temperature and variable‐gradient experiments. The molecular geometry information of the two isomers was established by molecular mechanics calculations, and good correlation between the two isomers' UV–vis spectra and their molecular geometry information was also discovered. In addition, it was revealed by molecular docking that the two isomers have different affinities to HCV NS3?4A protease, and the Z isomer, the minor form of telaprevir in solution, is the more effective inhibitor of HCV NS3?4A protease. The investigation can provide more structure information about telaprevir in solution and in the binding process of HCV NS3?4A protease.     

Serine protease acylation proceeds with a subtle re-orientation of the histidine ring at the tetrahedral intermediate

The acylation mechanism of a prototypical serine protease trypsin and its complete free energy reaction profile have been determined by Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling.     

Fast 3D molecular superposition and similarity search in databases of flexible molecules

We present a new method (fFLASH) for the virtual screening of compound databases that is based on explicit three-dimensional molecular superpositions. fFLASH takes the torsional flexibility of the database molecules fully into account, and can deal with an arbitrary number of conformation-dependent molecular features. The method utilizes a fragmentation-reassembly approach which allows for an efficient sampling of the conformational space. A fast clique-based pattern matching algorithm generates alignments of pairs of adjacent molecular fragments on the rigid query molecule that are subsequently reassembled to complete database molecules. Using conventional molecular features (hydrogen bond donors and acceptors, charges, and hydrophobic groups) we show that fFLASH is able to rapidly produce accurate alignments of medium-sized drug-like molecules. Experiments with a test database containing a diverse set of 1780 drug-like molecules (including all conformers) have shown that average query processing times of the order of 0.1 seconds per molecule can be achieved on a PC.