Biomedical subjects

Richard J Law

Publications and source records attributed to Richard J Law.

7 recordsLinked to original sources

PubMed2005

Channel opening motion of alpha7 nicotinic acetylcholine receptor as suggested by normal mode analysis.

The gating motion of the human nicotinic acetylcholine receptor (nAChR) alpha7 was investigated with normal mode analysis (NMA) of two homology models. The first model, referred to as model I, was built from both the Lymnaea stagnalis acetylcholine binding protein (AChBP) and the transmembrane (TM) domain of the Torpedo marmorata nAChR. The second model, referred to as model C, was based solely on the recent electron microscopy structure of the T. marmorata nAChR. Despite structural differences, both models exhibit nearly identical patterns of flexibility and correlated motions. In addition, both models show a similar global twisting motion that may represent channel gating. The similar results obtained for the two models indicate that NMA is most sensitive to the contact topology of the structure rather than its finer detail. The major difference between the low-frequency motions sampled for the two models is that a symmetrical pore-breathing motion, favoring channel opening, is present as the second most dominant motion in model I, whilst largely absent from model C. The absence of this mode in model C can be attributed to its less symmetrical architecture. Finally, as a further goal of the present study, an approximate open channel model, consistent with many experimental findings, has been produced.

Humans↗

PubMed2005

A gating mechanism proposed from a simulation of a human alpha7 nicotinic acetylcholine receptor.

The nicotinic acetylcholine receptor is a well characterized ligand-gated ion channel, yet a proper description of the mechanisms involved in gating, opening, closing, ligand binding, and desensitization does not exist. Until recently, atomic-resolution structural information on the protein was limited, but with the production of the x-ray crystal structure of the Lymnea stagnalis acetylcholine binding protein and the EM image of the transmembrane domain of the torpedo electric ray nicotinic channel, we were provided with a window to examine the mechanism by which this channel operates. A 15-ns all-atom simulation of a homology model of the homomeric human alpha7 form of the receptor was conducted in a solvated palmitoyl-2-oleoyl-sn-glycerol-phosphatidylcholine bilayer and examined in detail. The receptor was unliganded. The structure undergoes a twist-to-close motion that correlates movements of the C loop in the ligand binding domain, via the beta10-strand that connects the two, with the 10 degrees rotation and inward movement of two nonadjacent subunits. The Cys loop appears to act as a stator around which the alpha-helical transmembrane domain can pivot and rotate relative to the rigid beta-sheet binding domain. The M2-M3 loop may have a role in controlling the extent or kinetics of these relative movements. All of this motion, along with essential dynamics analysis, is suggestive of the direction of larger motions involved in gating of the channel.

Binding Sites↗

PubMed2005

Membrane protein structure quality in molecular dynamics simulation.

Our goal was to assess the relationship between membrane protein quality, output from protein quality checkers and output from molecular dynamics (MD) simulations. Membrane transport proteins are essential for a wide range of cellular processes. Structural features of integral membrane proteins are still under-explored due to experimental limitations in structure determination. Computational techniques can be used to exploit biochemical and medium resolution structural data, as well as sequence homology to known structures, and enable us to explore the structure-function relationships in several transmembrane proteins. The quality of the models produced is vitally important to obtain reliable predictions. An examination of the relationship between model stability in molecular dynamics (MD) simulations derived from RMSD (root mean squared deviation) and structure quality assessment from various protein quality checkers was undertaken. The results were compared to membrane protein structures, solved at various resolution, by either X-ray or electron diffraction techniques. The checking programs could predict the potential success of MD in making functional conclusions. MD stability was shown to be a good indicator for the quality of structures. The quality was also shown to be dependent on the resolution at which the structures were determined.

Computer Simulation↗

PubMed2005

Structure and function of two-pore-domain K+ channels: contributions from genetic model organisms.

K+ channels that possess two pore domains in each channel subunit are common in many animal tissues. Such channels are generated from large families of subunits and are implicated in several functions, including temperature sensation, responses to ischaemia, K+ homeostasis and setting the resting potential of the cell. Their activity can be modulated by polyunsaturated fatty acids, pH and oxygen, and some are candidate targets of volatile anaesthetics. However, despite their potential as targets for novel drugs for human health, comparatively little is known about the molecular basis of their diverse physiological and pharmacological properties. Genetic model organisms have considerable potential for improving our understanding of these channels. In this article, we review the contributions of some of these genetic model organisms to recent advances in our knowledge of two-pore-domain K+ channels.

Animals↗

PubMed2004

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1.

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.

Aquaporin 1↗

PubMed2002

Water transporters: how so fast yet so selective?

A high-resolution X-ray structure of an aquaporin has revealed water molecules bound within the transmembrane pore and provided new clues to the mechanisms of rapid water transport and high selectivity in this important class of membrane proteins.

Aquaporin 1↗

PubMed2002

Water in ion channels and pores--simulation studies.

The microscopic properties of water in narrow pores are relevant to the function of ion channels and related membrane transport proteins. The emergence of several high-resolution structures allows one to perform molecular dynamics simulation studies of water in such pores. Simulations of bundles of parallel alpha-helical peptides (e.g. alamethicin) have enabled development of methodologies and concepts appropriate to such investigations. In the narrow channels formed by such bundles, water molecules exhibit reduced rotational and translation motion. This reduction in water mobility may be a general property of narrow pores. We have used simplified channel models to explore the role of hydrophobicity/hydrophilicity in the entry of water into pores. Narrow pores with a hydrophobic lining, although physically open, may not admit water molecules, acting as a 'hydrophobic gate' that prevents water and ion permeation. Such a gate can be opened either by widening the pore or making its lining more polar. Simulations have been used to explore the behaviour of water in GlpF, a member of the aquaporin family of water pores, and OmpA, a bacterial outer membrane protein. Preliminary results suggest that a continuous water wire is not formed within the amphipathic GlpF pore. Simulations of OmpA, in which polar residues line the channel, indicate that a small conformational change in one of the channel lining side chains may open the channel. In summary, comparison of the behaviour of water in different narrow transmembrane pores suggests that an amphipathic pore is ideal for water permeation, and that either a highly hydrophobic pore lining or a charged pore-lining region can act as a gate.

Computer Simulation↗