Biomedical subjects

R C Wade

Publications and source records attributed to R C Wade.

At least 19 recordsLinked to original sources

PubMed1998

Hydration energy landscape of the active site cavity in cytochrome P450cam.

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein-ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol-1) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable.

Binding Sites↗

PubMed1998

Species dependence of enzyme-substrate encounter rates for triose phosphate isomerases.

Triose phosphate isomerase (TIM) is a diffusion-controlled enzyme whose rate is limited by the diffusional encounter of the negatively charged substrate glyceraldehyde 3-phosphate (GAP) with the homodimeric enzyme's active sites. Translational and orientational steering of GAP toward the active sites by the electrostatic field of chicken muscle TIM has been observed in previous Brownian dynamics (BD) simulations. Here we report simulations of the association of GAP with TIMs from four species with net charges at pH 7 varying from -12e to +12e. Computed second-order rate constants are in good agreement with experimental data. The BD simulations and computation of average Boltzmann factors of substrate-protein interaction energies show that the protein electrostatic potential enhances the rates for all the enzymes. There is much less variation in the computed rates than might be expected on the basis of the net charges. Comparison of the electrostatic potentials by means of similarity indices shows that this is due to conservation of the local electrostatic potentials around the active sites which are the primary determinants of electrostatic steering of the substrate.

Animals↗

PubMed1998

Electrostatic steering and ionic tethering in enzyme-ligand binding: insights from simulations.

To bind at an enzyme's active site, a ligand must diffuse or be transported to the enzyme's surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and beta-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as "ionic tethering." We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme's surroundings even when the substrate is nonpolar.

Computer Simulation↗

PubMed1998

Brownian dynamics simulation of protein-protein diffusional encounter.

Protein association events are ubiquitous in biological systems. Some protein associations and subsequent responses are diffusion controlled in vivo. Hence, it is important to be able to compute bimolecular diffusional association rates for proteins. The Brownian dynamics simulation methodology may be used to simulate protein-protein encounter, compute association rates, and examine their dependence on protein mutation and the nature of the physical environment (e.g., as a function of ionic strength or viscosity). Here, the theory for Brownian dynamics simulations is described, and important methodological aspects, particularly pertaining to the correct modeling of electrostatic forces and definition of encounter complex formation, are highlighted. To illustrate application of the method, simulations of the diffusional encounter of the extracellular ribonuclease, barnase, and its intracellular inhibitor, barstar, are described. This shows how experimental rates for a series of mutants and the dependence of rates on ionic strength can be reproduced well by Brownian dynamics simulations. Potential future uses of the Brownian dynamics method for investigating protein-protein association are discussed.

Bacterial Proteins↗

PubMed1998

Importance of explicit salt ions for protein stability in molecular dynamics simulation.

The accurate and efficient treatment of electrostatic interactions is one of the challenging problems of molecular dynamics simulation. Truncation procedures such as switching or shifting energies or forces lead to artifacts and significantly reduced accuracy. The particle mesh Ewald (PME) method is one approach to overcome these problems by providing a computationally efficient means of calculating all long-range electrostatic interactions in a periodic simulation box by use of fast Fourier transformation techniques. For the application of the PME method to the simulation of a protein with a net charge in aqueous solution, counterions are added to neutralize the system. The usual procedure is to add charge-balancing counterions close to charged residues to neutralize the protein surface. In the present article, we show that for MD simulation of a small protein of marginal stability, the YAP-WW domain, explicit modeling of 0.2 M ionic strength (in addition to the charge-balancing counterions) is necessary to maintain a stable protein structure. Without explicit ions throughout the periodic simulation box, the charge-balancing counterions on the protein surface diffuse away from the protein, resulting in destruction of the beta-sheet secondary structure of the WW domain.

Computer Simulation↗

PubMed1998

Classification of auxin plant hormones by interaction property similarity indices.

Although auxins were the first type of plant hormone to be identified, little is known about the molecular mechanism of this important class of plant hormones. We present a classification of a set of about 50 compounds with measured auxin activities, according to their interaction properties. Four classes of compounds were defined: strongly active, weakly active with weak antiauxin behaviour, inactive and inhibitory. All compounds were modeled in two low-energy conformations, 'P' and 'T', so as to obtain the best match to the 'planar' and 'tilted' conformations, respectively, of indole 3-acetic acid. Each set of conformers was superimposed separately using several different alignment schemes. Molecular interaction energy fields were computed for each molecule with five different chemical probes and then compared by computing similarity indices. Similarity analysis showed that the classes are on average distinguishable, with better differentiation achieved for the T conformers than the P conformers. This indicates that the T conformation might be the active one. Further, a screening was developed which could distinguish compounds with auxin activity from inactive compounds and most antiauxins using the T conformers. The classifications rationalize ambiguities in activity data found in the literature and should be of value in predicting the activities of new plant growth substances and herbicides.

Computer Simulation↗

PubMed1998

Prediction of protein hydration sites from sequence by modular neural networks.

The hydration properties of a protein are important determinants of its structure and function. Here, modular neural networks are employed to predict ordered hydration sites using protein sequence information. First, secondary structure and solvent accessibility are predicted from sequence with two separate neural networks. These predictions are used as input together with protein sequences for networks predicting hydration of residues, backbone atoms and sidechains. These networks are trained with protein crystal structures. The prediction of hydration is improved by adding information on secondary structure and solvent accessibility and, using actual values of these properties, residue hydration can be predicted to 77% accuracy with a Matthews coefficient of 0.43. However, predicted property data with an accuracy of 60-70% result in less than half the improvement in predictive performance observed using the actual values. The inclusion of property information allows a smaller sequence window to be used in the networks to predict hydration. It has a greater impact on the accuracy of hydration site prediction for backbone atoms than for sidechains and for non-polar than polar residues. The networks provide insight into the mutual interdependencies between the location of ordered water sites and the structural and chemical characteristics of the protein residues.

Amino Acid Sequence↗

PubMed1997

'Flu' and structure-based drug design.

The threat of a catastrophic outbreak of influenza is ever present. Vaccines are only partially effective and the two compounds, amantidine and rimantidine, used clinically against influenza A cause side-effects and rapid viral resistance. Recent advances bring hope that specific and potent drugs against influenza may soon be available in the clinic. These compounds were designed to inhibit influenza neuraminidase (NA), one of the viral coat glycoproteins, using the crystal structure of NA which was first published in 1983. In this review, the application of structure-based drug design approaches to the design of anti-influenza agents targeted at NA and haemagglutinin (HA), the other viral surface glycoprotein, is discussed.

Amines↗

PubMed1997

Exceptionally stable salt bridges in cytochrome P450cam have functional roles.

A long-standing puzzle in structure-function studies of cytochrome P450cam is how the substrate, camphor, reaches the buried active site. The crystal structure shows no channel from the surface to the active site large enough for substrate to pass through. Recent experiments indicate that access of the rather nonpolar substrate to the active site is controlled by electrostatic interactions and may involve rupture of the two salt links to Asp251 [Deprez, E., Gerber, N. C., Di Primo, C., Douzou, P., Sligar, S. G., & Hui Bon Hoa, G. (1994) Biochemistry 33, 14464-14468]. Consequently, we have computed the electrostatic strength of 53 ionic pairs, including 32 salt links, in cytochrome P450cam by numerical solution of the finite-difference linearized Poisson-Boltzmann equation. The calculated electrostatic free energies, delta Gtot, of the salt links range from -9 to +6 kcal/mol with approximately 60% of the salt links being energetically favorable and 40% being unfavorable with respect to mutation to their uncharged, nonpolar isosteres. Strikingly, of the four most stable salt links in the protein (delta Gtot < -6 kcal/mol), two involve the propionate groups of the heme and the other two involve Asp251. In the modeled D251N mutant, for which electrostatic effects on substrate binding are diminished, the latter two salt links lose their stability (delta Gtot > -2.4 kcal/mol). Thus it appears that cytochrome P450cam has evolved four unusually strong salt bridges, stabilized by surrounding charged and polar groups in the protein, to keep its heme cofactor in place and to regulate substrate binding.

Binding Sites↗

PubMed1997

Reliability of comparative molecular field analysis models: effects of data scaling and variable selection using a set of human synovial fluid phospholipase A2 inhibitors.

The effects of data pretreatment, data scaling, and variable selection on three-dimensional quantitative structure-activity relationships derived by comparative molecular field analysis (CoMFA) using the GRID energy function were studied in detail for a set of inhibitors of the human synovial fluid phospholipase A2 (HSF-PLA2). The quality of the models was evaluated for predictive power and ability to map the receptor binding site by (a) comparison of predicted and experimental activities using cross-validation and external validation sets and (b) comparison of the regions selected in space in the CoMFA models with a crystal structure of a HSF-PLA2-inhibitor complex, with optimized comparative binding energy analysis (COMBINE) models (Ortiz et al., 1995) and with structure-activity relationships derived previously for different sets of compounds. It is found that (1) data scaling and dielectric modeling strongly influence CoMFA results. Unscaled data and a uniform dielectric constant of 4 are well suited to GRID-CoMFA studies for the present compound set. (2) The GOLPE and Q2-GRS variable selection methods select variables in roughly the same regions in Cartesian space, but they produce different models in chemometric space and differ in their sensitivity to data scaling and pretreatment and their tendency to overfitting. (3) CoMFA models are consistent with COMBINE models in that they identify approximately the same intermolecular interactions as relevant for activity. Our study provides support for the qualitative receptor-mapping properties of CoMFA models and for the validity of variable selection when applied with care and also provides guidelines for how to evaluate the quality of CoMFA models.

Enzyme Inhibitors↗

PubMed1997

CASP2 molecular docking predictions with the LIGIN software.

Seven docking predictions were made with the LIGIN program. In six cases the location of the binding pocket was identified correctly by systematically docking everywhere within the protein structure. In two cases the ligand was docked to within 1.8 A RMSD of the experimentally determined structure. LIGIN has not been optimized to deal with highly flexible ligands that dock at the surface of proteins. Consequently, in three cases the exposed part of the ligand was docked poorly, although the buried parts were docked well, and made similar atomic contacts with the protein as in the experimentally determined structure.

Amiloride↗

PubMed1997

pKa calculations for class A beta-lactamases: methodological and mechanistic implications.

Beta-lactamases are responsible for resistance to penicillins and related beta-lactam compounds. Despite numerous studies, the identity of the general base involved in the acylation step is still unclear. It has been proposed, on the basis of a previous pKa calculation and analysis of structural data, that the unprotonated Lys73 in the active site could act as the general base. Using a continuum electrostatic model with an improved treatment of the multiple titration site problem, we calculated the pKa values of all titratable residues in the substrate-free TEM-1 and Bacillus licheniformis class A beta-lactamases. The pKa of Lys73 in both enzymes was computed to be above 10, in good agreement with recent experimental data on the TEM-1 beta-lactamase, but inconsistent with the proposal that Lys73 acts as the general base. Even when the closest titratable residue, Glu166, is mutated to a neutral residue, the predicted downward shift of the pKa of Lys73 shows that it is unlikely to act as a proton abstractor in either enzyme. These results support a mechanism in which the proton of the active Ser70 is transferred to the carboxylate group of Glu166.

Acylation↗

PubMed1997

Simulation of the diffusional association of barnase and barstar.

The rate of protein association places an upper limit on the response time due to protein interactions, which, under certain circumstances, can be diffusion-controlled. Simulations of model proteins show that diffusion-limited association rates are approximately 10(6)-10(7) M-1 s-1 in the absence of long-range forces (Northrup, S. H., and H. P. Erickson. 1992. Kinetics of protein-protein association explained by Brownian dynamics computer simulations. Proc. Natl. Acad. Sci. U.S.A. 89:3338-3342). The measured association rates of barnase and barstar are 10(8)-10(9) M-1 s-1 at 50 mM ionic strength, and depend on ionic strength (Schreiber, G., and A. R. Fersht. 1996. Rapid, electrostatically assisted association of proteins. Nat. Struct. Biol. 3:427-431), implying that their association is electrostatically facilitated. We report Brownian dynamics simulations of the diffusional association of barnase and barstar to compute association rates and their dependence on ionic strength and protein mutation. Crucial to the ability to reproduce experimental rates is the definition of encounter complex formation at the endpoint of diffusional motion. Simple definitions, such as a required root mean square (RMS) distance to the fully bound position, fail to explain the large influence of some mutations on association rates. Good agreement with experiments could be obtained if satisfaction of two intermolecular residue contacts was required for encounter complex formation. In the encounter complexes, barstar tends to be shifted from its position in the bound complex toward the guanine-binding loop on barnase.

Algorithms↗

PubMed1997

Local interactions of aromatic residues in short peptides in aqueous solution: a combined database and energetic analysis.

Although short peptides are usually structurally disordered in aqueous solution, particular peptide sequences display local structure. We performed database and conformational searches, along with molecular dynamics simulations, to study two local interactions detected by 1H-NMR in tetrapeptides excised from bovine pancreatic trypsin inhibitor: aromatic-(i+2)amide and (i-1)cisproline-aromatic. For both types of interaction, at least two major and distinct peptide conformations are identified in the folded state. The aromatic-(i+2)amide interaction can have parallel and perpendicular arrangements of the N-H bond and the aromatic ring. The (i-1)cisproline-aromatic interaction can have close packing of the aromatic ring to the (i-2)C alpha H or the (i-1)C gamma H but not both simultaneously. Although these local aromatic interactions are weak, they may influence folding and binding properties. The combination of search and simulation techniques provides a useful route towards obtaining an atomic-detail description of peptides exhibiting these interactions.

Amides↗

PubMed1996

Structural changes in cytochrome P-450cam effected by the binding of the enantiomers (1R)-camphor and (1S)-camphor.

A comparative study of the enantiomeric substrate [(1R)-camphor- and (1S)-camphor)-bound cytochrome P-450cam concerns the spin-state equilibrium, substrate dissociation, the thermal unfolding of the protein structure, and the subconformer equilibria observed in the infrared spectra of the carbon monoxide (CO) complex of cytochrome P-450cam. The behavior of the different conformational equilibria in dependence on temperature, pressure, pH-value, cosolvent, and cation binding led us to suggest that (1S)-camphor is more loosely and less optimally bound in the heme pocket, which facilitates the access of solvent molecules into the heme-iron environment. The spin reaction volume difference measured using the high pressure technique is smaller by 16 +/- 9 cm3/mol for (1S)-camphor-bound P-450cam compared to the (1R)-camphor-bound P-450cam, which might indicate a higher water content in the protein and in the heme environment in the (1S)-camphor complex. The half-transition temperature of the thermal unfolding of 53.8 degrees C for the (1S)-camphor-bound oxidized cytochrome P-450cam is one degree lower than the value for the (1R)-camphor-bound protein (54.8 degrees C). In the reduced, CO-bound form of cytochrome P-450cam at 290 K the (1S)-camphor complex reveals another CO stretch vibration population distribution with slightly higher frequencies [1940.2 cm-1 (major band) and 1946.3 cm-1 (minor band)] compared to the (1R)-camphor complex [1939.7 cm-1 (major band) and 1930 cm-1 (minor band)]. A loosening of the contact between the iron-bound CO ligand and amino acids of the I-helix, probably induced by compensating effects of the increased water content, is suggested. Assuming the carbon monoxide complex as a model for the dioxygen complex, the more loosened binding of (1S)-camphor, therefore the increased water accessibility, and the weaker contact of the iron ligand to the I-helix might explain the higher amount of uncoupling of the cytochrome P-450 reaction cycle compared to that when (1R)-camphor is used as substrate.

Binding Sites↗

PubMed1996

Improved binding of cytochrome P450cam substrate analogues designed to fill extra space in the substrate binding pocket.

Cytochrome P450cam catalyzes the 5-exo-hydroxylation of camphor. Camphor analogues were designed to fill an empty region of the substrate binding pocket with the expectation that they would bind more tightly than camphor itself due to increased van der Waals interactions with the protein and the displacement of any solvent occupying this site. A series of compounds (endo-borneol methyl ether, endo-borneol propyl ether, endo-borneol allyl ether and endo-borneol dimethyl allyl ether) were synthesized with substituents at the camphor carbonyl oxygen. The spin conversion and thermodynamic properties of this series of compounds were measured for wild type and Y96F mutant cytochrome P450cam and were interpreted in the context of molecular dynamics simulations of the camphor analogues in the P450 binding site and in solution. Compounds with a 3-carbon chain substituent were predicted to match the size of the unoccupied region most optimally and thus bind best. Consistent with this prediction, the borneol allyl ether binds to cytochrome P450cam with highest affinity with a Kd = 0.6 +/- 0.1 microM (compared to a Kd = 1.7 +/- 0.2 microM for camphor under the same experimental conditions). Binding of the camphor analogues to the Y96F mutant is much enhanced over the binding of camphor, indicating that hydrogen bonding plays a less important role in binding of these analogues. Binding enthalpies calculated from the simulations, taking all solvent contributions into account, agree very well with experimental binding enthalpies. Binding affinity is not however correlated with the calculated binding enthalpy because the binding of the substrate analogues is characterized by enthalpy-entropy compensation. The new compounds are useful probes for further studies of the mechanism of cytochrome P450cam due to their high binding affinities and high spin properties.

Binding Sites↗

PubMed1996

Molecular docking using surface complementarity.

A method is described to dock a ligand into a binding site in a protein on the basis of the complementarity of the intermolecular atomic contacts. Docking is performed by maximization of a complementarity function that is dependent on atomic contact surface area and the chemical properties of the contacting atoms. The generality and simplicity of the complementarity function ensure that a wide range of chemical structures can be handled. The ligand and the protein are treated as rigid bodies, but displacement of a small number of residues lining the ligand binding site can be taken into account. The method can assist in the design of improved ligands by indicating what changes in complementarity may occur as a result of the substitution of an atom in the ligand. The capabilities of the method are demonstrated by application to 14 protein-ligand complexes of known crystal structure.

Algorithms↗

PubMed1996

Theoretical investigation of the dynamics of the active site lid in Rhizomucor miehei lipase.

Interfacial activation of Rhizomucor miehei lipase is accompanied by a hinge-type motion of a single helix (residues 83-94) that acts as a lid over the active site. Activation of the enzyme involves the displacement of the lid to expose the active site, suggesting that the dynamics of the lid could be of mechanistic and kinetic importance. To investigate possible activation pathways and to elucidate the effect of a hydrophobic environment (as would be provided by a lipid membrane) on the lid opening, we have applied molecular dynamics and Brownian dynamics techniques. Our results indicate that the lipase activation is enhanced in a hydrophobic environment. In nonpolar low-dielectric surroundings, the lid opens in approximately 100 ns in the BD simulations. In polar high-dielectric (aqueous) surroundings, the lid does not always open up in simulations of up to 900 ns duration, but it does exhibit some gating motion, suggesting that the enzyme molecule may exist in a partially active form before the catalytic reaction. The activation is controlled by the charged residues ARG86 and ASP91. In the inactive conformation, ASP91 experiences repulsive forces and pushes the lid toward the open conformation. Upon activation ARG86 approaches ASP61, and in the active conformation, these residues form a salt bridge that stabilizes the open conformation.

Binding Sites↗