Molecular dynamics (MD) simulations offer the possibility to study biological processes at high spatial and temporal resolution often not reachable by experiments. Corresponding biomolecular force field parameters have been developed for a wide variety of molecules ranging from inorganic ligands and small organic molecules over proteins and lipids to nucleic acids. Force fields have typically been parametrized and validated on thermodynamic observables and structural characteristics of individual compounds, e.g. of soluble proteins or lipid bilayers. Less strictly, due to the added complexity and missing experimental data to compare to, force fields have hardly been tested on the properties of mixed systems, e.g. on protein-lipid systems. Their selection and combination for mixed systems is further complicated by the partially differing parametrization strategies. Additionally, the presence of other compounds in the system may shift the subtle balance of force field parameters. Here, we assessed the protein-lipid interactions as described in the four atomistic force fields GROMOS54a7, CHARMM36 and the two force field combinations Amber14sb/Slipids and Amber14sb/Lipid14. Four observables were compared, focusing on the membrane-water interface: the conservation of the secondary structure of transmembrane proteins, the positioning of transmembrane peptides relative to the lipid bilayer, the insertion depth of side chains of unfolded peptides absorbed at the membrane interface, and the ability to reproduce experimental insertion energies of Wimley-White peptides at the membrane interface. Significant differences between the force fields were observed that affect e.g. membrane insertion depths and tilting of transmembrane peptides.

Cysteine side chains can exist in distinct oxidation states depending on the pH and redox potential of the environment, and cysteine oxidation plays important yet complex regulatory roles. Compared with the effects of post-translational modifications such as phosphorylation, the effects of oxidation of cysteine to sulfenic, sulfinic, and sulfonic acid on protein structure and function remain relatively poorly characterized. We present an analysis of the role of cysteine reactivity as a regulatory factor in proteins, emphasizing the interplay between electrostatics and redox potential as key determinants of the resulting oxidation state. A review of current computational approaches suggests underdeveloped areas of research for studying cysteine reactivity through molecular simulations.

The metabotropic glutamate receptor (mGluR) 2 plays a key role in the central nervous system. mGluR2 has been shown to be regulated by its surrounding lipid environment, especially by cholesterol, by an unknown mechanism. Here, using a combination of biochemical approaches, photo-cross-linking experiments, and molecular dynamics simulations we show the interaction of cholesterol with at least two, but potentially five more, preferential sites on the mGluR2 transmembrane domain. Our simulations demonstrate that surface matching, rather than electrostatic interactions with specific amino acids, is the main factor defining cholesterol localization. Moreover, the cholesterol localization observed here is similar to the sterol-binding pattern previously described in silico for other members of the mGluR family. Biochemical assays suggest little influence of cholesterol on trafficking or dimerization of mGluR2. Nevertheless, simulations revealed a significant reduction of residue−residue contacts together with an alteration in the internal mechanical stress at the cytoplasmic side of the helical bundle when cholesterol was present in the membrane. These alterations may be related to destabilization of the basal state of mGluR2. Due to the high sequence conservation of the transmembrane domains of mGluRs, the molecular interaction of cholesterol and mGluR2 described here is also likely to be relevant for other members of the mGLuR family.

No abstract