Search PubMed⌕ Search

Biomedical subjects

R C Wade

Publications and source records attributed to R C Wade.

64 records · Page 4Linked to original sources
PubMed1992

Binding of an antiviral agent to a sensitive and a resistant human rhinovirus. Computer simulation studies with sampling of amino acid side-chain conformations. II. Calculation of free-energy differences by thermodynamic integration.

Thermodynamic-cycle perturbation theory and molecular dynamics simulations were used to calculate the difference in the free energy of binding of the antiviral compound WIN53338 to the wild-type human rhinovirus 14 and to a drug-resistant mutant of the virus in which valine 188 of the viral protein 1 is mutated to leucine. Because of the difficulty of achieving adequate sampling of all of the rotational isomers of amino acid side-chains in molecular dynamics simulations, an explicit treatment of the effects of the existence of multiple rotational isomers of residue 188 on the calculated free energies was used. The rotamers of residue 188 were first mapped by steric and energetic techniques as described in the accompanying article. Thermodynamic integration was then carried out during simulations of the virus, both with and without the antiviral compound bound, by mutating residue 188 while restraining its side-chain to one conformation. The contributions of the other rotamers of residue 188 to the free-energy changes for this mutation were then added to those calculated by thermodynamic integration as correction factors. Binding of WIN53338 to the wild-type virus was calculated to be favored over binding to the mutant virus by 1.7(+/- 3.0) kcal/mol. This is consistent with experimental data which, if differences in activity are assumed to be due to differences in binding, indicate that the binding affinity of WIN53338 for the wild-type virus is at least 0.15 to 1.7 kcal/mol greater than for the mutant virus. Thermodynamic integration was also performed in the conventional manner without restraints and was found to give less accurate results.

Antiviral Agents↗
PubMed1991

A molecular dynamics study of thermodynamic and structural aspects of the hydration of cavities in proteins.

The structure and activity of a protein molecule are strongly influenced by the extent of hydration of its cavities. This is, in turn, related to the free energy change on transfer of a water molecule from bulk solvent into a cavity. Such free energy changes have been calculated for two cavities in a sulfate-binding protein. One of these cavities contains a crystallographically observed water molecule while the other does not. Thermodynamic integration and perturbation methods were used to calculate free energies of hydration for each of the cavities from molecular dynamics simulations of two separate events: the removal of a water molecule from pure water, and the introduction of a water molecule into each protein cavity. From the simulations for the pure water system, the excess chemical potential of water was computed to be -6.4 +/- 0.4 kcal/mol, in accord with experiment and with other recent theoretical calculations. For the protein cavity containing an experimentally observed water molecule, the free energy change on hydrating it with one water molecule was calculated as -10.0 +/- 1.3 kcal/mol, indicating the high probability that this cavity is occupied by a water molecule. By contrast, for the cavity in which no water molecules were experimentally observed, the free energy change on hydrating it with one water molecule was calculated as 0.2 +/- 1.5 kcal/mol, indicating its low occupancy by water. The agreement of these results with experiment suggests that thermodynamic simulation methods may become useful for the prediction and analysis of internal hydration in proteins.

Hydrogen Bonding↗
PubMed1990

Solvation of the active site of cytochrome P450-cam.

Energetically favorable water binding sites in the substrate pocket of cytochrome P450-cam have been predicted by a molecular mechanics method. Binding sites corresponding to all the experimentally observed water sites in this region of the enzyme were located. The calculations also indicate the presence of two further water binding sites. One of these is located in a hydrophobic region of the protein where a water molecule would not bind tightly to the substrate-free enzyme. However, in the substrate-bound enzyme, a water molecule in this region could donate a hydrogen bond of optimum geometry to the carbonyl oxygen atom of the camphor substrate and could therefore contribute to the correct positioning of the camphor substrate for 5-exo-hydroxylation. These calculations also suggest that a steric analogue of camphor, containing an alkyl group which could prevent a water molecule from binding in this region, might inhibit cytochrome P450-cam by forming a more stable enzyme-ligand complex than camphor itself.

Binding Sites↗
PubMed1989

Chemical evolution and the origin of life: bibliography supplement 1985.

This bibliography is the seventeenth annual supplement to the comprehensive bibliography on the same subject which was published in Space Life Sci. 2(1970), 225-295; 3(1972), 293-304; 4(1973), 309-329 and in Origins of Life 5(1974), 505-527; 6(1975), 285-300; 7(1976), 75-85; 8(1977), 59-66; 9(1978), 67-74; 10(1980), 69-87; 10(1980), 379-404; 11(1981), 273-288; 12(1982), 93-118; 13(1983), 61-80; 15(1984), 55-69; 17(1987), 171-184; 17(1987), 185-206.

Chemical Phenomena↗
PubMed1989

Identifying targets for bioreductive agents: using GRID to predict selective binding regions of proteins.

A computational procedure is described for investigating potential binding sites of a target macromolecule for their ability to bind both a reduced probe molecule and an oxidized probe molecule. The interaction energies are obtained using a molecular mechanics method and can be displayed as three-dimensional (3D) energy contours, indicating regions of the target molecule that may have favorable interactions with the probe molecule. Differences in the interaction energies of the oxidized and reduced probe with the target can also be plotted as contours, indicating regions that are selective for the reduced probe. These selectivity contours can be used to show whether the macromolecule is a potential target for bioreductive agents. The method has been applied to the chicken liver dihydrofolate reductase enzyme and has indicated new binding regions that may be suitable binding sites for bioreductive agents.

Binding Sites↗
PubMed1989

New hydrogen-bond potentials for use in determining energetically favorable binding sites on molecules of known structure.

An empirical energy function designed to calculate the interaction energy of a chemical probe group, such as a carbonyl oxygen or an amine nitrogen atom, with a target molecule has been developed. This function is used to determine the sites where ligands, such as drugs, may bind to a chosen target molecule which may be a protein, a nucleic acid, a polysaccharide, or a small organic molecule. The energy function is composed of a Lennard-Jones, an electrostatic and a hydrogen-bonding term. The latter is dependent on the length and orientation of the hydrogen bond and also on the chemical nature of the hydrogen-bonding atoms. These terms have been formulated by fitting to experimental observations of hydrogen bonds in crystal structures. In the calculations, thermal motion of the hydrogen-bonding hydrogen atoms and lone-pair electrons may be taken into account. For example, in a alcoholic hydroxyl group, the hydrogen may rotate around the C-O bond at the observed tetrahedral angle. In a histidine residue, a hydrogen atom may be bonded to either of the two imidazole nitrogens and movement of this hydrogen will cause a redistribution of charge which is dependent on the nature of the probe group and the surrounding environment. The shape of some of the energy functions is demonstrated on molecules of pharmacological interest.

Binding Sites↗
PubMed1989

The role of hydrogen-bonds in drug binding.

Hydrogen-bonds play a crucial role in determining the specificity of ligand binding. Their important contribution is explicitly incorporated into a computational method, called GRID, which has been designed to detect energetically favourable ligand binding sites on a chosen target molecule of known structure. An empirical energy function consisting of a Lennard-Jones, an electrostatic and a hydrogen-bonding term is employed. The latter term is found to be necessary because spherically symmetric atom-centred forces alone may not adequately reproduce the geometry of two interacting molecules. The hydrogen-bonding term is dependent on the length and orientation of the hydrogen-bond. Its functional form also varies according to the chemical nature of both the hydrogen-bond donor and acceptor atoms, and has been modelled to fit experimental observations of crystal structures. The mobility of the hydrogen-bonding hydrogens is considered analytically in calculating the hydrogen-bond energy. The hydrogen-bonding energy functions will be described and their application will be demonstrated on molecules of pharmacological interest where hydrogen-bonds influence the binding of ligands.

Camphor 5-Monooxygenase↗
PubMedDate not supplied

On the protein-protein diffusional encounter complex.

When two proteins diffuse together to form a bound complex, an intermediate is formed at the end-point of diffusional association which is called the encounter complex. Its characteristics are important in determining association rates, yet its structure cannot be directly observed experimentally. Here, we address the problem of how to construct the ensemble of three-dimensional structures which constitute the protein-protein diffusional encounter complex using available experimental data describing the dependence of protein association rates on mutation and on solvent ionic strength and viscosity. The magnitude of the association rates is fitted well using a variety of definitions of encounter complexes in which the two proteins are located at up to about 17 A root-mean-squared distance from their relative arrangement in the bound complex. Analysis of the ionic strength dependence of bimolecular association rates shows that this is determined to a greater extent by the (protein charge) - (salt ion) separation distance than by the protein-protein charge separation distance. Consequently, ionic strength dependence of association rates provides little information about the geometry of the encounter complex. On the other hand, experimental data on electrostatic rate enhancement, mutation and viscosity dependence suggest a model of the encounter complex in which the two proteins form a subset of the contacts present in the bound complex and are significantly desolvated.

Amino Acid Substitution↗