Search

Binding of a T-cell receptor (TCR) to a peptide/major histocompatibility complex is the key interaction involved in antigen specificity of T cells. The recognition involves up to six complementarity determining regions (CDR) of the TCR. Efforts to examine the structural basis of these interactions and to exploit them in adoptive T-cell therapies has required the isolation of specific T-cell clones and their clonotypic TCRs. Here we describe a strategy using in vitro-directed evolution of a single TCR to change its peptide specificity, thereby avoiding the need to isolate T-cell clones. The human TCR A6, which recognizes the viral peptide Tax/HLA-A2, was converted to TCR variants that recognized the cancer peptide MART1/HLA-A2. Mutational studies and molecular dynamics simulations identified CDR residues that were predicted to be important in the specificity switch. Thus, in vitro engineering strategies alone can be used to discover TCRs with desired specificities.

The SNARE complex plays a vital role in vesicle fusion arising during neuronal exocytosis. Key components in the regulation of SNARE complex formation, and ultimately fusion, are the transmembrane and linker regions of the vesicle-associated protein, synaptobrevin. However, the membrane-embedded structure of synaptobrevin in its prefusion state, which determines its interaction with other SNARE proteins during fusion, is largely unknown. This study reports all-atom molecular-dynamics simulations of the prefusion configuration of synaptobrevin in a lipid bilayer, aimed at characterizing the insertion depth and the orientation of the protein in the membrane, as well as the nature of the amino acids involved in determining these properties. By characterizing the structural properties of both wild-type and mutant synaptobrevin, the effects of C-terminal additions on tilt and insertion depth of membrane-embedded synaptobrevin are determined. The simulations suggest a robust, highly tilted state for membrane-embedded synaptobrevin with a helical connection between the transmembrane and linker regions, leading to an apparently new characterization of structural elements in prefusion synaptobrevin and providing a framework for interpreting past mutation experiments.

Experimental challenges associated with characterization of the membrane-bound form of talin have prevented us from understanding the molecular mechanism of its membrane-dependent integrin activation. Here, utilizing what we believe to be a novel membrane mimetic model, we present a reproducible model of membrane-bound talin observed across multiple independent simulations. We characterize both local and global membrane-induced structural transitions that successfully reconcile discrepancies between biochemical and structural studies and provide insight into how talin might modulate integrin function. Membrane binding of talin, captured in unbiased simulations, proceeds through three distinct steps: initial electrostatic recruitment of the F2 subdomain to anionic lipids via several basic residues; insertion of an initially buried, conserved hydrophobic anchor into the membrane; and association of the F3 subdomain with the membrane surface through a large, interdomain conformational change. These latter two steps, to our knowledge, have not been observed or described previously. Electrostatic analysis shows talin F2F3 to be highly polarized, with a highly positive underside, which we attribute to the initial electrostatic recruitment, and a negative top face, which can help orient the protein optimally with respect to the membrane, thereby reducing the number of unproductive membrane collision events.

While small molecules have been used to induce anesthesia in a clinical setting for well over a century, a detailed understanding of the molecular mechanism remains elusive. In this study, we utilize ab initio calculations to develop a novel set of CHARMM-compatible parameters for the ubiquitous modern anesthetics desflurane, isoflurane, sevoflurane, and propofol for use in molecular dynamics (MD) simulations. The parameters generated were rigorously tested against known experimental physicochemical properties including dipole moment, density, enthalpy of vaporization, and free energy of solvation. In all cases, the anesthetic parameters were able to reproduce experimental measurements, signifying the robustness and accuracy of the atomistic models developed. The models were then used to study the interaction of anesthetics with the membrane. Calculation of the potential of mean force for inserting the molecules into a POPC bilayer revealed a distinct energetic minimum of 4-5 kcal/mol relative to aqueous solution at the level of the glycerol backbone in the membrane. The location of this minimum within the membrane suggests that anesthetics partition to the membrane prior to binding their ion channel targets, giving context to the Meyer-Overton correlation. Moreover, MD simulations of these drugs in the membrane give rise to computed membrane structural parameters, including atomic distribution, deuterium order parameters, dipole potential, and lateral stress profile, that indicate partitioning of anesthetics into the membrane at the concentration range studied here, which does not appear to perturb the structural integrity of the lipid bilayer. These results signify that an indirect, membrane-mediated mechanism of channel modulation is unlikely.

The respiratory chains of nearly all aerobic organisms are terminated by proton-pumping heme-copper oxygen reductases (HCOs). Previous studies have established that C-family HCOs contain a single channel for uptake from the bacterial cytoplasm of all chemical and pumped protons, and that the entrance of the K(C)-channel is a conserved glutamate in subunit III. However, the majority of the K(C)-channel is within subunit I, and the pathway from this conserved glutamate to subunit I is not evident. In the present study, molecular dynamics simulations were used to characterize a chain of water molecules leading from the cytoplasmic solution, passing the conserved glutamate in subunit III and extending into subunit I. Formation of the water chain, which controls the delivery of protons to the K(C)-channel, was found to depend on the conformation of Y241(Vc), located in subunit I at the interface with subunit III. Mutations of Y241(Vc) (to A/F/H/S) in the Vibrio cholerae cbb3 eliminate catalytic activity, but also cause perturbations that propagate over a 28-Å distance to the active site heme b3. The data suggest a linkage between residues lining the K(C)-channel and the active site of the enzyme, possibly mediated by transmembrane helix α7, which contains both Y241(Vc) and the active site cross-linked Y255(Vc), as well as two CuB histidine ligands. Other mutations of residues within or near helix α7 also perturb the active site, indicating that this helix is involved in modulation of the active site of the enzyme.

We report here a detailed thermodynamic description of water molecules inside a biological water channel. Taking advantage of high-resolution molecular dynamics trajectories calculated for an aquaporin (AQP) channel, we compute the spatial translational and rotational components of water diffusion and entropy in AQP. Our results reveal that the spontaneous filling and entry of water into the pore in AQPs are driven by an entropic gain. Specifically, water molecules exhibit an elevated degree of rotational motion inside the pore, while their translational motion is slow compared with bulk. The partial charges of the lining asparagine residues at the conserved signature Asn-Pro-Ala motifs play a key role in enhancing rotational diffusion and facilitating dipole flipping of water inside the pore. The frequencies of the translational and rotational motions in the power spectra overlap indicating a strong coupling of these motions in AQPs. A shooting mechanism with diffusive behavior is observed in the extracellular region which might be a key factor in the fast conduction of water in AQPs.

Characterizing large-scale structural transitions in biomolecular systems poses major technical challenges to both experimental and computational approaches. On the computational side, efficient sampling of the configuration space along the transition pathway remains the most daunting challenge. Recognizing this issue, we introduce a knowledge-based computational approach toward describing large-scale conformational transitions using (i) nonequilibrium, driven simulations combined with work measurements and (ii) free energy calculations using empirically optimized biasing protocols. The first part is based on designing mechanistically relevant, system-specific reaction coordinates whose usefulness and applicability in inducing the transition of interest are examined using knowledge-based, qualitative assessments along with nonequilirbrium work measurements which provide an empirical framework for optimizing the biasing protocol. The second part employs the optimized biasing protocol resulting from the first part to initiate free energy calculations and characterize the transition quantitatively. Using a biasing protocol fine-tuned to a particular transition not only improves the accuracy of the resulting free energies but also speeds up the convergence. The efficiency of the sampling will be assessed by employing dimensionality reduction techniques to help detect possible flaws and provide potential improvements in the design of the biasing protocol. Structural transition of a membrane transporter will be used as an example to illustrate the workings of the proposed approach.

Multidrug ATP binding cassette (ABC) exporters are ubiquitous ABC transporters that extrude cytotoxic molecules across cell membranes. Despite recent progress in structure determination of these transporters, the conformational motion that transduces the energy of ATP hydrolysis to the work of substrate translocation remains undefined. Here, we have investigated the conformational cycle of BmrCD, a representative of the heterodimer family of ABC exporters that have an intrinsically impaired nucleotide binding site. We measured distances between pairs of spin labels monitoring the movement of the nucleotide binding (NBD) and transmembrane domains (TMD). The results expose previously unobserved structural intermediates of the NBDs arising from asymmetric configuration of catalytically inequivalent nucleotide binding sites. The two-state transition of the TMD, from an inward- to an outward-facing conformation, is driven exclusively by ATP hydrolysis. These findings provide direct evidence of divergence in the mechanism of ABC exporters.DOI: http://dx.doi.org/10.7554/eLife.02740.001., (Copyright © 2014, Mishra et al.)

The potassium channel KcsA offers a unique opportunity to explicitly study the dynamics of the moving parts of ion channels, yet our understanding of the extent and dynamic behavior of the physiologically relevant structural changes at the inner gate in KcsA remains incomplete. Here, we use electron paramagnetic resonance, nuclear magnetic resonance, and molecular dynamics simulations to characterize the extent of pH-dependent conformational changes of the inner gate in lipid bilayers or detergent micelles. Our results show that under physiological conditions the inner gate experiences a maximal diagonal opening of ∼24 Å with the largest degree of dynamics near the pKa of activation (pH ∼3.9). These results extend the observation that the C-terminus is necessary to limit the extent of opening and imply that the inner gate regulates the extent of conformational change at the zone of allosteric coupling and at the selectivity filter.

Understanding the process of membrane insertion is an essential step in developing a detailed mechanism, for example, for peripheral membrane protein association and membrane fusion. The highly mobile membrane mimetic (HMMM) has been used to accelerate the membrane association and binding of peripheral membrane proteins in simulations by increasing the lateral diffusion of phospholipid headgroups while retaining an atomistic description of the interface. Through a comparative study, we assess the difference in insertion rate of a free phospholipid into an HMMM as well as into a conventional phospholipid bilayer and develop a detailed mechanistic model of free phospholipid insertion into biological membranes. The mechanistic insertion model shows that successful irreversible association of the free phospholipid to the membrane interface, which results in its insertion, is the rate-limiting step. Association is followed by independent, sequential insertion of the acyl tails of the free phospholipid into the membrane, with splayed acyl tail intermediates. Use of the HMMM is found to replicate the same intermediate insertion states as in the full phospholipid bilayer; however, it accelerates overall insertion by approximately a factor of 3, with the probability of successful association of phospholipid to the membrane being significantly enhanced.

Energetics of protein side chain partitioning between aqueous solution and cellular membranes is of fundamental importance for correctly capturing the membrane binding and specific protein-lipid interactions of peripheral membrane proteins. We recently reported a highly mobile membrane mimetic (HMMM) model that accelerates lipid dynamics by modeling the membrane interior partially as a fluid organic solvent while retaining a literal description of the lipid head groups and the beginning of the tails. While the HMMM has been successfully applied to study spontaneous insertion of a number of peripheral proteins into membranes, a quantitative characterization of the energetics of membrane-protein interactions in HMMM membranes has not been performed. We report here the free energy profiles for partitioning of 10 protein side chain analogues into a HMMM membrane. In the interfacial and headgroup regions of the membrane, the side chain free energy profiles show excellent agreement with profiles previously reported for conventional membranes with full-tail lipids. In regions where the organic solvent is prevalent, the increased dipole and fluidity of the solvent generally result in a less accurate description, most notably overstabilization of aromatic and polar amino acids. As an additional measure of the ability of the HMMM model to describe membrane-protein interactions, the water-to-membrane interface transfer energies were analyzed and found to be in agreement with the previously reported experimental and computational hydrophobicity scales. We discuss strengths and weaknesses of HMMM in describing protein-membrane interactions as well as further development of model membranes.

CLC transporters catalyze transmembrane exchange of chloride for protons. Although a putative pathway for Cl(-) has been established, the pathway of H(+) translocation remains obscure. Through a highly concerted computational and experimental approach, we characterize microscopic details essential to understanding H(+)-translocation. An extended (0.4 µs) equilibrium molecular dynamics simulation of membrane-embedded, dimeric ClC-ec1, a CLC from Escherichia coli, reveals transient but frequent hydration of the central hydrophobic region by water molecules from the intracellular bulk phase via the interface between the two subunits. We characterize a portal region lined by E202, E203, and A404 as the main gateway for hydration. Supporting this mechanism, site-specific mutagenesis experiments show that ClC-ec1 ion transport rates decrease as the size of the portal residue at position 404 is increased. Beyond the portal, water wires form spontaneously and repeatedly to span the 15-Å hydrophobic region between the two known H(+) transport sites [E148 (Glu(ex)) and E203 (Glu(in))]. Our finding that the formation of these water wires requires the presence of Cl(-) explains the previously mystifying fact that Cl(-) occupancy correlates with the ability to transport protons. To further validate the idea that these water wires are central to the H(+) transport mechanism, we identified I109 as the residue that exhibits the greatest conformational coupling to water wire formation and experimentally tested the effects of mutating this residue. The results, by providing a detailed microscopic view of the dynamics of water wire formation and confirming the involvement of specific protein residues, offer a mechanism for the coupled transport of H(+) and Cl(-) ions in CLC transporters.

The inability to rapidly generate accurate and robust parameters for novel chemical matter continues to severely limit the application of molecular dynamics simulations to many biological systems of interest, especially in fields such as drug discovery. Although the release of generalized versions of common classical force fields, for example, General Amber Force Field and CHARMM General Force Field, have posited guidelines for parameterization of small molecules, many technical challenges remain that have hampered their wide-scale extension. The Force Field Toolkit (ffTK), described herein, minimizes common barriers to ligand parameterization through algorithm and method development, automation of tedious and error-prone tasks, and graphical user interface design. Distributed as a VMD plugin, ffTK facilitates the traversal of a clear and organized workflow resulting in a complete set of CHARMM-compatible parameters. A variety of tools are provided to generate quantum mechanical target data, setup multidimensional optimization routines, and analyze parameter performance. Parameters developed for a small test set of molecules using ffTK were comparable to existing CGenFF parameters in their ability to reproduce experimentally measured values for pure-solvent properties (<15% error from experiment) and free energy of solvation (±0.5 kcal/mol from experiment)., (Copyright © 2013 Wiley Periodicals, Inc.)

During their transport cycle, ATP-binding cassette (ABC) transporters undergo large-scale conformational changes between inward- and outward-facing states. Using an approach based on designing system-specific reaction coordinates and using nonequilibrium work relations, we have performed extensive all-atom molecular dynamics simulations in the presence of explicit membrane/solvent to sample a large number of mechanistically distinct pathways for the conformational transition of MsbA, a bacterial ABC exporter whose structure has been solved in multiple functional states. The computational approach developed here is based on (i) extensive exploration of system-specific biasing protocols (e.g., using collective variables designed based on available low-resolution crystal structures) and (ii) using nonequilibrium work relations for comparing the relevance of the transition pathways. The most relevant transition pathway identified using this approach involves several distinct stages reflecting the complex nature of the structural changes associated with the function of the protein. The opening of the cytoplasmic gate during the outward- to inward-facing transition of apo MsbA is found to be disfavored when the periplasmic gate is open and facilitated by a twisting motion of the nucleotide-binding domains that involves a dramatic change in their relative orientation. These results highlight the cooperativity between the transmembrane and the nucleotide-binding domains in the conformational transition of ABC exporters. The approach introduced here provides a framework to study large-scale conformational changes of other membrane transporters whose computational investigation at an atomic resolution may not be currently feasible using conventional methods.

During the fusion of the influenza virus to the host cell, bending of the HA2 chain of hemagglutinin into a hairpin-shaped structure in a pH-dependent manner facilitates the fusion of the viral envelope and the endosomal membrane. To characterize the structural and dynamical responses of the hinge region of HA2 to pH changes and examine the role of a conserved histidine in this region (the hinge histidine), we have performed an extensive set of molecular dynamics (MD) simulations of 26-residue peptides encompassing the hinge regions of several hemagglutinin subtypes under both neutral and low pH conditions, modeled by the change of the protonation state of the hinge histidine. More than 70 sets of MD simulations (collectively amounting to 25.1 μs) were performed in both implicit and explicit solvents to study the effect of histidine protonation on structural dynamics of the hinge region. In both explicit and implicit solvent simulations, hinge bending was consistently observed upon the protonation of the histidine in all the simulations starting with an initial straight helical conformation, whereas the systems with a neutral histidine retained their primarily straight conformation throughout the simulations. Conversely, the MD simulations starting from an initially bent conformation resulted in the formation of a straight helical structure upon the neutralization of the hinge histidine, whereas the bent structure was maintained when the hinge histidine remained protonated. Finally, mutation of the hinge histidine to alanine abolishes the bending response of the peptide altogether. A molecular mechanism based on the interaction of the hinge histidine with neighboring acidic residues is proposed to be responsible for its role in controlling the conformation of the hinge. We propose that this might present a common mechanism for pH-controlled structural changes in helical structures when histidines act as the pH sensor., (Copyright © 2013 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

P-glycoprotein (Pgp) is one of the most biomedically relevant transporters in the ATP binding cassette (ABC) superfamily due to its involvement in developing multidrug resistance in cancer cells. Employing molecular dynamics simulations and double electron-electron resonance spectroscopy, we have investigated the structural dynamics of membrane-bound Pgp in the inward-facing state and found that Pgp adopts an unexpectedly wide range of conformations, highlighted by the degree of separation between the two nucleotide-binding domains (NBDs). The distance between the two NBDs in the equilibrium simulations covers a range of at least 20 Å, including, both, more open and more closed NBD configurations than the crystal structure. The double electron-electron resonance measurements on spin-labeled Pgp mutants also show wide distributions covering both longer and shorter distances than those observed in the crystal structure. Based on structural and sequence analyses, we propose that the transmembrane domains of Pgp might be more flexible than other structurally known ABC exporters. The structural flexibility of Pgp demonstrated here is not only in close agreement with, but also helps rationalize, the reported high NBD fluctuations in several ABC exporters and possibly represents a fundamental difference in the transport mechanism between ABC exporters and ABC importers. In addition, during the simulations we have captured partial entrance of a lipid molecule from the bilayer into the lumen of Pgp, reaching the putative drug binding site. The location of the protruding lipid suggests a putative pathway for direct drug recruitment from the membrane.

Aquaporins and Rh proteins can function as gas (CO₂ and NH₃) channels. The present study explores the urea, H₂O, CO₂, and NH₃ permeability of the human urea transporter B (UT-B) (SLC14A1), expressed in Xenopus oocytes. We monitored urea uptake using [¹⁴C]urea and measured osmotic water permeability (Pf) using video microscopy. To obtain a semiquantitative measure of gas permeability, we used microelectrodes to record the maximum transient change in surface pH (ΔpHS) caused by exposing oocytes to 5% CO₂/33 mM HCO₃⁻ (pHS increase) or 0.5 mM NH₃/NH₄⁺ (pHS decrease). UT-B expression increased oocyte permeability to urea by >20-fold, and Pf by 8-fold vs. H₂O-injected control oocytes. UT-B expression had no effect on the CO₂-induced ΔpHS but doubled the NH₃-induced ΔpHS. Phloretin reduced UT-B-dependent urea uptake (Jurea*) by 45%, Pf* by 50%, and (- ΔpHS*)NH₃ by 70%. p-Chloromercuribenzene sulfonate reduced Jurea* by 25%, Pf* by 30%, and (ΔpHS*)NH₃ by 100%. Molecular dynamics (MD) simulations of membrane-embedded models of UT-B identified the monomeric UT-B pores as the main conduction pathway for both H₂O and NH₃ and characterized the energetics associated with permeation of these species through the channel. Mutating each of two conserved threonines lining the monomeric urea pores reduced H₂O and NH₃ permeability. Our data confirm that UT-B has significant H₂O permeability and for the first time demonstrate significant NH₃ permeability. Thus the UTs become the third family of gas channels. Inhibitor and mutagenesis studies and results of MD simulations suggest that NH₃ and H₂O pass through the three monomeric urea channels in UT-B.

Aquaporins are membrane channels that facilitate the flow of water across biological membranes. Two conserved regions are central for selective function: the dual asparagine-proline-alanine (NPA) aquaporin signature motif and the aromatic and arginine selectivity filter (SF). Here, we present the crystal structure of a yeast aquaporin at 0.88 angstrom resolution. We visualize the H-bond donor interactions of the NPA motif's asparagine residues to passing water molecules; observe a polarized water-water H-bond configuration within the channel; assign the tautomeric states of the SF histidine and arginine residues; and observe four SF water positions too closely spaced to be simultaneously occupied. Strongly correlated movements break the connectivity of SF waters to other water molecules within the channel and prevent proton transport via a Grotthuss mechanism.

Cytochrome P450 3A4 (CYP3A4) is the most abundant membrane-associated isoform of the P450 family in humans and is responsible for biotransformation of more than 50% of drugs metabolized in the body. Despite the large number of crystallographic structures available for CYP3A4, no structural information for its membrane-bound state at an atomic level is available. In order to characterize binding, depth of insertion, membrane orientation, and lipid interactions of CYP3A4, we have employed a combined experimental and simulation approach in this study. Taking advantage of a novel membrane representation, highly mobile membrane mimetic (HMMM), with enhanced lipid mobility and dynamics, we have been able to capture spontaneous binding and insertion of the globular domain of the enzyme into the membrane in multiple independent, unbiased simulations. Despite different initial orientations and positions of the protein in solution, all the simulations converged into the same membrane-bound configuration with regard to both the depth of membrane insertion and the orientation of the enzyme on the surface of the membrane. In tandem, linear dichroism measurements performed on CYP3A4 bound to Nanodisc membranes were used to characterize the orientation of the enzyme in its membrane-bound form experimentally. The heme tilt angles measured experimentally are in close agreement with those calculated for the membrane-bound structures resulted from the simulations, thereby verifying the validity of the developed model. Membrane binding of the globular domain in CYP3A4, which appears to be independent of the presence of the transmembrane helix of the full-length enzyme, significantly reshapes the protein at the membrane interface, causing conformational changes relevant to access tunnels leading to the active site of the enzyme.

We present a novel free-energy calculation method that constructively integrates two distinct classes of nonequilibrium sampling techniques, namely, driven (e.g., steered molecular dynamics) and adaptive-bias (e.g., metadynamics) methods. By employing nonequilibrium work relations, we design a biasing protocol with an explicitly time- and history-dependent bias that uses on-the-fly work measurements to gradually flatten the free-energy surface. The asymptotic convergence of the method is discussed, and several relations are derived for free-energy reconstruction and error estimation. Isomerization reaction of an atomistic polyproline peptide model is used to numerically illustrate the superior efficiency and faster convergence of the method compared with its adaptive-bias and driven components in isolation.

Membrane transporters rely on highly coordinated structural transitions between major conformational states for their function, to prevent simultaneous access of the substrate binding site to both sides of the membrane--a mode of operation known as the alternating access model. Although this mechanism successfully accounts for the efficient exchange of the primary substrate across the membrane, accruing evidence on significant water transport and even uncoupled ion transport mediated by transporters has challenged the concept of perfect mechanical coupling and coordination of the gating mechanism in transporters, which might be expected from the alternating access model. Here, we present a large set of extended equilibrium molecular dynamics simulations performed on several classes of membrane transporters in different conformational states, to test the presence of the phenomenon in diverse transporter classes and to investigate the underlying molecular mechanism of water transport through membrane transporters. The simulations reveal spontaneous formation of transient water-conducting (channel-like) states allowing passive water diffusion through the lumen of the transporters. These channel-like states are permeable to water but occluded to substrate, thereby not hindering the uphill transport of the primary substrate, i.e., the alternating access model remains applicable to the substrate. The rise of such water-conducting states during the large-scale structural transitions of the transporter protein is indicative of imperfections in the coordinated closing and opening motions of the cytoplasmic and extracellular gates. We propose that the observed water-conducting states likely represent a universal phenomenon in membrane transporters, which is consistent with their reliance on large-scale motion for function.

Computational modeling and molecular simulation techniques have become an integral part of modern molecular research. Various areas of molecular sciences continue to benefit from, indeed rely on, the unparalleled spatial and temporal resolutions offered by these technologies, to provide a more complete picture of the molecular problems at hand. Because of the continuous development of more efficient algorithms harvesting ever-expanding computational resources, and the emergence of more advanced and novel theories and methodologies, the scope of computational studies has expanded significantly over the past decade, now including much larger molecular systems and far more complex molecular phenomena. Among the various computer modeling techniques, the application of molecular dynamics (MD) simulation and related techniques has particularly drawn attention in biomolecular research, because of the ability of the method to describe the dynamical nature of the molecular systems and thereby to provide a more realistic representation, which is often needed for understanding fundamental molecular properties. The method has proven to be remarkably successful in capturing molecular events and structural transitions highly relevant to the function and/or physicochemical properties of biomolecular systems. Herein, after a brief introduction to the method of MD, we use a number of membrane transport proteins studied in our laboratory as examples to showcase the scope and applicability of the method and its power in characterizing molecular motions of various magnitudes and time scales that are involved in the function of this important class of membrane proteins.

Membrane transporters facilitate active transport of their specific substrates, often against their electrochemical gradients across the membrane, through coupling the process to various sources of cellular energy, for example, ATP binding and hydrolysis in primary transporters, and pre-established electrochemical gradient of molecular species other than the substrate in the case of secondary transporters. In order to provide efficient energy-coupling mechanisms, membrane transporters have evolved into molecular machines in which stepwise binding, translocation, and transformation of various molecular species are closely coupled to protein conformational changes that take the transporter from one functional state to another during the transport cycle. Furthermore, in order to prevent the formation of leaky states and to be able to pump the substrate against its electrochemical gradient, all membrane transporters use the widely-accepted "alternating access mechanism," which ensures that the substrate is only accessible from one side of the membrane at a given time, but relies on complex and usually global protein conformational changes that differ for each family of membrane transporters. Describing the protein conformational changes of different natures and magnitudes is therefore at the heart of mechanistic studies of membrane transporters. Here, using a number of membrane transporters from diverse families, we present common protocols used in setting up and performing molecular dynamics simulations of membrane transporters and in analyzing the results, in order to characterize relevant motions of the system. The emphasis will be on highlighting how optimal design of molecular dynamics simulations combined with mechanistically oriented analysis can shed light onto key functionally relevant protein conformational changes in this family of membrane proteins.

As an adaptation to infrequent access to water, terrestrial mammals produce urine that is hyperosmotic to plasma. To prevent osmotic diuresis by the large quantity of urea generated by protein catabolism, the kidney epithelia contain facilitative urea transporters (UTs) that allow rapid equilibration between the urinary space and the hyperosmotic interstitium. Here we report the first X-ray crystal structure of a mammalian UT, UT-B, at a resolution of 2.36 Å. UT-B is a homotrimer and each protomer contains a urea conduction pore with a narrow selectivity filter. Structural analyses and molecular dynamics simulations showed that the selectivity filter has two urea binding sites separated by an approximately 5.0 kcal/mol energy barrier. Functional studies showed that the rate of urea conduction in UT-B is increased by hypoosmotic stress, and that the site of osmoregulation coincides with the location of the energy barrier.

Characterizing atomic details of membrane binding of peripheral membrane proteins by molecular dynamics (MD) has been significantly hindered by the slow dynamics of membrane reorganization associated with the phenomena. To expedite lateral diffusion of lipid molecules without sacrificing the atomic details of such interactions, we have developed a novel membrane representation, to our knowledge, termed the highly mobile membrane-mimetic (HMMM) model to study binding and insertion of various molecular species into the membrane. The HMMM model takes advantage of an organic solvent layer to represent the hydrophobic core of the membrane and short-tailed phospholipids for the headgroup region. We demonstrate that using these components, bilayer structures are formed spontaneously and rapidly, regardless of the initial position and orientation of the lipids. In the HMMM membrane, lipid molecules exhibit one to two orders of magnitude enhancement in lateral diffusion. At the same time, the membrane atomic density profile of the headgroup region produced by the HMMM model is essentially identical to those obtained for full-membrane models, indicating the faithful representation of the membrane surface by the model. We demonstrate the efficiency of the model in capturing spontaneous binding and insertion of peripheral proteins by using the membrane anchor (γ-carboxyglutamic-acid-rich domain; GLA domain) of human coagulation factor VII as a test model. Achieving full insertion of the GLA domain consistently in 10 independent unbiased simulations within short simulation times clearly indicates the robustness of the HMMM model in capturing membrane association of peripheral proteins very efficiently and reproducibly. The HMMM model will provide significant improvements to the current all-atom models by accelerating lipid dynamics to examine protein-membrane interactions more efficiently., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Blood clotting is triggered when the plasma serine protease factor VIIa binds to the cell-surface protein, tissue factor (TF); the resulting TF:FVIIa complex activates factors IX (FIX) and X (FX) by limited proteolysis. FVIIa, FIX and FX all bind reversibly to membranes via their gamma-carboxyglutamate-rich (GLA) domains, while TF is an integral membrane protein. Removing these proteases from the membrane surface is known to render them thousands of times less active, although the mechanisms by which blood clotting proteins bind to membranes-and the contributions of membranes to catalysis-remain very incompletely understood. Our recent and ongoing studies use a combination of nanoscale membrane bilayers (Nanodiscs), solid-state NMR and all-atom molecular dynamics (MD) simulations, enabling detailed insights into how GLA domains bind to phospholipid bilayers and how specific phospholipids enhance the catalytic activity of the TF:FVIIa complex., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Neisseria are obligate human pathogens causing bacterial meningitis, septicaemia and gonorrhoea. Neisseria require iron for survival and can extract it directly from human transferrin for transport across the outer membrane. The transport system consists of TbpA, an integral outer membrane protein, and TbpB, a co-receptor attached to the cell surface; both proteins are potentially important vaccine and therapeutic targets. Two key questions driving Neisseria research are how human transferrin is specifically targeted, and how the bacteria liberate iron from transferrin at neutral pH. To address these questions, we solved crystal structures of the TbpA-transferrin complex and of the corresponding co-receptor TbpB. We characterized the TbpB-transferrin complex by small-angle X-ray scattering and the TbpA-TbpB-transferrin complex by electron microscopy. Our studies provide a rational basis for the specificity of TbpA for human transferrin, show how TbpA promotes iron release from transferrin, and elucidate how TbpB facilitates this process.

Employing molecular dynamics (MD) simulations, the pathway and mechanism of substrate unbinding from the inward-facing state of the Na(+)-coupled galactose transporter, vSGLT, have been investigated. During a 200-ns equilibrium simulation, repeated spontaneous unbinding events of the substrate from its binding site have been observed. In contrast to the previously proposed gating role of a tyrosine residue (Y263), the unbinding mechanism captured in the present equilibrium simulation does not rely on the displacement and/or rotation of this side chain. Rather, the unbinding involves an initial lateral displacement of the substrate out of the binding site which allows the substrate to completely emerge from the region covered by the side chain of Y263 without any noticeable conformational changes of the latter. Starting with the snapshots taken from this equilibrium simulation with the substrate outside the binding site, steered MD (SMD) simulations were then used to probe the translocation of the substrate along the remaining of the release pathway within the protein's lumen and to characterize the nature of protein-substrate interactions involved in the process. Combining the results of the equilibrium and SMD simulations, we provide a description of the full translocation pathway for the substrate release from the binding site into the cytoplasm. Residues E68, N142, T431, and N267 facilitate the initial substrate's displacement out of the binding site, while the translocation of the substrate along the remainder of the exit pathway formed between TM6 and TM8 is facilitated by H-bond interactions between the substrate and a series of conserved, polar residues (Y138, N267, R273, S365, S368, N371, S372, and T375). The observed molecular events indicate that no gating is required for the release of the substrate from the crystallographically captured structure of the inward-facing state of SGLT, suggesting that this conformation might represent an open, rather than occluded, state of the transporter. This article is part of a Special Issue entitled: Membrane protein structure and function., (Copyright © 2011 Elsevier B.V. All rights reserved.)

Voltage sensor domains (VSD) are transmembrane proteins that respond to changes in membrane voltage and modulate the activity of ion channels, enzymes, or in the case of proton channels allow permeation of protons across the cell membrane. VSDs consist of four transmembrane segments, S1-S4, forming an antiparallel helical bundle. The S4 segment contains several positively charged residues, mainly arginines, located at every third position along the helix. In the voltage-gated Shaker K(+) channel, the mutation of the first arginine of S4 to a smaller uncharged amino acid allows permeation of cations through the VSD. These currents, known as ω-currents, pass through the VSD and are distinct from K(+) currents passing through the main ion conduction pore. Here we report molecular dynamics simulations of the ω-current in the resting-state conformation for Kv1.2 and for four of its mutants. The four tested mutants exhibit various degrees of conductivity for K(+) and Cl(-) ions, with a slight selectivity for K(+) over Cl(-). Analysis of the ion permeation pathway, in the case of a highly conductive mutant, reveals a negatively charged constriction region near the center of the membrane that might act as a selectivity filter to prevent permeation of anions through the pore. The residues R1 in S4 and E1 in S2 are located at the narrowest region of the ω-pore for the resting state conformation of the VSD, in agreement with experiments showing that the largest increase in current is produced by the double mutation E1D and R1S., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Sarco(endo)plasmic reticulum Ca(2+)-ATPase (SERCA) transports two Ca(2+) ions across the membrane of the sarco(endo)plasmic reticulum against the concentration gradient, harvesting the required energy by hydrolyzing one ATP molecule during each transport cycle. Although SERCA is one of the best structurally characterized membrane transporters, it is still largely unknown how the transported Ca(2+) ions reach their transmembrane binding sites in SERCA from the cytoplasmic side. Here, we performed extended all-atom molecular dynamics simulations of SERCA. The calculated electrostatic potential of the protein reveals a putative mechanism by which cations may be attracted to and bind to the Ca(2+)-free state of the transporter. Additional molecular dynamics simulations performed on a Ca(2+)-bound state of SERCA reveal a water-filled pathway that may be used by the Ca(2+) ions to reach their buried binding sites from the cytoplasm. Finally, several residues that are involved in attracting and guiding the cations toward the possible entry channel are identified. The results point to a single Ca(2+) entry site close to the kinked part of the first transmembrane helix, in a region loaded with negatively charged residues. From this point, a water pathway outlines a putative Ca(2+) translocation pathway toward the transmembrane ion-binding sites., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Autotransporters are secreted proteins produced by pathogenic Gram-negative bacteria. They consist of a membrane-embedded β-domain and an extracellular passenger domain that is sometimes cleaved and released from the cell surface. We solved the structures of three noncleavable mutants of the autotransporter EspP to examine how it promotes asparagine cyclization to cleave its passenger. We found that cyclization is facilitated by multiple factors. The active-site asparagine is sterically constrained to conformations favorable for cyclization, while electrostatic interactions correctly orient the carboxamide group for nucleophilic attack. During molecular dynamics simulations, water molecules were observed to enter the active site and to form hydrogen bonds favorable for increasing the nucleophilicity of the active-site asparagine. When the activated asparagine attacks its main-chain carbonyl carbon, the resulting oxyanion is stabilized by a protonated glutamate. Upon cleavage, this proton could be transferred to the leaving amine group, helping overcome a significant energy barrier. Together, these findings provide insight into factors important for asparagine cyclization, a mechanism broadly used for protein cleavage., (Published by Elsevier Ltd.)

Basic architecture of ABC transporters includes two transmembrane domains (TMDs) and two nucleotide-binding domains (NBDs). Although the transport process takes place in the TMDs, which provide the substrate translocation pathway across the cell membrane and control its accessibility between the two sides of the membrane, the energy required for the process is provided by conformational changes induced in the NBDs by binding and hydrolysis of ATP. Nucleotide-dependent conformational changes in the NBDs, therefore, need to be coupled to structural changes in the TMDs. Using molecular dynamics simulations, we have investigated the structural elements involved in the conformational coupling between the NBDs and the TMDs in the Escherichia coli maltose transporter, an ABC importer for which an intact structure is available both in inward-facing and outward-facing conformations. The prevailing model of coupling is primarily based on a single structural motif, known as the coupling helices, as the main structural element for the NBD-TMD coupling. Surprisingly, we find that in the absence of the NBDs the coupling helices can be conformationally decoupled from the rest of the TMDs, despite their covalent connection. That is, the structural integrity of the coupling helices and their tight coupling to the core of the TMDs rely on the contacts provided by the NBDs. Based on the conformational and dynamical analysis of the simulation trajectories, we propose that the core coupling elements in the maltose transporter involve contributions from several structural motifs located at the NBD-TMD interface, namely, the EAA loops from the TMDs, and the Q-loop and the ENI motifs from the NBDs. These three structural motifs in small ABC importers show a high degree of correlation in motion and mediate the necessary conformational coupling between the core of TMDs and the helical subdomains of NBDs. A comprehensive analysis of the structurally known ABC transporters shows a high degree of conservation of the identified 3-motif coupling elements only in the subfamily of small ABC importers, suggesting a distinct mode of NBD-TMD coupling from the other two major ABC transporter folds, namely large ABC importers and ABC exporters., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Many regulatory processes in biology involve reversible association of proteins with membranes. Clotting proteins bind to phosphatidylserine (PS) on cell surfaces, but a clear picture of this interaction has yet to emerge. We present a novel explanation for membrane binding by GLA domains of clotting proteins, supported by biochemical studies, solid-state NMR analyses, and molecular dynamics simulations. The model invokes a single "phospho-L-serine-specific" interaction and multiple "phosphate-specific" interactions. In the latter, the phosphates in phospholipids interact with tightly bound Ca(2+) in GLA domains. We show that phospholipids with any headgroup other than choline strongly synergize with PS to enhance factor X activation. We propose that phosphatidylcholine and sphingomyelin (the major external phospholipids of healthy cells) are anticoagulant primarily because their bulky choline headgroups sterically hinder access to their phosphates. Following cell damage or activation, exposed PS and phosphatidylethanolamine collaborate to bind GLA domains by providing phospho-L-serine-specific and phosphate-specific interactions, respectively.

Membrane binding of peripheral proteins, mediated by specialized anchoring domains, is a crucial step for their biological function. Computational studies of membrane insertion, however, have proven challenging and largely inaccessible, due to the time scales required for the complete description of the process, mainly caused by the slow diffusion of the lipid molecules composing the membrane. Furthermore, in many cases, the nature of the membrane "anchor", i.e., the part of the protein that inserts into the membrane, is also unknown. Here, we address some of these issues by developing and employing a simplified representation of the membrane by a biphasic solvent model which we demonstrate can be used efficiently to capture and describe the process of hydrophobic insertion of membrane anchoring domains in all-atom molecular dynamics simulations. Applying the model, we have studied the insertion of the anchoring domain of a coagulation protein (the GLA domain of human protein C), starting from multiple initial configurations varying with regard to the initial orientation and height of the protein with respect to the membrane. In addition to efficiently and consistently identifying the "keel" region as the hydrophobic membrane anchor, within a few nanoseconds each configuration simulated showed a convergent height (2.20 ± 1.04 Å) and angle with respect to the interface normal (23.37 ± 12.48°). We demonstrate that the model can produce the same results as those obtained from a full representation of a membrane, in terms of both the depth of penetration and the orientation of the protein in the final membrane-bound form with an order of magnitude decrease in the required computational time compared to previous models, allowing for a more exhaustive search for the correct membrane-bound configuration.

Membranes play key regulatory roles in biological processes, with bilayer composition exerting marked effects on binding affinities and catalytic activities of a number of membrane-associated proteins. In particular, proteins involved in diverse processes such as vesicle fusion, intracellular signaling cascades, and blood coagulation interact specifically with anionic lipids such as phosphatidylserine (PS) in the presence of Ca(2+) ions. While Ca(2+) is suspected to induce PS clustering in mixed phospholipid bilayers, the detailed structural effects of this ion on anionic lipids are not established. In this study, combining magic angle spinning (MAS) solid-state NMR (SSNMR) measurements of isotopically labeled serine headgroups in mixed lipid bilayers with molecular dynamics (MD) simulations of PS lipid bilayers in the presence of different counterions, we provide site-resolved insights into the effects of Ca(2+) on the structure and dynamics of lipid bilayers. Ca(2+)-induced conformational changes of PS in mixed bilayers are observed in both liposomes and Nanodiscs, a nanoscale membrane mimetic of bilayer patches. Site-resolved multidimensional correlation SSNMR spectra of bilayers containing (13)C,(15)N-labeled PS demonstrate that Ca(2+) ions promote two major PS headgroup conformations, which are well resolved in two-dimensional (13)C-(13)C, (15)N-(13)C, and (31)P-(13)C spectra. The results of MD simulations performed on PS lipid bilayers in the presence or absence of Ca(2+) provide an atomic view of the conformational effects underlying the observed spectra.

Conformational changes of proteins are involved in all aspects of protein function in biology. Almost all classes of proteins respond to changes in their environment, ligand binding, and interaction with other proteins and regulatory agents through undergoing conformational changes of various degrees and magnitudes. Membrane channels and transporters are the major classes of proteins that are responsible for mediating efficient and selective transport of materials across the cellular membrane. Similar to other proteins, they take advantage of conformational changes to make transitions between various functional states. In channels, large-scale conformational changes are mostly involved in the process of "gating", i.e., opening and closing of the pore of the channel protein in response to various signals. In transporters, conformational changes constitute various steps of the conduction process, and, thus, are more closely integrated in the transport process. Owing to significant progress in developing highly efficient parallel algorithms in molecular dynamics simulations and increased computational resources, and combined with the availability of high-resolution, atomic structures of membrane proteins, we are in an unprecedented position to use computer simulation and modeling methodologies to investigate the mechanism of function of membrane channels and transporters. While the entire transport cycle is still out of reach of current methodologies, many steps involved in the function of transport proteins have been characterized with molecular dynamics simulations. Here, we present several examples of such studies from our laboratory, in which functionally relevant conformational changes of membrane channels and transporters have been characterized using extended simulations.

The interaction of cytochrome c with ubiquinol-cytochrome c oxidoreductase (bc₁ complex) has been studied for >30 years, yet many aspects remain unclear or controversial. We report the first molecular dynamic simulations of the cyt c-bc₁ complex interaction. Contrary to the results of crystallographic studies, our results show that there are multiple dynamic hydrogen bonds and salt bridges in the cyt c-c₁ interface. These include most of the basic cyt c residues previously implicated in chemical modification studies. We suggest that the static nature of x-ray structures can obscure the quantitative significance of electrostatic interactions between highly mobile residues. This provides a clear resolution of the discrepancy between the structural data and functional studies. It also suggests a general need to consider dynamic interactions of charged residues in protein-protein interfaces. In addition, a novel structural change in cyt c is reported, involving residues 21-25, which may be responsible for cyt c destabilization upon binding. We also propose a mechanism of interaction between cyt c₁ monomers responsible for limiting the binding of cyt c to only one molecule per bc₁ dimer by altering the affinity of the cytochrome c binding site on the second cyt c₁ monomer., (Copyright © 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The transport cycle in the glutamate transporter (GlT) is catalyzed by the cotransport of three Na(+) ions. However, the positions of only two of these ions (Na1 and Na2 sites) along with the substrate have been captured in the crystal structures reported for both the outward-facing and the inward-facing states of Glt(ph). Characterizing the third ion binding site (Na3) is necessary for structure-function studies attempting to investigate the mechanism of transport in GlTs at an atomic level, particularly for the determination of the sequence of the binding events during the transport cycle. In this study, we report a series of molecular dynamics simulations performed on various bound states of Glt(ph) (the apo state, as well as in the presence of Na(+), the substrate, or both), which have been used to identify a putative Na3 site. The calculated trajectories have been used to determine the water accessibility of potential ion-binding residues in the protein, as a prerequisite for their ion binding. Combined with conformational analysis of the key regions in the protein in different bound states and several additional independent simulations in which a Na(+) ion was randomly introduced to the interior of the transporter, we have been able to characterize a putative Na3 site and propose a plausible binding sequence for the substrate and the three Na(+) ions to the transporter during the extracellular half of the transport cycle. The proposed Na3 site is formed by a set of highly conserved residues, namely, Asp(312), Thr(92), and Asn(310), along with a water molecule. Simulation of a fully bound state, including the substrate and the three Na(+) ions, reveals a stable structure--showing closer agreement to the crystal structure when compared to previous models lacking an ion in the putative Na3 site. The proposed sequence of binding events is in agreement with recent experimental models suggesting that two Na(+) ions bind before the substrate, and one after that. Our results, however, provide additional information about the sites involved in these binding events., (Copyright 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The leucine transporter (LeuT) has recently commanded exceptional attention due mainly to two distinctions; it provides the only crystal structures available for a protein homologous to the pharmacologically relevant neurotransmitter: sodium symporters (NSS), and, it exhibits a hallmark 5-TM inverted repeat ("LeuT-fold"), a fold recently discovered to also exist in several secondary transporter families, underscoring its general role in transporter function. Constructing the transport cycle of "LeuT-fold" transporters requires detailed structural and dynamic descriptions of the outward-facing (OF) and inward-facing (IF) states, as well as the intermediate states. To this end, we have modeled the structurally unknown IF state of LeuT, based on the known crystal structures of the OF state of LeuT and the IF state of vSGLT, a "LeuT-fold" transporter. The detailed methodology developed for the study combines structure-based alignment, threading, targeted MD and equilibrium MD, and can be applied to other proteins. The resulting IF-state models maintain the secondary structural features of LeuT. Water penetration and solvent accessibility calculations show that TM1, TM3, TM6 and TM8 line the substrate binding/unbinding pathway with TM10 and its pseudosymmetric partner, TM5, participating in the extracellular and intracellular halves of the lumen, respectively. We report conformational hotspots where notable changes in interactions occur between the IF and OF states. We observe Na2 exiting the LeuT-substrate- complex in the IF state, mainly due to TM1 bending. Inducing a transition in only one of the two pseudosymmetric domains, while allowing the second to respond dynamically, is found to be sufficient to induce the formation of the IF state. We also propose that TM2 and TM7 may be facilitators of TM1 and TM6 motion. Thus, this study not only presents a novel modeling methodology applied to obtain the IF state of LeuT, but also describes structural elements involved in a possibly general transport mechanism in transporters adopting the "LeuT-fold".

Transmembrane exchange of materials is a fundamental process in biology. Molecular dynamics provides a powerful method to investigate in great detail various aspects of the phenomenon, particularly the permeation of small uncharged molecules, which continues to pose a challenge to experimental studies. We will discuss some of the recent simulation studies investigating the role of lipid-mediated and protein-mediated mechanisms in permeation of water and gas molecules across the membrane.

The atomic models of the Kv1.2 potassium channel in the active and resting state, originally presented elsewhere, are here refined using molecular dynamics simulations in an explicit membrane-solvent environment. With a minor adjustment of the orientation of the first arginine along the S4 segment, the total gating charge of the channel determined from >0.5 mus of molecular dynamics simulation is approximately 12-12.7 e, in good accord with experimental estimates for the Shaker potassium channel, indicating that the final models offer a realistic depiction of voltage-gating. In the resting state of Kv1.2, the S4 segment in the voltage-sensing domain (VSD) spontaneously converts into a 3(10) helix over a stretch of 10 residues. The 3(10) helical conformation orients the gating arginines on S4 toward a water-filled crevice within the VSD and allows salt-bridge interactions with negatively charged residues along S2 and S3. Free energy calculations of the fractional transmembrane potential, acting upon key charged residues of the VSD, reveals that the applied field varies rapidly over a narrow region of 10-15 A corresponding to the outer leaflet of the bilayer. The focused field allows the transfer of a large gating charge without translocation of S4 across the membrane., (Copyright 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Fibronectins (FNs) are a major component of the extracellular matrix (ECM), and provide important binding sites for a variety of ligands outside and on the surface of the cell. Similar to other ECM proteins, FNs are consistently subject to mechanical stress in the ECM. Therefore, it is important to study their structure and binding properties under mechanical stress and understand how their binding and mechanical properties might affect each other. Although certain FN modules have been extensively investigated, no simulation studies have been reported for the FN type I (Fn1) domains, despite their prominent role in binding of various protein modules to FN polymers in the ECM. Using equilibrium and steered molecular dynamics simulations, we have studied mechanical properties of Fn1 modules in the presence or the absence of a specific FN-binding peptide (FnBP). We have also investigated how the binding of the FnBP peptide to Fn1 might be affected by tensile force. Despite the presence of disulfide bonds within individual Fn1 modules that are presumed to prevent their extension, it is found that significant internal structural changes within individual modules are induced by the forces applied in our simulations. These internal structural changes result in significant variations in the accessibility of different residues of the Fn1 modules, which affect their exposure, and, thus, the binding properties of the Fn1 modules. Binding of the FnBP appears to reduce the flexibility of the linker region connecting individual Fn1 modules (exhibited in the form of reduced fluctuation and motion of the linker region), both with regard to bending and stretching motions, and hence stabilizes the inter-domain configuration under force. Under large tensile forces, the FnBP peptide unbinds from Fn1. The results suggest that Fn1 modules in FN polymers do contribute to the overall extension caused by force-induced stretching of the polymer in the ECM, and that binding properties of Fn1 modules can be affected by mechanically induced internal protein conformational changes in spite of the presence of disulfide bonds which were presumed to completely abolish the capacity of Fn1 modules to undergo extension in response to external forces.

Posttranslational modifications of the RelA subunit of NF-kappaB, including acetylation and methylation, play a key role in controlling the strength and duration of its nuclear activity. Whether these modifications are functionally linked is largely unknown. Here, we show that the acetylation of lysine 310 of RelA impairs the Set9-mediated methylation of lysines 314 and 315, which is important for the ubiquitination and degradation of chromatin-associated RelA. Abolishing the acetylation of lysine 310 either by the deacetylase SIRT1 or by mutating lysine 310 to arginine enhances methylation. Conversely, enhancing the acetylation of lysine 310 by depleting SIRT1 or by replacing lysine 310 with acetyl-mimetic glutamine inhibits methylation, thereby decreasing ubiquitination, prolonging the stability of chromatin-associated RelA, and enhancing the transcriptional activity of NF-kappaB. The acetylation of lysine 310 of RelA interferes with its interaction with Set9. Based on structural modeling of the SET domain of Set9 with RelA, we propose that the positive charge of lysine 310 is critical for the binding of RelA to a negatively charged "exosite" within the SET domain of Set9. Together, these findings demonstrate for the first time an interplay between RelA acetylation and methylation and also provide a novel mechanism for the regulation of lysine methylation by acetylation.

Most steps of the blood clotting cascade require the assembly of a serine protease with its specific regulatory protein on a suitable phospholipid bilayer. Unfortunately, the molecular details of how blood clotting proteins bind to membrane surfaces remain poorly understood, owing to a dearth of techniques for studying protein-membrane interactions at high resolution. Our laboratories are tackling this question using a combination of approaches, including nanoscale membrane bilayers, solid-state NMR, and large-scale molecular dynamics simulations. These studies are now providing structural insights at atomic resolution into clotting protein-membrane interactions., (Copyright (c) 2010 Elsevier Ltd. All rights reserved.)

We here present a detailed study of the ligand-receptor interactions between single and triple-helical strands of collagen and the α2A domain of integrin (α2A), providing valuable new insights into the mechanisms and dynamics of collagen-integrin binding at a sub-molecular level. The occurrence of single and triple-helical strands of the collagen fragments was scrutinized with atom force microscopy (AFM) techniques. Strong interactions of the triple-stranded fragments comparable to those of collagen can only be detected for the 42mer triple-helical collagen-like peptide under study (which contains 42 amino acid residues per strand) by solid phase assays as well as by surface plasmon resonance (SPR) measurements. However, changes in NMR signals during titration and characteristic saturation transfer difference (STD) NMR signals are also detectable when α2A is added to a solution of the 21mer single-stranded collagen fragment. Molecular dynamics (MD) simulations employing different sets of force field parameters were applied to study the interaction between triple-helical or single-stranded collagen fragments with α2A. It is remarkable that even single-stranded collagen fragments can form various complexes with α2A showing significant differences in the complex stability with identical ligands. The results of MD simulations are in agreement with the signal alterations in our NMR experiments, which are indicative of the formation of weak complexes between single-stranded collagen and α2A in solution. These results provide useful information concerning possible interactions of α2A with small collagen fragments that are of relevance to the design of novel therapeutic A-domain inhibitors.

Phosphorylation of adenine nucleotide translocator 1 (ANT1) at residue Y194, which is part of the aromatic ladder located within the lumen of the carrier, critically regulates mitochondrial metabolism. Recent data support the concept that members of the Src family of nonreceptor tyrosine kinases are constitutively present in mitochondria and key to regulation of mitochondrial function. Herein, we demonstrate that site mutations of ANT1 (Y190-->F190, Y194-->F194) mimicking dephosphorylation of the aromatic ladder resulted in loss of oxidative growth and ADP/ATP exchange activity in respiration-incompetent yeast expressing mutant chimeric yN-hANT1. ANT1 is phosphorylated at Y194 by the Src family kinase members Src and Lck, and increased phosphorylation is tightly linked to reduced cell injury in preconditioned protected vs. unprotected cardiac mitochondria. Molecular dynamics simulations find the overall structure of the phosphorylated ANT1 stable, but with an increased steric flexibility in the region of the aromatic ladder, matrix loop m2, and four helix-linking regions. Combined with an analysis of the putative cytosolic salt bridge network, we reason that the effect of phosphorylation on transport is likely due to an accelerated transition between the main two conformational states (c<-->m) of the carrier during the transport cycle. Since "aromatic signatures" are typical for other mitochondrial carrier proteins with important biological functions, our results may be more general and applicable to these carriers.

Glycerol 3-phosphate transporter (GlpT) mediates the import of glycerol 3-phosphate (G3P) using the gradient of inorganic phosphate (P(i)). To study the process and mechanism of substrate binding and to investigate the protein's initial response, we performed equilibrium simulations of wild-type GlpT and several of its mutant forms in membranes in the presence of all physiologically relevant substrates (P(i)(-), P(i)(2-), G3P(-), and G3P(2-)). The simulations capture spontaneous substrate binding of GlpT, driven by the positive electrostatic potential of the lumen. K80 is found to act as a "hook" making the first encounter with the substrate and guiding it toward the binding site, where it binds tightly to R45, a key binding site residue that acts as a "fork" holding the substrate. R269 establishes no direct contact with the substrate during the simulations, a surprising behavior given its structural pseudosymmetry to R45. In all substrate-bound systems, partial closing of the cytoplasmic half of GlpT was observed. The substrate appears to stabilize the partially occluded state, as in the two apo simulations either no closing was observed or the protein reverted to its open form toward the end of the simulation, whereas in all substrate-bound systems, a stable partially closed state was produced. Along with the modulation of the periplasmic salt bridge network, these substrate-induced events destabilize the periplasmic half while inducing a closure in the cytoplasmic half, thus capturing the early stages of the proposed rocker-switch mechanism in GlpT.

Involvement of aquaporins in gas conduction across the membrane and the physiological significance of this process have attracted marked attention from both experimental and theoretical studies. Previous work demonstrated that AQP1 is permeable to both CO(2) and O(2). Here we employ various simulation techniques to examine the permeability of the brain aquaporin AQP4 to NO and O(2) and to describe energetics and pathways associated with these phenomena. The energy barrier to NO and O(2) permeation through AQP4 central pore is found to be only approximately 3 kcal mol(-1). The results suggest that the central pore of AQP4, similar to that of AQP1, can indeed conduct gas molecules. Interestingly, despite a longer and narrower central pore, AQP4 appears to provide an energetically more favorable permeation pathway for gas molecules than AQP1, mainly due to the different orientation of its charged residues near the pore entrance. Although the low barrier against gas permeation through AQP4 indicates that it can participate in gas conduction across the cellular membrane, physiological relevance of the phenomenon remains to be established experimentally, particularly since pure lipid bilayers appear to present a more favorable pathway for gas conduction across the membrane. With an energy well of -1.8 kcal mol(-1), the central pore of AQP4 may also act as a reservoir for NO molecules to accumulate in the membrane.

We present a new method to obtain information about the conformational dynamics of membrane-proteins using solid-state NMR experiments of oriented samples. By measuring the orientation-dependent (1)H-(15)N dipole-dipole coupling, (15)N anisotropic chemical shift, and (2)H quadrupole coupling parameters for a single residue, it is possible to obtain information about the local dynamics of each residue in the protein. This may be interpreted on an individual basis or through models extended to study conformational motion of membrane-protein segments. The method is demonstrated for the antimicrobial peptaibol alamethicin for which combined analysis of anisotropic interactions for the Aib(8) residue provides detailed information about helix-tilt angle, wobbling, and oscillatory rotation around the helix axis in the membrane bound state. This information is in very good agreement with coarse-grained MD simulations of the peptide in lipid bilayers.