Prediction of protein--ligand interactions: the complex of porcine pancreatic elastase with a valine-derived benzoxazinone

Prior to availability of the crystal structure of the complex, we evaluated models of the complex between porcine pancreatic elastase and a t-Boc–Val-derived benzoxazinone inhibitor. Models of the noncovalent and covalent complex were generated using computer graphics and each model was subjected to energy minimization using molecular mechanics. After the crystal structure became available, we found that the model with the lowest energy was in good agreement with the crystal structure, except for the position of the His57 side chain. Permissible conformations of the inhibitor were based on information from x-ray crystal structures and an earlier conformational energy investigation of t-Boc–amino acids. We did not, however, limit ourselves to these conformations. The conformation of the inhibitor in the lowest energy model and crystal structure, was not similar to any of the minimum-energy conformations of t-Boc–amino acids. This suggests that limiting proposed binding modes only to the lowest energy conformations of a ligand (prior to binding) may sometimes unfairly bias the procedure.     

Creation of a polyenzyme composition of proteases from Streptomyces griseus and Acremonium chrysogenum

Proteolytic enzymes of microorganisms have been studied for the possibility to create their polyenzymic composition in order to rise a degree of protein hydrolysis and to lower the process duration. Optimal action conditions are selected and a hydrolysis of a number of globular and fibrillar proteins is conducted by a polyenzymic system of Streptomyces griseus and Acremonium chrysogenum proteases.     

A hierarchical method for generating low-energy conformers of a protein-ligand complex

A novel dynamical protocol for finding the low-energy conformations of a protein-ligand complex is described. The energy functions examined consist of an empirical force field with four different dielectric screening models; the generalized Born/surface area model also is examined. Application of the method to three complexes of known crystal structure provides insights into the energy functions used for selecting low-energy docked conformations and into the structure of the binding-energy surface. Evidence is presented that the local energy minima of a ligand in a binding site are arranged in a hierarchical fashion. This observation motivates the construction of a hierarchical docking algorithm that substantially enriches the population of ligand conformations close to the crystal conformation. The algorithm is also adapted to permit docking into a flexible binding site and preliminary tests of this method are presented. Proteins 33:475–495, 1998. © 1998 Wiley-Liss, Inc.     

Crystal structure of casein kinase-1, a phosphate-directed protein kinase.

The structure of a truncated variant of casein kinase-1 from Schizosaccharomyces pombe, has been determined in complex with MgATP at 2.0 A resolution. The model resembles the 'closed', ATP-bound conformations of the cyclin-dependent kinase 2 and the cAMP-dependent protein kinase, with clear differences in the structure of surface loops that impart unique features to casein kinase-1. The structure is of unphosphorylated, active conformation of casein kinase-1 and the peptide-binding site is fully accessible to substrate.     

In Silico Vaccine Design Based on Molecular Simulations of Rhinovirus Chimeras Presenting HIV-1 gp41 Epitopes

A cluster of promising epitopes for the development of human immunodeficiency virus (HIV) vaccines is located in the membrane-proximal external region (MPER) of the gp41 subunit of the HIV envelope spike structure. The crystal structure of the peptide corresponding to the so-called ELDKWA epitope (HIV-1 HxB2 gp41 residues 662-668), in complex with the corresponding broadly neutralizing human monoclonal antibody 2F5, provides a target for structure-based vaccine design strategies aimed at finding macromolecular carriers that are able to present this MPER-derived epitope with optimal antigenic activity. To this end, a series of replica exchange molecular dynamics computer simulations was conducted to characterize the distributions of conformations of ELDKWA-based epitopes inserted into a rhinovirus carrier and to identify those with the highest fraction of conformations that are able to bind 2F5. The length, hydrophobic character, and precise site of insertion were found to be critical for achieving structural similarity to the target crystal structure. A construct with a high degree of complementarity to the corresponding determinant region of 2F5 was obtained. This construct was employed to build a high-resolution structural model of the complex between the 2F5 antibody and the chimeric human rhinovirus type 14:HIV-1 ELDKWA virus particle. Additional simulations, which were conducted to study the conformational propensities of the ELDKWA region in solution, confirm the hypothesis that the ELDKWA region of gp41 is highly flexible and capable of assuming helical conformations (as in the postfusion helical bundle structure) and β-turn conformations (as in the complex with the 2F5 antibody). These results also suggest that the ELDKWA epitope can be involved in intramolecular—and likely intermolecular—hydrophobic interactions. This tendency offers an explanation for the observation that mutations decreasing the hydrophobic character of the MPER in many cases result in conformational changes that increase the affinity of this region for the 2F5 antibody.     

Protein interactions and fluctuations in a proteomic network using an elastic network model

A set of protein conformations are analyzed by normal mode analysis. An elastic network model is used to obtain fluctuation and cooperativity of residues with low amplitude fluctuations across different species. Slow modes that are associated with the function of proteins have common features among different protein structures. We show that the degree of flexibility of the protein is important for proteins to interact with other proteins and as the species gets more complex its proteins become more flexible. In the complex organism, higher cooperativity arises due to protein structure and connectivity.     

Computational and experimental studies of the catalytic mechanism of Thermobifida fusca cellulase Cel6A (E2)

Mutagenesis experiments suggest that Asp79 in cellulase Cel6A (E2) from Thermobifida fusca has a catalytic role, in spite of the fact that this residue is more than 13 A from the scissile bond in models of the enzyme-substrate complex built upon the crystal structure of the protein. This suggests that there is a substantial conformational shift in the protein upon substrate binding. Molecular mechanics simulations were used to investigate possible alternate conformations of the protein bound to a tetrasaccharide substrate, primarily involving shifts of the loop containing Asp79, and to model the role of water in the active site complex for both the native conformation and alternative low-energy conformations. Several alternative conformations of reasonable energy have been identified, including one in which the overall energy of the enzyme-substrate complex in solution is lower than that of the conformation in the crystal structure. This conformation was found to be stable in molecular dynamics simulations with a cellotetraose substrate and water. In simulations of the substrate complexed with the native protein conformation, the sugar ring in the -1 binding site was observed to make a spontaneous transition from the (4)C(1) conformation to a twist-boat conformer, consistent with generally accepted glycosidase mechanisms. Also, from these simulations Tyr73 and Arg78 were found to have important roles in the active site. Based on the results of these various MD simulations, a new catalytic mechanism is proposed. Using this mechanism, predictions about the effects of changes in Arg78 were made which were confirmed by site-directed mutagenesis.     

Spatial structure of the rp142 peptide containing the immunodominant epitope of the HIV-1 gp120 protein. A theoretical study]

A model of spatial structure of the synthetic peptide rp142 (24 amino acid residues) containing the immunodominant epitope of the HIV-1 protein gp120 in the region Gly-10-Phe-15 was constructed by the method of "constrained" molecular mechanics, which uses the algorithms of theoretical conformational analysis, based on NMR spectroscopy data. A comparative analysis of calculated conformations revealed that the spatial structure of rp142 in solution can be described by a family of conformations to which nine different structural clusters involving the sets of topologically close conformers correspond. It is shown that the main chain of the peptide forms irregular but "structured" conformations in which the main portion of amino acid residues is incorporated into beta-turns and helix-like fragments, while Pro-11 and Gly-12 form in some structures inverse gamma-turns, which rarely occur in protein-peptide molecules. It was found that the spatial packing of the Gly-10-Phe-15 hexapeptide in different clusters is realized at different internal rotation angles, to which topologically close structures correspond. It is assumed that this invariant structural element describes the "conformation of complex formation" that is complementary to the antigen-binding center of antibodies and is responsible for their binding to the peptide.     

Mechanism of action of aspartyl proteinases. VI. Nonvalent enzyme-inhibitory and enzyme-substrate complexes of the aspartyl proteinase rhizopus pepsin

The structure of a complex of rhizopuspepsin, a fungal aspartyl protease, with Pro1-Phe2-His3-Phe4-psi[CH2-NH]-Phe5-Val6, its substrate-like inhibitor, was calculated by theoretical conformational analysis. The search for energetically favorable conformational variants of the ligand structure was based on the fragmental approach using the dynamic library of peptide fragments, which were successively extended in the potential field of the protein. The root-mean-square deviation of atom positions in the calculated and experimental inhibitor conformations was 0.56 A. A similar approach was used to model a noncovalent complex of rhizopuspepsin with Pro1-Phe2-His3-Lys4-Phe5-Val6, its specific substrate. As a result, two isoenergetic structures of the complex with different arrangements of the cleavable peptide group and a nucleophilic water molecule were calculated. The possibility of the achieving each of these conformations during the catalytic act is considered. It is shown that there are no structural prerequisites for the distortion of the cleavable bond in the active site of the enzyme. On the basis of the resulting structural data, the assumption was made that Asp35 may be protonated at a late stage of formation of the tetrahedral intermediate rather than at the basic state of the complex.     

SnugDock: Paratope Structural Optimization during Antibody-Antigen Docking Compensates for Errors in Antibody Homology Models

High resolution structures of antibody-antigen complexes are useful for analyzing the binding interface and to make rational choices for antibody engineering. When a crystallographic structure of a complex is unavailable, the structure must be predicted using computational tools. In this work, we illustrate a novel approach, named SnugDock, to predict high-resolution antibody-antigen complex structures by simultaneously structurally optimizing the antibody-antigen rigid-body positions, the relative orientation of the antibody light and heavy chains, and the conformations of the six complementarity determining region loops. This approach is especially useful when the crystal structure of the antibody is not available, requiring allowances for inaccuracies in an antibody homology model which would otherwise frustrate rigid-backbone docking predictions. Local docking using SnugDock with the lowest-energy RosettaAntibody homology model produced more accurate predictions than standard rigid-body docking. SnugDock can be combined with ensemble docking to mimic conformer selection and induced fit resulting in increased sampling of diverse antibody conformations. The combined algorithm produced four medium (Critical Assessment of PRediction of Interactions-CAPRI rating) and seven acceptable lowest-interface-energy predictions in a test set of fifteen complexes. Structural analysis shows that diverse paratope conformations are sampled, but docked paratope backbones are not necessarily closer to the crystal structure conformations than the starting homology models. The accuracy of SnugDock predictions suggests a new genre of general docking algorithms with flexible binding interfaces targeted towards making homology models useful for further high-resolution predictions.     

Protein Interactions and Fluctuations in a Proteomic Network using an Elastic Network Model

Abstract

A set of protein conformations are analyzed by normal mode analysis. An elastic network model is used to obtain fluctuation and cooperativity of residues with low amplitude fluctuations across different species. Slow modes that are associated with the function of proteins have common features among different protein structures. We show that the degree of flexibility of the protein is important for proteins to interact with other proteins and as the species gets more complex its proteins become more flexible. In the complex organism, higher cooperativity arises due to protein structure and connectivity.     

Crystal structure of the biologically active form of class Ib ribonucleotide reductase small subunit from Mycobacterium tuberculosis

Two nrdF genes of Mycobacterium tuberculosis code for different R2 subunits of the class Ib ribonucleotide reductase (RNR). The proteins are denoted R2F-1 and R2F-2 having 71% sequence identity. The R2F-2 subunit forms the biologically active RNR complex with the catalytic R1E-subunit. We present the structure of the reduced R2F-2 subunit to 2.2 A resolution. Comparison of the R2F-2 structure with a model of R2F-1 suggests that the important differences are located at the C-terminus. We found that within class Ib, the E-helix close to the iron diiron centre has two preferred conformations, which cannot be explained by the redox-state of the diiron centre. In the R2F-2 structure, we also could see a mobility of alphaE in between the two conformations.     

Apolipoprotein structure and dynamics

PURPOSE OF REVIEW: This review highlights recent advances in structural studies of exchangeable human apolipoproteins and the insights they provide into lipoprotein action in cardiovascular and amyloid diseases. RECENT FINDINGS: The high-resolution X-ray crystal structure of free apoA-II reveals a parallel helical array that may represent other lipid-poor apolipoproteins, and the structure in complex with detergent substantiates the belt model for the protein arrangement on lipoproteins. Nuclear magnetic resonance structures of apolipoprotein-detergent complexes show a repertoire of curved helical conformations, suggesting multiple helical arrangements on the lipid. Low-resolution spectroscopic analyses, interface studies and molecular modeling provide new insights into the 'hinge-domain' mechanism of apolipoprotein adaptation at variable lipoprotein surfaces. A kinetic mechanism for lipoprotein stabilization is proposed. SUMMARY: Cumulative evidence supports the belt model that provides a general structural basis for understanding the molecular mechanisms of functional apolipoprotein reactions, such as binding to lipoprotein receptors, lipid transporters, and the activation of lipophilic enzymes. However, the detailed protein and lipid conformations on lipoproteins and the underlying molecular interactions are unclear. New insights will hopefully emerge once the first detailed lipoprotein structure is solved.     

Molecular modeling of p-chlorophenoxyacetic acid binding to the CLC-0 channel

Molecular simulation techniques were applied to predict the interaction of the CLC-0 Cl(-) channel and the channel-blocking molecule p-chlorophenoxyacetic acid (CPA). A three-dimensional model of the CLC-0 channel was constructed on the basis of the homology with the bacterial Cl(-) channel StCLC, the structure of which has been solved by X-ray crystallography. Docking of the CPA molecule was obtained by using a geometric recognition algorithm, yielding 5000 possible conformations. By restraining the simulation to those conformations in which CPA is near the intracellular mouth of the channel, the CPA-protein complex models were reduced to three sets of conformations, which are interconvertible within 2 ns when molecular dynamics is applied to the system. Point mutations of CLC-0 at three different positions predicted to interact with CPA in these configurations did, however, not greatly alter CPA inhibition, suggesting a deeper final binding location. In the model, binding of CPA to a more internal position in the ionic pathway was obtained by applying a constant force vector to CPA, pushing it toward the center of the channel. This technique allowed us to outline the possible intrachannel pathway of CPA and to describe qualitatively the binding sites and energy barriers of this pathway. The consistency of the obtained models and the experimental data indicates that the CLC-0-CPA complex model is reasonable and can be used in further biological studies, such as rational design of blocking agents of and mutagenesis of CLC Cl(-) channels.     

Crystal structure of human interleukin-10 at 1.6 A resolution and a model of a complex with its soluble receptor.

The crystal structure of human interleukin-10 (IL-10) was refined at 1.6 A resolution against X-ray diffraction data collected at 100 K with the use of synchrotron radiation. Although similar to the IL-10 structure determined previously at room temperature, this low-temperature IL-10 structure contains, in addition, four N-terminal residues, three sulfate anions, and 175 extra water molecules. Whereas the main-chain conformation is preserved, about 30% of the side chains, most of them on the protein surface, assume different conformations. A computer model of a complex of IL-10 with its two soluble receptors was generated based on the topological similarity of IL-10 to interferon-gamma. The contact region between the cytokine and each receptor shows excellent complementarity of polar and hydrophobic interactions, suggesting that the model is generally correct and should be useful in guiding mutagenesis experiments.     

Application of a graphic method for analyzing the kinetics of multienzyme reactions

A new graphic method is proposed to solve kinetic equations for polyenzymic reactions. Each graph apex is corresponded by the transmitting function deduced from kinetic equations by means of Laplas transformation. Application of this procedure allows to simplify the solution of kinetic equations and its analysis. The procedure suggested makes it possible to use the methods of automatic control when solving theoretical problems of enzymology.     

Mechanism of action of aspartic proteases. III. Conformational characteristics of HIV-1 protease inhibitor JG-365]

A set of conformations was shown to be characteristic of the free-state spatial structure of substrate-like inhibitor JG-365 for aspartic protease from HIV-1. Among them, the lowest-energy conformations have a folded form of the peptide backbone. The inhibitor has a noncleavable hydroxyethylamine group with an additional chiral center in its structure. Our calculations showed that only the S-isomer of the inhibitor displays conformational characteristics that practically coincide with those of the native substrate for HIV-1 protease. One of the calculated conformations with a completely extended main chain and a relative energy of 9.5 kcal/mol very closely mimics the experimentally observed structure of the inhibitor in the enzyme-inhibitor complex. The realization of this structure is unlikely for a free inhibitor, because it has only a small number of interresidual noncovalent interactions in the extended conformation; these are presumably compensated for by intermolecular interactions at the active site of the enzyme.     

X-ray structure of HIV-1 protease tethered dimer complexed to ritonavir

We report here the 2.5A structure of HIV-1 protease tethered-dimer ritonavir complex. The inhibitor bound in the active site has different conformations in the two orientations. There is only one hydrogen bond between the inhibitor and the enzyme. The conserved flap-water is not found in the present complex.