Search

The insertion and folding nascent or fully synthesized polypeptides into membranes is assisted by insertases. Here, the authors use a range of biophysical approaches to provide molecular details of how the transmembrane insertase YidC facilitates the insertion a protein into a phospholipid membrane.

Abstract Background Lipid-protein interactions stabilize protein oligomers, shape their structure, and modulate their function. Whereas in vitro experiments already account for the functional importance of lipids by using natural lipid extracts, in silico methods lack behind by embedding proteins in single component lipid bilayers. However, to accurately complement in vitro experiments with molecular details at very high spatio-temporal resolution, molecular dynamics simulations have to be performed in natural(-like) lipid environments. Results To enable more accurate MD simulations, we have prepared four membrane models of E. coli polar lipid extract, a typical model organism, each at all-atom (CHARMM36) and coarse-grained (Martini3) representations. These models contain all main lipid headgroup types of the E. coli inner membrane, i.e., phosphatidylethanolamines, phosphatidylglycerols, and cardiolipins, symmetrically distributed between the membrane leaflets. The lipid tail (un)saturation and propanylation stereochemistry represent the bacterial lipid tail composition of E. coli grown at 37∘C until 3/4 of the log growth phase. The comparison of the Simple three lipid component models to the complex 14-lipid component model Avanti over a broad range of physiologically relevant temperatures revealed that the balance of lipid tail unsaturation and propanylation in different positions and inclusion of lipid tails of various length maintain realistic values for lipid mobility, membrane area compressibility, lipid ordering, lipid volume and area, and the bilayer thickness. The only Simple model that was able to satisfactory reproduce most of the structural properties of the complex Avanti model showed worse agreement of the activation energy of basal water permeation with the here performed measurements. The Martini3 models reflect extremely well both experimental and atomistic behavior of the E. coli polar lipid extract membranes. Aquaporin-1 embedded in our native(-like) membranes causes partial lipid ordering and membrane thinning in its vicinity. Moreover, aquaporin-1 attracts and temporarily binds negatively charged lipids, mainly cardiolipins, with a distinct cardiolipin binding site in the crevice at the contact site between two monomers, most probably stabilizing the tetrameric protein assembly. Conclusions The here prepared and validated membrane models of E. coli polar lipids extract revealed that lipid tail complexity, in terms of double bond and cyclopropane location and varying lipid tail length, is key to stabilize membrane properties over a broad temperature range. In addition, they build a solid basis for manifold future simulation studies on more realistic lipid membranes bridging the gap between simulations and experiments.

G protein-coupled receptors (GPCRs) are the largest class of human membrane proteins that bind extracellular ligands at their orthosteric binding pocket to transmit signals to the cell interior. Ligand binding evokes conformational changes in GPCRs that trigger the binding of intracellular interaction partners (G proteins, G protein kinases, and arrestins), which initiate diverse cellular responses. It has become increasingly evident that the preference of a GPCR for a certain intracellular interaction partner is modulated by a diverse range of factors, e.g., ligands or lipids embedding the transmembrane receptor. Here, by means of molecular dynamics simulations of the β2-adrenergic receptor and β-arrestin2, we study how membrane lipids and receptor phosphorylation regulate GPCR-arrestin complex conformation and dynamics. We find that phosphorylation drives the receptor’s intracellular loop 3 (ICL3) away from a native negatively charged membrane surface to interact with arrestin. If the receptor is embedded in a neutral membrane, the phosphorylated ICL3 attaches to the membrane surface, which widely opens the receptor core. This opening, which is similar to the opening in the G protein-bound state, weakens the binding of arrestin. The loss of binding specificity is manifested by shallower arrestin insertion into the receptor core and higher dynamics of the receptor-arrestin complex. Our results show that receptor phosphorylation and the local membrane composition cooperatively fine-tune GPCR-mediated signal transduction. Moreover, the results suggest that deeper understanding of complex GPCR regulation mechanisms is necessary to discover novel pathways of pharmacological intervention.

Abstract Facilitated water permeation through narrow biological channels is fundamental for all forms of life. Despite its significance in health and disease as well as for biotechnological applications, the energetics of water permeation are still elusive. Gibbs free energy of activation is composed of an enthalpic and an entropic component. Whereas the enthalpic contribution is readily accessible via temperature dependent water permeability measurements, estimation of the entropic contribution requires information on the temperature dependence of the rate of water permeation. Here, we estimate, by means of accurate activation energy measurements of water permeation through Aquaporin-1 and by determining the accurate single channel permeability, the entropic barrier of water permeation through a narrow biological channel. Thereby the calculated value for $$\triangle {S}^{{{\ddagger}} }$$ △ S ‡ = 2.01 ± 0.82 J/(mol·K) links the activation energy of 3.75 ± 0.16 kcal/mol with its efficient water conduction rate of ~1010 water molecules/second. This is a first step in understanding the energetic contributions in various biological and artificial channels exhibiting vastly different pore geometries.

Gasdermin-A3 pore formation propagates along diverse pathways. It begins with membrane attachment and oligomeric pre-assembly. Once inserted in the membrane, the oligomers re-assemble into various shapes and sizes, which open their lytic pores.

Gasdermins execute programmatory cell death, known as pyroptosis, by forming medium-sized membrane pores. Recently, the molecular structure of those pores as well as the diversity in their shape and size have been revealed by cryoTEM and atomic force microscopy, respectively. Even though a growth of smaller to larger oligomers and reshaping from slits to rings could be documented, the initiation of the gasdermin pore formation remains a mystery. In one hypothesis, gasdermin monomers insert into membranes before associating into oligomeric pores. In the other hypothesis, gasdermin oligomers preassemble on the membrane surface prior to membrane insertion. Here, by studying the behavior of monomeric membrane-inserted gasdermin-A3 (GSDMA3), we unveil that a monomeric gasdermin prefers the membrane-adsorbed over the membrane-inserted state. Our results thus support the hypothesis of oligomers preassembling on the membrane surface before membrane penetration. At the same time, our simulations of small membrane-inserted arcs of GSDMA3 suggest that the inserting oligomer can be small and does not have to comprise a full ring of approximately 26–30 subunits. Moreover, our simulations have revealed an astonishingly large impact of salt-bridge formation and protein surroundings on the transmembrane passage of charged residues, reducing the energetic cost by up to 53% as compared to their free forms. The here observed free energy barrier of mere 5.6 kcal/mol for the membrane insertion of monomeric GSDMA3 explains the surprising ability of gasdermins to spontaneously self-insert into cellular membranes.

Lipids and proteins, as essential components of biological cell membranes, exhibit a significant degree of freedom for different kinds of motions including lateral long-range mobility. Due to their interactions, they not only preserve the cellular membrane but also contribute to many important cellular functions as e.g., signal transport or molecular exchange of the cell with its surrounding. Many of these processes take place on a short time (up to some nanoseconds) and length scale (up to some nanometers) which is perfectly accessible by quasielastic neutron scattering (QENS) experiments and molecular dynamics (MD) simulations. In order to probe the influence of a peptide, a transmembrane sequence of the transferrin receptor (TFRC) protein, on the dynamics of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) large unilamellar vesicles (LUVs) on a nanosecond time scale, high-resolution QENS experiments and complementary MD simulations have been utilized. By using different scattering contrasts in the experiment (chain-deuterated lipids and protonated lipids, respectively), a model could be developed which allows to examine the lipid and peptide dynamics separately. The experimental results revealed a restricted lipid lateral mobility in the presence of the TFRC transmembrane peptides. Also the apparent self-diffusion coefficient of the lateral movement of the peptide molecules could be determined quantitatively for the probed short-time regime. The findings could be confirmed very precisely by MD simulations. Furthermore, the article presents an estimation for the radius of influence of the peptides on the lipid long-range dynamics which could be determined by consistently combining results from experiment and simulation.

The structure of dithienobicyclo[4.4.1]undeca-3,8-diene-11-one ethylene glycol ketal (database code RESVAN) was determined using the wave function theory (WFT) as well as density functional theory (DFT) methods combined with various Gaussian AO basis sets. The apparently most accurate procedure, employing the CCSD(T)/complete basis set (CBS), provides an S−S distance and an angle between the two thiophene rings which differ considerably from experimental values. The best agreement with the experimental data among all WFT methods was surprisingly obtained at the MP3/aug-cc-pVDZ and MP3/CBS(B) levels (the correction term to CBS was obtained by the aug-cc-pVDZ basis set). The very good results obtained by the CCSD(T)/6-31G* method are clearly a consequence of fortunate error compensation. MP2 calculations, even with a small basis set, overestimate the attraction between the thiophene rings, and the worst agreement with experimental data was found in full MP2/QZVP method optimizations (i.e., a strong distortion of the thiophene rings was observed). The SCS(MI)-MP2 and SCS-MP2 methods exhibit improvement over the MP2 procedure. All standard DFT approaches fail to predict reasonable S−S distances. The lack of intramolecular London dispersion energy results in too great distance between the thiophene rings. Much better agreement with experiment was obtained if advanced DFT methods, covering dispersion effects, were used. The best results were obtained at the TPSS–D/TZVP, M06-L/TZVP and B2PLYP-D/def2-TZVP levels. When a larger basis (LP in the case of TPSS functional) or more advanced versions of the new Truhlar functionals (M06-2X) was used, the agreement with experiment deteriorated. The accurate description of this molecule is highly functional/basis dependent and this dependence is hardly predictable. To estimate effects of neighboring molecules in the experimental crystal structure, an optimization in the electric field of the 26 closest RESVAN molecules was performed, which, however, leads to only moderate (<0.05 Å) changes of the S−S distance. [ABSTRACT FROM AUTHOR]

A detailed quantum chemical study on five peptides (WG, WGG, FGG, GGF and GFA) containing the residues phenylalanyl (F), glycyl (G), tryptophyl (W) and alanyl (A)—where F and W are of aromatic character—is presented. When investigating isolated small peptides, the dispersion interaction is the dominant attractive force in the peptide backbone–aromatic side chain intramolecular interaction. Consequently, an accurate theoretical study of these systems requires the use of a methodology covering properly the London dispersion forces. For this reason we have assessed the performance of the MP2, SCS-MP2, MP3, TPSS-D, PBE-D, M06-2X, BH&H, TPSS, B3LYP, tight-binding DFT-D methods and ff99 empirical force field compared to CCSD(T)/complete basis set (CBS) limit benchmark data. All the DFT techniques with a ‘-D’ symbol have been augmented by empirical dispersion energy while the M06-2X functional was parameterized to cover the London dispersion energy. For the systems here studied we have concluded that the use of the ff99 force field is not recommended mainly due to problems concerning the assignment of reliable atomic charges. Tight-binding DFT-D is efficient as a screening tool providing reliable geometries. Among the DFT functionals, the M06-2X and TPSS-D show the best performance what is explained by the fact that both procedures cover the dispersion energy. The B3LYP and TPSS functionals—not covering this energy—fail systematically. Both, electronic energies and geometries obtained by means of the wave-function theory methods compare satisfactorily with the CCSD(T)/CBS benchmark data. [ABSTRACT FROM AUTHOR]

Stabilisation energies of stacked structures of C6H6⋯C6X6 (X = F, Cl, Br, CN) complexes were determined at the CCSD(T) complete basis set (CBS) limit level. These energies were constructed from MP2/CBS stabilisation energies and a CCSD(T) correction term determined with a medium basis set (6-31G**). The former energies were extrapolated using the two-point formula of Helgaker et al. from aug-cc-pVDZ and aug-cc-pVTZ Hartree–Fock energies and MP2 correlation energies. The CCSD(T) correction term is systematically repulsive. The final CCSD(T)/CBS stabilisation energies are large, considerably larger than previously calculated and increase in the series as follows: hexafluorobenzene (6.3 kcal mol−1), hexachlorobenzene (8.8 kcal mol−1), hexabromobenzene (8.1 kcal mol−1) and hexacyanobenzene (11.0 kcal mol−1). MP2/SDD** relativistic calculations performed for all complexes mentioned and also for benzene⋯hexaiodobenzene have clearly shown that due to relativistic effects the stabilisation energy of the hexaiodobenzene complex is lower than that of hexabromobenzene complex. The decomposition of the total interaction energy to physically defined energy components was made by using the symmetry adapted perturbation treatment (SAPT). The main stabilisation contribution for all complexes investigated is due to London dispersion energy, with the induction term being smaller. Electrostatic and induction terms which are attractive are compensated by their exchange counterparts. The stacked motif in the complexes studied is very stable and might thus be valuable as a supramolecular synthon. [ABSTRACT FROM AUTHOR]

Antigen binding to B-cell antigen receptors (BCRs) followed by signaling initiates the humoral immune response. The signaling is intimately coupled to nanoclustering of BCRs and their sorting to specific membrane domains, a process that is ruled by interactions between the BCR transmembrane domain and lipids. While the structure of the extracellular domains of BCRs has been resolved, little is known about the configuration of the constituting four immunoglobulin domains spanning the membrane. Here, we modeled the structure of the transmembrane (TM) domain of the IgM B-cell receptor using self-assembly coarse-grained molecular dynamics simulations. The obtained quaternary structure was validated against available experimental data and atomistic simulations. The IgM-BCR-TM domain configuration shows a 1:1 stoichiometry between the homodimeric membrane-bound domain of IgM (mIgM) and a Ig-α/Ig-β heterodimer. The mIgM homodimer is based on an asymmetric association of two mIgM domains. We show that a specific site of the Ig-α/Ig-β heterodimer is responsible for the association of IgM-BCRs with lipid rafts. Our results further suggest that this site is blocked in small-sized IgM-BCR clusters. The BCR TM structure provides a molecular basis for the previously suggested dissociation activation model of B-cell receptors. Self-assembly molecular dynamics simulations at the coarse-grained scale here proved as a versatile tool in the study of receptor complexes.

Chemokine receptors, a subclass of G protein coupled receptors (GPCRs), play essential roles in the human immune system, they are involved in cancer metastasis as well as in HIV-infection. A plethora of studies show that homo- and heterodimers or even higher order oligomers of the chemokine receptors CXCR4, CCR5, and CCR2 modulate receptor function. In addition, membrane cholesterol affects chemokine receptor activity. However, structural information about homo- and heterodimers formed by chemokine receptors and their interplay with cholesterol is limited. Here, we report homo- and heterodimer configurations of the chemokine receptors CXCR4, CCR5, and CCR2 at atomistic detail, as obtained from thousands of molecular dynamics simulations. The observed homodimerization patterns were similar for the closely related CC chemokine receptors, yet they differed significantly between the CC receptors and CXCR4. Despite their high sequence identity, cholesterol modulated the CC homodimer interfaces in a subtype-specific manner. Chemokine receptor heterodimers display distinct dimerization patterns for CXCR4/CCR5 and CXCR4/CCR2. Furthermore, associations between CXCR4 and CCR5 reveal an increased cholesterol-sensitivity as compared to CXCR4/CCR2 heterodimerization patterns. This work provides a first comprehensive structural overview over the complex interaction network between chemokine receptors and indicates how heterodimerization and the interaction with the membrane environment diversifies the function of closely related GPCRs.

G protein coupled receptors (GPCRs) allow for the transmission of signals across biological membranes. For a number of GPCRs, this signaling was shown to be coupled to prior dimerization of the receptor. The chemokine receptor type 4 (CXCR4) was reported before to form dimers and their functionality was shown to depend on membrane cholesterol. Here, we address the dimerization pattern of CXCR4 in pure phospholipid bilayers and in cholesterol-rich membranes. Using ensembles of molecular dynamics simulations, we show that CXCR4 dimerizes promiscuously in phospholipid membranes. Addition of cholesterol dramatically affects the dimerization pattern: cholesterol binding largely abolishes the preferred dimer motif observed for pure phospholipid bilayers formed mainly by transmembrane helices 1 and 7 (TM1/TM5-7) at the dimer interface. In turn, the symmetric TM3,4/TM3,4 interface is enabled first by intercalating cholesterol molecules. These data provide a molecular basis for the modulation of GPCR activity by its lipid environment.