The Martini coarse-grained force field has been successfully used for simulating a wide range of (bio)molecular systems. Recent progress in our ability to test the model against fully atomistic force fields, however, has revealed some shortcomings. Most notable, phenylalanine and proline were too hydrophobic, and dimers formed by polar residues in apolar solvents did not bind strongly enough. Here, we reparametrize these residues either through reassignment of particle types or by introducing embedded charges. The new parameters are tested with respect to partitioning across a lipid bilayer, membrane binding of Wimley–White peptides, and dimerization free energy in solvents of different polarity. In addition, we improve some of the bonded terms in the Martini protein force field that lead to a more realistic length of α-helices and to improved numerical stability for polyalanine and glycine repeats. The new parameter set is denoted Martini version 2.2.

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. 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.

Biological responses to mechanical stress require strain-sensing molecules, whose mechanically induced conformational changes are relayed to signaling cascades mediating changes in cell and tissue properties. In vertebrate muscle, the giant elastic protein titin is involved in strain sensing via its C-terminal kinase domain (TK) at the sarcomeric M-band and contributes to the adaptation of muscle in response to changes in mechanical strain. TK is regulated in a unique dual autoinhibition mechanism by a C-terminal regulatory tail, blocking the ATP binding site, and tyrosine autoinhibition of the catalytic base. For access to the ATP binding site and phosphorylation of the autoinhibitory tyrosine, the C-terminal autoinhibitory tail needs to be removed. Here, we use AFM-based single-molecule force spectroscopy, molecular dynamics simulations, and enzymatics to study the conformational changes during strain-induced activation of human TK. We show that mechanical strain activates ATP binding before unfolding of the structural titin domains, and that TK can thus act as a biological force sensor. Furthermore, we identify the steps in which the autoinhibition of TK is mechanically relieved at low forces, leading to binding of the cosubstrate ATP and priming the enzyme for subsequent autophosphorylation and substrate turnover.

Proteins that can be reversibly photoswitched between a fluorescent and a nonfluorescent state bear enormous potential in diverse fields, such as data storage, in vivo protein tracking, and subdiffraction resolution light microscopy. However, these proteins could hitherto not live up to their full potential because the molecular switching mechanism is not resolved. Here, we clarify the molecular photoswitching mechanism of asFP595, a green fluorescent protein (GFP)-like protein that can be transferred from a nonfluorescent “off” to a fluorescent “on” state and back again, by green and blue light, respectively. To this end, we establish reversible photoswitching of fluorescence in whole protein crystals and show that the switching kinetics in the crystal is identical with that in solution. Subsequent x-ray analysis demonstrated that upon the absorption of a green photon, the chromophore isomerizes from a trans (off) to a cis (on) state. Molecular dynamics calculations suggest that isomerization occurs through a bottom hula twist mechanism with concomitant rotation of both bonds of the chromophoric methine ring bridge. This insight into the switching mechanism should facilitate the targeted design of photoswitchable proteins. Reversible photoswitching of the protein chromophore system within intact crystals also constitutes a step toward the use of fluorescent proteins in three-dimensional data recording.

Cell membranes are comprised of multicomponent lipid and protein mixtures that exhibit a complex partitioning behavior. Regions of structural and compositional heterogeneity play a major role in the sorting and self-assembly of proteins, and their clustering into higher-order oligomers. Here, we use computer simulations and optical microscopy to study the sorting of transmembrane helices into the liquid-disordered domains of phase-separated model membranes, irrespective of peptide–lipid hydrophobic mismatch. Free energy calculations show that the enthalpic contribution due to the packing of the lipids drives the lateral sorting of the helices. Hydrophobic mismatch regulates the clustering into either small dynamic or large static aggregates. These results reveal important molecular driving forces for the lateral organization and self-assembly of transmembrane helices in heterogeneous model membranes, with implications for the formation of functional protein complexes in real cells.