Search PubMed⌕ Search

Biomedical subjects

R E Hubbard

Publications and source records attributed to R E Hubbard.

At least 37 records · Page 2Linked to original sources

Conformational change in the activation of lipase: an analysis in terms of low-frequency normal modes.

The interfacial activation of Rhizomucor miehei lipase (RmL) involves the motion of an alpha-helical region (residues 82-96) which acts as a "lid" over the active site of the enzyme, undergoing a displacement from a "closed" to an "open" conformation upon binding of substrate. Normal mode analyses performed in both low and high dielectric media reveal that low-frequency vibrational modes contribute significantly to the conformational transition between the closed and open conformations. In these modes, the lid displacement is coupled to local motions of active site loops as well as global breathing motions. Atomic fluctuations of the first hinge of the lid (residues 83-84) are substantially larger in the low dielectric medium than in the high dielectric medium. Our results also suggest that electrostatic interactions of Arg86 play an important role in terms of both the intrinsic stability of the lid and its displacement, through enhancement of hinge mobility in a high dielectric medium. Additional calculations demonstrate that the observed patterns of atomic fluctuations are an intrinsic feature of the protein structure and not dependent on the nature of specific energy minima.

Binding Sites↗

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↗

Towards structural models of molecular recognition in olfactory receptors.

The G protein coupled receptors (GPCR) are an important class of proteins that act as signal transducers through the cytoplasmic membrane. Understanding the structure and activation mechanism of these proteins is crucial for understanding many different aspects of cellular signalling. The olfactory receptors correspond to the largest family of GPCRs. Very little is known about how the structures of the receptors govern the specificity of interaction which enables identification of particular odorant molecules. In this paper, we review recent developments in two areas of molecular modelling: methods for modelling the configuration of trans-membrane helices and methods for automatic docking of ligands into receptor structures. We then show how a subset of these methods can be combined to construct a model of a rat odorant receptor interacting with lyral for which experimental data are available. This modelling can help us make progress towards elucidating the specificity of interactions between receptors and odorant molecules.

Aldehydes↗

Molecular basis of agonism and antagonism in the oestrogen receptor.

Oestrogens are involved in the growth, development and homeostasis of a number of tissues. The physiological effects of these steroids are mediated by a ligand-inducible nuclear transcription factor, the oestrogen receptor (ER). Hormone binding to the ligand-binding domain (LBD) of the ER initiates a series of molecular events culminating in the activation or repression of target genes. Transcriptional regulation arises from the direct interaction of the ER with components of the cellular transcription machinery. Here we report the crystal structures of the LBD of ER in complex with the endogenous oestrogen, 17beta-oestradiol, and the selective antagonist raloxifene, at resolutions of 3.1 and 2.6 A, respectively. The structures provide a molecular basis for the distinctive pharmacophore of the ER and its catholic binding properties. Agonist and antagonist bind at the same site within the core of the LBD but demonstrate different binding modes. In addition, each class of ligand induces a distinct conformation in the transactivation domain of the LBD, providing structural evidence of the mechanism of antagonism.

Binding Sites↗

Domain motions in dihydrofolate reductase: a molecular dynamics study.

Molecular dynamics simulations have been carried out on the enzyme dihydrofolate reductase from Lactobacillus casei complexed with methotrexate, NADPH and 264 crystallographic water molecules. Analysis of correlations in atomic fluctuations reveal the presence of highly correlated motion (correlation coefficient > 0.6) in the region between residues 30 to 35 and 85 to 90 leading to the identification of two domains, an "adenosine-binding domain" and a "large domain", which rotate by 3 to 4 degrees with respect to each other. The strongest correlation (> 0.6) within the large domain involves a coupling between the motions of the "teen-loop", and the spatially contiguous loops linking beta 6-beta 7 and beta 7-beta 8. Moreover, there is a significant correlation (approximately 0.5) between the adenosine fragment of NADPH and the pteridine and p-aminobenzoyl fragments of methotrexate, which are separated by approximately 17 A, and is lost on removal of "rigid-body" motion from the original trajectory. This provides support for the idea that the relative motion of the two domains is a means by which the occupation of the binding site for the adenosine end of the coenzyme can affect methotrexate binding and vice versa. Quasiharmonic vibrational analysis of the trajectory reveals that the overall dynamics of the system are governed by domain motions whose contributions are dominant at low frequencies. In addition, different low-frequency modes are responsible for separately coupling the adenosine-binding site and parts of methotrexate.

Binding Sites↗

Can drugs be designed?

Many examples are now emerging of the successful use of rational, structure-based methods in drug discovery. Of particular note is the development of imaginative NMR-based methods for rapid routes to ligand design. Our understanding of the chemistry underlying protein-ligand interactions, however, remains relatively poor and a major limitation in our ability to truly design drugs.

Drug Design↗

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↗

The x-ray crystal structure of phosphomannose isomerase from Candida albicans at 1.7 angstrom resolution.

Phosphomannose isomerase (PMI) catalyses the reversible isomerization of fructose-6-phosphate (F6P) and mannose-6-phosphate (M6P). Absence of PMI activity in yeasts causes cell lysis and thus the enzyme is a potential target for inhibition and may be a route to antifungal drugs. The 1.7 A crystal structure of PMI from Candida albicans shows that the enzyme has three distinct domains. The active site lies in the central domain, contains a single essential zinc atom, and forms a deep, open cavity of suitable dimensions to contain M6P or F6P The central domain is flanked by a helical domain on one side and a jelly-roll like domain on the other.

Binding Sites↗

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↗

Characterising the geometric diversity of functional groups in chemical databases.

We have developed a program, HookSpace, which provides a simplistic approach to assessing the diversity of molecular databases. The spatial relationship between pairs of intramolecular functional groups can be analysed in a variety of ways to provide both qualitative and quantitative measures of diversity. Results are described and contrasted for two commercially available databases and a combinatorial library of benzodiazepam derivatives. HookSpace highlights the main differences in molecular content of these data sets.

Crystallography, X-Ray↗

EXTRACT: a program to extract three-dimensional coordinates from stereo diagrams of proteins.

The program EXTRACT has been developed to extract accurate three-dimensional coordinates from published stereo alpha-carbon diagrams of protein structures. The approach is based on the display of scanned images of the left and right eye views of the diagram on a stereo-equipped workstation, allowing construction of a molecular model using the diagram as a guide. A number of structural checks assess the building, including probability maps derived for alpha-carbon geometry in protein structures. The procedure has also been extended to produce less accurate models from mono images.

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↗

Investigating the high affinity and low sequence specificity of calmodulin binding to its targets.

Calmodulin (CaM) is a calcium binding protein that regulates a wide range of enzymes. Recently the structures of a number of complexes between CaM and synthetic target peptides have been determined. The peptides correspond to the CaM-binding domain of skeletal and smooth muscle myosin light-chain kinase (MLCK) and calmodulin-dependent protein kinase II alpha. Comparison of the peptide-free and peptide-bound structures reveals that CaM undergoes a large conformational change when forming a complex, resulting in the formation of a binding surface that provides for an optimal interaction with its target. In this work, the available co-ordinates of the NMR solution structure of CaM-skeletal MLCK peptide are used as a basis upon which several molecular models of binding are built. The detailed features of the protein's peptide binding surface are revealed through two-dimensional topographical projections. Negatively charged margins at the binding surface extremities interact strongly with basic peptide residues separated by nine or ten positions. The binding surface core is hydrophobic and displays a groove with four deep pockets, which can accommodate bulky peptide residues at relative positions 4 and 8 (pocket A), 11 (pocket B), 13 (pocket C), 14 and 17 (pocket D). Therefore, both electrostatic and van der Waals' features contribute to the high affinity binding. A search for alternative peptide placements in the binding tunnel reveals the dominant role of specific electrostatic interactions in the binding energy. Apolar interactions are more permissive, such that the hydrophobic side-chains that line the binding tunnel adapt in order to maintain favourable van der Waals' contacts. The model suggests that the structure can accommodate large peptide translations (up to 5 A) and a reversed peptide binding mode, with a little loss in binding interaction energy. These calculations are compared with available experimental data, providing a structural rationale for the low sequence specificity of the CaM target recognition.

Amino Acid Sequence↗

Crystallization and preliminary X-ray analysis of Candida albicans phosphomannose isomerase.

Crystals of recombinant phosphomannose isomerase from Candida albicans have been obtained in a form suitable for X-ray diffraction analysis. The enzyme plays a key role in the biosynthesis of the mannan component of the fungal cell wall. It crystallizes in monoclinic space group C2, with cell dimensions a = 124.9 A, b = 52.9 A, c = 85.9 A and beta = 127.4 degrees. The crystals diffract to Bragg spacings beyond 1.7 A, native data have been collected to 2.4 A and a search for heavy-metal derivatives is in progress. The asymmetric unit contains one molecule of the enzyme (M(r) approximately 49,000) with a Vm of 2.3 A3/Da.

Candida albicans↗

Analysis of C alpha geometry in protein structures.

The polypeptide of a protein molecule can be considered as a chain of C alpha atoms linked by pseudobonds between the C alpha atoms of successive amino acid residues. This paper presents an analysis of the angle and dihedral angles made by these pseudobonds in protein structures determined at high resolution by X-ray crystallography. This analysis reveals a strong correlation between C alpha geometry and the protein fold. The regular features of protein secondary structure such as alpha-helix and beta-sheet are very clearly defined. In addition, it is possible to identify with some confidence the discrete populations of particular conformations of beta-turn. Comparison with the traditional Ramachandran type of plot demonstrates that an analysis of protein structure on the basis of C alpha geometry provides a richer description of protein conformation. In addition, the characteristics of this geometry could be a useful guide in model building of protein structure.

Crystallography, X-Ray↗