PubMed Health⌕ Search

Biomedical subjects

Emad Tajkhorshid

Publications and source records attributed to Emad Tajkhorshid.

At least 19 recordsLinked to original sources
PubMed2006

Sugar binding and protein conformational changes in lactose permease.

Lactose permease is an integral membrane protein that uses the cell membrane's proton gradient for import of lactose. Based on extensive biochemical data and a substrate-bound crystal structure, intermediates involved in lactose/H(+) co-transport have been suggested. Yet, the transport mechanism, especially the coupling of protonation states of essential residues and protein conformational changes involved in the transport, is not understood. Here we report molecular-dynamics simulations of membrane-embedded lactose permease in different protonation states, both in the presence and in the absence of lactose. The results analyzed in terms of pore diameter, salt-bridge formation, and substrate motion, strongly implicate Glu(269) as one of the main proton translocation sites, whose protonation state controls several key steps of the transport process. A critical ion pair (Glu(269) and Arg(144)) was found to keep the cytoplasmic entrance open, but via a different mechanism than the currently accepted model. After protonation of Glu(269), the salt bridge between Glu(269) and Arg(144) was found to break, and Arg(144) to move away from Glu(269), establishing a new salt bridge with Glu(126); furthermore, neutralization of Glu(269) and the displacement of Arg(144) and consequently of water molecules from the interdomain region was seen to initiate the closing of the cytoplasmic half channel (2.6-4.0 A reduction in diameter in the cytoplasmic constriction region in 10 ns) by allowing hydrophobic surfaces of the N- and C-domains to fuse. Charged Glu(269) was found to strongly bind the lactose permeant, indicating that proton transfer from water or another residue to Glu(269) is a prerequisite for unbinding of lactose from the binding pocket.

Arginine↗
PubMed2006

Color tuning in rhodopsins: the mechanism for the spectral shift between bacteriorhodopsin and sensory rhodopsin II.

The mechanism of color tuning in the rhodopsin family of proteins has been studied by comparing the optical properties of the light-driven proton pump bacteriorhodopsin (bR) and the light detector sensory rhodopsin II (sRII). Despite a high structural similarity, the maximal absorption is blue-shifted from 568 nm in bR to 497 nm in sRII. The molecular mechanism of this shift is still a matter of debate, and its clarification sheds light onto the general mechanisms of color tuning in retinal proteins. The calculations employ a combined quantum mechanical/molecular mechanical (QM/MM) technique, using a DFT-based method for ground state properties and the semiempirical OM2/MRCI method and ab initio SORCI method for excited state calculations. The high efficiency of the methodology has allowed us to study a wide variety of aspects including dynamical effects. The absorption shift as well as various mutation experiments and vibrational properties have been successfully reproduced. Our results indicate that several sources contribute to the spectral shift between bR and sRII. The main factors are the counterion region at the extracellular side of retinal and the amino acid composition of the binding pocket. Our analysis allows a distinction and identification of the different effects in detail and leads to a clear picture of the mechanism of color tuning, which is in good agreement with available experimental data.

Bacteriorhodopsins↗
PubMed2006

Dynamics of K+ ion conduction through Kv1.2.

The crystallographic structure of a potassium channel, Kv1.2, in an open state makes it feasible to simulate entire K(+) ion permeation events driven by a voltage bias and, thereby, elucidate the mechanism underlying ion conduction and selectivity of this type of channel. This Letter demonstrates that molecular dynamics simulations can provide movies of the overall conduction of K(+) ions through Kv1.2. As suggested earlier, the conduction is concerted in the selectivity filter, involving 2-3 ions residing mainly at sites identified previously by crystallography and modeling. The simulations reveal, however, the jumps of ions between these sites and identify the sequence of multi-ion configurations involved in permeation.

Computer Simulation↗
PubMed2006

Mechanism of gating and ion conductivity of a possible tetrameric pore in aquaporin-1.

While substrate permeation through monomeric pores of aquaporins is well characterized, little is known about the possible tetrameric pore. AQP1 has been suggested to function as an ion channel upon cGMP activation, although this idea has been controversial. Taking a theoretical and experimental approach, we demonstrate that the current might arise through the tetrameric pore and propose a plausible mechanism for conduction and gating. In response to simulated ion permeation, immediate hydration of the putative central pore was facilitated by moderate conformational changes of pore-lining residues. cGMP is found to interact with an unusually arginine-rich, cytoplasmic loop (loop D) facilitating its outward motion, which is hypothesized to trigger the opening of a cytoplasmic gate. Physiological analyses of wild-type AQP1 and a designed mutant in which two arginines of the gating loop are replaced by alanine provide experimental support for identifying a key component of the proposed mechanism.

Amino Acid Sequence↗
PubMed2005

Structural mechanism of plant aquaporin gating.

Plants counteract fluctuations in water supply by regulating all aquaporins in the cell plasma membrane. Channel closure results either from the dephosphorylation of two conserved serine residues under conditions of drought stress, or from the protonation of a conserved histidine residue following a drop in cytoplasmic pH due to anoxia during flooding. Here we report the X-ray structure of the spinach plasma membrane aquaporin SoPIP2;1 in its closed conformation at 2.1 A resolution and in its open conformation at 3.9 A resolution, and molecular dynamics simulations of the initial events governing gating. In the closed conformation loop D caps the channel from the cytoplasm and thereby occludes the pore. In the open conformation loop D is displaced up to 16 A and this movement opens a hydrophobic gate blocking the channel entrance from the cytoplasm. These results reveal a molecular gating mechanism which appears conserved throughout all plant plasma membrane aquaporins.

Aquaporins↗
PubMed2005

Determination of structural and functional overlap/divergence of five proto-type galectins by analysis of the growth-regulatory interaction with ganglioside GM1 in silico and in vitro on human neuroblastoma cells.

The growth-regulatory interplay between ganglioside GM1 on human SK-N-MC neuroblastoma cells and an endogenous lectin provides a telling example for glycan (polysaccharide) functionality. Galectin-1 is the essential link between the sugar signal and the intracellular response. The emerging intrafamily complexity of galectins raises the question on defining extent of their structural and functional overlap/divergence. We address this problem for proto-type galectins in this system: ganglioside GM1 as ligand, neuroblastoma cells as target. Using the way human galectin-1 interacts with this complex natural ligand as template, we first defined equivalent positioning for distinct substitutions in the other tested proto-type galectins, e.g., Lys63 vs. Leu60/Gln72 in galectins-2 and -5. As predicted from our in silico work, the tested proto-type galectins have affinity for the pentasaccharide of ganglioside GM1. In contrast to solid-phase assays, cell surface presentation of the ganglioside did not support binding of galectin-5, revealing the first level of regulation. Next, a monomeric proto-type galectin (CG-14) can impair galectin-1-dependent negative growth control by competitively blocking access to the shared ligand without acting as effector. Thus, the quaternary structure of proto-type galectins is an efficient means to give rise to functional divergence. The identification of this second level of regulation is relevant for diagnostic monitoring. It might be exploited therapeutically by producing galectin variants tailored to interfere with galectin activities associated with the malignant phenotype. Moreover, the given strategy for comparative computational analysis of extended binding sites has implications for the rational design of galectin-type-specific ligands.

Amino Acid Sequence↗
PubMed2005

Scalable molecular dynamics with NAMD.

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

Algorithms↗
PubMed2005

Molecular dynamics simulations of proteins in lipid bilayers.

With recent advances in X-ray crystallography of membrane proteins promising many new high-resolution structures, molecular dynamics simulations will become increasingly valuable for understanding membrane protein function, as they can reveal the dynamic behavior concealed in the static structures. Dramatic increases in computational power, in synergy with more efficient computational methodologies, now allow us to carry out molecular dynamics simulations of any structurally known membrane protein in its native environment, covering timescales of up to 0.1 micros. At the frontiers of membrane protein simulations are ion channels, aquaporins, passive and active transporters, and bioenergetic proteins.

Computer Simulation↗
PubMed2005

What makes an aquaporin a glycerol channel? A comparative study of AqpZ and GlpF.

The recent availability of high-resolution structures of two structurally highly homologous, but functionally distinct aquaporins from the same species, namely Escherichia coli AqpZ, a pure water channel, and GlpF, a glycerol channel, presents a unique opportunity to understand the mechanism of substrate selectivity in these channels. Comparison of the free energy profile of glycerol conduction through AqpZ and GlpF reveals a much larger barrier in AqpZ (22.8 kcal/mol) than in GlpF (7.3 kcal/mol). In either channel, the highest barrier is located at the selectivity filter. Analysis of substrate-protein interactions suggests that steric restriction of AqpZ is the main contribution to this large barrier. Another important difference is the presence of a deep energy well at the periplasmic vestibule of GlpF, which was not found in AqpZ. The latter difference can be attributed to the more pronounced structural asymmetry of GlpF, which may play a role in attracting glycerol.

Aquaporins↗
PubMed2004

Collective diffusion model for water permeation through microscopic channels.

Water permeation through nanometric channels, generally a coupled many-body process, is described in this study by a single collective coordinate. A collective diffusion model is proposed in which water movement at equilibrium is characterized as an unbiased diffusion along this coordinate and water transport in the presence of a chemical potential difference is described as one-dimensional diffusion in a linear potential. The model allows one to determine the osmotic permeability of a water channel from equilibrium simulations.

Aquaporins↗
PubMed2004

Role of hydrogen-bond network in energy storage of bacteriorhodopsin's light-driven proton pump revealed by ab initio normal-mode analysis.

Vibrational modes of the hydrogen-bond network in the binding site of bacteriorhodopsin (bR), a protein in halobacteria functioning as a light-driven proton pump, were investigated by an ab initio quantum mechanical/molecular mechanical (QM/MM) method. Normal-mode analysis calculations for O-D and N-D stretching modes of internal water molecules and the Schiff base of the retinal chromophore in the early intermediate state, K, reproduced well experimentally observed vibrational spectra. Supported by agreement with observed spectra, the QM/MM calculation suggests that weakened hydrogen bonds upon photoisomerization of the chromophore are an important means of energy storage in bR.

Bacteriorhodopsins↗
PubMed2004

The mechanism of proton exclusion in aquaporin channels.

The mechanism of proton exclusion in aquaporin channels is elucidated through free energy calculations of the pathway of proton transport. The second generation multistate empirical valence bond (MS-EVB2) model was applied to simulate the interaction of an excess proton with the channel environment. Jarzynski's equality was employed for rapid convergence of the free energy profile. A barrier sufficiently high to block proton transport is located near the channel center at the NPA motif-a site involved in bi-orientational ordering of the embedded water-wire in absence of the excess proton. A second and lower barrier is observed at the selectivity filter near the periplasmic outlet where the channel is narrowest. This secondary barrier may be essential in filtering other large solutes and cations.

Aquaporins↗
PubMed2004

Classical force field parameters for the heme prosthetic group of cytochrome c.

Accurate force fields are essential for describing biological systems in a molecular dynamics simulation. To analyze the docking of the small redox protein cytochrome c (cyt c) requires simulation parameters for the heme in both the reduced and oxidized states. This work presents parameters for the partial charges and geometries for the heme in both redox states with ligands appropriate to cyt c. The parameters are based on both protein X-ray structures and ab initio density functional theory (DFT) geometry optimizations at the B3LYP/6-31G* level. The simulations with the new parameter set reproduce the geometries of the X-ray structures and the interaction energies between water and heme prosthetic group obtained from B3LYP/6-31G* calculations. The parameter set developed here will provide new insights into docking processes of heme containing redox proteins.

Algorithms↗
PubMed2004

Molecular basis of proton blockage in aquaporins.

Water transport channels in membrane proteins of the aquaporin superfamily are impermeable to ions, including H+ and OH-. We examine the molecular basis for the blockage of proton translocation through the single-file water chain in the pore of a bacterial aquaporin, GlpF. We compute the reversible thermodynamic work for the two complementary steps of the Grotthuss "hop-and-turn" relay mechanism: consecutive transfers of H+ along the hydrogen-bonded chain (hop) and conformational reorganization of the chain (turn). In the absence of H+, the strong preference for the bipolar orientation of water around the two Asn-Pro-Ala (NPA) motifs lining the pore over both unidirectional polarization states of the chain precludes the reorganization of the hydrogen-bonded network. Inversely, translocation of an excess proton in either direction is opposed by a free-energy barrier centered at the NPA region. Both hop and turn steps of proton translocation are opposed by the electrostatic field of the channel.

Aquaporins↗
PubMed2004

Theory and simulation of water permeation in aquaporin-1.

We discuss the difference between osmotic permeability pf and diffusion permeability pd of single-file water channels and demonstrate that the pf/pd ratio corresponds to the number of effective steps a water molecule needs to take to permeate a channel. While pd can be directly obtained from equilibrium molecular dynamics simulations, pf can be best determined from simulations in which a chemical potential difference of water has been established on the two sides of the channel. In light of this, we suggest a method to induce in molecular dynamics simulations a hydrostatic pressure difference across the membrane, from which pf can be measured. Simulations using this method are performed on aquaporin-1 channels in a lipid bilayer, resulting in a calculated pf of 7.1 x 10(-14) cm(3)/s, which is in close agreement with observation. Using a previously determined pd value, we conclude that pf/pd for aquaporin-1 measures approximately 12. This number is explained in terms of channel architecture and conduction mechanism.

Aquaporin 1↗
PubMed2003

Unique conformer selection of human growth-regulatory lectin galectin-1 for ganglioside GM1 versus bacterial toxins.

Endogenous lectins induce effects on cell growth by binding to antennae of natural glycoconjugates. These complex carbohydrates often present more than one potential lectin-binding site in a single chain. Using the growth-regulatory interaction of the pentasaccharide of ganglioside GM(1) with homodimeric galectin-1 on neuroblastoma cell surfaces as a model, we present a suitable strategy for addressing this issue. The approach combines NMR spectroscopic and computational methods and does not require isotope-labeled glycans. It involves conformational analysis of the two building blocks of the GM(1) glycan, i.e., the disaccharide Galbeta1-3GalNAc and the trisaccharide Neu5Acalpha2-3Galbeta1-4Glc. Their bound-state conformations were determined by transferred nuclear Overhauser enhancement spectroscopy. Next, measurements on the lectin-pentasaccharide complex revealed differential conformer selection regarding the sialylgalactose linkage in the tri- versus pentasaccharide (Phi and Psi value of -70 degrees and 15 degrees vs 70 degrees and 15 degrees, respectively). To proceed in the structural analysis, the characteristic experimentally detected spatial vicinity of a galactose unit and Trp68 in the galectin's binding site offered a means, exploiting saturation transfer from protein to carbohydrate protons. Indeed, we detected two signals unambiguously assigned to the terminal Gal and the GalNAc residues. Computational docking and interaction energy analyses of the entire set of ligands supported and added to experimental results. The finding that the ganglioside's carbohydrate chain is subject to differential conformer selection at the sialylgalactose linkage by galectin-1 and GM(1)-binding cholera toxin (Phi and Psi values of -172 degrees and -26 degrees, respectively) is relevant for toxin-directed drug design. In principle, our methodology can be applied in studies aimed at blocking galectin functionality in malignancy and beyond glycosciences.

Acetylgalactosamine↗
PubMed2003

Developing an energy landscape for the novel function of a (beta/alpha)8 barrel: ammonia conduction through HisF.

HisH-hisF is a multidomain globular protein complex; hisH is a class I glutamine amidotransferase that hydrolyzes glutamine to form ammonia, and hisF is a (beta/alpha)8 barrel cyclase that completes the ring formation of imidizole glycerol phosphate synthase. Together, hisH and hisF form a glutamine amidotransferase that carries out the fifth step of the histidine biosynthetic pathway. Recently, it has been suggested that the (beta/alpha)8 barrel participates in a novel function: to channel ammonia from the active site of hisH to the active site of hisF. The present study presents a series of molecular dynamic simulations that investigate the channeling function of hisF. This article reconstructs potentials of mean force for the conduction of ammonia through the channel, and the entrance of ammonia through the strictly conserved channel gate, in both a closed and a hypothetical open conformation. The resulting energy landscape within the channel supports the idea that ammonia does indeed pass through the barrel, interacting with conserved hydrophilic residues along the way. The proposed open conformation, which involves an alternate rotamer state of one of the gate residues, presents only an approximately 2.5-kcal energy barrier to ammonia entry. Another alternate open-gate conformation, which may play a role in non-nitrogen-fixing organisms, is deduced through bioinformatics.

Algorithms↗