Revisiting Hydrophobic Mismatch with Free Energy Simulation Studies of Transmembrane Helix Tilt and Rotation

This is the publisher's version. Copyright 2010 by Elsevier.Protein-lipid interaction and bilayer regulation of membrane protein functions are largely controlled by the hydrophobic match between the transmembrane (TM) domain of membrane proteins and the surrounding lipid bilayer. To systematically characterize responses of a TM helix and lipid adaptations to a hydrophobic mismatch, we have performed a total of 5.8-μs umbrella sampling simulations and calculated the potentials of mean force (PMFs) as a function of TM helix tilt angle under various mismatch conditions. Single-pass TM peptides called WALPn (n = 16, 19, 23, and 27) were used in two lipid bilayers with different hydrophobic thicknesses to consider hydrophobic mismatch caused by either the TM length or the bilayer thickness. In addition, different flanking residues, such as alanine, lysine, and arginine, instead of tryptophan in WALP23 were used to examine their influence. The PMFs, their decomposition, and trajectory analysis demonstrate that 1), tilting of a single-pass TM helix is the major response to a hydrophobic mismatch; 2), TM helix tilting up to ∼10° is inherent due to the intrinsic entropic contribution arising from helix precession around the membrane normal even under a negative mismatch; 3), the favorable helix-lipid interaction provides additional driving forces for TM helix tilting under a positive mismatch; 4), the minimum-PMF tilt angle is generally located where there is the hydrophobic match and little lipid perturbation; 5), TM helix rotation is dependent on the specific helix-lipid interaction; and 6), anchoring residues at the hydrophilic/hydrophobic interface can be an important determinant of TM helix orientation

Quantification of Drive-Response Relationships Between Residues During Protein Folding

Mutual correlation and cooperativity are commonly used to describe residue-residue interactions in protein folding/function. However, these metrics do not provide any information on the causality relationships between residues. Such drive-response relationships are poorly studied in protein folding/function and difficult to measure experimentally due to technical limitations. In this study, using the information theory transfer entropy (TE) that provides a direct measurement of causality between two times series, we have quantified the drive-response relationships between residues in the folding/unfolding processes of four small proteins generated by molecular dynamics simulations. Instead of using a time-averaged single TE value, the time-dependent TE is measured with the Q-scores based on residue-residue contacts and with the statistical significance analysis along the folding/unfolding processes. The TE analysis is able to identify the driving and responding residues that are different from the highly correlated residues revealed by the mutual information analysis. In general, the driving residues have more regular secondary structures, are more buried, and show greater effects on the protein stability as well as folding and unfolding rates. In addition, the dominant driving and responding residues from the TE analysis on the whole trajectory agree with those on a single folding event, demonstrating that the drive-response relationships are preserved in the non-equilibrium process. Our study provides detailed insights into the protein folding process and has potential applications in protein engineering and interpretation of time-dependent residue-based experimental observables for protein function

Brownian dynamics simulations of ions channels: A general treatment of electrostatic reaction fields for molecular pores of arbitrary geometry

Copyright 2001 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in the Journal of Chemical Physics and may be found at http://dx.doi.org/10.1063/1.1390507.A general method has been developed to include the electrostatic reaction field in Brownian dynamics (BD) simulations of ions diffusing through complex molecular channels of arbitrary geometry. Assuming that the solvent is represented as a featureless continuum dielectric medium, a multipolar basis-set expansion is developed to express the reaction field Green’s function. A reaction field matrix, which provides the coupling between generalized multipoles, is calculated only once and stored before the BD simulations. The electrostatic energy and forces are calculated at each time step by updating the generalized multipole moments. The method is closely related to the generalized solvent boundary potential [Im et al., J. Chem. Phys. 114, 2924 (2001)] which was recently developed to include the influence of distant atoms on a small region part of a large macromolecular system in molecular dynamics simulations. It is shown that the basis-set expansion is accurate and computationally inexpensive for three simple models such as a spherical ionic system, an impermeable membrane system, and a cylindrical pore system as well as a realistic system such as OmpF porin with all atomic details. The influence of the static field and the reaction field on the ion distribution and conductance in the OmpF channel is studied and discussed

Glycan fragment database: a database of PDB-based glycan 3D structures

The glycan fragment database (GFDB), freely available at http://www.glycanstructure.org, is a database of the glycosidic torsion angles derived from the glycan structures in the Protein Data Bank (PDB). Analogous to protein structure, the structure of an oligosaccharide chain in a glycoprotein, referred to as a glycan, can be characterized by the torsion angles of glycosidic linkages between relatively rigid carbohydrate monomeric units. Knowledge of accessible conformations of biologically relevant glycans is essential in understanding their biological roles. The GFDB provides an intuitive glycan sequence search tool that allows the user to search complex glycan structures. After a glycan search is complete, each glycosidic torsion angle distribution is displayed in terms of the exact match and the fragment match. The exact match results are from the PDB entries that contain the glycan sequence identical to the query sequence. The fragment match results are from the entries with the glycan sequence whose substructure (fragment) or entire sequence is matched to the query sequence, such that the fragment results implicitly include the influences from the nearby carbohydrate residues. In addition, clustering analysis based on the torsion angle distribution can be performed to obtain the representative structures among the searched glycan structures.University of Kansas General Research Fund [2301388-003]; Kanas-COBRE NIH P20 [GM103420]; NSF [MCB-0918374]; TeraGrid/XSEDE resources [TG-MCB070009]. Funding for open access charge: NSF [MCB-0918374]

NMR Observable-Based Structure Refinement of DAP12-NKG2C Activating Immunoreceptor Complex in Explicit Membranes

This is the published version. Copyright 2012 by Elsevier.NMR observables, such as NOE-based distance measurements, are increasingly being used to characterize membrane protein structures. However, challenges in membrane protein NMR studies often yield a relatively small number of such restraints that can create ambiguities in defining critical side chain-side chain interactions. In the recent solution NMR structure of the DAP12-NKG2C immunoreceptor transmembrane helix complex, five functionally required interfacial residues (two Asps and two Thrs in the DAP12 dimer and one Lys in NKG2C) display a surprising arrangement in which one Asp side chain faces the membrane hydrophobic core. To explore whether these side-chain interactions are energetically optimal, we used the published distance restraints for molecular dynamics simulations in explicit micelles and bilayers. The structures refined by this protocol are globally similar to the published structure, but the side chains of those five residues form a stable network of salt bridges and hydrogen bonds, leaving the Asp side chain shielded from the hydrophobic core, which is also consistent with available experimental observations. Moreover, the simulations show similar short-range interactions between the transmembrane complex and lipid/detergent molecules in micelles and bilayers, respectively. This study illustrates the efficacy of NMR membrane protein structure refinements in explicit membrane systems

Transmembrane Helix Orientation and Dynamics: Insights from Ensemble Dynamics with Solid-State NMR Observables

This is the publisher's version. Copyright 2011 by Elsevier.As the major component of membrane proteins, transmembrane helices embedded in anisotropic bilayer environments adopt preferential orientations that are characteristic or related to their functional states. Recent developments in solid-state nuclear magnetic resonance (SSNMR) spectroscopy have made it possible to measure NMR observables that can be used to determine such orientations in a native bilayer environment. A quasistatic single conformer model is frequently used to interpret the SSNMR observables, but important motional information can be missing or misinterpreted in the model. In this work, we have investigated the orientation of the single-pass transmembrane domain of viral protein ”u“ (VpuTM) from HIV-1 by determining an ensemble of structures using multiple conformer models based on the SSNMR ensemble dynamics technique. The resulting structure ensemble shows significantly larger orientational fluctuations while the ensemble-averaged orientation is compatible with the orientation based on the quasistatic model. This observation is further corroborated by comparison with the VpuTM orientation from comparative molecular dynamics simulations in explicit bilayer membranes. SSNMR ensemble dynamics not only reveals the importance of transmembrane helix dynamics in interpretation of SSNMR observables, but also provides a means to simultaneously extract both transmembrane helix orientation and dynamics information from the SSNMR measurements

NMR Observable-Based Structure Refinement of DAP12-NKG2C Activating Immunoreceptor Complex in Explicit Membranes

This is the published version. Copyright 2012 by Elsevier.NMR observables, such as NOE-based distance measurements, are increasingly being used to characterize membrane protein structures. However, challenges in membrane protein NMR studies often yield a relatively small number of such restraints that can create ambiguities in defining critical side chain-side chain interactions. In the recent solution NMR structure of the DAP12-NKG2C immunoreceptor transmembrane helix complex, five functionally required interfacial residues (two Asps and two Thrs in the DAP12 dimer and one Lys in NKG2C) display a surprising arrangement in which one Asp side chain faces the membrane hydrophobic core. To explore whether these side-chain interactions are energetically optimal, we used the published distance restraints for molecular dynamics simulations in explicit micelles and bilayers. The structures refined by this protocol are globally similar to the published structure, but the side chains of those five residues form a stable network of salt bridges and hydrogen bonds, leaving the Asp side chain shielded from the hydrophobic core, which is also consistent with available experimental observations. Moreover, the simulations show similar short-range interactions between the transmembrane complex and lipid/detergent molecules in micelles and bilayers, respectively. This study illustrates the efficacy of NMR membrane protein structure refinements in explicit membrane systems

Ligand Binding Site Detection b Local Structure Alignment and Its Performance Complementarity

Accurate determination of potential ligand binding sites (BS) is a key step for protein function characterization and structure-based drug design. Despite promising results of template-based BS prediction methods using global structure alignment (GSA), there is a room to improve the performance by properly incorporating local structure alignment (LSA) because BS are local structures and often similar for proteins with dissimilar global folds. We present a template-based ligand BS prediction method using G-LoSA, our LSA tool. A large benchmark set validation shows that G-LoSA predicts drug-like ligands’ positions in single-chain protein targets more precisely than TM-align, a GSA-based method, while the overall success rate of TM-align is better. G-LoSA is particularly efficient for accurate detection of local structures conserved across proteins with diverse global topologies. Recognizing the performance complementarity of G-LoSA to TM-align and a non-template geometry-based method, fpocket, a robust consensus scoring method, CMCS-BSP (Complementary Methods and Consensus Scoring for ligand Binding Site Prediction), is developed and shows improvement on prediction accuracy. The G-LoSA source code is freely available at http://im.bioinformatics.ku.edu/GLoSA