Correlated motions are a fundamental property of ß-sheets.

Correlated motions in proteins can mediate fundamental biochemical processes such as signal transduction and allostery. The mechanisms that underlie these processes remain largely unknown due mainly to limitations in their direct detection. Here, based on a detailed analysis of protein structures deposited in the protein data bank, as well as on state-of-the art molecular simulations, we provide general evidence for the transfer of structural information by correlated backbone motions, mediated by hydrogen bonds, across ß-sheets. We also show that the observed local and long-range correlated motions are mediated by the collective motions of ß-sheets and investigate their role in large-scale conformational changes. Correlated motions represent a fundamental property of ß-sheets that contributes to protein function.

Hydrophobic mismatch of mobile transmembrane helices: Merging theory and experiments.

Hydrophobic mismatch still represents a puzzle for transmembrane peptides, despite the apparent simplicity of this concept and its demonstrated validity in natural membranes. Using a wealth of available experimental ((2))H NMR data, we provide here a comprehensive explanation of the orientation and dynamics of model peptides in lipid bilayers, which shows how they can adapt to membranes of different thickness. The orientational adjustment of transmembrane α-helices can be understood as the result of a competition between the thermodynamically unfavorable lipid repacking associated with peptide tilting and the optimization of peptide/membrane hydrophobic coupling. In the positive mismatch regime (long-peptide/thin-membrane) the helices adapt mainly via changing their tilt angle, as expected from simple geometrical predictions. However, the adaptation mechanism varies with the peptide sequence in the flanking regions, suggesting additional effects that modulate hydrophobic coupling. These originate from re-adjustments of the peptide hydrophobic length and they depend on the hydrophobicity of the flanking region, the strength of interfacial anchoring, the structural flexibility of anchoring side-chains and the presence of alternative anchoring residues.

Understanding biomolecular motion, recognition, and allostery by use of conformational ensembles.

We review the role conformational ensembles can play in the analysis of biomolecular dynamics, molecular recognition, and allostery. We introduce currently available methods for generating ensembles of biomolecules and illustrate their application with relevant examples from the literature. We show how, for binding, conformational ensembles provide a way of distinguishing the competing models of induced fit and conformational selection. For allostery we review the classic models and show how conformational ensembles can play a role in unravelling the intricate pathways of communication that enable allostery to occur. Finally, we discuss the limitations of conformational ensembles and highlight some potential applications for the future.

Weak long-range correlated motions in a surface patch of ubiquitin involved in molecular recognition.

Long-range correlated motions in proteins are candidate mechanisms for processes that require information transfer across protein structures, such as allostery and signal transduction. However, the observation of backbone correlations between distant residues has remained elusive, and only local correlations have been revealed using residual dipolar couplings measured by NMR spectroscopy. In this work, we experimentally identified and characterized collective motions spanning four ß-strands separated by up to 15 Å in ubiquitin. The observed correlations link molecular recognition sites and result from concerted conformational changes that are in part mediated by the hydrogen-bonding network.

A lipocentric view of peptide-induced pores.

Although lipid membranes serve as effective sealing barriers for the passage of most polar solutes, nonmediated leakage is not completely improbable. A high activation energy normally keeps unassisted bilayer permeation at a very low frequency, but lipids are able to self-organize as pores even in peptide-free and protein-free membranes. The probability of leakage phenomena increases under conditions such as phase coexistence, external stress or perturbation associated to binding of nonlipidic molecules. Here, we argue that pore formation can be viewed as an intrinsic property of lipid bilayers, with strong similarities in the structure and mechanism between pores formed with participation of peptides, lipidic pores induced by different types of stress, and spontaneous transient bilayer defects driven by thermal fluctuations. Within such a lipocentric framework, amphipathic peptides are best described as pore-inducing rather than pore-forming elements. Active peptides bound to membranes can be understood as a source of internal surface tension which facilitates pore formation by diminishing the high activation energy barrier. This first or immediate action of the peptide has some resemblance to catalysis. However, the presence of membrane-active peptides has the additional effect of displacing the equilibrium towards the pore-open state, which is then maintained over long times, and reducing the size of initial individual pores. Thus, pore-inducing peptides, regardless of their sequence and oligomeric organization, can be assigned a double role of increasing the probability of pore formation in membranes to high levels as well as stabilizing these pores after they appear.

Pores formed by Baxα5 relax to a smaller size and keep at equilibrium.

Pores made by amphipathic cationic peptides (e.g., antimicrobials and fragments of pore-forming proteins) are typically studied by examining the kinetics of vesicle leakage after peptide addition or obtaining structural measurements in reconstituted peptide-lipid systems. In the first case, the pores have been considered transient phenomena that allow the relaxation of the peptide-membrane system. In the second, they correspond to equilibrium structures at minimum free energy. Here we reconcile both approaches by investigating the pore activity of the α5 fragment from the proapoptotic protein Bax (Baxα5) before and after equilibrium of peptide/vesicle complexes. Quenching assays on suspensions of large unilamellar vesicles suggest that in the presence of Baxα5, the vesicles maintain a leaky state for hours under equilibrium conditions. We proved and analyzed stable pores on single giant unilamellar vesicles (GUVs) in detail by monitoring the entrance of dyes added at different times after incubation with the peptide. When the GUVs came in contact with Baxα5, leakage started stochastically, was delayed for various periods of time, and in the majority of cases proceeded rapidly to completion. After hours in the presence of the peptide, the same individual GUVs that refilled completely at first instance maintained a porated state, which could be observed in subsequent leak-in events for serially added dyes. However, these long-term pores were smaller in size than the initial equilibration pores. Stable pores were also detected in GUVs made in the presence of Baxα5. The latter pores can be considered equilibrium states and may correspond to structures measured previously in bilayer stacks. Although pore formation may occur as a kinetic process, equilibrium pores may also be functionally relevant structures, especially in highly regulated systems such as the apoptotic mitochondrial pores induced by Bax.

Role of membrane lipids for the activity of pore forming peptides and proteins.

Bilayer lipids, far from being passive elements, have multiple roles in polypeptide-dependent pore formation. Lipids participate at all stages of the formation of pores by providing the binding site for proteins and peptides, conditioning their active structure and modulating the molecular reorganization of the membrane complex. Such general functions of lipids superimpose to other particular roles, from electrostatic and curvature effects to more specific actions in cases like cholesterol, sphingolipids or cardiolipin. Pores are natural phenomena in lipid membranes. Driven by membrane fluctuations and packing defects, transient water pores are related to spontaneous lipid flip-flop and non-assisted ion permeation. In the absence ofproteins or peptides, these are rare short living events, with properties dependent on the lipid composition of the membrane. Their frequency increases under conditions of internal membrane disturbance of the lipid packing, like in the presence of membrane-bound proteins or peptides. These latter molecules, in fact, form dynamic supramolecular assemblies together with the lipids and transmembrane pores are one of the possible structures of the complex. Active peptides and proteins can thus be considered inducers or enhancers of pores which increase their probability and lifetime by modifying the thermodynamic membrane balance. This includes destabilizing the membrane lamellar structure, lowering the activation energy for pore formation and stabilizing the open pore structure.

Refinement of ensembles describing unstructured proteins using NMR residual dipolar couplings.

Residual dipolar couplings (RDCs) are unique probes of the structural and dynamical properties of biomolecules on the sub-millisecond time scale that can be used as restraints in ensemble molecular dynamics simulations to study the relationship between macromolecular motion and biological function. To date, however, this powerful strategy is applicable only to molecules that do not undergo shape changes on the time scale sampled by RDCs, thus preventing the study of key biological macromolecules such as multidomain and unstructured proteins. In this work, we circumvent this limitation by using an algorithm that explicitly computes the individual alignment tensors of the different ensemble members from their coordinates at each step in the simulation. As a first application, we determine an ensemble of conformations that accurately describes the structure and dynamics of chemically denatured ubiquitin. In analogy to dynamic refinement of folded, globular proteins, where simulations are initiated from average structures, we use statistical coil models as starting configuration because they represent the best available descriptions of unstructured proteins. We find that refinement with RDCs causes significant structural corrections and yields an ensemble that is in complete agreement with the measured RDCs and presents transient mid-range inter-residue interactions between strands beta1 and beta2 of the native protein, also observed in other studies based on trans-hydrogen bond (3)J(NC') scalar couplings and paramagnetic relaxation enhancements. Finally, and in spite of the high structural heterogeneity of the refined ensemble, we find that it can be cross-validated against RDCs not used to restrain the simulation. This method increases the range of systems that can be studied using ensemble simulations restrained by RDCs and is likely to yield new insights into how the large-scale motions of macromolecules relate to biological function.

Solid state NMR analysis of peptides in membranes: Influence of dynamics and labeling scheme.

The functional state of a membrane-active peptide is often defined by its conformation, molecular orientation, and its oligomeric state in the lipid bilayer. These "static" structural properties can be routinely studied by solid state NMR using isotope-labeled peptides. In the highly dynamic environment of a liquid crystalline biomembrane, however, the whole-body fluctuations of a peptide are also of paramount importance, although difficult to address and most often ignored. Yet it turns out that disregarding such motional averaging in calculating the molecular alignment from orientational NMR-constraints may give a misleading, if not false picture of the system. Here, we demonstrate that the reliability of a simplified static or an advanced dynamic data analysis depends critically on the choice of isotope labeling scheme used. Two distinctly different scenarios have to be considered. When the labels are placed on the side chains of a helical peptide (such as a CD(3)- or CF(3)-group attached to the C(alpha)C(beta) bond), their nuclear spin interaction tensors are very sensitive to motional averaging. If this effect is not properly accounted for, the helix tilt angle tends to be severely underestimated. At the same time, the analysis of labels in the side chains allows to extract valuable dynamical information about whole-body fluctuations of the peptide helix in the membrane. On the other hand, the alternative labeling scheme where (15)N-labels are accommodated within the peptide backbone, will yield nearly correct helix tilt angles, irrespective as to whether dynamics are taken into account or not.

Stability of asymmetric lipid bilayers assessed by molecular dynamics simulations.

The asymmetric insertion of amphiphiles into biological membranes compromises the balance between the inner and outer monolayers. As a result, area expansion of the receiving leaflet and curvature strain may lead to membrane permeation, shape changes, or membrane fusion events. We have conducted both atomistic and coarse-grained molecular dynamics simulations of dipalmitoyl-phosphatidylcholine (DPPC) bilayers to study the effect of an asymmetric distribution of lipids between the two monolayers on membrane stability. Highly asymmetric lipid bilayers were found to be surprisingly stable within the submicrosecond time span of the simulations. Even the limiting case of a monolayer immersed in water ruptured spontaneously only after at least 20 ns simulation. A thermal shock could destabilize these kinetically trapped states. We also studied mixed systems composed of DPPC and short tail diC(8)PC lipids, showing that the presence of the cone-shaped short tail lipid facilitates the release of tension in the asymmetric systems via formation of a transmembrane pore. Thus, asymmetric area expansion and curvature stress cooperate to yield bilayer disruption. It appears that, although asymmetric area expansion destabilizes the bilayer structure, the activation energy for transmonolayer lipid re-equilibration is increased. Such a large kinetic barrier can be reduced by lipids with positive spontaneous curvature. These effects are important at the onset of bilayer destabilization phenomena, such as lipid pore formation and membrane fusion, and should be considered for the mechanism of induction of such processes by peptides and proteins.

Orientational landscapes of peptides in membranes: prediction of (2)H NMR couplings in a dynamic context.

Unlike soluble proteins, membrane polypeptides face an anisotropic milieu. This imposes restraints on their orientation and provides a reference that makes structure prediction tractable by minimalistic thermodynamic models. Here we use this framework to build orientational distributions of monomeric membrane-bound peptides and to predict their expected solid-state (2)H NMR quadrupolar couplings when labeled at specific side chain positions. Using a complete rigid-body sampling of configurations relative to an implicit lipid membrane, peptide free energy landscapes are calculated. This allows us to obtain probability distributions of the peptide tilt, azimuthal rotation, and depth of membrane insertion. The orientational distributions are broad and originate from an interplay among the three relevant rigid-body degrees of freedom, which allows population of multiple states in shallow free energy minima. Remarkably, only when the orientational distributions are taken into account do we obtain a close correlation between predicted (2)H NMR splittings and values measured in experiments. Such a good correlation is not seen with splittings calculated from single configurations, being either the averaged or the lowest free energy state, showing there are distributions, rather than single structures, that best define the peptide-membrane systems. Moreover, we propose that these distributions contribute to the understanding of the rigid-body dynamics of the system.

Orientation and dynamics of peptides in membranes calculated from 2H-NMR data.

Solid-state (2)H-NMR is routinely used to determine the alignment of membrane-bound peptides. Here we demonstrate that it can also provide a quantitative measure of the fluctuations around the distinct molecular axes. Using several dynamic models with increasing complexity, we reanalyzed published (2)H-NMR data on two representative alpha-helical peptides: 1), the amphiphilic antimicrobial peptide PGLa, which permeabilizes membranes by going from a monomeric surface-bound to a dimeric tilted state and finally inserting as an oligomeric pore; and 2), the hydrophobic WALP23, which is a typical transmembrane segment, although previous analysis had yielded helix tilt angles much smaller than expected from hydrophobic mismatch and molecular dynamics simulations. Their (2)H-NMR data were deconvoluted in terms of the two main helix orientation angles (representing the time-averaged peptide tilt and azimuthal rotation), as well as the amplitudes of fluctuation about the corresponding molecular axes (providing the dynamic picture). The mobility of PGLa is found to be moderate and to correlate well with the respective oligomeric states. WALP23 fluctuates more vigorously, now in better agreement with the molecular dynamics simulations and mismatch predictions. The analysis demonstrates that when (2)H-NMR data are fitted to extract peptide orientation angles, an explicit representation of the peptide rigid-body angular fluctuations should be included.

Influence of whole-body dynamics on 15N PISEMA NMR spectra of membrane proteins: a theoretical analysis.

Membrane proteins and peptides exhibit a preferred orientation in the lipid bilayer while fluctuating in an anisotropic manner. Both the orientation and the dynamics have direct functional implications, but motions are usually not accessible, and structural descriptions are generally static. Using simulated data, we analyze systematically the impact of whole-body motions on the peptide orientations calculated from two-dimensional polarization inversion spin exchange at the magic angle (PISEMA) NMR. Fluctuations are found to have a significant effect on the observed spectra. Nevertheless, wheel-like patterns are still preserved, and it is possible to determine the average peptide tilt and azimuthal rotation angles using simple static models for the spectral fitting. For helical peptides undergoing large-amplitude fluctuations, as in the case of transmembrane monomers, improved fits can be achieved using an explicit dynamics model that includes Gaussian distributions of the orientational parameters. This method allows extracting the amplitudes of fluctuations of the tilt and azimuthal rotation angles. The analysis is further demonstrated by generating first a virtual PISEMA spectrum from a molecular dynamics trajectory of the model peptide, WLP23, in a lipid membrane. That way, the dynamics of the system from which the input spectrum originates is completely known at atomic detail and can thus be directly compared with the dynamic output obtained from the fit. We find that fitting our dynamics model to the polar index slant angles wheel gives an accurate description of the amplitude of underlying motions, together with the average peptide orientation.

The dynamic orientation of membrane-bound peptides: bridging simulations and experiments.

The structural organization in a peptide/membrane supramolecular complex is best described by knowledge of the peptide orientation plus its time-dependent and spatial fluctuations. The static orientation, defined by the peptide tilt and a rotation about its molecular axis, is accessible through a number of spectroscopic methods. However, peptide dynamics, although relevant to understand the functionality of these systems, remains largely unexplored. Here, we describe the orientation and dynamics of Trp-flanked and Lys-flanked hydrophobic peptides in a lipid bilayer from molecular dynamics simulations. A novel view is revealed, where collective nontrivial distributions of time-evolving and ensemble peptide orientations closely represent the systems as studied experimentally. Such global distributions are broad and unveil the existence of orientational states, which depend on the anchoring mode of interfacial residues. We show that this dynamics modulates (2)H quadrupolar splittings and introduces ambiguity in the analysis of NMR data. These findings demonstrate that structural descriptions of peptide/membrane complexes are incomplete, and in cases even imprecise, without knowledge of dynamics.

Self-assembling of peptide/membrane complexes by atomistic molecular dynamics simulations.

Model biological membranes consisting of peptide/lipid-bilayer complexes can nowadays be studied by classical molecular dynamics (MD) simulations at atomic detail. In most cases, the simulation starts with an assumed state of a peptide in a preformed bilayer, from which equilibrium configurations are difficult to obtain due to a relatively slow molecular diffusion. As an alternative, we propose an extension of reported work on the self-organization of unordered lipids into bilayers, consisting of including a peptide molecule in the initial random configuration to obtain a membrane-bound peptide simultaneous to the formation of the lipid bilayer. This strategy takes advantage of the fast reorganization of lipids, among themselves and around the peptide, in an aqueous environment. Model peptides of different hydrophobicity, CH3-CO-W2L18W2-NH2 (WL22) and CH3-CO-W2A18W2-NH2 (WA22), in dipalmitoyl-phosphatidylcholine (DPPC) are used as test cases. In the equilibrium states of the peptide/membrane complexes, achieved in time ranges of 50-100 ns, the two peptides behave as expected from experimental and theoretical studies. The strongly hydrophobic WL22 is inserted in a transmembrane configuration and the marginally apolar, alanine-based WA22 is found in two alternative states: transmembrane inserted or parallel to the membrane plane, embedded close to the bilayer interface, with similar stability. This shows that the spontaneous assembly of peptides and lipids is an unbiased and reliable strategy to produce and study models of equilibrated peptide/lipid complexes of unknown membrane-binding mode and topology.