Search

Nuclear magnetic resonance (NMR) spin relaxation is the most informative approach to experimentally probe the internal dynamics of proteins on the picosecond to nanosecond timescale. At the same time, molecular dynamics (MD) simulations of biological macromolecules are steadily improving through better physical models, enhanced sampling algorithms, and increasing computational power, and they provide exquisite information about flexibility and its role in protein stability and molecular interactions. Many examples have shown that MD is now adept in probing protein backbone motion, but improvements are still required towards a quantitative description of the dynamics of side chains, for example probed by the dynamics of methyl groups. Thus far, the comparison of computation with experiment for side chains has primarily focused on the relaxation of 13C and 2H nuclei induced by auto-correlated variation of spin interactions. However, the cross-correlation of 13C-1H dipolar interactions in methyl groups offers an attractive alternative. Here, we establish a methodological framework to extract cross-correlation relaxation parameters of methyl groups in proteins from all-atom MD simulations. To demonstrate the utility of the approach, cross-correlation relaxation rates of ubiquitin are computed from MD simulations performed with the AMBER99SB*-ILDN and CHARMM36 force fields. The simulation results were found to agree well with those obtained by experiment. Moreover, the data obtained with the two force fields are highly consistent.

Molecular simulations of intrinsically disordered proteins (IDPs) are challenging because they require sampling a very large number of relevant conformations, corresponding to a multitude of shallow minima in a flat free energy landscape. However, in the presence of a binding partner, the free energy landscape of an IDP can be dominated by few deep minima. This characteristic imposes high demands on the accuracy of the force field used to describe the molecular interactions. Here, as a model system for an IDP that is unstructured in solution but folds upon binding to a structured interaction partner, the transactivation domain of c-Myb was studied both in the unbound (free) form and when bound to the KIX domain. Six modern biomolecular force fields were systematically tested and compared in terms of their ability to describe the structural ensemble of the IDP. The protein force field/water model combinations included in this study are AMBER ff99SB-disp with its corresponding water model that was derived from TIP4P-D, CHARMM36m with TIP3P, ff15ipq with SPC/Eb, ff99SB*-ILDNP with TIP3P and TIP4P-D, and FB15 with TIP3P-FB water. Comparing the results from REST2-enhanced sampling simulations with experimental CD spectra and secondary chemical shifts reveals that the ff99SB-disp force field can realistically capture the broad and mildly helical structural ensemble of free c-Myb. The structural ensembles yielded by CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 are also mildly helical; however, each of these force fields can be assigned a specific subset of c-Myb residues for which the simulations could not reproduce the experimental secondary chemical shifts. In addition, microsecond-timescale MD simulations of the KIX/c-Myb complex show that most force fields used preserve a stable helix fold of c-Myb in the complex. Still, all force fields predict a KIX/c-Myb complex interface that differs slightly from the structures provided by NMR because several NOE-derived distances between KIX and c-Myb were exceeded in the simulations. Taken together, the ff99SB-disp force field in the first place but also CHARMM36m, ff99SB*-ILDNP together with TIP4P-D water, and FB15 can be suitable choices for future simulation studies of the coupled folding and binding mechanism of the KIX/c-Myb complex and potentially also other IDPs.

Drug discovery, in particular optimization of candidates using medicinal chemistry, is generally guided by structural biology. However, for optimizing binding kinetics, relevant for efficacy and off‐target effects, information on protein motion is important. Herein, we demonstrate for the prototypical textbook example of an allegedly “rigid protein” that substantial active‐site dynamics have generally remained unrecognized, despite thousands of medicinal‐chemistry studies on this model over decades. Comparing cryogenic X‐ray structures, solid‐state NMR on micro‐crystalline protein at room temperature, and solution NMR structure and dynamics, supported by MD simulations, we show that under physiologically relevant conditions the pocket is in fact shaped by pronounced open/close conformational‐exchange dynamics. The study, which is of general significance for pharmacological research, evinces a generic pitfall in drug discovery routines., Human carbonic anhydrase II (hCAII) plays an important role as a rigid model protein in medicinal chemistry. Under physiological conditions in solution, NMR spectroscopy combined with molecular dynamics simulations demonstrates that the active site shows pronounced open/close conformational exchange, undetected in crystallography‐based work.

ATP-binding cassette (ABC) transporters chemomechanically couple ATP binding and hydrolysis to large-scale conformational changes, ultimately leading to substrate translocation across biological membranes. Despite recent progress in the structure determination of substrate-bound ABC exporters, the inherently dynamic mechanism of substrate transport remains unclear at the atomic level. In this work, we capture substrate translocation in the heterodimeric ABC exporter TM287/288 from the hyperthermophilic bacterium Thermotoga maritima using all-atom molecular dynamics (MD) simulations. Unguided multimicrosecond simulations at 375 K show how the drugs daunorubicin and verapamil, which were initially docked into the ABC transporter, get translocated through the exporter by following its large-scale alternating access conformational transitions between an inward-facing (IF) and an outward-facing (OF) conformation. Triggered by the affinity difference due to differential solvation of the binding cavity in the IF and OF conformations, the substrates unbind from the OF transporter and partition into the lipid bilayer. While daunorubicin is stably inserted into the outer leaflet of the bilayer, verapamil dynamically flip flops between the bilayer leaflets, possibly rendering its net transport futile.

Solution-state NMR relaxation experiments are the cornerstone to study internal protein dynamics at an atomic resolution on time scales that are faster than the overall rotational tumbling time τR. Since the motions described by NMR relaxation parameters are connected to thermodynamic quantities like conformational entropies, the question arises how much of the total entropy is contained within this tumbling time. Using all-atom molecular dynamics simulations of the T4 lysozyme, we found that entropy buildup is rather fast for the backbone, such that the majority of the entropy is indeed contained in the short-time dynamics. In contrast, the contribution of the slow dynamics of side chains on time scales beyond τR on the side-chain conformational entropy is significant and should be taken into account for the extraction of accurate thermodynamic properties.

ABC exporters harness the energy of ATP to pump substrates across membranes. Extracellular gate opening and closure are key steps of the transport cycle, but the underlying mechanism is poorly understood. Here, we generated a synthetic single domain antibody (sybody) that recognizes the heterodimeric ABC exporter TM287/288 exclusively in the presence of ATP, which was essential to solve a 3.2 Å crystal structure of the outward-facing transporter. The sybody binds to an extracellular wing and strongly inhibits ATPase activity by shifting the transporter’s conformational equilibrium towards the outward-facing state, as shown by double electron-electron resonance (DEER). Mutations that facilitate extracellular gate opening result in a comparable equilibrium shift and strongly reduce ATPase activity and drug transport. Using the sybody as conformational probe, we demonstrate that efficient extracellular gate closure is required to dissociate the NBD dimer after ATP hydrolysis to reset the transporter back to its inward-facing state., ABC exporters hydrolyze ATP to pump substrates across membranes, but critical steps of the transport mechanism remain poorly understood. Here, the authors solve the crystal structure of outward-facing TM287/288 with the help of a state-specific sybody and gain insights into the role of the extracellular gate.

Solvent fluctuations have been explored in detail for idealized and rigid hydrophobic model systems, but so far it has remained unclear how internal protein motions and their coupling to the surrounding solvent affect the dynamics of ligand binding to biomolecular surfaces. Here, molecular dynamics simulations were used to elucidate the solvent-mediated binding of a model ligand to the hydrophobic surface patch of ubiquitin. The ligand's friction profiles reveal pronounced long-time correlations and enhanced friction in the vicinity of the protein, similar to idealized hydrophobic surfaces. Interestingly, these effects are shaped by internal protein motions. Protein flexibility modulates water density fluctuations near the hydrophobic surface patch and smooths out the friction profile of ligand binding.

Enhanced sampling techniques are a promising approach to obtain reliable binding free energy profiles for flexible protein-ligand complexes from molecular dynamics (MD) simulations. To put four popular enhanced sampling techniques to a biologically relevant and challenging test, we studied the partial dissociation of an antigenic peptide from the Major Histocompatibility Complex I (MHC I) HLA-B*35:01 to systematically investigate the performance of Umbrella Sampling (US), Replica Exchange with Solute Tempering 2 (REST2), Bias Exchange Umbrella Sampling (BEUS, or replica-exchange umbrella sampling), and well-tempered Metadynamics (MTD). With regard to the speed of sampling and convergence, the peptide-MHC I complex (pMHC I) under study showcases intrinsic strengths and weaknesses of the four enhanced sampling techniques used. We found that BEUS can handle best the sampling challenges that arise from the coexistence of an enthalpically and an entropically stabilized free energy minimum in the pMHC I under study. These findings might be relevant also for other flexible biomolecular systems with competing enthalpically and entropically stabilized minima.

An approach is presented to directly simulate the dynamics of methyl groups in protein side-chains, as accessible via NMR spin relaxation measurements, by all-atom MD simulations. The method, which does not rely on NMR information or any system-specific adjustable parameters, is based on calculating the time-correlation functions (TCFs) of the C-H bonds in methyl groups and explicitly takes the truncation of the TCFs due to overall tumbling of the molecule into account. Using ubiquitin as a model protein, we show (i) that an accurate description of the methyl dynamics requires reparametrization of the potential energy barriers of methyl group rotation in the AMBER ff99SB*-ILDN force field (and related parameter sets), which was done with CCSD(T) coupled cluster calculations of frequency isolated dipeptides as reference, and (ii) that the TIP4P/2005 solvation model yields overall tumbling correlation times that are in close agreement with experimental data. The methyl axis squared order parameters S-axis(2) and associated correlation times tau(f), obtained within the Lipari-Szabo formalism, are in good agreement with the values derived from NMR deuterium relaxation experiments. Importantly, the relaxation rates and spectral densities derived from MD and NMR agree as well, enabling a direct comparison without assumptions inherent to simplified motional models.

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side chain molecular motions and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations even now display very good agreement with the experiment. Results herein showcase the performance of present-day MD force fields and manifest their refined ability to accurately describe internal protein dynamics.

The collective behaviour of individual lipid molecules determines the properties of phospholipid membranes. However, the collective molecular motions often remain challenging to characterise at the desired spatial and temporal resolution. Here we study collective vibrational motion on picosecond time scales in dioleoylphosphatidylcholine lipid bilayers with varying cholesterol content using all-atom molecular dynamics simulations. Cholesterol is found to not only laterally compact the lipid bilayer, but also to change the velocity of longitudinal density fluctuations propagating in the plane of the membrane. Cholesterol-induced reduction of the area per lipid alters the collective dynamics of the lipid headgroups, but not of the lipid tails. The introduction of cholesterol reduces the number of water molecules interacting with the lipid headgroups, leading to a decrease in the velocity of the laterally-propagating sound mode. Thus, the stiffening effect of cholesterol is found to be indirect: decreasing the area per lipid weakens the interactions between the lipid headgroups and water. The collective modes characterised in this work can enable the membrane to dissipate excess energy and thus maintain its structural integrity, e.g., under mechanical stress.

Correlated vibrational motion on the sub-picosecond timescale and associated collective dynamics in a protein–membrane environment are characterized using molecular dynamics simulations. We specifically analyze correlated motion of a membrane-associated protein and a lipid bilayer for distinct separation distances. Correlated vibrations persist up to distances of 25 A between both biomolecular surfaces. These correlations are mediated by separating layers of water molecules, whose collective properties are altered by the simultaneous presence of protein and lipid bilayer interfaces.

Water dynamics in the hydration shell of the peripheral membrane protein annexin B12 were studied using MD simulations and Overhauser DNP-enhanced NMR. We show that retardation of water motions near phospholipid bilayers is extended by the presence of a membrane-bound protein, up to around 10 Å above that protein. Near the membrane surface, electrostatic interactions with the lipid head groups strongly slow down water dynamics, whereas protein-induced water retardation is weaker and dominates only at distances beyond 10 Å from the membrane surface. The results can be understood from a simple model based on additive contributions from the membrane and the protein to the activation free energy barriers of water diffusion next to the biomolecular surfaces. Furthermore, analysis of the intermolecular vibrations of the water network reveals that retarded water motions near the membrane shift the vibrational modes to higher frequencies, which we used to identify an entropy gradient from the membrane surface towards the bulk water. Our results have implications for processes that take place at lipid membrane surfaces, including molecular recognition, binding, and protein-protein interactions.

Proteins depend on defined molecular plasticity for their functionality. How to comprehensively capture dynamics correctly is of ubiquitous biological importance. Approaches commonly used to probe protein dynamics include model-free elucidation of site-specific motion by NMR relaxation, molecular dynamics (MD)-based approaches, and capturing the substates within a dynamic ensemble by recent eNOE-based multiple-structure approaches. Even though MD is sometimes combined with ensemble-averaged NMR restraints, these approaches have largely been developed and used individually. Owing to the different underlying concepts and practical requirements, it has remained unclear how they compare, and how they cross-validate and complement each other. Here, we extract and compare the differential information contents of MD simulations, NMR relaxation measurements, and eNOE-based multi-state structures for the SH3 domain of chicken α-spectrin. The data show that a validated, consistent, and detailed picture is feasible both for timescales and actual conformational states sampled in the dynamic ensemble. This includes the biologically important side-chain plasticity, for which experimentally cross-validated assessment is a significant challenge.

Nuclear magnetic resonance (NMR) spin relaxation has become the mainstay technique to sample protein dynamics at atomic resolution, expanding its repertoire from backbone 15N to side-chain 2H probes. At the same time, molecular dynamics (MD) simulations have become increasingly powerful to study protein dynamics due to steady improvements of physical models, algorithms, and computational power. Good agreement between generalized Lipari-Szabo order parameters derived from experiment and MD simulation has been observed for the backbone dynamics of a number of proteins. However, the agreement for the more dynamic side-chains, as probed by methyl group relaxation, was much worse. Here, we use T4 lysozyme (T4L), a protein with moderate tumbling anisotropy, to showcase a number of improvements that reduce this gap by a combined evaluation of NMR relaxation experiments and MD simulations. By applying a protein force field with accurate methyl group rotation barriers in combination with a solvation model that yields correct protein rotational diffusion times, we find that properly accounting for anisotropic protein tumbling is an important factor to improve the match between NMR and MD in terms of methyl axis order parameters, spectral densities, and relaxation rates. The best agreement with the experimentally measured relaxation rates is obtained by a posteriori fitting the appropriate internal time correlation functions, truncated by anisotropic overall tumbling. In addition, MD simulations led us to account for a hitherto unrealized artifact in deuterium relaxation experiments arising from strong coupling for leucine residues in uniformly 13C-enriched proteins. For T4L, the improved analysis reduced the RMSD between MD and NMR derived methyl axis order parameters from 0.19 to 0.11. At the level of the spectral density functions, the improvements allow us to extract the most accurate parameters that describe protein side-chain dynamics. Further improvement is challenging not only due to force field and sampling limitations in MD, but also due to inherent limitations of the Lipari-Szabo model to capture complex dynamics.

In bundled SPC water models, the relative motion of groups of four water molecules is restrained by distance-dependent potentials. Bundled SPC models have been used in hybrid all-atom/coarse-grained (AA/CG) multiscale simulations, since they enable to couple atomistic SPC water with supra-molecular CG water models that effectively represent more than a single water molecule. In the present work, we systematically validated and critically tested bundled SPC water models as solvent for biomolecular simulations. To that aim, we investigated both thermodynamic and structural properties of various biomolecular systems through molecular dynamics (MD) simulations. Potentials of mean force of dimerization of pairs of amino acid side chains as well as hydration free energies of single side chains obtained with bundled SPC and standard (unrestrained) SPC water agree closely with each other and with experimental data. Decomposition of the hydration free energies into enthalpic and entropic contributions reveals that in bundled SPC, this favorable agreement of the free energies is due to a larger degree of error compensation between hydration enthalpy and entropy. The Ramachandran maps of Ala3, Ala5, and Ala7 peptides are similar in bundled and unrestrained SPC, whereas for the (GS)2 peptide, bundled water leads to a slight overpopulation of extended conformations. Analysis of the end-to-end distance autocorrelation times of the Ala5 and (GS)2 peptides shows that sampling in more viscous bundled SPC water is about two times slower. Pronounced differences between the water models were found for the structure of a coiled-coil dimer, which is instable in bundled SPC but not in standard SPC. In addition, the hydration of the active site of the serine protease α-chymotrypsin depends on the water model. Bundled SPC leads to an increased hydration of the active site region, more hydrogen bonds between water and catalytic triad residues, and a significantly slower exchange of water molecules between the active site and the bulk. Our results form a basis for assessing the accuracy that can be expected from bundled SPC water models. At the same time, this study also highlights the importance of evaluating beforehand the effects of water bundling on the biomolecular system of interest for a particular multiscale simulation application.

Monoclonal antibody (mAb)-based therapeutics often require high-concentration formulations. Unfortunately, highly concentrated antibody solutions often have biophysical properties that are disadvantageous for therapeutic development, such as high viscosity, solubility limitations, precipitation issues, or liquid-liquid phase separation. In this work, we present a computational rational design principle for improving the thermodynamic stability of mAb solutions through targeted point mutations. Two publicly available IgG1 monoclonal antibodies that exhibit high viscosity at high concentrations were used as model systems. Guided by a computationally efficient approach that combines molecular dynamics simulations with three-dimensional reference interaction site model theory, point mutations of charged residues were introduced in the variable Fv regions in such a manner that the hydration free energy was optimized. Two selected point mutants were then produced by transient expression and characterized experimentally. Both engineered mAbs have reduced viscosity at high concentration, less negative second virial coefficient, and improved solubility compared to the respective wild-types. The results obtained with the suggested straightforward design principle underline the relevance of solvation effects for understanding, and ultimately optimizing, the properties of highly concentrated mAb solutions, with possible implications also for other biomolecular systems.

Solvation plays an important role in virtually all biomolecular recognition and binding processes. However, the consequences of changes in solvation conditions often remain elusive. In this work, we combined isothermal titration calorimetry (ITC) and molecular dynamics (MD) simulations to investigate the effect of solvent composition on the thermodynamics of protein-ligand binding. We studied the binding of p-aminobenzamidine (PAB) to trypsin in various water/methanol mixtures as a model system for a biomolecular complex. Our ITC experiments show that the free energy of binding changes only very modestly with methanol concentration, and that this small change is due to strong enthalpy-entropy compensation. The MD and free energy simulations not only reproduce the experimental binding free energies, but also provide atomic-level insights into the mechanisms behind the thermodynamic observations. The more favorable binding enthalpy at increased methanol concentrations (when compared to pure water) is attributed to stronger protein-ligand and intramolecular protein-protein interactions. The stronger protein-ligand interaction is linked to a small-scale conformational rearrangement of the L2 binding pocket loop, which senses the solvent environment. Remarkably, the stronger interactions do not substantially reduce the configurational entropy of the protein. Instead, the more disfavorable entropy contribution to the binding free energy at increased methanol concentrations is due to the desolvation of the ligand from the bulk, which is more favorable in pure aqueous solution than in the presence of methanol. Our work thus underpins the importance of including conformational flexibility, even for an overall rather rigid complex, since even small-amplitude motions can significantly alter the binding energetics. Furthermore, the ability of our combined ITC/MD approach to assign different thermodynamic contributions to distinct conformational states might contribute to an enhanced understanding of biomolecular binding processes in general.

We investigated the effect of amino acid composition and hydrophobic length of α-helical transmembrane peptides and the role of electrostatic interactions on the lateral diffusion of the peptides in lipid membranes. Model peptides of varying length and composition, and either tryptophans or lysines as flanking residues, were synthesized. The peptides were labeled with the fluorescent label Alexa Fluor 488 and incorporated into phospholipid bilayers of different hydrophobic thickness and composition. Giant unilamellar vesicles were formed by electroformation, and the lateral diffusion of the transmembrane peptides (and lipids) was determined by fluorescence correlation spectroscopy. In addition, we performed coarse-grained molecular-dynamics simulations of single peptides of different hydrophobic lengths embedded in planar membranes of different thicknesses. Both the experimental and simulation results indicate that lateral diffusion is sensitive to membrane thickness between the peptides and surrounding lipids. We did not observe a difference in the lateral diffusion of the peptides with respect to the presence of tryptophans or lysines as flanking residues. The specific lipid headgroup composition of the membrane has a much less pronounced impact on the diffusion of the peptides than does the hydrophobic thickness.

Light-driven molecular rotary motors derived from chiral overcrowded alkenes represent a broad class of compounds for which photochemical rearrangements lead to large scale motion of one part of the molecule with respect to another. It is this motion/change in molecular shape that is employed in many of their applications. A key group in this class are the molecular rotary motors that undergo unidirectional light-driven rotation about a double bond through a series of photochemical and thermal steps. In the present contribution we report a combined quantum chemical and molecular dynamics study of the mechanism of the rotational cycle of the fluorene-based molecular rotary motor 9-(2,4,7-trimethy1-2,3-dihydro-1H-inden-1-ylidene)-9H-fluorene (1). The potential energy surfaces of the ground and excited singlet states of I were calculated, and it was found that conical intersections play a central role in the mechanism of photo conversion between the stable conformer of 1 and its metastable conformer. Molecular dynamics simulations indicate that the average lifetime of the fluorene motor in the excited state is 1.40 +/- 0.10 ps when starting from the stable conformer, which increases to 1.77 +/- 0.13 ps for the reverse photoisomerization. These simulations indicate that the quantum yield of photoisomerization of the stable conformer is 0.92, whereas it is only 0.40 for the reverse photoisomerization. For the first time, a theoretical understanding of the experimentally observed photostationary state of 1 is reported that provides a detailed picture of the photoisomerization dynamics in overcrowded alkene-based molecular motor 1. The analysis of the electronic structure of the fluorene molecular motor holds considerable implications for the design of molecular motors. Importantly, the role of pyramidalization and conical intersections offer new insight into the factors that dominate the photostationary state achieved in these systems.

Substrate translocation by ATP-binding cassette (ABC) transporters involves coupling of ATP binding and hydrolysis in the nucleotide-binding domains (NBDs) to conformational changes in the transmembrane domains. We used molecular dynamics simulations to investigate the atomic-level mechanism of conformational coupling in the ABC transporter BtuCD-F, which imports vitamin B12 across the inner membrane of Escherichia coli. Our simulations show how an engineered disulfide bond across the NBD dimer interface reduces conformational fluctuations and hence configurational entropy. As a result, the disulfide bond is under substantial mechanical stress. Releasing this entropic spring, as is the case in the wild-type transporter, combined with analyzing the pairwise forces between individual residues, unravels the coupling mechanism. The identified pathways along which force is propagated from the NBDs via the coupling helix to the transmembrane domains are composed of highly conserved residues, underlining their functional relevance. This study not only reveals the details of conformational coupling in BtuCD-F, it also provides a promising approach to other long-range conformational couplings, e.g., in ABC exporters or other ATP-driven molecular machines.

Immune recognition of infected or malignantly transformed cells relies on antigenic peptides exposed at the cell surface by major histocompatibility complex class I (MHC I) molecules. Selection and loading of peptides onto MHC I is orchestrated by the peptide-loading complex (PLC), a multiprotein assembly whose structure has not yet been resolved. Tapasin, a central component of the PLC, stabilises MHC I and catalyses the exchange of low-affinity against high-affinity, immunodominant peptides. Up to now, the molecular basis of this peptide editing mechanism remained elusive. Here, using all-atom molecular dynamics (MD) simulations, we unravel the atomic details of how tapasin and antigen peptides act on the MHC I binding groove. Force distribution analysis reveals an intriguing molecular tug-of-war mechanism: only high-affinity peptides can exert sufficiently large forces to close the binding groove, thus overcoming the opposite forces exerted by tapasin to open it. Tapasin therefore accelerates the release of low-affinity peptides until a high-affinity antigen binds, promoting subsequent PLC break-down. Fluctuation and entropy analyses show how tapasin chaperones MHC I by stabilising it in a peptide-receptive conformation. Our results explain previous experiments and mark a key step towards a better understanding of adaptive immunity.

Salt bridges in lipid bilayers play a decisive role in the dynamic assembly and downstream signaling of the natural killer and T-cell receptors. Here, we describe the identification of an inter-subunit salt bridge in the membrane within yet another key component of the immune system, the peptide-loading complex (PLC). The PLC regulates cell surface presentation of self-antigens and antigenic peptides via molecules of the major histocompatibility complex class I. We demonstrate that a single salt bridge in the membrane between the transporter associated with antigen processing TAP and the MHC I-specific chaperone tapasin is essential for the assembly of the PLC and for efficient MHC I antigen presentation. Molecular modeling and all-atom molecular dynamics simulations suggest an ionic lock-switch mechanism for the binding of TAP to tapasin, in which an unfavorable uncompensated charge in the ER-membrane is prevented through complex formation. Our findings not only deepen the understanding of the interaction network within the PLC, but also provide evidence for a general interaction principle of dynamic multiprotein membrane complexes in immunity.

The peptide-loading complex plays a pivotal role in Ag processing and is thus central to the efficient immune recognition of virally and malignantly transformed cells. The underlying mechanism by which MHC class I (MHC I) molecules sample immunodominant peptide epitopes, however, remains poorly understood. In this article, we delineate the interaction between tapasin (Tsn) and MHC I molecules. We followed the process of peptide editing in real time after ultra-fast photoconversion to pseudoempty MHC I molecules. Tsn discriminates between MHC I loaded with optimal and MHC I bound to suboptimal cargo. This differential interaction is key to understanding the kinetics of epitope proofreading. To elucidate the underlying mechanism at the atomic level, we modeled the Tsn/MHC I complex using all-atom molecular dynamics simulations. We present a catalytic working cycle, in which Tsn binds to MHC I with suboptimal cargo and thereby adjusts the energy landscape in favor of MHC I complexes with immunodominant epitopes.

We investigated structural and dynamical properties of nanodiscs comprising dimyristoylphosphatidylcholine (DMPC) lipids and major scaffold protein MSP1Δ(1-22) from human apolipoprotein A-1 using combined all-atom and coarse-grained (CG) molecular dynamics (MD) simulations. The computational efficiency of the Martini-CG force field enables the spontaneous self-assembly of lipids and scaffold proteins into stable nanodisc structures on time scales up to tens of microseconds. Subsequent all-atom and CG-MD simulations reveal that the lipids in the nanodisc have lower configurational entropy and higher acyl tail order than in a lamellar bilayer phase. These altered average properties arise from rather differential behavior of lipids, depending on their location in the nanodisc. Since the scaffold proteins exert constrictive forces from the outer rim of the disc toward its center, lipids at the center of the nanodisc are highly ordered, whereas annular lipids that are in contact with the MSP proteins are remarkably disordered due to perturbed packing. Although specific differences between all-atom and CG simulations are also evident, the results obtained at both levels of resolution are in overall good agreement with each other and provide atomic level interpretations of recent experiments. Thus, the present study highlights the applicability of multiscale simulation approaches for nanodisc systems and opens the way for future applications, including the study of nanodisc-embedded membrane proteins.

We studied compositionally heterogeneous multi-component model membranes comprised of saturated lipids, unsaturated lipids, cholesterol, and a-helical TM protein models using coarse-grained molecular dynamics simulations. Reducing the mismatch between the length of the saturated and unsaturated lipid tails reduced the driving force for segregation into liquid-ordered (l(o)) and liquid-disordered (l(d)) lipid domains. Cholesterol depletion had a similar effect, and binary lipid mixtures without cholesterol did not undergo large-scale phase separation under the simulation conditions. The phase-separating ternary dipalmitoyl-phosphatidylcholine (DPPC)/dilinoleoyl-PC (DLiPC)/cholesterol bilayer was found to segregate into l(o) and l(d) domains also in the presence of a high concentration of TM helices. The l(d) domain was highly crowded with TM helices (protein-to-lipid ratio -1:5), slowing down lateral diffusion by a factor of 5-10 as compared to the dilute case, with anomalous (sub)-diffusion on the mu s time scale. The membrane with the less strongly unsaturated palmitoyl-linoleoyl-PC instead of DLiPC, which in the absence of TM alpha-helices less strongly deviated from ideal mixing, could be brought closer to a miscibility critical point by introducing a high concentration of TM helices. Finally, the 7-TM protein bacteriorhodopsin was found to partition into the l(d) domains irrespective of hydrophobic matching. These results show that it is possible to directly study the lateral reorganization of lipids and proteins in compositionally heterogeneous and crowded model biomembranes with coarse-grained molecular dynamics simulations, a step toward simulations of realistic, compositionally complex cellular membranes. This article is part of a Special Issue entitled: Protein Folding in Membranes. (C) 2011 Elsevier B.V. All rights reserved.

With today's available computer power, free energy calculations from equilibrium molecular dynamics simulations "via counting" become feasible for an increasing number of reactions. An example is the dimerization reaction of transmembrane alpha-helices. If an extended simulation of the two helices covers sufficiently many dimerization and dissociation events, their binding free energy is readily derived from the fraction of time during which the two helices are observed in dimeric form. Exactly how the correct value for the free energy is to be calculated, however, is unclear, and indeed several different and contradictory approaches have been used. In particular, results obtained via Boltzmann statistics differ from those determined via the law of mass action. Here, we develop a theory that resolves this discrepancy. We show that for simulation systems containing two molecules, the dimerization free energy is given by a formula of the form Delta G alpha ln(P(1)/P(0)). Our theory is also applicable to high concentrations that typically have to be used in molecular dynamics simulations to keep the simulation system small, where the textbook dilute approximations fail. It also covers simulations with an arbitrary number of monomers and dimers and provides rigorous error estimates. Comparison with test simulations of a simple Lennard Jones system with various particle numbers as well as with reference free energy values obtained from radial distribution functions show full agreement for both binding free energies and dimerization statistics. (C) 2011 Wiley Periodicals, Inc. J Comput Chem 32: 1919-1928, 2011

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models., Author Summary Many biomolecular processes involve charged species moving between regions of high polarity, such as the water phase, and regions of lower polarity, such as the lipid membrane. Due to the change in electrostatic screening between these two environments, the strength of the interactions between the moving charge and the surrounding molecules also changes. This has important consequences for the way biological activity is controlled. To help understand the forces driving the movement of biomolecules, we developed a computational model which is capable of describing these processes at near-atomic detail. To do so efficiently, we use a coarse-grained description of the molecules, in which some of the atomistic detail is averaged out. To capture the inhomogeneous nature of the dielectric response, we re-introduce some detail in the water model; the new model effectively mimics the orientational polarizability of real water molecules, and screens electrostatic interactions realistically. This enables the study of a number of important biological processes that were hitherto considered challenging for coarse-grained models, such as the permeation of ions across a lipid membrane and the rupture of membranes due to an electrostatic field, at relatively low computational cost.

We present an algorithm to reconstruct atomistic structures from their corresponding coarse-grained (CG) representations and its implementation into the freely available molecular dynamics (MD) program package GROMACS. The central part of the algorithm is a simulated annealing MD simulation in which the CG and atomistic structures are coupled via restraints. A number of examples demonstrate the application of the reconstruction procedure to obtain low-energy atomistic structural ensembles from their CO counterparts. We reconstructed individual molecules in vacuo (NCQ tripeptide, dipalmitoylphosphatidylcholine, and cholesterol), bulk water, and a WALP transmembrane peptide embedded in a solvated lipid bilayer. The first examples serve to optimize the parameters for the reconstruction procedure, whereas the latter examples illustrate the applicability to condensed-phase biomolecular systems. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 31 : 1333-1343, 2010

Protein-water interactions have widespread effects on protein structure and dynamics. As such, the function of many biomacromolecules can be directly related to the presence and exchange of water molecules. While the presence of structural water sites can be easily detected by X-ray crystallography, the dynamics within functional water-protein network architectures is largely elusive. Here we use solid-state NMR relaxation dispersion measurements with a focus on those active-site residues in the enzyme human carbonic anhydrase II (hCAII) that constitute the evolutionarily conserved water pocket, key for CAs' enzymatic catalysis. Together with chemical shifts, peak broadening, and results of molecular dynamics (MD) and DFT shift calculations, the relaxation dispersion data suggest the presence of a widespread fast μs-time-scale dynamics in the pocket throughout the protein-water network. This process is abrogated in the presence of an inhibitor which partially disrupts the network. The time scale of the protein-water pocket motion coincides both with the estimated residence time of Zn-bound water/OH- in the pocket showing the longest lifetimes in earlier magnetic relaxation dispersion experiments as well as with the rate-limiting step of catalytic turnover. As such, the reorganization of the water pocket:enzyme architecture might constitute an element of importance for enzymatic activity of this and possibly other proteins.