Research Topics - Membrane Biophysics
Exchange of materials and information between a living cell and its environment is mediated and regulated by membrane proteins. These proteins are involved in the regulation of electrical activity of the cell, transport of water and water soluble materials across the membrane, and production of ATP. Membrane receptors are the sites for detection informational signals, such as neurotransmitters and hormones, light, and even mechanical stress. The atomic-resolution structures of a few membrane proteins have been solved, and recent advances in structure determination of membrane proteins promise more structure to be solved soon. Our group studies the mechanism of function of membrane proteins with various computational methodologies. The proteins are simulated in their natural environment, i.e., embedded in fully hydrated patches of lipid bilayers and under constant pressure and temperature conditions. The main objective is to understand how specific structural motifs and/or chemical interactions in a protein play a role in its function.
image size:
82.5KB
made with VMD
Aquaporins are membrane water channels that play critical roles in controlling the water contents of cells. These channels are widely distributed in all kingdoms of life, including bacteria, plants, and mammals. More than ten different aquaporins have been found in human body, and several diseases, such as congenital cataracts and nephrogenic diabetes insipidus, are connected to the impaired function of these channels. They form tetramers in the cell membrane, and facilitate the transport of water and, in some cases, other small solutes across the membrane. However, the water pores are completely impermeable to charged species, such as protons, a remarkable property that is critical for the conservation of membrane's electrochemical potential, but paradoxical at the same time, since protons can usually be transfered readily through water molecules. The results of our simulations have now provided new insight into the mechanism underlying this fascinating property. Water molecules passing the channel are forced, by the protein's electrostatic forces, to flip at the center of the channel (see the animation), thereby breaking the alternative donor-acceptor arrangement that is necessary for proton translocation (read the complete story in our Science paper).
image size:
153.6KB
made with VMD
High-density lipoproteins (HDL), often called "good cholesterol", are protein-lipid particles which circulate in the blood collecting cholesterol from peripheral tissues and transporting them to the liver for degradation. Native HDL are heterogeneous particles which exhibit a variety of shapes and sizes, thus making structural studies of the major protein component apolipoprotein A-I (apo A-I) difficult. However, nanodiscs which are reconstituted discoidal HDL mimics being developed as platforms in which to embed membrane proteins, can be assembled into homogeneous particles. Thus, we utilized the extensive characterization of nanodiscs in furthering our understanding of the structure of apo A-I as well as the assembly of lipoprotein particles. Whereas, the structure of lipoprotein particles can be studied using all-atom molecular dynamics, the long time scales needed for assembly simulations required the development of a coarse grained protein-lipid model.
image size:
215.8KB
made with VMD
image size:
40.4KB
made with VMD
Anyone who has attempted to fit a long piece of thread through a needle's eye realizes how difficult fitting something so small and flexible into such a small hole can be. Yet this action is carried out every second in every living cell. Flexible polypeptides, proteins, often have to cross a cellular membrane to get to their correct location, whether that location is an organelle within the cell or even outside of it. To accomplish this, they are pushed through a protein pore in an unfolded conformation much like a long string. The channel that accepts the string-like proteins, the protein translocon, allows only certain proteins to pass, while restricting access to molecules even much smaller than the macromolecular proteins. As reported in a recent publication, computer simulations using the molecular dynamics program NAMD helped to answer the question of how such a small channel could achieve this feat, demonstrating how the channel itself can be flexible yet resilient during a protein-crossing event and also elucidating in part how it can maintain such tight control over what is permitted to cross. For more information, see our Protein Translocation website.
In a biological cell, membrane channels act like miniature valves regulating the flow of ions and other solutes between intracellular compartments and across the cell's boundary. Assembled in complex circuits, they generate, transmit, and amplify signals orchestrating cell function. To investigate how membrane channels work, life scientists, using an extremely fine pipette, isolate a tiny patch of a cell membrane and, in so-called patch clamp measurements, determine electric currents in response to applied electric potentials. Dramatic increase in computational power and its efficient utilization by NAMD allows one today to reproduce such studies computationally, calculating the permeability of a membrane channel to ions and water directly from its atomic structure. In what is one of the largest molecular dynamics simulation to date, described in a recent paper as well as on our web site (here), one copy of the membrane channel alpha-hemolysin, submerged in a lipid membrane and water, was subject to an external electric field that drove ions and water through the channel. The calculations produced also an image of the electrostatic potential across the channel (see figure).
One of the best tools for making the connection between the structure and function of ion channels is molecular dynamics~(MD) simulation, which allows one to foll ow conformational changes in the structure, and movement of K+ ions across the potassium channel. The difficulty in simulating the passing of K+ or Na+ ions through the channel can be overcome by using SMD.
Responding to a great variety of ligands including hormones, neurotransmitters, and peptides, G-protein-coupled receptors (GPCRs) are the most important drug targets. Yet there is still little known about the first step of signal transduction. Rhodopsin is a member of this largest group of transmembrane receptors. The high resolution structure of rhodopsin determined by Palzcewski et al. (2000) gives us the opportunity to take a close look at the activation mechanism. Rhodopsin's ligand, retinal, which is covalently bound to the binding pocket, absorbs light at 500nm and undergoes a cis-trans isomerization triggering conformational changes in the protein. In order to understand the effects of the isomerization on the protein structure, we are performing molecular dynamics of a 40,000 atom model of rhodopsin in membrane.
Adenosine triphosphate (ATP) is the primary energy "currency" in most living organisms. ATP synthase is a large (about 100,000 atoms) protein, which includes a transmembrane F0 unit coupled to a solvent-exposed F1 unit via a central stalk gamma. The F0 unit utilizes a transmembrane electrochemical potential (proton motive force), converting it into the mechanical energy of the stalk rotation. The rotation leads to cyclic conformational changes in the catalytic sites in the F1 unit, thereby driving ATP synthesis. We seek to identify and explore the chain of the elementary chemical (proton transfer) and mechanical (domain motion) events involved in the process of converting the electrochemical energy of the transmembrane proton gradient into the mechanical energy of the c subunit oligomer rotation.
The purple membrane (PM) of halobacterium salinarium consists of bacteriorhodopsin trimers, lipids and water molecules. Upon light absorption by the chromophore inside bacteriorhodopsin, protons are pumped across the PM to the extra-cellular side of the membrane. With molecular dynamics simulations of the whole PM, we are studying the influence of the native environment on the protein dynamics.
image size:
139.9KB
made with VMD
We are applying the steered molecular dynamics method to investigate the action of human synovial protein phospholipase A2 (PLA2) at the lipid water interface. Our hypothesis is that prior to extruding the phospholipid, PLA2 must form the tightly bound complex, while the loosely bound complex should not lead to catalysis.
image size:
301.5KB
made with VMD
Formation of bacteriorhodopsin (bR) from the apoprotein and retinal has been studied experimentally, but the actual pathway, including the site of retinal entry, is little understood. Molecular dynamics simulations provide a surprisingly clear prediction.




