Search PubMed⌕ Search

Biomedical subjects

R C Wade

Publications and source records attributed to R C Wade.

At least 55 records · Page 3Linked to original sources
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↗
PubMed1996

TRAJAN: a tool for analyzing trajectories from molecular simulations.

Molecular dynamics simulations of biological systems are notoriously difficult to analyze because of the complexity of the information that contain. We describe a new method for analyzing trajectories from simulations in order to extract important features of the motion. The trajectory of each atom is partitioned into conformation wells, in each of which its motion is assumed to be predominantly harmonic oscillation around an average position. Thus each atom moves anharmonically through a sequence of wells during the trajectory. The movement of atoms between wells, their ellipsoids of motion within each well, and correlations in the motion of atoms are quantified and can be visualized with molecular graphics. The TRAJAN analysis procedure is applicable to trajectories from equilibrium and nonequilibrium simulations and is not restricted to molecular dynamics simulations. Its application is demonstrated for a range of model systems.

Algorithms↗
PubMed1996

Analytically defined surfaces to analyze molecular interaction properties.

Molecular surfaces are widely used for characterizing molecules and displaying and quantifying their interaction properties. Here we consider molecular surfaces defined as isocontours of a function (a sum of exponential functions centered on each atom) that approximately represents electron density. The smoothness is advantageous for surface mapping of molecular properties (e.g., electrostatic potential). By varying parameters, these surfaces can be constructed to represent the van der Waals or solvent-accessible surface of a molecular with any accuracy. We describe numerical algorithms to operate on the analytically defined surfaces. Two applications are considered: (1) We define and locate extremal points of molecular properties on the surfaces. The extremal points provide a compact representation of a property on a surface, obviating the necessity to compute values of the property on an array of surface points as is usually done; (2) a molecular surface patch or interface is projected onto a flat surface (by introducing curvilinear coordinates) with approximate conservation of area for analysis purposes. Applications to studies of protein-protein interactions are described.

Computer Graphics↗
PubMed1996

Double-blind parallel study of a combination of chlorthalidone 50 mg and triamterene 50 mg in patients with mild and moderate hypertension.

In a double-blind, partial cross-over study of 141 patients with mild to moderate hypertension, chlorthalidone 50 mg was compared with chlorthalidone 50 mg/triamterene 50 mg combination treatment. The study assessed clinical efficacy in the reduction of blood pressure, and safety with regard to serum potassium concentration and adverse reactions. Both treatment regimens were found to offer equivalent therapeutic success, but the decrease in serum potassium concentration was significantly less with the combination treatment than with chlorthalidone alone.

Adult↗
PubMed1995

Prediction of drug binding affinities by comparative binding energy analysis.

A new computational method for deducing quantitative structure-activity relationships (QSARs) using structural data from ligand-macromolecule complexes is presented. First, the ligand-macromolecule interaction energy is computed for a set of ligands using molecular mechanics calculations. Then, by selecting and scaling components of the ligand-macromolecule interaction energy that show good predictive ability, a regression equation is obtained in which activity is correlated with the interaction energies of parts of the ligands and key regions of the macromolecule. Application to the interaction of the human synovial fluid phospholipase A2 with 26 inhibitors indicates that the derived QSAR has good predictive ability and provides insight into the mechanism of enzyme inhibition. The method, which we term comparative binding energy (COMBINE) analysis, is expected to be applicable to ligand-receptor interactions in a range of contexts including rational drug design, host-guest systems, and protein engineering.

Binding Sites↗
PubMed1995

Thermodynamics of water mediating protein-ligand interactions in cytochrome P450cam: a molecular dynamics study.

Ordered water molecules are observed by crystallography and nuclear magnetic resonance to mediate protein-ligand interactions. Here, we examine the energetics of hydrating cavities formed at protein-ligand interfaces using molecular dynamics simulations. The free energies of hydrating two cavities in the active site of two liganded complexes of cytochrome P450cam were calculated by multiconfigurational thermodynamic integration. The complex of cytochrome P450cam with 2-phenyl-imidazole contains a crystallographically well defined water molecule mediating hydrogen bonds between the protein and the inhibitor. We calculate that this water molecule is stabilized by a binding free energy of -11.6 +/- kJ/mol. The complex of cytochrome P450cam with its natural substrate, camphor, contains a cavity that is empty in the crystal structure although a water molecule in it could make a hydrogen bond to camphor. Here, solvation of this cavity is calculated to be unfavorable by +15.8 +/- 5.0 kJ/mol. The molecular dynamics simulations can thus distinguish a hydrated interfacial cavity from an empty one. They also provide support for the notion that protein-ligand complexes can accommodate empty interfacial cavities and that such cavities are likely to be unhydrated unless more than one hydrogen bond can be made to a water molecule in the cavity.

Binding Sites↗
PubMed1995

Simulation of enzyme-substrate interactions: the diffusional encounter step.

Enzymes are targets for structure-based redesign and drug design. One potentially important design strategy is to create molecules--enzymes and ligands--with desired diffusional encounter properties. Rates of diffusive bimolecular encounter can be calculated by Brownian dynamics simulation. This methodology and its application to study the factors influencing the rates of diffusion-influenced enzymes are described here.

Diffusion↗
PubMed1994

Receptor binding properties of four-helix-bundle growth factors deduced from electrostatic analysis.

Hormones of the hematopoietin class mediate signal transduction by binding to specific transmembrane receptors. Structural data show that the human growth hormone (hGH) forms a complex with a homodimeric receptor and that hGH is a member of a class of hematopoietins possessing an antiparallel 4-alpha-helix bundle fold. Mutagenesis experiments suggest that electrostatic interactions may have an important influence on hormone-receptor recognition. In order to examine the specificity of hormone-receptor complexation, an analysis was made of the electrostatic potentials of hGH, interleukin-2 (IL-2), interleukin-4 (IL-4), granulocyte colony-stimulating factor (G-CSF), granulocyte-macrophage colony-stimulating factor (GM-CSF), and the hGH and IL-4 receptors. The binding surfaces of hGH and its receptor, and of IL-4 and its receptor, show complementary electrostatic potentials. The potentials of the hGH and its receptor display approximately 2-fold rotational symmetry because the receptor subunits are identical. In contrast, the potentials of GM-CSF and IL-2 lack such symmetry, consistent with their known high affinity for hetero-oligomeric receptors. Analysis of the electrostatic potentials supports a recently proposed hetero-oligomeric model for a high-affinity IL-4 receptor and suggests a possible new receptor binding mode for G-CSF; it also provides valuable information for guiding structural and mutagenesis studies of signal-transducing proteins and their receptors.

Electrochemistry↗
PubMed1994

Homology modeling of the Abl-SH3 domain.

A tertiary structure model of the Abl-SH3 domain is predicted by using homology modeling techniques coupled to molecular dynamics simulations. Two template proteins were used, Fyn-SH3 and Spc-SH3. The refined model was extensively checked for errors using criteria based on stereochemistry, packing, solvation free-energy, accessible surface areas, and contact analyses. The different checking methods do not totally agree, as each one evaluates a different characteristic of protein structures. Several zones of the protein are more susceptible to incorporating errors. These include residues 13, 15, 35, 39, 45, 46, 50, and 60. An interesting finding is that the measurement of the C alpha chirality correlated well with the rest of the criteria, suggesting that this parameter might be a good indicator of correct local conformation. Deviations of more than 4 degrees may be indicative of poor local structure.

Computer Simulation↗
PubMed1994

Simulation of enzyme-substrate encounter with gated active sites.

We describe a brownian dynamics simulation method that allows investigation of the effects of receptor flexibility on ligand binding rates. The method is applied to the encounter of substrate, glyceraldehyde 3-phosphate, with triose phosphate isomerase, a diffusion-controlled enzyme with flexible peptide loops at its active sites. The simulations show that while the electrostatic field surrounding the enzyme steers the substrate into its active sites, the flexible loops appear to have little influence on the substrate binding rate. The dynamics of the loops may therefore have been optimized during evolution to minimize their interference with the substrate's access to the active sites. The calculated and experimental rate constants are in good agreement.

Animals↗
PubMed1993

Further development of hydrogen bond functions for use in determining energetically favorable binding sites on molecules of known structure. 1. Ligand probe groups with the ability to form two hydrogen bonds.

The directional properties of hydrogen bonds play a major role in determining the specificity of intermolecular interactions. An energy function which takes explicit account of these properties has been developed for use in the determination of energetically favorable ligand binding sites on molecules of known structure by the GRID method (Goodford, P.J.J. Med. Chem. 1985, 28, 849. Boobbyer, D.N.A.; Goodford, P.J.; McWhinnie, P.M.; Wade, R.C.J. Med. Chem. 1989, 32, 1083). In this method, the interaction energy between a target molecule and a small chemical group (a probe), which may be part of a larger ligand, was calculated using an energy function consisting of Lennard-Jones, electrostatic, and hydrogen bond terms. The latter term was a function of the length of the hydrogen bond, its orientation at the hydrogen-bonding atoms, and their chemical nature. We now describe hydrogen bond energy functions which take account of the spatial distribution of the hydrogen bonds made by probes with the ability to form two hydrogen bonds. These functions were designed so as to model the experimentally observed angular dependence of the hydrogen bonds. We also describe the procedure to locate the position and orientation of the probe at which the interaction energy is optimized. The use of this procedure is demonstrated by examples of biological and pharmacological interest which show that it can produce results that are consistent with other theoretical approaches and with experimental observations.

Binding Sites↗
PubMed1993

Further development of hydrogen bond functions for use in determining energetically favorable binding sites on molecules of known structure. 2. Ligand probe groups with the ability to form more than two hydrogen bonds.

The specificity of interactions between biological macromolecules and their ligands may be partially attributed to the directional properties of hydrogen bonds. We have now extended the GRID method (Goodford, P. J. J. Med. Chem. 1985, 28, 849. Boobbyer, D. N. A.; Goodford, P. J.; McWhinnie, P. M.; Wade, R. C. J. Med. Chem. 1989, 32, 1083), of determining energetically favorable ligand binding sites on molecules of known structure, in order to improve the treatment of groups which can make multiple hydrogen bonds. In this method, the interaction energy between a probe (a small chemical group that may be part of a larger ligand) and a target molecule is calculated using an energy function which includes a hydrogen bond term which is dependent on the length of the hydrogen bond, its orientation at the hydrogen-bonding atoms, and their chemical character. The methods described in the preceding paper (Wade, R. C.; Clark, K. J.; Goodford, P. J. J. Med. Chem., preceding paper in this issue) for probes capable of making two hydrogen bonds are here extended to the following probes which have the ability to make more than two hydrogen bonds: ammonium-NH3+, amine-NH2, sp3-hybridized hydroxyl, and water. Use of the improved GRID procedure is demonstrated by the determination of the conformation of an amino acid side chain at the subunit interface in hemoglobin and of the location of water binding sites in human lysozyme.

Binding Sites↗
PubMed1993

Gating of the active site of triose phosphate isomerase: Brownian dynamics simulations of flexible peptide loops in the enzyme.

The enzyme triose phosphate isomerase has flexible peptide loops at its active sites. The loops close over these sites upon substrate binding, suggesting that the dynamics of the loops could be of mechanistic and kinetic importance. To investigate these issues, the loop motions in the dimeric enzyme were simulated by Brownian dynamics. The two loops, one on each monomer, were represented by linear chains of appropriately parameterized spheres, each sphere corresponding to an amino acid residue. The loops moved in the electrostatic field of the rest of the enzyme, which was held rigid in its crystallographically observed conformation. In the absence of substrate, the loops exhibited gating of the active site with a period of about 1 ns and occupied "closed" conformations for about half of the time. As the period of gating is much shorter than the enzyme-substrate relaxation time, the motion of the loops does not reduce the rate constant for the approach of substrate from its simple diffusion-controlled value. This suggests that the flexible loops may have evolved to create the appropriate environment for catalysis while, at the same time, minimizing the kinetic penalty for gating the active site.

Amino Acid Sequence↗
PubMed1992

Binding of an antiviral agent to a sensitive and a resistant human rhinovirus. Computer simulation studies with sampling of amino acid side-chain conformation. I. Mapping the rotamers of residue 188 of viral protein 1.

The mutation of valine 188 to leucine in the viral protein 1 of human rhinovirus 14 renders the virus resistant to certain antiviral compounds. Thermodynamic-cycle perturbation theory provides a means of calculating the difference in the binding free energies of an antiviral compound to the wild-type virus and to the mutant virus. In calculating the relevant free-energy differences in molecular dynamics simulations, it is important to sample the multiple rotational isomers of residue 188 correctly. In general, these rotamers will not be fully sampled during a single molecular dynamics simulation. However, the contributions of all the rotamers to the free-energy differences associated with mutation of residue 188 may be considered explicitly once they have been identified and their relative free energies determined. Therefore, we describe here the mapping of the rotamers of residue 188 by steric-bump search and energy minimization techniques, and by the computation of potentials of mean force (p.m.f.s.) using umbrella sampling. The usefulness, validity and efficiency of these methods of examining rotameric states is discussed. Adiabatic mapping by energy minimization was found to be unreliable for this residue due to the small magnitude of its interactions with the surrounding protein atoms. Ambiguities in the adiabatic maps were resolved by computing p.m.f.s. The p.m.f. for valine 188 in the unliganded wild-type virus shows a minimum corresponding to the crystallographically observed conformation of valine 188. The p.m.f.s. for valine 188 in the liganded virus and for leucine 188 in the unliganded mutant virus suggest that the experimentally observed conformations may be interpreted as averages of a number of conformations corresponding to those at the minima in the p.m.f.s. The calculations suggest also that the conformation of leucine 188 may change when the ligand binds. The use of the calculated p.m.f.s. to compute the difference in the free energy of binding of an antiviral compound to the wild-type and mutant rhinoviruses is described in the accompanying article.

Antiviral Agents↗