PubMed HealthSearch

Biomedical subjects

P Herzyk

Publications and source records attributed to P Herzyk.

13 recordsLinked to original sources

Combined biophysical and biochemical information confirms arrangement of transmembrane helices visible from the three-dimensional map of frog rhodopsin.

The electron density projection map of frog rhodopsin at 6 A resolution had been until recently the most direct evidence for the three-dimensional structure of a transmembrane domain of any G-protein-coupled receptor. Only three out of seven transmembrane helices are clearly defined, whilst the other four are hidden in a patch of unresolved electron density. A model of the seven-helix bundle has been created by generating positions and orientations for the four unresolved helices through performing a conformational search directed by structural restraints derived from other experimental data. These four helices are significantly tilted with respect to the membrane normal, and the cytosolic end of helix C is inserted between helices D and E. These calculations produce positions and orientations for these additional helices that are consistent with the recently published low-resolution three-dimensional map, and provide a template for more detailed modelling of rhodopsin structure and function.

Amino Acid Sequence

Using experimental information to produce a model of the transmembrane domain of the ion channel phospholamban.

Molecular models of the transmembrane domain of the phospholamban pentamer have been generated by a computational method that uses the experimentally measured effects of systematic single-site mutations as a guiding force in the modeling procedure. This method makes the assumptions that 1) the phospholamban transmembrane domain is a parallel five-helix bundle, and 2) nondisruptive mutation positions are lipid exposed, whereas 3) disruptive or partially disruptive mutations are not. Our procedure requires substantially less computer time than systematic search methods, allowing rapid assessment of the effects of different experimental results on the helix arrangement. The effectiveness of the approach is investigated in test calculations on two helix-dimer systems of known structure. Two independently derived sets of mutagenesis data were used to define the restraints for generating models of phospholamban. Both resulting models are left-handed, highly symmetrical pentamers. Although the overall bundle geometry is very similar in the two models, the orientation of individual helices differs by approximately 50 degrees, resulting in different sets of residues facing the pore. This demonstrates how differences in restraints can have an effect on the model structures generated, and how the violation of these restraints can identify inconsistent experimental data.

Amino Acid Sequence

Modeling of the three-dimensional structure of the human melanocortin 1 receptor, using an automated method and docking of a rigid cyclic melanocyte-stimulating hormone core peptide.

A model is presented of the melanocortin 1 receptor (MC1R), constructed by use of an unbiased, objective method. The model is created directly from data derived from multiple sequence analysis, a low-resolution EM-projection map of rhodopsin, and the approximate membrane thickness. The model agrees well with available data concerning natural mutations of MC1Rs occurring in different species. A model is also presented of the most rigid ligand for this receptor, the cyclic pentapeptide cHFRWG, shown docked in the receptor model. The receptor-ligand complex model agrees well with available experimental data. The ligand is located between transmembrane region 1 (TM1), TM2, TM3, TM6, and TM7 of the receptor. Multiple interactions occur between ligand and receptor, including interactions with Leu-48 (TM1), Ser-52 (TM1), Glu-55 (TM1), Asn-91 (TM2), Glu-94 (TM2), Thr-95 (TM2) Ile-98 (TM2), Asp-121 (TM3), Thr-124 (TM3), Phe-257 (TM6), Phe-283 (TM7), Asn-290 (TM7), and Asp-294 (TM7) of the receptor.

Amino Acid Sequence

Automated modelling of the transmembrane region of G-protein coupled receptor by Swiss-model.

Molecular modelling of the transmembrane helices of G-protein coupled receptors is an increasingly used method to identify the possible three-dimensional environment of key residues. Thereby site-directed mutagenesis experiments, aimed at the understanding of the receptor-ligand interactions, can be designed in a rational way. The modelling methods are however not generally available to experimentalists, and often require expensive software and hardware. To overcome these limitations, we have constructed a World Wide Web server for the automated protein modelling of user-defined transmembrane helices. The service is freely available at this address: http:/(/)expasy.hcuge.ch/swissmod/SWISS-MODEL.++ +html.

Animals

Automated method for modeling seven-helix transmembrane receptors from experimental data.

A rule-based automated method is presented for modeling the structures of the seven transmembrane helices of G-protein-coupled receptors. The structures are generated by using a simulated annealing Monte Carlo procedure that positions and orients rigid helices to satisfy structural restraints. The restraints are derived from analysis of experimental information from biophysical studies on native and mutant proteins, from analysis of the sequences of related proteins, and from theoretical considerations of protein structure. Calculations are presented for two systems. The method was validated through calculations using appropriate experimental information for bacteriorhodopsin, which produced a model structure with a root mean square (rms) deviation of 1.87 A from the structure determined by electron microscopy. Calculations are also presented using experimental and theoretical information available for bovine rhodopsin to assign the helices to a projection density map and to produce a model of bovine rhodopsin that can be used as a template for modeling other G-protein-coupled receptors.

Amino Acid Sequence

A reduced representation of proteins for use in restraint satisfaction calculations.

A reduced representation of proteins has been developed for use in restraint satisfaction calculations with dynamic simulated annealing. Each amino acid residue is represented by up to four spherical virtual atoms. The virtual bonds and excluded volume of these atoms has been parameterized by analysis of 83 protein structures determined at high resolution by X-ray crystallography. The use of the new representation in NOE distance restraint satisfaction has been compared with the standard all-atom representation for the determination of the structures of crambin, echistatin, and protein G. Using the reduced representation, there is a 30-fold decrease in the computer time needed for generating a single structure, and up to a 20-fold decrease in the time taken to produce an acceptable structure compared to using the all-atom representation. The root mean square deviation between the mean structure obtained with all-atom and reduced representations is between 1.5 and 1.7 A for C alpha atoms. The new representation is adequate for describing the "low-resolution" features of protein structure such as the general fold and the positions of secondary structure elements. It can also provide an initial structure for more detailed refinement with the full all-atom representation.

Computer Simulation

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

The solution structure of echistatin: evidence for disulphide bond rearrangement in homologous snake toxins.

The solution structure of the fibrinogen antagonist, echistatin, has been determined by a combination of NMR and simulated annealing methods. While the structure of the disulphide-linked core is well-defined by the NMR data, the N- and C-termini and the loop bearing the RGD sequence (which is responsible for the fibrinogen antagonist properties) are poorly defined. The pattern of disulphide bridges, which could not be determined by classical methods, was predicted by a statistical analysis of the simulated annealing structures. This pattern is distinct from that for the homologous protein kistrin, leading to the novel suggestion that homologous proteins possess non-conserved patterns of disulphide bridges.

Amino Acid Sequence

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

Molecular-mechanics modelling of drug-DNA structures; the effects of differing dielectric treatment on helix parameters and comparison with a fully solvated structural model.

This study analyses the influence that the nature of the dielectric constant has on the final structures obtained from in vacuo molecular mechanics calculations on a drug-DNA complex and compares these structures with the energy minimised complex including explicit solvent molecules. Minimisations have been performed on a proflavine-decanucleotide structure, where the drug was intercalated at the d(CpG) site of the d(GpApTpApCpGpApTpApC) decamer duplex, using two expressions for the dielectric constant: a distance-independent, epsilon ij = EPS, and a distance-dependent, epsilon ij = EPS*Rij, form and for values of EPS from 1 to 8. Significantly different structures are obtained for the distance-independent and the distance-dependent expressions of the dielectric constant. The use of a distance-independent dielectric constant leads to distorted structures, which are very sensitive to slight changes in the value of EPS. The use of a distance-dependent dielectric constant leads to less distorted and more stable structures. The effects on helical parameters are analysed in detail. The structures obtained for different values of EPS (within the distance-dependent formalism) seem to converge for values of EPS equal to 4 or greater. Based on these results a distance-dependent form of the dielectric with an EPS value of 4 is recommended in order to produce reliable refined nucleic acid structures by molecular mechanics. These conclusions have been supported by molecular-mechanics minimisation of the same structure with the inclusion of explicit water molecules and counter-ions.

Base Composition

A molecular model for proflavine-DNA intercalation.

A molecular model has been derived for the intercalation of proflavine into the CpG site of the decamer duplex of d(GATACGATAC). The starting geometry of the intercalation site was taken from previous crystallographic studies on the d(CpG)-proflavine complex, and molecular mechanics used to obtain a stereochemically acceptable structure. This has widened grooves compared to standard A- or B- double helices, as well as distinct conformational, roll, twist and tilt features.

Acridines

Conformational properties of 3'-azido-3'deoxy-thymidine (AZT), an inhibitor of HIV reverse transcriptase.

The low-energy conformations of 3'-azido-3'-deoxy-thymidine, (AZT), an inhibitor of retroviral reverse transcriptase, have been studied by molecular mechanics techniques. A force-field has been developed for the azido group by quantum-mechanical methods, and used in the analysis. The global low-energy structure of AZT has C3'-endo sugar pucker, an anti glycosidic angle, and a g+ C4'-C5' conformation. It is concluded that the AZT molecule has conformational properties that are very similar to those of standard deoxypyrimidines.

Antiviral Agents