Explicit treatment of force contribution from alignment tensor using overdetermined linear equations and its application in NMR structure determination

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.

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.