Transmembrane helix orientation and dynamics: insights from ensemble dynamics with solid-state NMR observables

Hydrophobic matching of HIV-1 Vpu transmembrane helix-helix interactions is optimized for subcellular location

The HIV-1 accessory protein Vpu mediates the downregulation of several host cell proteins, an activity that is critical for viral replication in vivo. As the first step in directing cell-surface proteins to internal cellular compartments, and in many cases degradation, Vpu binds a subset of its target proteins through their transmembrane domains. Each of the known targets of Vpu are synthesized in the ER, and must traverse the different membrane environments found along the secretory pathway, thus it is important to consider how membrane composition might influence the interactions between Vpu and its targets. We have used Förster resonance energy transfer (FRET) to measure the oligomerization of Vpu with the transmembrane domains of target proteins in model membranes of varying lipid composition. Our data show that both lipid bilayer thickness and acyl chain order can significantly influence monomer-oligomer equilibria within the Vpu-target system. Changes in oligomerization levels were found to be non-specific with no single Vpu-target interaction being favored under any condition. Our analysis of the influence of the membrane environment on the strength of helix-helix interactions between Vpu and its targets in vitro suggests that the strength of Vpu-target interactions in vivo will be partially dependent on the membrane environment found in specific membrane compartments.

The conjugation of indolicidin to polyethylenimine for enhanced gene delivery with reduced cytotoxicity

The hydrophobic domains of conjugated peptides can stabilize their insertion into the cell membrane to promote transportation.

What controls open-pore and residual currents in the first sensing zone of alpha-hemolysin nanopore? Combined experimental and theoretical study

The electrophoretic transport of single-stranded DNA through biological nanopores such as alpha-hemolysin (αHL) is a promising and cost-effective technology with the potential to revolutionize genomics. The rational design of pores with the controlled polymer translocation rates and high contrast between different nucleotides could improve significantly nanopore sequencing applications. Here, we apply a combination of theoretical and experimental methods in an attempt to elucidate several selective modifications in the pore which were proposed to be central for the effective discrimination between purines and pyrimidines. Our nanopore test set includes the wild type αHL and six mutants (E111N/M113X/K147N) in which the cross-section and chemical functionality of the first constriction zone of the pore are modified. Electrophysiological recordings were combined with all-atom Molecular Dynamics simulations (MD) and a recently developed Brownian Dynamics (BROMOC) protocol to investigate residual ion currents and pore-DNA interactions for two homo-polymers e.g. poly(dA)40 or poly(dC)40 blocking the pore. The calculated residual currents and contrast in the poly(dA)40/poly(dC)40 blocked pore are in qualitative agreement with the experimental recordings. We showed that a simple structural metric allows rationalization of key elements in the emergent contrast between purines and pyrimidines in the modified αHL mutants. The shape of the pore and its capacity for hydrogen bonding to a translocated polynucleotide are two essential parameters for contrast optimization. To further probe the impact of these two factors in the ssDNA sensing, we eliminated the effect of the primary constriction using serine substitutions (i.e. E111S/M113S/T145S/K147S) and increased the hydrophobic volume of the central residue in the secondary constriction (L135I). This pore modification sharply increased the contrast between Adenine (A) and Cytosine (C).

Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes

Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys(40), residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.

Structural Refinement of Proteins by Restrained Molecular Dynamics Simulations with Non-interacting Molecular Fragments

The knowledge of multiple conformational states is a prerequisite to understand the function of membrane transport proteins. Unfortunately, the determination of detailed atomic structures for all these functionally important conformational states with conventional high-resolution approaches is often difficult and unsuccessful. In some cases, biophysical and biochemical approaches can provide important complementary structural information that can be exploited with the help of advanced computational methods to derive structural models of specific conformational states. In particular, functional and spectroscopic measurements in combination with site-directed mutations constitute one important source of information to obtain these mixed-resolution structural models. A very common problem with this strategy, however, is the difficulty to simultaneously integrate all the information from multiple independent experiments involving different mutations or chemical labels to derive a unique structural model consistent with the data. To resolve this issue, a novel restrained molecular dynamics structural refinement method is developed to simultaneously incorporate multiple experimentally determined constraints (e.g., engineered metal bridges or spin-labels), each treated as an individual molecular fragment with all atomic details. The internal structure of each of the molecular fragments is treated realistically, while there is no interaction between different molecular fragments to avoid unphysical steric clashes. The information from all the molecular fragments is exploited simultaneously to constrain the backbone to refine a three-dimensional model of the conformational state of the protein. The method is illustrated by refining the structure of the voltage-sensing domain (VSD) of the Kv1.2 potassium channel in the resting state and by exploring the distance histograms between spin-labels attached to T4 lysozyme. The resulting VSD structures are in good agreement with the consensus model of the resting state VSD and the spin-spin distance histograms from ESR/DEER experiments on T4 lysozyme are accurately reproduced.

Knowledge of multiple conformational states of membrane transport proteins is a prerequisite to understand their function. However, the determination of atomic structures for all these states with conventional high-resolution approaches can be very challenging due to inherent difficulties in high yield purification of functional membrane transport proteins. Various complementary structural information of proteins in their native states can be obtained by a variety of biophysical and biochemical methods with site-directed mutations. Here, a novel restrained molecular dynamics structural refinement method is developed to help derive a structural model that is consistent with experimental data by incorporating all the experimental constraints simultaneously through the use of non-interacting all-atom molecular fragments. The method can be easily and effectively extended to incorporate many kinds of structural constraints from a variety of biophysical and biochemical experiments, and should be very useful in generating and refining models of proteins in specific functional states.

Structural determination of virus protein U from HIV-1 by NMR in membrane environments

Virus protein U (Vpu) from HIV-1, a small membrane protein composed of a transmembrane helical domain and two α-helices in an amphipathic cytoplasmic domain, down modulates several cellular proteins, including CD4, BST-2/CD317/tetherin, NTB-A, and CCR7. The interactions of Vpu with these proteins interfere with the immune system and enhance the release of newly synthesized virus particles. It is essential to characterize the structure and dynamics of Vpu in order to understand the mechanisms of the protein-protein interactions, and potentially to discover antiviral drugs. In this article, we describe investigations of the cytoplasmic domain of Vpu as well as full-length Vpu by NMR spectroscopy. These studies are complementary to earlier analysis of the transmembrane domain of Vpu. The results suggest that the two helices in the cytoplasmic domain form a U-shape. The length of the inter-helical loop in the cytoplasmic domain and the orientation of the third helix vary with the lipid composition, which demonstrate that the C-terminal helix is relatively flexible, providing accessibility for interaction partners.

Preferred orientations of phosphoinositides in bilayers and their implications in protein recognition mechanisms

Phosphoinositides (PIPs), phosphorylated derivatives of phosphatidylinositol (PI), are essential regulatory lipids involved in various cellular processes, including signal transduction, membrane trafficking, and cytoskeletal remodeling. To gain insight into the protein-PIPs recognition process, it is necessary to study the inositol ring orientation (with respect to the membrane) of PIPs with different phosphorylation states. In this study, 8 PIPs (3 PIP, 2 PIP2, and 3 PIP3) with different phosphorylation and protonation sites have been separately simulated in two mixed bilayers (one with 20% phosphatidylserine (PS) lipids and another with PS lipids switched to phosphatidylcholine (PC) lipids), which roughly correspond to yeast membranes. Uniformity of the bilayer properties including hydrophobic thickness, acyl chain order parameters, and heavy atom density profiles is observed in both PS-contained and PC-enriched membranes due to the same hydrophobic core composition. The relationship between the inositol ring orientation (tilt and rotation angles) and its solvent-accessible surface area indicates that the orientation is mainly determined by its solvation energy. Different PIPs exhibit a clear preference in the inositol ring rotation angle. Surprisingly, a larger proportion of PIPs inositol rings stay closer to the surface of PS-contained membranes compared to PC-enriched ones. Such a difference is rationalized with the formation of more hydrogen bonds between the PS/PI headgroups and the PIPs inositol rings in PS-contained membranes. This hydrogen bond network could be functionally important; thus, the present results can potentially add important and detailed features into the existing protein-PIPs recognition mechanism.

NMR-based simulation studies of Pf1 coat protein in explicit membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.

ST-analyzer: a web-based user interface for simulation trajectory analysis

Molecular dynamics (MD) simulation has become one of the key tools to obtain deeper insights into biological systems using various levels of descriptions such as all-atom, united-atom, and coarse-grained models. Recent advances in computing resources and MD programs have significantly accelerated the simulation time and thus increased the amount of trajectory data. Although many laboratories routinely perform MD simulations, analyzing MD trajectories is still time consuming and often a difficult task. ST-analyzer, http://im.bioinformatics.ku.edu/st-analyzer, is a standalone graphical user interface (GUI) toolset to perform various trajectory analyses. ST-analyzer has several outstanding features compared to other existing analysis tools: (i) handling various formats of trajectory files from MD programs, such as CHARMM, NAMD, GROMACS, and Amber, (ii) intuitive web-based GUI environment--minimizing administrative load and reducing burdens on the user from adapting new software environments, (iii) platform independent design--working with any existing operating system, (iv) easy integration into job queuing systems--providing options of batch processing either on the cluster or in an interactive mode, and (v) providing independence between foreground GUI and background modules--making it easier to add personal modules or to recycle/integrate pre-existing scripts utilizing other analysis tools. The current ST-analyzer contains nine main analysis modules that together contain 18 options, including density profile, lipid deuterium order parameters, surface area per lipid, and membrane hydrophobic thickness. This article introduces ST-analyzer with its design, implementation, and features, and also illustrates practical analysis of lipid bilayer simulations.

Combining experiments and simulations using the maximum entropy principle

A key component of computational biology is to compare the results of computer modelling with experimental measurements. Despite substantial progress in the models and algorithms used in many areas of computational biology, such comparisons sometimes reveal that the computations are not in quantitative agreement with experimental data. The principle of maximum entropy is a general procedure for constructing probability distributions in the light of new data, making it a natural tool in cases when an initial model provides results that are at odds with experiments. The number of maximum entropy applications in our field has grown steadily in recent years, in areas as diverse as sequence analysis, structural modelling, and neurobiology. In this Perspectives article, we give a broad introduction to the method, in an attempt to encourage its further adoption. The general procedure is explained in the context of a simple example, after which we proceed with a real-world application in the field of molecular simulations, where the maximum entropy procedure has recently provided new insight. Given the limited accuracy of force fields, macromolecular simulations sometimes produce results that are at not in complete and quantitative accordance with experiments. A common solution to this problem is to explicitly ensure agreement between the two by perturbing the potential energy function towards the experimental data. So far, a general consensus for how such perturbations should be implemented has been lacking. Three very recent papers have explored this problem using the maximum entropy approach, providing both new theoretical and practical insights to the problem. We highlight each of these contributions in turn and conclude with a discussion on remaining challenges.

Ensemble-based interpretations of NMR structural data to describe protein internal dynamics

NMR spectroscopy is the leading technique to characterize protein internal dynamics at the atomic level and on multiple time scales. However, the structural interpretation of the observables obtained by various measurements is not always straightforward and in many cases dynamics-related parameters are only used to “decorate” static structural models without offering explicit description of conformational heterogeneity. To overcome such limitations, several computational techniques have been developed to generate ensemble-based representations of protein structure and dynamics with the use of NMR-derived data. An important common aspect of the methods is that NMR observables and derived parameters are interpreted as properties of the ensemble instead of individual conformers. The resulting ensembles reflect the experimentally determined internal mobility of proteins at a given time scale and can be used to understand the role of internal motions in biological processes at atomic detail. In this review we provide an overview of the calculation methods currently available and examples of biological insights obtained by the ensemble-based models of the proteins investigated.

On the statistical equivalence of restrained-ensemble simulations with the maximum entropy method

An issue of general interest in computer simulations is to incorporate information from experiments into a structural model. An important caveat in pursuing this goal is to avoid corrupting the resulting model with spurious and arbitrary biases. While the problem of biasing thermodynamic ensembles can be formulated rigorously using the maximum entropy method introduced by Jaynes, the approach can be cumbersome in practical applications with the need to determine multiple unknown coefficients iteratively. A popular alternative strategy to incorporate the information from experiments is to rely on restrained-ensemble molecular dynamics simulations. However, the fundamental validity of this computational strategy remains in question. Here, it is demonstrated that the statistical distribution produced by restrained-ensemble simulations is formally consistent with the maximum entropy method of Jaynes. This clarifies the underlying conditions under which restrained-ensemble simulations will yield results that are consistent with the maximum entropy method.

Structural refinement from restrained-ensemble simulations based on EPR/DEER data: application to T4 lysozyme

DEER (double electron-electron resonance) is a powerful pulsed ESR (electron spin resonance) technique allowing the determination of distance histograms between pairs of nitroxide spin-labels linked to a protein in a native-like solution environment. However, exploiting the huge amount of information provided by ESR/DEER histograms to refine structural models is extremely challenging. In this study, a restrained ensemble (RE) molecular dynamics (MD) simulation methodology is developed to address this issue. In RE simulation, the spin-spin distance distribution histograms calculated from a multiple-copy MD simulation are enforced, via a global ensemble-based energy restraint, to match those obtained from ESR/DEER experiments. The RE simulation is applied to 51 ESR/DEER distance histogram data from spin-labels inserted at 37 different positions in T4 lysozyme (T4L). The rotamer population distribution along the five dihedral angles connecting the nitroxide ring to the protein backbone is determined and shown to be consistent with available information from X-ray crystallography. For the purpose of structural refinement, the concept of a simplified nitroxide dummy spin-label is designed and parametrized on the basis of these all-atom RE simulations with explicit solvent. It is demonstrated that RE simulations with the dummy nitroxide spin-labels imposing the ESR/DEER experimental distance distribution data are able to systematically correct and refine a series of distorted T4L structures, while simple harmonic distance restraints are unsuccessful. This computationally efficient approach allows experimental restraints from DEER experiments to be incorporated into RE simulations for efficient structural refinement.

Restrained-ensemble molecular dynamics simulations based on distance histograms from double electron-electron resonance spectroscopy

DEER (double electron-electron resonance) spectroscopy is a powerful pulsed ESR (electron spin resonance) technique allowing the determination of spin-spin distance histograms between site-directed nitroxide label sites on a protein in their native environment. However, incorporating ESR/DEER data in structural refinement is challenging because the information from the large number of distance histograms is complex and highly coupled. Here, a novel restrained-ensemble molecular dynamics simulation method is developed to incorporate the information from multiple ESR/DEER distance histograms simultaneously. Illustrative tests on three coupled spin-labels inserted in T4 lysozyme show that the method efficiently imposes the experimental distance distribution in this system. Different rotameric states of the χ1 and χ2 dihedrals in the spin-labels are also explored by restrained ensemble simulations. Using this method, it is hoped that experimental restraints from ESR/DEER experiments can be used to refine structural properties of biological systems.

Dynamic Heterogeneous Dielectric Generalized Born (DHDGB): An implicit membrane model with a dynamically varying bilayer thickness

An extension to the heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is presented to allow dynamic membrane deformations in response to membrane-inserted biomolecules during molecular dynamic simulations. The flexible membrane is implemented through additional degrees of freedom that represent the membrane deformation at the contact points of a membrane-inserted solute with the membrane. The extra degrees of freedom determine the dielectric and non-polar solvation free energy profiles that are used to obtain the solvation free energy in the presence of the membrane and are used to calculate membrane deformation free energies according to an elastic membrane model. With the dynamic HDGB (DHDGB) model the membrane is able to deform in response to the insertion of charged molecules thereby avoiding the overestimation of insertion free energies with static membrane models. The DHDGB model also allows the membrane to respond to the insertion of membrane-spanning solutes with hydrophobic mismatch. The model is tested with the membrane insertion of amino acid side chain analogs, arginine-containing helices, the WALP23 peptide, and the gramicidin A channel.

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4-5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations-except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.

Solid-state NMR ensemble dynamics as a mediator between experiment and simulation

Solid-state NMR (SSNMR) is a powerful technique to describe the orientations of membrane proteins and peptides in their native membrane bilayer environments. The deuterium ((2)H) quadrupolar splitting (DQS), one of the SSNMR observables, has been used to characterize the orientations of various single-pass transmembrane (TM) helices using a semistatic rigid-body model such as the geometric analysis of labeled alanine (GALA) method. However, dynamic information of these TM helices, which could be related to important biological function, can be missing or misinterpreted with the semistatic model. We have investigated the orientation of WALP23 in an implicit membrane of dimyristoylglycerophosphocholine by determining an ensemble of structures using multiple conformer models with a DQS restraint potential. When a single conformer is used, the resulting helix orientation (tilt angle (τ) of 5.6 ± 3.2° and rotation angle (ρ) of 141.8 ± 40.6°) is similar to that determined by the GALA method. However, as the number of conformers is increased, the tilt angles of WALP23 ensemble structures become larger (26.9 ± 6.7°), which agrees well with previous molecular dynamics simulation results. In addition, the ensemble structure distribution shows excellent agreement with the two-dimensional free energy surface as a function of WALP23's τ and ρ. These results demonstrate that SSNMR ensemble dynamics provides a means to extract orientational and dynamic information of TM helices from their SSNMR observables and to explain the discrepancy between molecular dynamics simulation and GALA-based interpretation of DQS data.

Orientational and motional narrowing of solid-state NMR lineshapes of uniaxially aligned membrane proteins

A unified theory for the NMR line shapes of aligned membrane proteins arising from uniaxial disorder (mosaic spread) and global rotational diffusion about the director axis is presented. A superoperator formalism allows one to take into account the effects of continuous radiofrequency irradiation and frequency offsets in the presence of dynamics. A general method based on the Stochastic Liouville Equation makes it possible to bridge the static and dynamic limits in a single model. Simulations of solid-state NMR spectra are performed for a uniform α helix by considering orientational disorder and diffusion of the helix as a whole relative to the alignment axis. The motional narrowing of the resonance lines is highly inhomogeneous and can be used as an additional angular restraint in structure calculations. Experimental solid-state NMR spectra of Pf1 coat protein support the conclusions of the theory for two limiting cases. The static disorder dominates the (15)N NMR spectra of Pf1 aligned on a phage, while fast uniaxial diffusion provides a line narrowing mechanism for the Pf1 protein reconstituted in magnetically aligned bicelles.

An ensemble dynamics approach to decipher solid-state NMR observables of membrane proteins

Solid-state NMR (SSNMR) is an invaluable tool for determining orientations of membrane proteins and peptides in lipid bilayers. Such orientational descriptions provide essential information about membrane protein functions. However, when a semi-static single conformer model is used to interpret various SSNMR observables, important dynamics information can be missing, and, sometimes, even orientational information can be misinterpreted. In addition, over the last decade, molecular dynamics (MD) simulation and semi-static SSNMR interpretation have shown certain levels of discrepancies in terms of transmembrane helix orientation and dynamics. Dynamic fitting models have recently been proposed to resolve these discrepancies by taking into account transmembrane helix whole body motions using additional parameters. As an alternative approach, we have developed SSNMR ensemble dynamics (SSNMR-ED) using multiple conformer models, which generates an ensemble of structures that satisfies the experimental observables without any fitting parameters. In this review, various computational methods for determining transmembrane helix orientations are discussed, and the distributions of VpuTM (from HIV-1) and WALP23 (a synthetic peptide) orientations from SSNMR-ED simulations are compared with those from MD simulations and semi-static/dynamic fitting models. Such comparisons illustrate that SSNMR-ED can be used as a general means to extract both membrane protein structure and dynamics from the SSNMR measurements. This article is part of a Special Issue entitled: Membrane protein structure and function.