PubMed Health⌕ Search

Biomedical subjects

D J Osguthorpe

Publications and source records attributed to D J Osguthorpe.

At least 19 recordsLinked to original sources
PubMed1992

Modeling of substrate and inhibitor binding to phospholipase A2.

Molecular graphics and molecular mechanics techniques have been used to study the mode of ligand binding and mechanism of action of the enzyme phospholipase A2. A substrate-enzyme complex was constructed based on the crystal structure of the apoenzyme. The complex was minimized to relieve initial strain, and the structural and energetic features of the resultant complex analyzed in detail, at the molecular and residue level. The minimized complex was then used as a basis for examining the action of the enzyme on modified substrates, binding of inhibitors to the enzyme, and possible reaction intermediate complexes. The model is compatible with the suggested mechanism of hydrolysis and with experimental data about stereoselectivity, efficiency of hydrolysis of modified substrates, and inhibitor potency. In conclusion, the model can be used as a tool in evaluating new ligands as possible substrates and in the rational design of inhibitors, for the therapeutic treatment of diseases such as rheumatoid arthritis, atherosclerosis, and asthma.

Amino Acid Sequence↗
PubMed1992

FOCUS: a program for analyzing molecular dynamics simulations, featuring digital signal-processing techniques.

FOCUS is a program for analyzing molecular dynamics simulations. It enables the researcher to monitor structural and energetic properties during the trajectory, and to calculate the corresponding statistical averages, correlation functions and Fourier transforms. In addition to these conventional methods, the program also utilizes novel methods based on digital signal-processing techniques to characterize the various motions. The characteristic frequencies in the system are revealed by the frequency distribution function g(v), which is calculated from the Fourier transform of the atomic coordinates. A filtering technique is employed to remove uninteresting motion (e.g., high-frequency bond stretching) while retaining and focusing on important motion (e.g., low-frequency conformational motion). The filtering technique enables fast display of slow events without getting a blurry or jittery picture due to the high-frequency motions. Another new way for analyzing the motion is by extracting "characteristic modes" and associated frequencies. This yields a pictorial description of the oscillatory motions in a manner analogous to normal mode analysis.

Acetamides↗
PubMed1992

BIOSITE: a program for the interactive comparison of aligned homologous protein sequences.

A program, BIOSITE, providing for the interactive visual comparison of aligned homologous amino-acid sequences is presented, including an example of its application. The program allows for two types of comparison sequence to be generated: an 'identity' sequence and a 'difference' sequence. These may be used on subsets of sequences and in further comparisons to identify candidate sites involved in a distinct functional property. The program should prove a useful tool for biologists engaged in understanding sequence--function relationships.

Amino Acid Sequence↗
PubMed1991

Relative-residue surface-accessibility patterns reveal myoglobin and catalase similarity.

A novel sliding-window search method using relative-residue surface-accessibility patterns identified extensive, but unsuspected, structural similarity over a 3-helix region in the C-terminus of the evolutionarily unrelated proteins sperm-whale myoglobin and beef liver catalase. This clear example of structural similarity between non-homologous proteins highlights the importance of relative-residue surface-accessibility patterns in understanding the local folded structure in proteins.

Amino Acid Sequence↗
PubMed1991

Conformational analysis of peptide surrogates. Reduced and retro-amide links in blocked alanine and in secondary structures.

We have investigated the conformational effects of modifying the amide link of peptides. We studied a reverse amide bond psi [NHCO], a reduced amide bond psi [CH2NH] and a retro-reduced bond psi [NHCH2] as surrogates for the amide link [CONH] in native peptides. A complete search of the conformational space available to residues with these modified links was carried out. The local minima and the rotational barriers were described and compared to the minima of the native residue. The results are compatible with the available observed structural data. These modified links have been incorporated in secondary structure units such as beta turns, alpha helices, and parallel and anti-parallel beta sheets. It was found that a reduced amide link can lead to stabilised beta turns, while the retro modification can be incorporated in stable beta sheets. A significant reduction in the stability of alpha helices is caused by the retro links, while a reduced amide link results in only a small destabilisation.

Alanine↗
PubMed1990

Extraction of the energetics of selected types of motion from molecular dynamics trajectories by filtering.

A novel method for analyzing molecular dynamics trajectories has been developed which enables the study of selected motions and the corresponding energetics. In particular, it is possible to filter out the high-frequency motions and focus on the structural and energetic features of low-frequency collective motions. The trajectories of the properties of interest are Fourier transformed to the frequency domain, a filtering function is applied, and then an inverse transformation back to the time domain yields the filtered trajectory. The method is demonstrated for harmonic fluctuations and conformational transitions of acetamide and N-acetylalanine N-methylamide, as models for peptides and proteins.

Alanine↗
PubMed1990

Melanin-concentrating hormone: a structural and conformational study based on synthesis, biological activity, high-field NMR, and molecular modeling techniques.

A series of Melanin-concentrating hormone (MCH) fragments have been synthesized and their biological activities compared with the parent peptide. The substructural units, 5-14 linear and 5-14 cyclic, have been used as models for MCH-- H-Asp1-Thr-Met-Arg-Cys-Met-Val-Gly-Arg HO-Val17-Glu-Trp-Cys-Pro-Arg-Tyr-Val in 1H-nmr conformational studies. Conformational features predicted by molecular dynamics analyses find support in the nmr experiments.

Amino Acid Sequence↗
PubMed1990

Accessible conformations of melanin-concentrating hormone: a molecular dynamics approach.

Molecular dynamics simulations have been used to search for the accessible conformations of the melanin-concentrating hormone (MCH). The studies have been performed on native MCH and two of its peptide fragments, a cyclic MCH(5-14) fragment and a linear MCH(5-14) fragment. An analysis of the molecular dynamics trajectories of the three peptides indicates that two regions of the peptide have characteristic conformational properties that may be important for the biological activity. One is a region around Gly8, which is conformationally mobile, and the other is around Pro13, which shows unusual rigidity. The molecular dynamics simulation results are discussed in terms of backbone structural features like beta turns, side-chain interactions, and orientations of the disulfide bridge. The results of this analysis are used to suggest new analogues that will modify the conformational features of the peptide and further define the conformational requirements for activity. Finally, the results are related to nmr studies of the peptide and reveal agreements between the experimental nuclear Overhauser effect constraints and some of the accessible conformations obtained from the simulation.

Amino Acid Sequence↗
PubMed1990

Modeling of agonist binding to the ligand-gated ion channel superfamily of receptors.

A generalized model is presented of agonist binding to ligand-gated ion channels (LGICs). Broad similarity in the structure of agonists suggests that the binding sites of LGICs may have evolved from a protobinding site. Aligned sequence data identified as a candidate for such a site a highly conserved 15 residue stretch of primary structure in the N-terminal extracellular region of all known LGIC subunits. We modeled this subregion, termed the cys-loop, as a rigid, amphiphilic beta-hairpin and propose that it may form a major determinant of a conserved structural binding cleft. In the model of the binding complex (1) an invariant aspartate residue at position 11 of the cys-loop is the anionic site interacting with the positively charged amine group of agonists, (2) a local dipole within the pi-electron system of agonists is favorably oriented in the electrostatic field of the invariant aspartate, (3) the epsilon ring-proton of a conserved aromatic residue at the turn of the cys-loop interacts orthogonally with the agonist pi-electron density at its electronegative center, and (4) selective recognition is partly a result of the type of amino acid residue at position 6 of the cys-loop. Additionally, formation of a hydrogen bond between the electronegative atom of the pi-electron system of agonist and a complementary group in the receptor may be important in the high-affinity binding of agonists.

Amino Acid Sequence↗
PubMed1990

The conformational preferences of gamma-lactam and its role in constraining peptide structure.

The conformational constraints imposed by gamma-lactams in peptides have been studied using valence force field energy calculations and flexible geometry maps. It has been found that while cyclisation restrains the psi of the lactam, non-bonded interactions contribute to the constraints on psi of the lactam. The gamma-lactam also affects the (psi, psi) of the residue after it in a peptide sequence. For an L-lactam, the ring geometry restricts psi to about -120 degrees, and psi has two minima, the lowest energy around -140 degrees and a higher minimum (5 kcal/mol higher) at 60 degrees, making an L-gamma-lactam more favourably accommodated in a near extended conformation than in position 2 of a type II' beta-turn. The energy of the psi approximately +60 degrees minimum can be lowered substantially until it is more favoured than the -140 degrees minimum by progressive substitution of bulkier groups on the amide N of the L-gamma-lactam. The (psi, psi) maps of the residue succeeding a gamma-lactam show subtle differences from those of standard N-methylated residues. The dependence of the constraints on the chirality of gamma-lactams and N-substituted gamma-lactams, in terms of the formation of secondary structures like beta-turns is discussed and the comparison of the theoretical conformations with experimental results is highlighted.

Computer Simulation↗
PubMed1990

Modelling of binding sites of the nicotinic acetylcholine receptor and their relation to models of the whole receptor.

Models for the acetylcholine (ACh)-binding site of the nicotinic acetylcholine receptor (nAChR) are proposed. These models have been developed by using the concept of the ligand-gated ion-channel (LGIC) superfamily of receptors that have evolved from a common ancestor. An initial component of the binding site was identified as a highly conserved 15-residue stretch of primary structure in the N-terminal extracellular region of all known LGIC subunits, based on aligned sequence data of LGICs. This subregion, termed the Cys-loop, was modelled as an amphiphilic beta-hairpin and we propose that it forms a major determinant of the binding cleft for agonists. This initial, partial binding-site model has been extended to include residues biochemically identified as spatially adjacent to the binding cleft. A recently developed technique for rapidly scanning the known protein structural database for 'non-homologous similarity' using just sequence information identified the known structure of the enzyme pyrophosphatase (PPase) as a candidate scaffold for the N-terminal domain of the nAChR. This similarity was investigated further using sequence alignments. A framework model of the full N-terminal domain in which the position of the Cys-loop and other binding-site determinants, as well as the main immunogenic region (MIR), have been mapped on to the PPase structure.

Amino Acid Sequence↗
PubMed1989

A novel beta-turn location in an LHRH antagonist: a combined conformational search and molecular dynamics study.

A 50 pico-second molecular dynamics simulation on a cyclic LHRH antagonist analogue Ac-D-Phe1-D-Phe2-D-Trp3-Ser4-Glu5-D-Arg6-Leu7-Lys8+ ++-Pro9-D-Ala10-NH2 (where the cyclisation is via an amide linkage between the Glu5 and Lys8 side chains), reveals some hitherto unseen conformational features. The LHRH analogue is found to adopt a near beta-sheet type of conformation with the reversal in the chain being brought about by a D-Trp3-Ser4-Glu5-D-Arg6 beta turn. The N- and C-terminal ends of the peptide come close together and interact through a network of hydrogen bonds. Additional hydrogen bonds expected of a sheet type of conformation stabilise the lowest energy minima. A conformational search of all possible cyclic structures of a model system c(Glu-D-Ala-Ala-Lys) which was used to determine the starting structure for the simulation studies of the cyclic LHRH antagonist analogue is also highlighted. The influence of the cyclic part on the conformation of this LHRH analogue is discussed.

Amino Acid Sequence↗
PubMed1989

Filtering molecular dynamics trajectories to reveal low-frequency collective motions: phospholipase A2.

A novel method for analysing molecular dynamics trajectories has been developed, which filters out high frequencies using digital signal processing techniques and facilitates focusing on the low-frequency collective motions of proteins. These motions involve low energy slow motions, which lead to important biological phenomena such as domain closure and allosteric effects in enzymes. The filtering method treats each of the atomic trajectories obtained from the molecular dynamics simulation as a "signal". The trajectories of each of the atoms in the system (or any subset of interest) are Fourier transformed to the frequency domain, a filtering function is applied and then an inverse transformation back to the time domain yields the filtered trajectory. The filtering method has been used to study the dynamics of the enzyme phospholipase A2. In the filtered trajectory, all the high frequency bond and valence angle vibrations were eliminated, leaving only low-frequency motion, mainly fluctuations in torsions and conformational transitions. Analysis of this trajectory revealed interesting motions of the protein, including concerted movements of helices, and changes in shape of the active site cavity. Unlike normal mode analysis, which has been used to study the motion of proteins, this method does not require converged minimizations or diagonalization of a matrix of second derivatives. In addition, anharmonicity, multiple minima and conformational transitions are treated explicitly. Thus, the filtering method avoids most of the approximations implicit in other investigations of the dynamic behaviour of large systems.

Binding Sites↗
PubMed1989

Sequence alignment of citrate synthase proteins using a multiple sequence alignment algorithm and multiple scoring matrices.

The alignment of Escherichia coli citrate synthase to pig heart citrate synthase and the multiple alignment of the known sequences of the citrate synthase family of enzymes have been performed using six different amino acid similarity scoring matrices and a large range of gap penalty ratios for insertions and deletions of amino acids. The alignment studies have been performed as the first step in a project aimed at homology modelling E. coli citrate synthase (a hexamer) from pig heart citrate synthase (a dimer) in a molecular modelling approach to the study of multi-subunit enzymes. The effects of several important variables in producing realistic alignments have been investigated. The difference between multiple alignment of the family of enzymes versus simple pairwise alignment of the pig heart and E. coli proteins was explored. The effects of initial separate multiple alignments of the most highly related or most homologous species of the family of enzymes upon a subsequent pairwise alignment between species was evaluated. The value of 'fingerprinting' certain residues to bias the alignment in favour of matching those residues, as well as the worth of the computerized approach compared to an intuitive alignment technique, were assessed.

Algorithms↗
PubMed1988

Structure and energetics of ligand binding to proteins: Escherichia coli dihydrofolate reductase-trimethoprim, a drug-receptor system.

A study of the binding of the antibacterial agent trimethoprim to Escherichia coli dihydrofolate reductase was carried out using energy minimization techniques with both a full, all-atom valence force field and a united atom force field. Convergence criteria ensured that no significant structural or energetic changes would occur with further minimization. Root-mean-square (RMS) deviations of both minimized structures with the experimental structure were calculated for selected regions of the protein. In the active site, the all-atom minimized structure fit the experimental structure much better than did the united atom structure. To ascertain what constitutes a good fit, the RMS deviations between crystal structures of the same enzyme either from different species or in different crystal environments were compared. The differences between the active site of the all-atom minimized structure and the experimental structure are similar to differences observed between crystal structures of the same protein. Finally, the energetics of ligand binding were analyzed for the all-atom minimized coordinates. Strain energy induced in the ligand, the corresponding entropy loss due to shifts in harmonic frequencies, and the role of specific residues in ligand binding were examined. Water molecules, even those not in direct contact with the ligand, were found to have significant interaction energies with the ligand. Thus, the inclusion of at least one shell of waters may be vital for accurate simulations of enzyme complexes.

Computer Simulation↗