Predicting structural effects in HIV-1 protease mutant complexes with flexible ligand docking and protein side-chain optimization

We present a computational approach for predicting structures of ligand-protein complexes and analyzing binding energy landscapes that combines Monte Carlo simulated annealing technique to determine the ligand bound conformation with the dead-end elimination algorithm for side-chain optimization of the protein active site residues. Flexible ligand docking and optimization of mobile protein side-chains have been performed to predict structural effects in the V32I/I47V/V82I HIV-1 protease mutant bound with the SB203386 ligand and in the V82A HIV-1 protease mutant bound with the A77003 ligand. The computational structure predictions are consistent with the crystal structures of these ligand-protein complexes. The emerging relationships between ligand docking and side-chain optimization of the active site residues are rationalized based on the analysis of the ligand-protein binding energy landscape. Proteins 33:295–310, 1998. © 1998 Wiley-Liss, Inc.     

Stabilization of an alpha-helical conformation in an isolated hexapeptide inhibitor of calmodulin.

The conformational properties of two hexapeptides, Ac-LWRILW-NH(2) and its D-amino acid counterpart Ac-lwrilw-NH(2), identified as calmodulin inhibitors using mixture-based synthetic combinatorial library approaches, have been characterised by NMR and CD spectroscopy. The peptides fold into an alpha-helical conformation in aqueous solution. The observed short- and medium-range nuclear Overhauser effects were consistent with the formation of an alpha-helical structure and a reasonably well-defined set of structures was obtained by using restraints from the NMR data in simulated annealing calculations. Analysis of glycine-substitution analogues demonstrated that all the amino acids that make up the peptide sequence are important for the stabilization of the alpha-helical conformation. The results suggest that a well-defined set of interactions is indispensable to allow alpha-helix formation in this short hexapeptide.     

The Bowman-Birk inhibitor reactive site loop sequence represents an independent structural beta-hairpin motif

We have determined the NMR structure in aqueous solution of a disulphide-cyclised 11-residue peptide that forms a stable beta-hairpin, incorporating a type VIb beta-turn. The structure is found to be extremely well ordered for a short peptide, with the 30 lowest energy simulated annealing structures having an average pairwise r.m.s. deviation of only 0.36 A over the backbone. All but three side-chains adopt distinct conformations, allowing a detailed analysis of their involvement in cross-strand interactions. The peptide sequence analysed originates from a previously reported study, which identified potent inhibitors of human leukocyte elastase from screening a combinatorial peptide library based on the short protein beta-sheet segment that forms the reactive site loop of Bowman-Birk inhibitors. A detailed comparison of the peptide's solution structure with the corresponding region in the whole protein structure reveals a very good correspondence not only for the backbone (r.m.s. deviation approximately 0.7 A) but also for the side-chains. This isolated beta-hairpin retains the biologically active "canonical conformation" typical of small serine proteinase inhibitor proteins, which explains why it retains inhibitory activity. Since the structural integrity is sequence-inherent and does not depend upon the presence of the remaining protein, this beta-hairpin represents an independent structural motif and so provides a useful model of this type of protein architecture and its relation to biological function. The relationship between the conformation of this beta-hairpin and its biological activity is discussed.     

High-resolution three-dimensional structure of reduced recombinant human thioredoxin in solution

The solution structure of recombinant human thioredoxin (105 residues) has been determined by nuclear magnetic resonance (NMR) spectroscopy combined with hybrid distance geometry-dynamical simulated annealing calculations. Approximate interproton distance restraints were derived from nuclear Overhauser effect (NOE) measurements. In addition, a large number of stereospecific assignments for beta-methylene protons and torsion angle restraints for phi, psi, and chi 1 were obtained by using a conformational grid search on the basis of the intraresidue and sequential NOE data in conjunction with 3JHN alpha and 3J alpha beta coupling constants. The structure calculations were based on 1983 approximate interproton distance restraints, 52 hydrogen-bonding restraints for 26 hydrogen bonds, and 98 phi, 71 psi, and 72 chi 1 torsion angle restraints. The 33 final simulated annealing structures obtained had an average atomic rms distribution of the individual structures about the mean coordinate positions of 0.40 +/- 0.06 A for the backbone atoms and 0.78 +/- 0.05 A for all atoms. The solution structure of human thioredoxin consists of a five-stranded beta-sheet surrounded by four alpha-helices, with an active site protrusion containing the two redox-active cysteines. The overall structure is similar to the crystal and NMR structures of oxidized [Katti, S. K., LeMaster, D. M., & Eklund, H. (1990) J. Mol. Biol. 212, 167-184] and reduced [Dyson, J. H., Gippert, G. P., Case, D. A., Holmgren, A., & Wright, P. (1990) Biochemistry 29, 4129-4136] Escherichia coli thioredoxin, respectively, despite the moderate 25% amino acid sequence homology. Several differences, however, can be noted. The human alpha 1 helix is a full turn longer than the corresponding helix in E. coli thioredoxin and is characterized by a more regular helical geometry. The helix labeled alpha 3 in human thioredoxin has its counterpart in the 3(10) helix of the E. coli protein and is also longer in the human protein. In contrast to these structural differences, the conformation of the active site loop in both proteins is very similar, reflecting the perfect sequence identity for a stretch of eight amino acid residues around the redox-active cysteines.     

Substrate recognition by the Lyn protein-tyrosine kinase. NMR structure of the immunoreceptor tyrosine-based activation motif signaling region of the B cell antigen receptor

The immunoreceptor tyrosine-based activation motif (ITAM) plays a central role in transmembrane signal transduction in hematopoietic cells by mediating responses leading to proliferation and differentiation. An initial signaling event following activation of the B cell antigen receptor is phosphorylation of the CD79a (Ig-alpha) ITAM by Lyn, a Src family protein-tyrosine kinase. To elucidate the structural basis for recognition between the ITAM substrate and activated Lyn kinase, the structure of an ITAM-derived peptide bound to Lyn was determined using exchange-transferred nuclear Overhauser NMR spectroscopy. The bound substrate structure has an irregular helix-like character. Docking based on the NMR data into the active site of the closely related Lck kinase strongly favors ITAM binding in an orientation similar to binding of cyclic AMP-dependent protein kinase rather than that of insulin receptor tyrosine kinase. The model of the complex provides a rationale for conserved ITAM residues, substrate specificity, and suggests that substrate binds only the active conformation of the Src family tyrosine kinase, unlike the ATP cofactor, which can bind the inactive form.     

Modeling of the three-dimensional structure of luffin-α and its simulated reaction with the substrate oligoribonucleotide GAGA

A fundamental problem in biochemistry and molecular biology is understanding the spatial structure of macromolecules and then analyzing their functions. In this study, the three-dimensional structure of a ribosome-inactivating protein luffin-α was predicted using a neural network method and molecular dynamics simulation. A feedforward neural network with the backpropagation learning algorithm were trained on model class of homologous proteins including trichosanthin andα-momorcharin. The distance constraints for the Cα atoms in the protein backbone were utilized to generate a folded crude conformation of luffin-α by model building and the steepest descent minimization approach. The crude conformation was refined by molecular dynamics techniques and a simulated annealing procedure. The interaction between luffin-α and its analogous substrate GAGA was also simulated to understand its action mechanism.     

Biologically active conformations of thymopentin. Studies with conformationally restricted analogs

Four cyclic analogs of thymopentin were synthesized and evaluated for biological activity on the human T cell line CEM. Three of these conformationally restricted analogs were biologically active. The one analog which most closely mimicked the conformation predicted from NMR and theoretical energy minimization calculations proved to be inactive. These studies establish that the biologically active conformations of thymopentin differ from the most probable conformation predicted from solution NMR and theoretical energy minimization studies.     

Protein Folding: A Few Random Thoughts

Abstract

Both the aqueous and lipid-induced structure of Kassinin, a dodecapeptide of amphibian origin, has been studied by two-dimensional proton nuclear magnetic resonance (2D ¹H-NMR) spectroscopy and distance geometry calculations. Unambiguous NMR assignments of protons have been made with the aid of correlation spectroscopy (DQF-COSY and TOCSY) experiments and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The distance constraints obtained from the NMR data have been utilized in a distance geometry algorithm to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that, while in water Kassinin prefers to be in an extended chain conformation, in the presence of perdeuterated dodecylphosphocholine (DPC) micelles, a membrane model system, helical conformation is induced in the central core and C-terminal region (K4-M12) of the peptide. N-terminus though less defined also displays some degree of order and a possible turn structure. The conformation adopted by Kassinin in the presence of DPC micelles is consistent with the structural motif typical of neurokinin-1 selective agonists and with that reported for Eledoisin in hydrophobic environment.     

Relaxation matrix refinement of the solution structure of squash trypsin inhibitor

The structure of the small squash trypsin inhibitor CMTI-I is refined by directly minimizing the difference between the observed two-dimensional nuclear Overhauser enhancement (NOE) intensities and those calculated by the full relaxation matrix approach. To achieve this, a term proportional to this difference was added to the potential energy function of the molecular dynamics program X-PLOR. Derivatives with respect to atomic co-ordinates are calculated analytically. Spin diffusion effects are thus accounted for fully during the refinement. Initial structures for the refinement were those determined recently by solution nuclear magnetic resonance using the isolated two-spin approximation to derive distance range estimates. The fits to the nuclear magnetic resonance data improve significantly with only small shifts in the refined structures during a few cycles of conjugate gradient minimization. However, larger changes (approximately 1 A) in the conformation occur during simulated annealing, which is accompanied by a further reduction of the difference between experimental and calculated two-dimensional NOE intensities. The refined structures are closer to the X-ray structure of the inhibitor complexed with trypsin than the initial structures. The root-mean-square difference for backbone atoms between the initial structures and the X-ray structure is 0.96 A, and that between the refined structures and the X-ray structure 0.61 A.     

NMR three-dimensional solution structure of the serine protease inhibitor cyclotheonamide A

The nmr solution conformation of cyclotheonamide A (CtA) was determined in aqueous media. The data produced 15 distance and 10 torsional constraints which were used to generate conformations using restrained simulated annealing (SA) and distance geometry/simulated annealing (DG/SA) calculations. Two different calculation protocols were performed to ensure proper sampling of conformational space and even though the torsional restraints were input differently, both calculation methods yielded the same conformation of CtA. In the structure calculations, all solutions of the Karplus equation were sampled simultaneously using the restrained SA protocol and large ranges were used for the dihedral restraints in the DG/SA protocol because all solutions to the Karplus equation could not be sampled simultaneously. The solution conformation was also compared to the solid state x-ray conformations of CtA bound to thrombin and trypsin. The conformation of the residues important for active site binding (d -Phe, h-Arg, and Pro) are nearly identical in aqueous solution and solid state with largest differences at the a-Ala and v-Tyr residues. CtA appears to be preordered in structure and does not undergo a significant conformational change upon binding to the enzyme active site. © 1997 John Wiley & Sons, Inc.     

Conformation of a novel tetrapeptide inhibitor NH2-D-Trp-D-Met-Phe(pCl)-Gla-NH2 bound to farnesyl-protein transferase.

Farnesyl-protein transferase (FPTase) catalyzes the posttranslational farnesylation of the cysteine residue located in the C-terminal tetrapeptide of the Ras oncoprotein. Prenylation of this residue is essential for membrane association and cell-transforming activities of ras. Inhibitors of FPTase have been demonstrated to display antitumor activity in both tissue culture and animal models, and thus represent a potential therapeutic strategy for the treatment of human cancers. A synthetic tetrapeptide library, which included an expanded set of 68 L-, D- and noncoded amino acids, has been screened for inhibitors of FPTase activity. The tetrapeptide, NH2-D-Trp-D-Met-L-Phe(pCl)-L-Gla-NH2 was shown to be competitive with the isoprenyl cosubstrate, farnesyl diphosphate (FPP) but not with the peptide substrate, the C-terminal tetrapeptide of the Ras protein. The FPTase-bound conformation of the inhibitor, NH2-D-Trp-D-Met-L-Phe(pCl)-L-Gla-NH2 was determined by NMR spectroscopy. Distance constraints were derived from two-dimensional transferred nuclear Overhauser effect (TRNOE) experiments. Ligand competition experiments identified the NOEs that originate from the active-site conformation of the inhibitor. Structures were calculated using a combination of distance geometry and restrained energy minimization. The peptide backbone is shown to adopt a reverse-turn conformation most closely approximating a type II' beta-turn. The resolved conformation of the inhibitor represents a distinctly different structural motif from that determined for Ras-competitive inhibitors. Knowledge of the bound conformation of this novel inhibitor provides a template and future direction for the design of new classes of FPTase antagonists.     

Three-dimensional structure of the mammalian tachykinin peptide neurokinin A bound to lipid micelles

The solution structure of NKA, a decapeptide of mammalian origin, has been characterized by CD spectropolarimetry and 2D proton nuclear magnetic resonance (2D 1H-NMR) spectroscopy in both aqueous and membrane mimetic solvents. Unambiguous NMR assignments of protons have been made with the aid of correlation spectroscopy (DQF-COSY and TOCSY) experiments and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The distance constraints obtained from the NMR data have been utilized to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that in water NKA prefers to be in an extended chain conformation whereas a helical conformation is induced in the central core and the C-terminal region (D4-M10) of the peptide in the presence of perdeuterated dodecylphosphocholine (DPC) micelles, a membrane model system. Though less defined the N-terminus also displays some degree of order and a possible turn structure. The conformation adopted by NKA in the presence of DPC micelles represents a structural motif typical of neurokinin-2 selective agonists and is similar to that reported for eledoisin in hydrophobic environment.     

Substrate binding and catalytic mechanism in phospholipase C from Bacillus Cereus: A molecular mechanics and molecular dynamics study

For the first time a consistent catalytic mechanism of phospholipase C from Bacillus cereus is reported based on molecular mechanics calculations. We have identified the position of the nucleophilic water molecule, which is directly involved in the hydrolysis of the natural substrate, phosphatidylcholine, in phospholipase C. This catalytically essential water molecule, after being activated by an acidic residue (Asp55), performs the nucleophilic attack on the phosphorus atom in the substrate, leading to a trigonal bipyramidal pentacoordinated intermediate (and structurally similar transition state). The subsequent collapse of the intermediate, regeneration of the enzyme, and release of the products has to involve a not yet identified second water molecule. The catalytic mechanism reported here is based on a series of molecular mechanics calculations. First, the x-ray structure of phospholipase C from B. cereus including a docked substrate molecule was subjected to a stepwise molecular mechanics energy minimization. Second, the location of the nucleophilic water molecule in the active site of the fully relaxed enzyme–substrate complex was determined by evaluation of nonbonded interaction energies between the complex and a water molecule. The nucleophilic water molecule is positioned at a distance (3.8 Å) from the phosphorus atom in the substrate, which is in good agreement with experimentally observed distances. Finally, the stability of the complex between phospholipase C, the substrate, and the nucleophilic water molecule was verified during a 100 ps molecular dynamics simulation. During the simulation the substrate undergoes a conformational change, but retains its localization in the active site. The contacts between the enzyme, the substrate, and the nucleophilic water molecule display some fluctuations, but remain within reasonable limits, thereby confirming the stability of the enzyme–substrate–water complex. The protocol developed for energy minimization of phospholipase C containing three zinc ions located closely together at the bottom of the active site cleft is reported in detail. In order to handle the strong electrostatic interactions in the active site realistically during energy minimization, delocalization of the charges from the three zinc ions was considered. Therefore, quantum mechanics calculations on the zinc ions and the zinc-coordinating residues were carried out prior to the molecular mechanics calculations, and two different sets of partial atomic charges (MNDO-Mulliken and AM1-ESP) were applied. After careful assignment of partial atomic charges, a complete energy minimization of the protein was carried out by a stepwise procedure without explicit solvent molecules. Energy minimization with either set of charges yielded structures, which were very similar both to the x-ray structure and to each other, although using AM1-ESP partial atomic charges and a dielectric constant of 4, yielded the best protein structure. © 1997 John Wiley and Sons, Inc. Biopoly 42: 319–336, 1997     

Searching for the Lowest Energy Conformation of Substrates in the Carboxypeptidase A Active Site Using Monte Carlo/Energy Minimization Techniques

More than fifteen substrates of carboxypeptidase A (CPA) have been “docked” to the active site of the enzyme and searched for the lowest energy conformation of the substrates bound to the active site. The method employed combines Monte Carlo procedures with energy minimization (MC/EM procedures). The distances of P_1′–P₃to S_1′–S₃were measured. The computational results are consistent with the proposed binding interactions of CPA and its substrates and also with the promoted-water pathway of CPA catalysis. This study has demonstrated that MC/EM procedures are very useful for searching the conformations of substrates available in the active site of the enzyme.     

The active site and substrates binding mode of malonyl-CoA synthetase determined by transferred nuclear Overhauser effect spectroscopy, site-directed mutagenesis, and comparative modeling studies

The active sites and substrate bindings of Rhizobium trifolii molonyl-CoA synthetase (MCS) catalyzing the malonyl-CoA formation from malonate and CoA have been determined based on NMR spectroscopy, site-directed mutagenesis, and comparative modeling methods. The MCS-bound conformation of malonyl-CoA was determined from two-dimensional-transferred nuclear Overhauser effect spectroscopy data. MCS protein folds into two structural domains and consists of 16 alpha-helices, 24 beta-strands, and several long loops. The core active site was determined as a wide cleft close to the end of the small C-terminal domain. The catalytic substrate malonate is placed between ATP and His206 in the MCS enzyme, supporting His206 in its catalytic role as it generates reaction intermediate, malonyl-AMP. These findings are strongly supported by previous biochemical data, as well as by the site-directed mutagenesis data reported here. This structure reveals the biochemical role as well as the substrate specificity that conservative residues of adenylate-forming enzymes have.     

Studies of the mechanism of action and regulation of cAMP-dependent protein kinase

Magnetic resonance and kinetic studies of the catalytic subunit of a Type II cAMP-dependent protein kinase from bovine heart have established the active complex to be an enzyme-ATP-metal bridge. The metal ion is β,γ coordinated with Δ chirality at the β-phosphorous atom. The binding of a second metal ion at the active site which bridges the enzyme to the three phosphoryl groups of ATP, partially inhibits the reaction. Binding of the metal-ATP substrate to the enzyme occurs in a diffusion-controlled reaction followed by a 40 ° change in the glycosidic torsional angle. This conformational change results from strong interaction of the nucleotide base with the enzyme. NMR studies of four ATP-utilizing enzymes show a correlation between such conformational changes and high nucleotide base specificity. Heptapeptide substrates and substrate analogs bind to the active site of the catalytic subunit at a rate significantly lower than collision frequency indicating conformational selection by the enzyme or a subsequent slow conformational change. NMR studies of the conformation of the enzyme-bound peptide substrates have ruled out α-helical and β-pleated sheet structures. The results of kinetic studies of peptide substrates in which the amino acid sequence was systematically varied were used to rule out the obligatory requirement for all possible β-turn conformations within the heptapeptide although an enzymatic preference for a β_2–5 or β_3–6 turn could not be excluded. Hence if protein kinase has an absolute requirement for a specific secondary structure, then this structure must be a coil. In the enzyme-substrate complex the distance along the reaction coordinate between the γ-P of ATP and the serine oxygen of the peptide substrate (5.3 ± 0.7 Å) allows room for a metaphosphate intermediate. This finding together with kinetic observations as well as the location of the inhibitory metal suggest a dissociative mechanism for protein kinase, although a mechanism with some associative character remains possible. Regulation of protein kinase is accomplished by competition between the regulatory subunit and peptide or protein substrates at the active site of the catalytic subunit. Thus, the regulatory subunit is found by NMR to block the binding of the peptide substrate to the active site of protein kinase but allows the binding of the nucleotide substrate and divalent cations. The dissociation constant of the regulatory subunit from the active site (10^?10m) is increased ~10-fold by phosphorylation and ~10⁴-fold by the binding of cAMP, to a value (10^?5m) which exceeds the intracellular concentration of the R₂C₂ holoenzyme complex (10^?6m). The resulting dissociation of the holoenzyme releases the catalytic subunit, permitting the active site binding of peptide or protein substrates.     

Structural evidence for a functionally relevant second camphor binding site in P450cam: model for substrate entry into a P450 active site

P450cam has long served as a prototype for the cytochrome P450 (CYP) gene family. But, little is known about how substrate enters its active site pocket, and how access is achieved in a way that minimizes exposure of the reactive heme. We hypothesize that P450cam may first bind substrate transiently near the mobile F-G helix that covers the active site pocket. Such a two-step binding process is kinetically required if P450cam rarely populates an open conformation-as suggested by previous literature and the inability to obtain a crystal structure of P450cam in an open conformation. Such a mechanism would minimize exposure of the heme by allowing P450cam to stay in a closed conformation as long as possible, since only brief flexing into an open conformation would be required to allow substrate entry. To test this model, we have attempted to dock a second camphor molecule into the crystal structure of camphor-bound P450cam. The docking identified only one potential entry site pocket, a well-defined cavity on the F-helix side of the F-G flap, 16 A from the heme iron. Location of this entry site pocket is consistent with our NMR T1 relaxation-based measurements of distances for a camphor that binds in fast exchange (active site camphor is known to bind in slow exchange). Presence of a second camphor binding site is also confirmed with [(1)H-(13)C] HSQC titrations of (13)CH3-threonine labeled P450cam. To confirm that camphor can bind outside of the active site pocket, (13)CH3-S-pyridine was bound to the heme iron to physically block the active site, and to serve as an NMR chemical shift probe. Titration of this P450cam-pyridine complex confirms that camphor can bind to a site outside the active site pocket, with an estimated Kd of 43 microM. The two-site binding model that is proposed based on these data is analogous to that recently proposed for CYP3A4, and is consistent with recent crystal structures of P450cam bound to tethered-substrates, which force a partially opened conformation.     

Engineering an Mg2+ site to replace a structurally conserved arginine in the catalytic center of histidyl-tRNA synthetase by computer experiments

Histidyl-tRNA synthetase (HisRS) differs from other class II aminoacyl-tRNA synthetases (aaRS) in that it harbors an arginine at a position where the others bind a catalytic Mg²⁺ ion. In computer experiments, four mutants of HisRS from Escherichia coli were engineered by removing the arginine and introducing a Mg²⁺ ion and residues from seryl-tRNA synthetase (SerRS) that are involved in Mg²⁺ binding. The mutants recreate an active site carboxylate pair conserved in other class II aaRSs, in two possible orders: Glu-Asp or Asp-Glu, replacing Glu-Thr in native HisRS. The mutants were simulated by molecular dynamics in complex with histidyl-adenylate. As controls, the native HisRS was simulated in complexes with histidine, histidyl-adenylate, and histidinol. The native structures sampled were in good agreement with experimental structures and biochemical data. The two mutants with the Glu-Asp sequence showed significant differences in active site structure and Mg²⁺ coordination from SerRS. The others were more similar to SerRS, and one of them was analyzed further through simulations in complex with histidine, and His+ATP. The latter complex sampled two Mg²⁺ positions, depending on the conformation of a loop anchoring the second carboxylate. The lowest energy conformation led to an active site geometry very similar to SerRS, with the principal Mg²⁺ bridging the α- and β-phosphates, the first carboxylate (Asp) coordinating the ion through a water molecule, and the second (Glu) coordinating it directly. This mutant is expected to be catalytically active and suggests a basis for the previously unexplained conservation of the active site Asp-Glu pair in class II aaRSs other than HisRS. Proteins 32:362–380, 1998. © 1998 Wiley-Liss, Inc.     

Analysis of the activation of acetylcholinesterase by carbon nanoparticles using a monolithic immobilized enzyme microreactor: role of the water molecules in the active site gorge

A biochromatographic system was used to study the direct effect of carbon nanoparticles (CNPs) on the acetylcholinesterase (AChE) activity. The AChE enzyme was covalently immobilized on a monolithic CIM-disk via its NH₂ residues. Our results showed an increase in the AChE activity in presence of CNPs. The catalytic constant (k_cat) was increased while the Michaelis constant (K_m) was slightly decreased. This indicated an increase in the enzyme efficiency with increase of the substrate affinity to the active site. The thermodynamic data of the activation mechanism of the enzyme, i.e. ΔH* and ΔS*, showed no change in the substrate interaction mechanism with the anionic binding site. The increase of the enthalpy (ΔH*) and the entropy (ΔS*) with decrease in the free energy of activation (Ea) was related to structural conformation change in the active site gorge. This affected the stability of water molecules in the active site gorge and facilitated water displacement by substrate for entering to the active site of the enzyme.