Lars V. Schäfer

Polarizable Water Model for the Coarse-Grained MARTINI Force Field

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.

Improved Parameters for the Martini Coarse-Grained Protein Force Field

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.

Mechanoenzymatics of titin kinase

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.

Lipid packing drives the segregation of transmembrane helices into disordered lipid domains in model membranes

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.

Reconstruction of atomistic details from coarse-grained structures.

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

Transmembrane helices can induce domain formation in crowded model membranes.

We studied compositionally heterogeneous multi-component model membranes comprised of saturated lipids, unsaturated lipids, cholesterol, and α-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 ΤΜ 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 μs time scale. The membrane with the less strongly unsaturated palmitoyl-linoleoyl-PC instead of DLiPC, which in the absence of TM α-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.

Surface Hopping Excited-State Dynamics Study of the Photoisomerization of a Light-Driven Fluorene Molecular Rotary Motor

We report a theoretical study of the photoisomerization step in the operating cycle of a prototypical fluorene-based molecular rotary motor (1). The potential energy surfaces of the ground electronic state (S0) and the first singlet excited state (S1) are explored by semiempirical quantum-chemical calculations using the orthogonalization-corrected OM2 Hamiltonian in combination with a multireference configuration interaction (MRCI) treatment. The OM2/MRCI results for the S0 and S1 minima of the relevant 1-P and 1-M isomers and for the corresponding S0 transition state are in good agreement with higher-level calculations, both with regard to geometries and energetics. The S1 surface is characterized at the OM2/MRCI level by locating two S0-S1 minimum-energy conical intersections and nearby points on the intersection seam and by computing energy profiles for pathways toward these MECIs. Semiclassical Tully-type trajectory surface hopping (TSH) simulations with on-the-fly OM2/MRCI calculations are carried out to study the excited-state dynamics after photoexcitation to the S1 state. Fast relaxation to the ground state is observed through the conical intersection regions, predominantly through the lowest-energy one with a strongly twisted central C═C double bond and pyramidalized central carbon atom. The excited-state lifetimes for the direct and inverse photoisomerization reactions (1.40 and 1.79 ps) and the photostationary state ratio (2.7:1) from the TSH-OM2 simulations are in good agreement with the available experimental data (ca. 1.7 ps and 3:1). Excited-state lifetimes, photostationary state ratio, and dynamical details of the TSH-OM2 simulations also agree with classical molecular dynamics simulations using a reparametrized optimized potentials for liquid simulations (OPLS) all-atom force field with ad-hoc surface hops at predefined conical intersection points. The latter approach allows for a more extensive statistical sampling.

Chromophore Protonation State Controls Photoswitching of the Fluoroprotein asFP595

Fluorescent proteins have been widely used as genetically encodable fusion tags for biological imaging. Recently, a new class of fluorescent proteins was discovered that can be reversibly light-switched between a fluorescent and a non-fluorescent state. Such proteins can not only provide nanoscale resolution in far-field fluorescence optical microscopy much below the diffraction limit, but also hold promise for other nanotechnological applications, such as optical data storage. To systematically exploit the potential of such photoswitchable proteins and to enable rational improvements to their properties requires a detailed understanding of the molecular switching mechanism, which is currently unknown. Here, we have studied the photoswitching mechanism of the reversibly switchable fluoroprotein asFP595 at the atomic level by multiconfigurational ab initio (CASSCF) calculations and QM/MM excited state molecular dynamics simulations with explicit surface hopping. Our simulations explain measured quantum yields and excited state lifetimes, and also predict the structures of the hitherto unknown intermediates and of the irreversibly fluorescent state. Further, we find that the proton distribution in the active site of the asFP595 controls the photochemical conversion pathways of the chromophore in the protein matrix. Accordingly, changes in the protonation state of the chromophore and some proximal amino acids lead to different photochemical states, which all turn out to be essential for the photoswitching mechanism. These photochemical states are (i) a neutral chromophore, which can trans-cis photoisomerize, (ii) an anionic chromophore, which rapidly undergoes radiationless decay after excitation, and (iii) a putative fluorescent zwitterionic chromophore. The overall stability of the different protonation states is controlled by the isomeric state of the chromophore. We finally propose that radiation-induced decarboxylation of the glutamic acid Glu215 blocks the proton transfer pathways that enable the deactivation of the zwitterionic chromophore and thus leads to irreversible fluorescence. We have identified the tight coupling of trans-cis isomerization and proton transfers in photoswitchable proteins to be essential for their function and propose a detailed underlying mechanism, which provides a comprehensive picture that explains the available experimental data. The structural similarity between asFP595 and other fluoroproteins of interest for imaging suggests that this coupling is a quite general mechanism for photoswitchable proteins. These insights can guide the rational design and optimization of photoswitchable proteins.

Understanding the Dynamics Behind the Photoisomerization of a Light-Driven Fluorene Molecular Rotary Motor

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-trimethyl-2,3-dihydro-1H-inden-1-ylidene)-9H-fluorene (1). The potential energy surfaces of the ground and excited singlet states of 1 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.

Influence of Hydrophobic Mismatch and Amino Acid Composition on the Lateral Diffusion of Transmembrane Peptides

We investigated the effect of amino acid composition and hydrophobic length of alpha-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.