Search PubMed⌕ Search

Biomedical subjects

J M Goodfellow

Publications and source records attributed to J M Goodfellow.

At least 37 records · Page 2Linked to original sources
PubMed1993

The pore dimensions of gramicidin A.

The ion channel forming peptide gramicidin A adopts a number of distinct conformations in different environments. We have developed a new method to analyze and display the pore dimensions of ion channels. The procedure is applied to two x-ray crystal structures of gramicidin that adopt distinct antiparallel double helical dimer conformations and a nuclear magnetic resonance (NMR) structure for the beta6.3 NH2-terminal to NH2-terminal dimer. The results are discussed with reference to ion conductance properties and dependence of pore dimensions on the environment.

Amino Acid Sequence↗
PubMed1992

Conformation and dynamics of drug-DNA intercalation.

Molecular dynamics simulations have been undertaken for a B-form dodecanucleotide duplex in solution with and without an intercalated proflavine molecule between the central C.G base pairs. The introduction of this simple intercalator affects both the conformational features and dynamic properties of the oligonucleotide double helix. Changes are seen in the rms atomic fluctuations and anisotropy of phosphate, sugar and base atoms. The backbone conformation is slightly changed on average and more sugars adopt the C3' endo conformation in the simulation of the complex compared with the simulation of the oligonucleotide alone. Both major and minor grooves becomes wider on average with the addition of the intercalating drug. Flanking A.T base pairs on both sides of the intercalation site have undergone an increase in flexibility, with the base pairs, especially at the 5' side, having the N1...N3 hydrogen bonds being broken.

Base Composition↗
PubMed1992

Molecular dynamics of alkylated DNA.

The effect of methylation of the O4 atom of thymine in two oligonucleotide sequences is investigated by molecular dynamics simulations. Three types of environments are considered including: (i) in vacuo calculation, with a distance-dependent dielectric function and unhydrated counter-ions; (ii) in vacuo calculation, with a distance-dependent dielectric constant and hydrated counter-ions; and (iii) with a 9 A thick explicit water layer and counter-ions. In all environments, the oligonucleotide sequence containing the chemically modified thymine paired with guanine is more stable than the oligonucleotide sequence in which the modified thymine is paired with adenine. The methyl group attached to the O4 atom of thymine is found in a syn configuration with respect to the N3 atom. The best fit between the experimental NMR results and the molecular dynamics simulations is obtained using the environment with hydrated counter-ions.

Alkylation↗
PubMed1992

Hydration of amino acid side chains: dependence on secondary structure.

Energy calculations have been used to study the hydration sites around the polar groups of serine, threonine and tyrosine side chains. These hydration sites depend not only on the hybridization of the polar group but also on the local secondary structure, the chi 1 side chain torsion angle and the position of the hydroxyl hydrogen atom. For tyrosine side chains, two solvent sites are found approximately in the plane of the ring. Even for serine and threonine side chains only two minimum energy sites are found in general of which one is in an expected position within hydrogen bonding of the hydroxyl hydrogen atom (unless this is blocked from interaction with solvent molecules by, for example, Oi-4 or Oi-3. The position of the second of these sites depends not only on the position of the hydroxyl oxygen but also on neighbouring main chain atoms to which it can also hydrogen bond. There is good agreement with the solvent distributions obtained from crystallographic data.

Amino Acids↗
PubMed1991

Analysis of protein main-chain solvation as a function of secondary structure.

We have analysed the hydration of main-chain carbonyl and amide groups in 24 high-resolution well-refined protein structures as a function of the secondary structure in which these polar groups occur. We find that main-chain atoms in beta-sheets are as hydrated as those in alpha-helices, with most interactions involving "free" amide and carbonyl groups that do not participate in secondary structure hydrogen bonds. The distributions of water molecules around these non-bonded carbonyl groups reflect specific steric interactions due to the local secondary structure. Approximately 20% and 4%, respectively of bonded carbonyl and amide groups interact with solvent. These include interactions with carbonyl groups on the exposed faces of alpha-helices that have been correlated previously with bending of the helix. Water molecules interacting with alpha-helices occur mainly at the amino and carbonyl termini of the helices, in which case the solvent sites maintain the hydrogen bonding by bridging between residues i and i-3 or i-4 at the amino terminus and between i and i+3 or i+4 at the carbonyl terminus. We also see a number of solvent-mediated Ncap and Ccap interactions. The water molecules interacting with beta-sheets occur mainly at the edges, in which case they extend the sheet structure, or at the ends of strands, in which case they extend the beta-ladder. In summary, the solvent networks appear to extend the hydrogen-bonding structure of the secondary structures. In beta-turns, which usually occur at the surface of a protein, exposed amide and carbonyl groups are often hydrated, especially close to glycine residues. Occasionally water molecules form a bridge between residues i and i+3 in the turn and this may provide extra stabilization.

Amides↗
PubMed1991

Molecular dynamics simulations of dinucleoside and dinucleoside-drug crystal hydrates.

Molecular dynamics simulations have been performed on the dinucleoside monophosphates rGpC and dCpG, the latter in its intercalation complex with the acridine drug proflavine. The simulations were performed on the crystal structures, with crystallographically-located solvent molecules. It was found that satisfactory results were best obtained with restraints placed on the movements of the water molecules. Motions of individual atoms have been examined in terms of rms fluctuations and anisotropy and correlation functions. Relative motions of groups (phosphates, sugars, bases and proflavine molecules) have been analysed.

Computer Simulation↗
PubMed1991

Modelling of solvent positions around polar groups in proteins.

Previous analysis of the distribution of experimental solvent molecule positions around amino acid side chains showed that distinct clustering occurred close to polar or charged atoms in proteins. We have used those data to predict likely solvent positions around proteins not used in our initial analysis. We envisage that this algorithm, AQUARIUS, will be useful for finding solvent positions in electron density maps generated by protein crystallography and as useful starting positions for solvent molecules in computer simulation studies of macromolecules.

Algorithms↗
PubMed1990

Sequence dependent hydration of DNA.

The transitions between the different helical conformations of DNA depend on the base sequence and the ambient conditions such as humidity and counter-ion concentration. In this study energy minimization techniques have been used to locate water molecule sites around nucleotides especially those which form hydrogen bonds between two or more nucleotide atoms and thus form solvent mediated bridges. We have studied several sequences and find that those which are known not to exist in the low hydration 'A' form have very similar number of bridging sites in both 'A' and 'B' conformations. Those sequences which are found in the 'A' conformation have considerably more bridging sites in this low hydration form than in the 'B' conformation. Sequence related solvent effects for a given conformation have also been analysed.

Base Sequence↗
PubMed1990

Free energy changes associated with amino acid substitution in proteins.

The estimation of free energy differences from computer simulation of macromolecular systems is important for rational strategies for drug design and for protein engineering. As an example of one mutation, we have studied the free energy change resulting from the conversion of a polar group (OH) to an apolar group (CH3) in aqueous solution. We have estimated the effect of various local environments on the magnitude of the free energy difference and find that significant environmental effects are found. We have also studied the reliability of the results in detail.

Amino Acid Sequence↗
PubMed1990

Influence of secondary structure on the hydration of serine, threonine and tyrosine residues in proteins.

Previous analysis of experimental data on the solvation of high resolution protein structures has shown that preferred interaction sites for water molecules exist around most amino acid side chains. We have extended this analysis to look in more detail at the distributions around serine, threonine and tyrosine. We find that for serine and threonine side chains the preferred interaction sites of solvent molecules with the hydroxyl group depends on secondary structure and the chi 1 torsion angle of the side chain. For tyrosine side chains the hydroxyl group is too far from the main chain to reflect secondary structure influences. Specific patterns of hydration are observed in which water molecules 'bridge' between the hydroxyl side-chain atom and another main chain or side-chain atom.

Amino Acids↗
PubMed1989

Solvent interactions stabilising nucleic acid conformers.

The transition of oligonucleotides from the B to the A conformation has been studied by the use of simple geometric calculations aimed at finding possible hydration sites which could stabilize these conformations. The method involves the classification of equally spaced grid points, surrounding the oligonucleotide, into groups depending on whether a water molecule, so placed, could form single, multiple or nil contacts to polar oligonucleotide atoms. The occurrence of the multiple, and therefore bridging, sites is more extensive for the 'A' than the 'B' conformation. Thus, more general evidence is presented in support of the economy of hydration hypothesis in which phosphate groups, in the low humidity 'A' form, have been seen to be bridged by single water molecules. Similar calculations for the 'Z' DNA conformation show a different preference for multiple bridging sites.

DNA↗
PubMed1989

Solvent bridging sites in A- and B-DNA helices.

Nucleotide hydration is important for the understanding of the stability of and the transitions between the different helical conformations of DNA. We have used energy minimization and geometric criteria in order to look for possible sites for solvent which can bridge more than one polar or charged atomic group on a nucleotide. Such bridging sites between phosphate groups have been seen experimentally and used to explain the A to B transition. We show that these phosphate bridging sites occur at energy minima around A-DNA but do not occur around B-DNA. We also find that there are further low energy bridging sites which depend on sequence and which enable the more economical hydration of the A form.

Base Composition↗
PubMed1988

Distributions of water around amino acid residues in proteins.

The atomic co-ordinates from 16 high-resolution (less than or equal to 1.7 A = 0.1 nm), non-homologous proteins have been used to study the distributions of water molecule sites around the 20 different amino acid residues. The proportion of residues whose main-chain atoms are in contact with water molecules was fairly constant (between 40% and 60%), irrespective of the nature of the side-chain. However, the proportion of residues whose side-chain atoms were in contact with water molecules showed a clear (inverse) correlation with the hydrophobicity of the residue, being as low as 14% for leucine and isoleucine but greater than 80% for asparagine and arginine. Despite the problems in determining accurate water molecule sites from X-ray diffraction data and the complexity of the protein surface, distinct non-random distributions of water molecules were found. These hydration patterns are consistent with the expected stereochemistry of the potential hydrogen-bonding sites on the polar side-chains. The water molecules around apolar side-chains lie predominantly at van der Waals' contact distances, but most of these have a primary, shorter contact with a neighbouring polar atom. Further analysis of these distributions, combined with energy minimization techniques, should lead to improved modelling of protein structures, including their primary shells of hydration.

Amino Acid Sequence↗
PubMed1987

Ammonium ion representation in Monte Carlo simulations of biomolecular solutions.

Monte Carlo computer simulation techniques may be used to predict structural properties of solvent networks in helical fragments of nucleic acids, provided that suitable potential functions are available to describe the interactions between nucleic acid atoms, water and counterions. Previous studies have shown that simple non-bonded and point charge parameters are adequate for mononuclear ions such as sodium and calcium. In this study a model interaction potential for the polynuclear ammonium ion is evaluated. The parameters used take account of the distribution of charge over the constituent atoms in the ion. Simulations are carried out on the ammonium salt of a small nucleic acid crystal hydrate and a comparison is made between the predicted and experimental results. It is shown that the simulated structure is in reasonable agreement with experiment. It is therefore feasible to use this potential in studies of ammonium-containing bimolecular systems.

Computer Simulation↗
PubMed1987

Monte Carlo simulations of nucleotide crystal hydrates and their counter-ions.

A knowledge of structural and energetic aspects of water- and ion-nucleic acid interactions is essential for the understanding of the role of solvent and counterions in stabilising the various helical forms of nucleic acids. In this study, Monte Carlo computer simulation techniques have been used to predict structural properties of solvent networks in small nucleic acid crystal hydrates containing the ions sodium, ammonium and calcium. Appropriate parameters to describe the interaction potentials of the ions are added to those previously developed for water and nucleic acid atoms. A comparison is made between the predicted and experimental results and it is concluded that the potential functions used lead to simulated solvent structure in reasonable agreement with experimental data, at least in the cases of sodium and calcium. It is now feasible to use these functions in studies of hydration of larger helical fragments of nucleic acids of more direct biological interest.

Calcium↗
PubMed1985

Computer simulation of aqueous biomolecular systems.

Computer simulation techniques are increasingly being used to predict structural and thermodynamic properties of large heterogeneous macromolecule and solvent assemblies. We discuss, with examples from our own studies, some problems we and others have experienced in using these techniques, which were originally devised for simple liquids. In particular, we consider the problems which arise from the large size and heterogeneity of macromolecule water systems, comparisons with experimental data and equilibrium and sampling procedures.

Computer Simulation↗
PubMed1984

Solvent interactions in nucleic acid crystal hydrates.

A detailed knowledge of structural and energetic aspects of water-nucleic acid interactions is essential for understanding the role of solvent in stabilizing the various helical forms of nucleic acids. In this study, computer simulation techniques have been used to predict structural properties of solvent networks in small nucleic acid crystal hydrates. A detailed comparison of predicted and experimental results on the structure of the solvent networks is presented and includes an analysis of both the local environment and hydrogen bond pattern of each water molecule. A correlation between the environment of each unique water molecule and its energetic properties (such a dipole moment and binding energy) is seen. As in the previous studies on small amino acid hydrate crystals, non-pair additive (cooperative) effects are found to be non-negligible. It is concluded that the potential functions used in this initial study lead to simulated solvent networks in reasonable agreement with experimental data. Thus, it is now feasible to use them in studies of hydration of larger helical fragments of nucleic acids of more direct biological interest.

Crystallization↗