Computational Design of Novel Allosteric Inhibitors for Plasmodium falciparum DegP

The serine protease, DegP exhibits proteolytic and chaperone activities, essential for cellular protein quality control and normal cell development in eukaryotes. The P. falciparum DegP is essential for the parasite survival and required to combat the oscillating thermal stress conditions during the infection, protein quality checks and protein homeostasis in the extra-cytoplasmic compartments, thereby establishing it as a potential target for drug development against malaria. Previous studies have shown that diisopropyl fluorophosphate (DFP) and the peptide SPMFKGV inhibit E. coli DegP protease activity. To identify novel potential inhibitors specific to PfDegP allosteric and the catalytic binding sites, we performed a high throughput in silico screening using Malaria Box, Pathogen Box, Maybridge library, ChEMBL library and the library of FDA approved compounds. The screening helped identify five best binders that showed high affinity to PfDegP allosteric (T0873, T2823, T2801, {"type":"entrez-protein","attrs":{"text":"RJC02337","term_id":"1480409465","term_text":"RJC02337"}}RJC02337, CD00811) and the catalytic binding site (T0078L, T1524, T2328, BTB11534 and 552691). Further, molecular dynamics simulation analysis revealed {"type":"entrez-protein","attrs":{"text":"RJC02337","term_id":"1480409465","term_text":"RJC02337"}}RJC02337, BTB11534 as the best hits forming a stable complex. WaterMap and electrostatic complementarity were used to evaluate the novel bio-isosteric chemotypes of {"type":"entrez-protein","attrs":{"text":"RJC02337","term_id":"1480409465","term_text":"RJC02337"}}RJC02337, that led to the identification of 231 chemotypes that exhibited better binding affinity. Further analysis of the top 5 chemotypes, based on better binding affinity, revealed that the addition of electron donors like nitrogen and sulphur to the side chains of butanoate group are more favoured than the backbone of butanoate group. In a nutshell, the present study helps identify novel, potent and Plasmodium specific inhibitors, using high throughput in silico screening and bio-isosteric replacement, which may be experimentally validated.     

Ratio of Protein Fractions in Cereal Grains Cultivated in Uzbekistan

The glutelin fraction was used as the food value index of local corn. The base-soluble fraction was used as a test of the glassiness of local strains of rice grain. The content of gliadins and glutelins did not change over 24 h after wheat grain cultivated in Uzbekistan was moistened.     

The search for a new model structure of β-Factor XIIa

We present the search for a new model of -factor XIIa, a blood coagulation enzyme, with an unknown experimental 3D-structure. We decided to build not one but three different models using different homologous proteins as well as different techniques and different modellers. Additional studies, including extensive molecular dynamics simulations on the solvated state, allowed us to draw several conclusions concerning homology modelling, in general, and -factor XIIa, in particular.     

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.     

Modeling of serine protease prototype reactions with the flexible effective fragment potential quantum mechanical/molecular mechanical method

A complete cycle of chemical transformations for the serine protease prototype reaction is modeled following calculations with the flexible effective fragment quantum mechanical/molecular mechanical (QM/MM) method. The initial molecular model is based on the crystal structure of the trypsin–bovine pancreatic trypsin inhibitor complex including all atoms of the enzyme within approximately 15–18 Å of the oxygen center O of the catalytic serine residue. Several selections of the QM/MM partitioning are considered. Fractions of the side chains of the residues from the catalytic triad (serine, histidine and aspartic acid) and a central part of a model substrate around the C–N bond to be cleaved are included into the QM subsystem. The remaining part, or the MM subsystem, is represented by flexible chains of small effective fragments, whose potentials explicitly contribute to the Hamiltonian of the QM part, but the corresponding fragment–fragment interactions are described by the MM force fields. The QM/MM boundaries are extended over the C–C bonds of the peptides assigned to the QM subsystem in the enzyme, C–C and C–N bonds in model substrates. Multiple geometry optimizations have been performed by using the RHF/6-31G method in the QM part and OPLSAA or AMBER sets of MM parameters, resulting in a series of stationary points on the complex potential-energy surfaces. All structures generally accepted for the serine protease catalytic cycle have been located. Energies at the stationary points found have been recomputed at the MP2/6-31+G^* level for the QM part in the protein environment. Structural changes along the reaction path are analyzed with special attention to hydrogen-bonding networks. In the case of a model substrate selected as a short peptide CH₃(NHCO-CH₂)₂ – HN–CO–(CH₂–NHCO)CH₃ the computed energy profile for the acylation step shows too high activation energy barriers. The energetics of this rate-limiting step is considerably improved, if more realistic model for the substrate is considered, following the motifs of the ThrI11–GlyI12–ProI13-–CysI14–LysI15–AlaI16–ArgI17–IleI18–IleI19 sequence of the bovine pancreatic trypsin inhibitor.     

Rapid quantification of HIV protease inhibitors in human plasma by high-performance liquid chromatography coupled with electrospray ionization tandem mass spectrometry

HIV protease inhibitors are important antiretroviral drugs which have substantially reduced the morbidity and mortality associated with HIV-1 infection. Recent data have shown relationships between plasma concentrations of the protease inhibitors and clinical response, which makes therapeutic drug monitoring valuable. We have developed and validated an assay, using liquid chromatography coupled with electrospray tandem mass spectrometry (LC/MS/MS), for the routine quantification of the six licensed protease inhibitors (amprenavir, indinavir, lopinavir, nelfinavir, ritonavir and saquinavir) and the pharmacologically active nelfinavir metabolite M8 in plasma. The sample pretreatment consisted of protein precipitation with a mixture of methanol and acetronitrile using only 100 microl of plasma. Chromatographic separation was performed on an Inertsil ODS3 column (50 x 2.0 mm i.d., particle size 5 microm), with a quick stepwise gradient using an acetate buffer (pH 5) and methanol, at a flow rate of 0.5 ml min(-1). The analytical run time was 5.5 min. The use of a 96-well plate autosampler allowed batch sizes up to 150 patient samples. The triple-quadrupole mass spectrometer was operated in the positive ion mode and multiple reaction monitoring was used for drug quantification. The method was validated over the concentration ranges 0.01-10 microg ml(-1) for indinavir and saquinavir, 0.1-10 microg ml(-1) for amprenavir, 0.05-10 microg ml(-1) for nelfinavir and ritonavir, 0.1-20 microg ml(-1) for lopinavir and 0.01-5 microg ml(-1) for M8. Saquinavir-d(5) and indinavir-d(6) were used as internal standards. The coefficients of variation were always <10% for both intra-day and inter-day precisions for each compound. Mean accuracies were also between the designated limits (+/-15%). The validated concentration ranges proved to be adequate in daily practice. This robust and fast LC/MS/MS assay is now successfully applied for routine therapeutic drug monitoring and pharmacokinetic studies in our hospital.     

有机溶剂性质及其水含量对酶催化合成含D－氨基酸残基的二肽衍生物的影响

有机溶剂性质及其水含量对酶催化合成保护的寡肽及其衍生物有显著影响。系 统地研究了有机溶剂性质及在不同有机溶剂中水含量对α－胰凝乳蛋白酶催化合成 含D－氨基酸残基的二肽衍生物产率的影响。用α－胰凝乳蛋白酶、枯草杆菌蛋白 酶和嗜热杆菌蛋白酶在有机溶剂中催化合成了一系列保护的二肽及其衍生物，并研 究了水含量对反应产率的影响，得到了十分有意义的结果。     

Studies on alkaline protease produced by Bacillus sp. NG312

An alkalophilic hyperproducer of alkaline protease, Bacillus sp. NG312, was isolated, and the enzyme showed maximum activity at pH 11.0 and 60°C. The temperature optimum was increased by 10°C in presence of Ca²⁺. The crudeenzyme was found to have half-life of 11 d at 37°C and maximum stability at pH 9.0–10.0. It also exhibited very good stability in presence of detergent components and some locally available commericial detergent powders.     

An approach to sequence-specific antibody proteases

We describe here a novel strategy for the isolation of antibodies with sequence-specific protease activity: the synthesis of dipeptide haptens in which the targeted peptide bond has been replaced by a ring-strained or torsionally strained hydroxyethylene transition-state analog. Thus, the analogs mimic both a peptide bond in a distorted, reactive conformation and the transition state for peptide bond hydrolysis. In order to obtain sequence-specific antibody proteases, these analogs have been flanked with additional amino acid residues in preparation for immunization. In particular, we have synthesized peptides containing analogs such as 2-cis-amino-3-cis-hydroxycyclobutane carboxylic acid andendo-(3-amino-2-hydroxy)bicyclo[2.2.1]heptane-7-anti-carboxylic acid. We have also prepared a series of peptide derivatives containing analogs, such as 2-[3-amino-2-oxo-1-azetidinyl]-3-methylbutanoic acid, in which the targeted peptide bond has been incorporated into a β-lactam ring. Since the “peptide bond” has been left intact, these species mimic only a distorted ground state. At present, antibodies are being elicited against a number of the above peptide derivatives.     

NMR chemical shifts in the low–pH form of a-chymotrypsin. A QM/MM and ONIOM–NMR study

 The relationship between hydrogen bonding and NMR chemical shifts in the catalytic triad of low-pH α-chymotrypsin is investigated by combined use of the effective fragment potential [(2001) J Phys Chem A 105:293] and ONIOM–NMR [(2000) Chem Phys Lett 317:589] methods. Our study shows that while the His57 N^δ1−H bond is stretched by a relatively modest amount (to about 1.060 ?) this lengthening, combined with the polarization due to the molecular environment, is sufficient to explain the experimentally observed chemical shifts of 18.2 ppm. Furthermore, the unusual down-field shift of H^ɛ1 (9.2 ppm) observed experimentally is reproduced and shown to be induced by interactions with the C=O group of Ser214 as previously postulated. The free-energy cost of moving H^δ1 from His57 to Asp102 is predicted to be 5.5 kcal/mol. Received: 26 September 2001 / Accepted: 6 September 2002 / Published online: 21 January 2003 Contribution to the Proceedings of the Symposium on Combined QM/MM Methods at the 222nd National Meeting of the American Chemical Society, 2001 Correspondence to: J. H. Jensen e-mail: jan-jensen@uiowa.edu Acknowledgements. This work was supported by a Research Innovation Award from the Research Corporation and a type G starter grant from the Petroleum Research Fund. The calculations were performed on IBM RS/6000 workstations obtained through a CRIF grant from the NSF (CHE-9974502) and on supercomputers at the National Center for Supercomputer Applications at Urbana-Champaign. The authors are indebted to Visvaldas Kairys for help with the CHARMM program, and to Daniel Quinn for many helpful discussions.