Modeling of substrate and inhibitory complexes of histidine-aspartic protease

A three-dimensional structure of histo-aspartic protease (HAP), a pepsin-like enzyme from the causative agent of malaria Plasmodium falciparum, is suggested on the basis of homologous modeling followed by equilibration by the method of molecular dynamics. The presence of a His residue in the catalytic site instead of an Asp residue, which is characteristic of pepsin-like enzymes, and replacement of some other conserved residues in the active site make it possible for the enzyme to function by the covalent mechanism inherent in serine proteases. The detailed structures of HAP complexes with pepstatin, a noncovalent inhibitor of aspartic proteases, and phenylmethylsulfonyl fluoride, a covalent inhibitor of serine proteases, as well as with a pentapeptide substrate are discussed.     

Mechanistic role of NS4A and substrate in the activation of HCV NS3 protease

Hepatitis C virus (HCV) NS3 protease is the key enzyme for its maturation. Three hypotheses have been advanced in the literature to demonstrate the mechanism of the activation of the HCV NS3 protease. A virus-encoded protein NS4A and substrate are proposed to be involved in the activation of the HCV NS3 protease. However, the three hypotheses are not completely consistent with one another. Multiple molecular dynamics simulations were performed on various NS3 protease systems: free NS3 protease, NS3/4A, NS3/inhibitor, and NS3/4A/inhibitor complexes, to further unravel the mechanism of the activation of the NS3 protease. Simulation results suggest that the binding of NS4A induces a classic serine protease conformation of the catalytic triad of the NS3 protease. NS4A rearranges the secondary structure of both the N-terminus and catalytic site of the NS3 protease, reduces the mobility of the global structure of the NS3 protease, especially the catalytic site, and provides a rigid and tight structure, except for the S1 pocket, for the binding and hydrolysis of substrates. The binding of substrate also contributes to the activation of the NS3 protease by an induced-fit of the classic serine protease catalytic triad. However, the global structure of the NS3 protease is still loose and highly flexible without stable secondary structural elements, such as helix α0 at the N-terminus and helix α1 and β-sheet E1-F1 at the catalytic site. The structure of the NS3 protease without NS4A is not suitable for the binding and hydrolysis of substrates.     

Analysis of the structure of HIV-1 protease complexed with a hexapeptide inhibitor. Part II: Molecular dynamic studies of the active site region

Six models of the catalytic site of HIV-1 protease complexed with a reduced peptide inhibitor, MVT-101, were investigated. These studies focused on the details of protonation of the active site, its total net charge and hydrogen bonding pattern, which was consistent with both the observed coplanar configuration of the acidic groups of the catalytic aspartates (Asp-25 and Asp-125) and the observed binding mode of the inhibitor. Molecular dynamic simulations using AMBER 4.0 indicated that the active site should be neutral. The planarity of the aspartate dyad may be due to the formation of two hydrogen bonds: one between the inner O_δ1oxygen atoms of the two catalytic aspartates and another between the O_δ2atom of Asp-125 and the nitrogen atom of the reduced peptide bond of the bound inhibitor. This would require two additional protonations, either of both aspartates, or of one Asp and the amido nitrogen atom of Nle-204. Our results favor the Asp-inhibitor protonation but the other one is not excluded. Implications of these findings for the mechanism of enzymatic catalysis are discussed. Dynamic properties of the hydrogen bond network in the active site and an analysis of the interaction energy between the inhibitor and the protease are presented. © 1997 Wiley-Liss, Inc.     

蛋白质结构与功能研究中的分子模拟技术

分子模拟技术为蛋白质的研究提供了一种崭新的手段,在理论上解决了结构预测和功能分析以及蛋白质工程实施方面所面临的难题。它在蛋白质的结构预测和模建工作中占有举足轻重的地位,实现了生物技术与计算机技术的完美结合。本文简要阐述了该技术的基本步骤和工作原理,并以目前应用最广的生物大分子领域的商品化分子模拟软件Accelrys公司基于Linux系统开发的InsightII为例,介绍了相关程序模块的功能和作用,同时结合该技术在蛋白质的结构预测和模建、结构与功能关系分析、分子设计等过程中的开发与应用,加以具体说明和展望。     

Internal dynamics and ionization states of the macrophage migration inhibitory factor: comparison between wild-type and mutant forms

The macrophage migration inhibitory factor (MIF) is a cytokine that shares a common structural architecture and catalytic strategy with three isomerases: 4-oxalocrotonate tautomerase, 5-carboxymethyl-2-hydroxymuconate isomerase, and D-dopachrome tautomerase. A highly conserved N-terminal proline acts as a base-acid during the proton transfer reaction catalyzed by these enzymes. Such unusual catalytic strategy appears to be possible only due to the N-terminal proline pK(a) shifted to 5.0-6.0 units. Mutations of this residue result in a significant decrease of the catalytic activity of MIF. Two hypotheses have been proposed to explain the catalytic inefficiency of MIF: the lower basicity of primary amines with regard to secondary ones and the increased flexibility resulting from the replacement of a proline by residues like glycine. To investigate that, we have performed molecular dynamics simulations of MIF wild-type and its mutant P1G, as well as calculated the protonation properties of several mutant forms. It was found that the N-terminal glycine does not show larger fluctuations compared to proline, but the former residue is more exposed to the solvent throughout the simulations. The apparent pK(a) of these residues displays very little change (as expected from the structural rigidity of MIF) and is not significantly affected by the surrounding ionizable residues. Instead, the hydrophobic character of the active site seems to be the main factor in determining the pKa of the N-terminal residue and the catalytic efficiency of MIF.     

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.     

Navigation inside a protease: substrate selection and product exit in the tricorn protease from Thermoplasma acidophilum

Chaperonins are multi-subunit double-ring complexes that mediate the folding of nascent or denatured proteins. Gene duplication has been a potent force in the evolution of chaperonins in Archaea. Here we show that gene conversion has also been an important factor. We utilized a novel maximum likehood-based phylogenetic method for scanning DNA sequence alignments for regions of anomalous phylogenetic signal, such as those affected by gene conversion. Our results suggest that in crenarchaeotes, where an ancient gene duplication producing alpha and beta subunits took place in the common ancestor of the Pyrodictium, Aeropyrum, Pyrobaculum and Sulfolobus lineages, multiple independent gene conversions have occurred between the alpha and beta genes independently in each of these groups. Significantly, the conversions have repeatedly homogenized the region of the gene encoding the substrate-binding domain. This suggests that while the alpha and beta subunits in crenarchaeotes share only 50-60% overall amino acid sequence identity, they do not possess distinct roles in the binding of substrate. Cryptic gene conversion between distantly related paralogs may be more common than is currently appreciated, and could be a significant factor in slowing the functional differentiation of proteins encoded by duplicate genes long after their duplication.     

Computational site‐directed mutagenesis studies of the role of the hydrophobic triad on substrate binding in cholesterol oxidase

Cholesterol oxidase (ChOx) is a flavoenzyme that oxidizes and isomerizes cholesterol (CHL) to form cholest‐4‐en‐3‐one. Molecular docking and molecular dynamics simulations were conducted to predict the binding interactions of CHL in the active site. Several key interactions (E361‐CHL, N485‐FAD, and H447‐CHL) were identified and which are likely to determine the correct positioning of CHL relative to flavin‐adenine dinucleotide (FAD). Binding of CHL also induced changes in key residues of the active site leading to the closure of the oxygen channel. A group of residues, Y107, F444, and Y446, known as the hydrophobic triad, are believed to affect the binding of CHL in the active site. Computational site‐directed mutagenesis of these residues revealed that their mutation affects the conformations of key residues in the active site, leading to non‐optimal binding of CHL and to changes in the structure of the oxygen channel, all of which are likely to reduce the catalytic efficiency of ChOx. Proteins 2017; 85:1645–1655. © 2017 Wiley Periodicals, Inc.     

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.     

Molecular dynamics of HIV-1 protease.

Molecular dynamics simulations have been carried out based on the GROMOS force field on the aspartyl protease (PR) of the human immunodeficiency virus HIV-1. The principal simulation treats the HIV-1 PR dimer and 6990 water molecules in a hexagonal prism cell under periodic boundary conditions and was carried out for a trajectory of 100 psec. Corresponding in vacuo simulations, i.e., treating the isolated protein without solvent, were carried out to study the influence of solvent on the simulation. The results indicate that including waters explicitly in the simulation results in a model considerably closer to the crystal structure than when solvent is neglected. Detailed conformational and helicoidal analysis was performed on the solvated form to determine the exact nature of the dynamical model and the exact points of agreement and disagreement with the crystal structure. The calculated dynamical model was further elucidated by means of studies of the time evolution of the cross-correlation coefficients for atomic displacements of the atoms comprising the protein backbone. The cross-correlation analysis revealed significant aspects of structure originating uniquely in the dynamical motions of the molecule. In particular, an unanticipated through-space, domain-domain correlation was found between the mobile flap region covering the active site and a remote regions of the structure, which collectively act somewhat like a molecular cantilever. The significance of these results is discussed with respect to the inactivation of the protease by site-specific mutagenesis, and in the design of inhibitors.     

Is histoaspartic protease a serine protease with a pepsin‐like fold?

The primary structure of the so-called histoaspartic protease from Plasmodium falciparum has a very high percentage of identity and homology with the pepsin-like enzyme plasmepsin II. A homology modeling approach was used to calculate the three-dimensional structure of the enzyme. Molecular dynamics (MD) simulations were applied to find those structural properties of the histoaspartic protease that had a tendency to remain stable during all runs. The results have shown that hydrogen-bonded residues Ser37-His34-Asp214 are arranged without any strain, in a manner that resembles the active site of a serine protease, while Ser38 and Asn39 take up positions appropriate to formation of an oxyanion hole. Although there are several important differences between the enzyme and plasmepsin II, all of the structural features associated with a typical pepsin-like aspartic protease are present in the final model of the histoaspartic protease. A possibility that this enzyme may function as a serine protease is discussed.     

Crystal structure of a complex of HIV-1 protease with a dihydroxyethylene-containing inhibitor: comparisons with molecular modeling.

The structure of a crystal complex of recombinant human immunodeficiency virus type 1 (HIV-1) protease with a peptide-mimetic inhibitor containing a dihydroxyethylene isostere insert replacing the scissile bond has been determined. The inhibitor is Noa-His-Hch psi [CH(OH)CH(OH)]Vam-Ile-Amp (U-75875), and its Ki for inhibition of the HIV-1 protease is < 1.0 nM (Noa = 1-naphthoxyacetyl, Hch = a hydroxy-modified form of cyclohexylalanine, Vam = a hydroxy-modified form of valine, Amp = 2-pyridylmethylamine). The structure of the complex has been refined to a crystallographic R factor of 0.169 at 2.0 A resolution by using restrained least-squares procedures. Root mean square deviations from ideality are 0.02 A and 2.4 degrees, for bond lengths and angles, respectively. The bound inhibitor diastereomer has the R configurations at both of the hydroxyl chiral carbon atoms. One of the diol hydroxyl groups is positioned such that it forms hydrogen bonds with both the active site aspartates, whereas the other interacts with only one of them. Comparison of this X-ray structure with a model-built structure of the inhibitor, published earlier, reveals similar positioning of the backbone atoms and of the side-chain atoms in the P2-P2' region, where the interaction with the protein is strongest. However, the X-ray structure and the model differ considerably in the location of the P3 and P3' end groups, and also in the positioning of the second of the two central hydroxyl groups. Reconstruction of the central portion of the model revealed the source of the hydroxyl discrepancy, which, when corrected, provided a P1-P1' geometry very close to that seen in the X-ray structure.     

Characterization of the mechanism of bile salt hydrolase substrate specificity by experimental and computational analyses

1. Download : Download high-res image (282KB)
2. Download : Download full-size image

     

The structural basis for catalysis and substrate specificity of a rhomboid protease

Rhomboids are intramembrane proteases that use a catalytic dyad of serine and histidine for proteolysis. They are conserved in both prokaryotes and eukaryotes and regulate cellular processes as diverse as intercellular signalling, parasitic invasion of host cells, and mitochondrial morphology. Their widespread biological significance and consequent medical potential provides a strong incentive to understand the mechanism of these unusual enzymes for identification of specific inhibitors. In this study, we describe the structure of Escherichia coli rhomboid GlpG covalently bound to a mechanism‐based isocoumarin inhibitor. We identify the position of the oxyanion hole, and the S₁‐ and S₂′‐binding subsites of GlpG, which are the key determinants of substrate specificity. The inhibitor‐bound structure suggests that subtle structural change is sufficient for catalysis, as opposed to large changes proposed from previous structures of unliganded GlpG. Using bound inhibitor as a template, we present a model for substrate binding at the active site and biochemically test its validity. This study provides a foundation for a structural explanation of rhomboid specificity and mechanism, and for inhibitor design.     

Clp protease complexes and their diversity in chloroplasts

The Clp proteases represent a large, ancient ATP-dependent protease family which in higher plants is known to be located in chloroplasts. The soluble, presumably multisubunit, enzyme of the organelle stroma is of dual genetic origin. It consists of a nuclear-encoded, regulatory subunit ClpC, which is an ATPase, and a plastid-encoded proteolytic subunit ClpP, which is a serine protease. An additional, nuclear-encoded proteolytic subunit resembling ClpP has been recently reported from tomato (Schaller and Ryan, 1995 plant gene Register 95–00). We demonstrate that in both tomato Lycopersicon esculentum Mill. and Arabidopsis thaliana, (L.) Heynh. the nuclear-encoded ClpP (nClpP) is made as a precursor molecule that can be imported into isolated intact chloroplasts of spinach (Spinacia oleracea L.) and processed in two or three steps, respectively, to the size of the authentic protein. Furthermore, both gel electrophoresis under non-denaturing conditions and size-exclusion chromatography verified that the three proteins can form distinct heteromeric supramolecular complexes of approximately 860, 1380 and 1700 kDa (probably also of 600 kDa) molecular mass. The size ranges of the former two are reminiscent of those of Clp complexes described from Escherichia coli. In addition, various complexes between 160 and 560 kDa are detectable with the individual components. Both the processing “intermediates” and the mature nClpP are found in assembled form. Received: 11 March 1998 / Accepted: 8 July 1998     

Structure-based phenotyping predicts HIV-1 protease inhibitor resistance

Mutations in HIV-1 drug targets lead to resistance and consequent therapeutic failure of antiretroviral drugs. Phenotypic resistance assays are time-consuming and costly, and genotypic rules-based interpretations may fail to predict the effects of multiple mutations. We have developed a computational procedure that rapidly evaluates changes in the binding energy of inhibitors to mutant HIV-1 PR variants. Models of WT complexes were produced from crystal structures. Mutant complexes were built by amino acid substitutions in the WT complexes with subsequent energy minimization of the ligand and PR binding site residues. Accuracy of the models was confirmed by comparison with available crystal structures and by prediction of known resistance-related mutations. PR variants from clinical isolates were modeled in complex with six FDA-approved PIs, and changes in the binding energy (DeltaE(bind)) of mutant versus WT complexes were correlated with the ratios of phenotypic 50% inhibitory concentration (IC(50)) values. The calculated DeltaE(bind) of five PIs showed significant correlations (R(2) = 0.7-0.8) with IC(50) ratios from the Virco Antivirogram assay, and the DeltaE(bind) of six PIs showed good correlation (R(2) = 0.76-0.85) with IC(50) ratios from the Virologic PhenoSense assay. DeltaE(bind) cutoffs corresponding to a four-fold increase in IC(50) were used to define the structure-based phenotype as susceptible, resistant, or equivocal. Blind predictions for 78 PR variants gave overall agreement of 92% (kappa = 0.756) and 86% (kappa = 0.666) with PhenoSense and Antivirogram phenotypes, respectively. The structural phenotyping predicted drug resistance of clinical HIV-1 PR variants with an accuracy approaching that of frequently used cell-based phenotypic assays.     

Modeling conformational changes in cyclosporin A.

NMR and X-ray structures for the immunosuppressant cyclosporin A (CsA) reveal a remarkable difference between the unbound (free) conformation in organic solvents and the conformation bound to cyclophilin. We have performed computer simulations of the molecular dynamics of CsA under a variety of conditions and confirmed the stability of these two conformations at room temperature in water and in vacuum. However, when the free conformation was modeled in vacuum at 600 K, a transition pathway leading to the bound conformation was observed. This involved a change in the cis MeLeu-9 peptide bond to a trans conformation and the movement of the side chains forming the dominant hydrophobic cluster (residues MeBmt-1, MeLeu-4, MeLeu-6, and MeLeu-10) to the opposite side of the plane formed by the backbone atoms in the molecular ring. The final conformation had a backbone RMS deviation from the bound conformation of 0.53 A and was as stable in dynamics simulations as the bound conformation. Our calculations allowed us to make a detailed analysis of a transition pathway between the free and the bound conformations of CsA and to identify two distinct regions of coordinated movement in CsA, both of which underwent transitions independently.     

Modeling Tetrahymena thermophila growth and protease production

Oxygen concentrations stimulated growth (maximum number of cells) and protease secretion by Tetrahymena thermophila. Agitation and aeration conditions for growth and protease secretion were optimised by a central composite design. The best optimised combination was a stirrer speed of 338 rpm and an aeration of 1 vvm. Journal of Industrial Microbiology & Biotechnology (2000) 25, 58–61. Received 24 September 1999/ Accepted in revised form 06 March 2000