Denaturation of proteins in methanol/water mixtures

The solvophobic theory developed earlier by Sinanoglu introducing the use of molecular surface areas and microthermodynamic surface and interfacial tensions at molecular dimensions is applied to the interpretation of calorimetric data on denaturation of lysozyme in a wide range of methanol/water mixtures. The experimental values of standard unitary free energies of denaturation correlate well with our predictions. The molecular surface area change of the protein upon denaturation is evaluated using the solvophobic theory. The maximum in the stability of the native form of the protein is predicted to occur at 8% (v/v) methanol. This is found to be in agreement with the experimental results.     

Probing the ultra-high resolution structure of aldose reductase with molecular modelling and noncovalent mass spectrometry

Aldose reductase, the first and rate-limiting enzyme of the polyol pathway, is a target for drug design for the treatment of diabetes complications. The structures of aldose reductase in complex with the cyclic imide inhibitors Fidarestat and Minalrestat were recently determined at ultra-high resolution (Proteins 2004, 55, 805). We have used the detailed structural information revealed at atomic resolution, including the assignment of protonation states for the inhibitors and active site residues, together with molecular modelling and noncovalent mass spectrometry to characterise the type and strength of the interactions between the enzyme and the inhibitors, and to attempt the design of novel potential inhibitors with enhanced binding energies of the complexes. The VC(50) values measured by mass spectrometry (accelerated voltage of ions needed to dissociate 50% of a noncovalent complex in the gas phase) for the aldose reductase inhibitors correlate with the IC(50) values (concentration of inhibitor giving 50% inhibition in solution) and with the electrostatic binding energies calculated between the active site residues Tyr48, His110 and Trp111 and the inhibitors, suggesting that electrostatic interactions play a major role in inhibitor binding. Our molecular modelling and design studies suggest that the replacement of the fluorine atom in Minalrestat's bromo-fluorobenzyl group with nitro, amide and carboxylate functional groups enhanced the predicted net binding energies of the complexes by 16%, 31% and 68%, respectively. When the carbamoyl group of Fidarestat was replaced with a nitro, 4-hydroxyl phenyl and carboxylate functional groups, the predicted net binding energies of the complexes were enhanced by 13%, 34% and 46%, respectively.     

TAP association influences the conformation of nascent MHC class I molecules

The influence of TAP-MHC class I interactions on peptide binding to the class I heavy chain is assessed during TAP-dependent assembly using Kb-specific Abs that recognize conformational changes induced by assembly with beta2-microglobulin (beta2m) and by peptide binding. A significant portion (45%) of Kb molecules in TAP+, RMA-derived microsomes are associated with the TAP complex as measured by coimmunoisolation of Kb using anti-TAP1 Abs, while only 20% of the Kb heavy chain molecules are isolated as Kbbeta2m complexes with the alpha-Kb-specific Abs, Y-3 or K-10-56. The amount of Kb isolated with Y-3 and K-10-56 increases in proportion to transport and binding of peptide to the Kb molecules within the RMA microsomes. In contrast, less than 5% of the Kb within TAP2-RMA-S microsomes associated with the remaining TAP1 subunit. However, greater than 60% of Kb heavy chain is isolated as K-10-56- and Y-3-reactive Kbbeta2m complexes. We propose that a TAP-MHC class I interaction serves to stabilize the MHC class I:beta2m complex in an immature conformation (Y-3 and K-10-56 nonreactive) prior to high affinity peptide binding, preventing the export of class I molecules complexed with low affinity peptide ligands from the ER.     

Interaction energies between beta-lactam antibiotics and E. coli penicillin-binding protein 5 by reversible thermal denaturation

Penicillin-binding proteins (PBPs) catalyze the final stages of bacterial cell wall biosynthesis. PBPs form stable covalent complexes with beta-lactam antibiotics, leading to PBP inactivation and ultimately cell death. To understand more clearly how PBPs recognize beta-lactam antibiotics, it is important to know their energies of interaction. Because beta-lactam antibiotics bind covalently to PBPs, these energies are difficult to measure through binding equilibria. However, the noncovalent interaction energies between beta-lactam antibiotics and a PBP can be determined through reversible denaturation of enzyme-antibiotic complexes. Escherichia coli PBP 5, a D-alanine carboxypeptidase, was reversibly denatured by temperature in an apparently two-state manner with a temperature of melting (T(m)) of 48.5 degrees C and a van't Hoff enthalpy of unfolding (H(VH)) of 193 kcal/mole. The binding of the beta-lactam antibiotics cefoxitin, cloxacillin, moxalactam, and imipenem all stabilized the enzyme significantly, with T(m) values as high as +4.6 degrees C (a noncovalent interaction energy of +2.7 kcal/mole). Interestingly, the noncovalent interaction energies of these ligands did not correlate with their second-order acylation rate constants (k(2)/K'). These rate constants indicate the potency of a covalent inhibitor, but they appear to have little to do with interactions within covalent complexes, which is the state of the enzyme often used for structure-based inhibitor design.     

Magnetic resonance studies on the copper site of dopamine beta-monooxygenase in the presence of cyanide and azide anions

In order to elucidate the coordination state of water molecules in the Cu(II) site of dopamine [( 3,4-dihydroxyphenyl)ethylamine] beta-monooxygenase, measurements of the paramagnetic 1H nuclear magnetic relaxation rate of solvent water in the enzyme solution containing cyanide or azide as an exogenous ligand were carried out to obtain the values of intrinsic paramagnetic relaxation rate decrements Rp1 and Rp2 for the ligand-enzyme 1:1 and 2:1 complexes, respectively. Rp1 (percent) values were 53 (pH 5.5) and 52 (pH 7.0) for cyanide and 38 (pH 5.5) and 32 (pH 7.0) for azide, while Rp2 (percent) values were 98 (pH 5.5) and 96 (pH 7.0) for azide. Although no Rp2 values for cyanide were obtained because of its reducing power at the Cu(II) site, the Rp1 and Rp2 values obtained above prove that the Cu(II) center has two coordinated water molecules that are exchangeable for exogenous ligands at either pH. Supporting evidence was provided by electron paramagnetic resonance (EPR) titration, in which the enzyme solution containing cyanide-enzyme (1:1) complex in an equal proportion to uncomplexed enzyme gave an observed paramagnetic relaxation rate decrement, Rp, of 23%. Another characteristic of the Rp1 and Rp2 values was their invariability with respect to pH, indicating that the three-dimensional structure of the Cu(II) site is pH-invariant within the range examined. Binding constants of ligand to enzyme Kb1 and Kb2 for 1:1 and 2:1 complex formation, respectively, were also determined through an analysis of the Rp values; it was found that Kb1 was larger than Kb2 irrespective of pH. (ABSTRACT TRUNCATED AT 250 WORDS)     

Tapasin-mediated retention and optimization of peptide ligands during the assembly of class I molecules

The murine class I H-2Kb molecule achieves high level surface expression in tapasin-deficient 721.220 human cells. Compared with their behavior in wild-type cells, Kb molecules expressed on 721.220 cells are more receptive to exogenous peptide, undergo more rapid surface decay, and fail to form macromolecular peptide loading complexes. As a result, they are rapidly transported to the cell surface, reflecting a failure of endoplasmic reticulum retention mechanisms in the absence of loading complex formation. Despite the failure of Kb molecules to colocalize to the TAP and their rapid egress to the cell surface, Kb is still capable of presenting TAP-dependent peptides in the absence of tapasin. Furthermore, pool sequencing of peptides eluted from these molecules revealed strict conservation of their canonical H-2Kb-binding motif. There was a reduction in the total recovery of peptides associated with Kb molecules purified from the surface of tapasin-deficient cells. Comparison of the peptides bound to Kb in the presence and absence of tapasin revealed considerable overlap in peptide repertoire. These results indicate that in the absence of an interaction with tapasin, Kb molecules fail to assemble with calreticulin and TAP, yet they are still capable of acquiring a diverse array of peptides. However, a significant proportion of these peptides appear to be suboptimal, resulting in reduced cell surface stability of Kb complexes. Taken together, the findings indicate that tapasin plays an essential role in the formation of the class I loading complex, which retains class I heterodimers in the endoplasmic reticulum until optimal ligand selection is completed.     

Temperature effects upon the time course of contraction and re-extension in Stentor coeruleus

A double-flash microphotographic technique has been used to follow the variation with temperature of the following kinetic parameters related to the contraction and re-extension of the ciliate Stentor coeruleus, namely the rate of contraction, the initiation time before contraction, the rate of re-extension and the initiation time before re-extension, all described by first order kinetics. Activation enthalpies, entropies and free energies related to the above mentioned parameters were calculated from the variation of the rate constants with temperature. The enthalpies and entropies appear to be of minor interest compared to the free energies. For the contraction and the initiation of contraction the delta G transition state values obtained were 14 and 15 kcal/mole, respectively, while the re-extension and the initiation of re-extension both were represented by a value of delta G transition state about 19 kcal/mole. These results are compared to activation parameters for different motile systems and for the formation and breakdown of ATP-myosin complexes. A model for the contraction and re-extension processes is proposed in accordance with the results measured.     

An evaluation of Poisson-Boltzmann electrostatic free energy calculations through comparison with experimental mutagenesis data

For systems involving highly and oppositely charged proteins, electrostatic forces dominate association and contribute to biomolecular complex stability. Using experimental or theoretical alanine-scanning mutagenesis, it is possible to elucidate the contribution of individual ionizable amino acids to protein association. We evaluated our electrostatic free energy calculations by comparing calculated and experimental data for alanine mutants of five protein complexes. We calculated Poisson-Boltzmann electrostatic free energies based on a thermodynamic cycle, which incorporates association in a reference (Coulombic) and solvated (solution) state, as well as solvation effects. We observe that Coulombic and solvation free energy values correlate with experimental data in highly and oppositely charged systems, but not in systems comprised of similarly charged proteins. We also observe that correlation between solution and experimental free energies is dependent on dielectric coefficient selection for the protein interior. Free energy correlations improve as protein dielectric coefficient increases, suggesting that the protein interior experiences moderate dielectric screening, despite being shielded from solvent. We propose that higher dielectric coefficients may be necessary to more accurately predict protein-protein association. Additionally, our data suggest that Coulombic potential calculations alone may be sufficient to predict relative binding of protein mutants.     

Somatic cell variants express altered H-2Kb allodeterminants recognized by cytolytic T cell clones

The present studies have made use of in vitro derived H-2Kb mutants to analyze the fine specificity of alloreactive cytotoxic T lymphocytes (CTL). The variants were derived by negatively selecting mutagenized tumor cells with a monoclonal anti-H-2Kb antibody and positively selecting for residual cells expressing serologically altered H-2Kb molecules. Details of this procedure are described in the companion paper. Selected populations of bulk alloreactive and cloned CTL were examined for recognition of the variants. In contrast to the serologic findings presented in the companion paper, there does not appear to be a correlation between the monoclonal antibody used to select the R8 variant and the CTL specificities recognized. In several instances, CTL clones could discriminate between variants having identical serologic profiles. Therefore, it would appear that the CTL have a large repertoire of allorecognition, even when generated across a mutant anti-Kb combination reflecting only a few amino acid differences. In addition, a diverse set of epitopes can be recognized on the Kb molecule. Finally, in some instances a change in what would appear to be a single amino acid resulted in a profound alteration of CTL recognition even though the Kb mutant molecule expressed limited serologic changes. These results support the idea that small changes in the H-2Kb molecule can have dramatic effects on CTL even though there are relatively little effects on serologic recognition of the target molecule.     

Computational method for relative binding energies of enzyme-substrate complexes.

A computational method for estimating the relative binding free energies of enzyme-substrate complexes is described that combines electrostatic and solvation models and X-ray crystallographic data. The polar contribution is evaluated by the Poisson-Boltzman equation. The nonpolar contribution is evaluated by solvent transfer data and surface area calculations. This algorithm was used to calculate the relative binding energies of 63 pairs of nine different mutant proteins with seven different substituted R-malate substrates of Escherichia coli isocitrate dehydrogenase. Comparison of calculated values with the experimentally observed values shows a high degree of correlation.     

In vitro peptide binding to soluble empty class I major histocompatibility complex molecules isolated from transfected Drosophila melanogaster cells.

A soluble form of a mouse class I major histocompatibility antigen (H-2Kb) has been expressed in transfected Drosophila melanogaster cells. These molecules were efficiently secreted (up to 4 mg/liter) as noncovalent heterodimers and purified to homogeneity from cell supernatants. The isolated soluble Kb molecules were devoid of endogenous peptides. Using these molecules, we have characterized the Kb heavy chain-beta 2-microglobulin (beta 2m) assembly as well as peptide binding in vitro. In detergent-free solution the heavy chains readily re-assembled with beta 2m even in the absence of peptides. Kinetic analyses showed that the peptide binding is rapid and reversible and dependent on the heavy chains being assembled with beta 2m. Likewise, peptide dissociated from Kb molecules without the displacement of beta 2m. Equilibrium binding experiments using various peptides confirmed that octapeptides bind to Kb molecules with the highest affinity and form the most stable complexes. However, in contrast to earlier studies, the amino-terminal positioning of peptide to Kb molecules was more crucial than the carboxyl-terminal positioning and amidation of the peptide carboxylate did not affect the binding. Soluble Kb molecules could selectively bind allele-specific peptides among a mixture of randomly synthesized octapeptides in vitro; however, no dominant residue was observed at the carboxyl terminus of bound peptides. This suggests that the previously observed hydrophobic residues at the carboxyl terminus of peptides may reflect the specificity of enzyme(s) or protein(s) involved in peptide processing in vivo.     

Lanthanide octaethylprophyrins: preparation, association, and interaction with axial ligands

The synthesis of virtually all the lanthanide octaethylporphyrin complexes have been achieved by heating appropriate anhydrous lanthanide halide and octaethylporphyrin in imidazole melt at 210 degrees C for two hours. The lighter lanthanide porphyrin complexes are very susceptible to hydrolysis, the middle lanthanide porphyrin complexes are moderately stable, and the heavier lanthanide porphyrin complexes are relatively more stable to hydrolysis. Two out of four lanthanide porphyrin complexes studied in detail, namely ytterbium and lutetium octaethylporphyrins, aggregate in benzene and the Soret bands in their absorption spectra are about 6 nm shifted to higher energies upon a hundred-fold increase in their concentrations. The aggregations of these lanthanide porphyrin complexes in non-coordinating solvents have been further verified by 1H NMR spectral studies. This spectral behavior can be interpreted qualitatively in terms of the model of the molecular exciton interactions with stacking of at least two prophyrins. A dimeric structure of these lanthanide porphyrin complexes has been proposed on the basis of geometrical considerations. On the contrary, the europium and gadolinium octaethylporphyrins associate very weakly in benzene in the concentration range studied. All four lanthanide porphyrin complexes interact with pyridine and piperidine, and the Soret bands in their absorption spectra are about 8 nm shifted to low energies as compared with their values in pure benzene.     

The solvation contribution to the binding energy of DNA with non-intercalating antibiotics.

The influence of the solvent on the binding energies to DNA of six non-intercalating antibiotics - netropsin, distamycin-3, distamycin-2, SN 18071, berenil and stilbamidine - is evaluated by combining the effect of the first hydration shell with that of bulk water. The first effect is computed by a methodology based on a spherical/point dipole model of water and limited to electrostatic interaction energies. Hydration shells are obtained which are energy optimized with respect to both water-solute and water-water interactions for the complexes and for the isolated DNA oligomers and ligands. The method allows even very large complexes to be studied in reasonable computation times. The second effect is introduced via a cavity treatment. It is shown that if the vacuum interaction energies already predict correctly the preference of the ligands for the minor groove of AT sequences of B-DNA, the introduction of the solvation effect is indispensable for reproducing the order of affinity of the ligands and for bringing the values of the complexation energies into close agreement with experimental data.     

A kinetic study of glucose-6-phosphate dehydrogenase.

The steady state kinetics of pig liver glucose-6-phosphate dehydrogenase is consistent with an ordered, sequential mechanism in which NADP is bound first and NADPH released last. Kia is 9.0 muM, Ka is 4.8 muM, and Kb is 36 muM. Glucosamine 6-phosphate, a substrate analogue and competitive inhibitor, is used to help rule out a possible random mechanism. ADP is seen to form a complex with the free form of the enzyme whereas ATP forms a complex with both the free and E-NADP forms of the enzyme. The KI for the E-ADP complex is 1.9 mM, while the Ki values for the E-ATP and E-NADP-ATP complexes are 7.2 and 4.5 mM, respectively.     

Analysis of cyclosporin A and a set of analogs as inhibitors of a T. cruzi cyclophilin by docking and molecular dynamics

Cyclophilins (CyPs) are enzymes involved in protein folding. In Trypanosoma cruzi (T. cruzi), the most abundantly expressed CyP is the isoform TcCyP19. It has been shown that TcCyP19 is inhibited by the immunosuppressive drug cyclosporin A (CsA) and analogs, which also proved to have potent trypanosomicidal activity in vitro. In this work, we continue and expand a previous study on the molecular interactions of CsA, and a set of analogs modeled in complexes with TcCyP19. The modeled complexes were used to evaluate binding free energies by molecular dynamics (MD), applying the Linear Interaction Energy (LIE) method. In addition, putative binding sites were identified by molecular docking. In our analysis, the binding free energy calculations did not correlate with experimental data. The heterogeneity of the non-bonded energies and the variation in the pattern of hydrogen bonds suggest that the systems may not be suitable for the application of the LIE method. Further, the docking calculations identified two other putative binding sites with comparable scoring energies to the active site, a fact that may also explain the lack of correlation found. Kinetic experiments are needed to confirm or reject the multiple binding sites hypothesis. In the meantime, MD simulations at the alternative sites, employing other methods to compute binding free energies, might be successful at finding good correlations with the experimental data.     

Influence of temperature on the in vitro activity of GDP-mannose dolicholphosphate mannosyltransferase in rat and trout liver microsomes

Temperature optimum of mannosyltransferase activity in liver microsomes is higher in trout than in rat, but this enzymatic activity for rat is higher than trout. Activation energies calculated for mannosyltransferase activity for trout and rat do not correlate with environmental temperature. For a given incubation temperature, Vm values for rat are higher than trout, whereas Km values for trout are lower than rat.     

Correlation of binding constants and molecular modelling of inhibitors in the active sites of aldose reductase and aldehyde reductase

Aldose reductase (ALR2) belongs to the aldo–keto reductase (AKR) superfamily of enzymes, is the first enzyme involved in the polyol pathway of glucose metabolism and has been linked to the pathologies associated with diabetes. Molecular modelling studies together with binding constant measurements for the four inhibitors Tolrestat, Minalrestat, quercetin and 3,5-dichlorosalicylic acid (DCL) were used to determine the type of inhibition, and correlate inhibitor potency and binding energies of the complexes with ALR2 and the homologous aldehyde reductase (ALR1), another member of the AKR superfamily. Our results show that the four inhibitors follow either uncompetitive or non-competitive inhibition pattern of substrate reduction for ALR1 and ALR2. Overall, there is correlation between the IC₅₀ (concentration giving 50% inhibition) values of the inhibitors for the two enzymes and the binding energies (ΔH) of the enzyme–inhibitor complexes. Additionally, the results agree with the detailed structural information obtained by X-ray crystallography suggesting that the difference in inhibitor binding for the two enzymes is predominantly mediated by non-conserved residues. In particular, Arg312 in ALR1 (missing in ALR2) contributes favourably to the binding of DCL through an electrostatic interaction with the inhibitor’s electronegative halide atom and undergoes a conformational change upon Tolrestat binding. In ALR2, Thr113 (Tyr116 in ALR1) forms electrostatic interactions with the fluorobenzyl moiety of Minalrestat and the 3- and 4-hydroxy groups on the phenyl ring of quercetin. Our modelling studies suggest that Minalrestat’s binding to ALR1 is accompanied by a conformational change including the side chain of Tyr116 to achieve the selectivity for ALR1 over ALR2.     

Empirical calculation of the relative free energies of peptide binding to the molecular chaperone DnaK

We describe a methodology to calculate the relative free energies of protein-peptide complex formation. The interaction energy was decomposed into nonpolar, electrostatic and entropic contributions. A free energy-surface area relationship served to calculate the nonpolar free energy term. The electrostatic free energy was calculated with the finite difference Poisson-Boltzmann method and the entropic contribution was estimated from the loss in the conformational entropy of the peptide side chains. We applied this methodology to a series of DnaK*peptide complexes. On the basis of the single known crystal structure of the peptide-binding domain of DnaK with a bound heptapeptide, we modeled ten other DnaK*heptapeptide complexes with experimentally measured K(d) values from 0.06 microM to 11 microM, using molecular dynamics to refine the structures of the complexes. Molecular dynamic trajectories, after equilibration, were used for calculating the energies with greater accuracy. The calculated relative binding free energies were compared with the experimentally determined free energies. Linear scaling of the calculated terms was applied to fit them to the experimental values. The calculated binding free energies were between -7.1 kcal/mol and - 9.4 kcal/mol with a correlation coefficient of 0.86. The calculated nonpolar contributions are mainly due to the central hydrophobic binding pocket of DnaK for three amino acid residues. Negative electrostatic fields generated by the protein increase the binding affinity for basic residues flanking the hydrophobic core of the peptide ligand. Analysis of the individual energy contributions indicated that the nonpolar contributions are predominant compared to the other energy terms even for peptides with low affinity and that inclusion of the change in conformational entropy of the peptide side chains does not improve the discriminative power of the calculation. The method seems to be useful for predicting relative binding energies of peptide ligands of DnaK and might be applicable to other protein-peptide systems, particularly if only the structure of one protein-ligand complex is available.     

Structural prediction of peptides bound to MHC class I

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values < or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values < or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.