Insights into lipid-protein interactions from computer simulations

Utilizing Molecular Dynamics Simulations, Machine Learning, Cryo-EM, and NMR Spectroscopy to Predict and Validate Protein Dynamics

Protein dynamics play a crucial role in biological function, encompassing motions ranging from atomic vibrations to large-scale conformational changes. Recent advancements in experimental techniques, computational methods, and artificial intelligence have revolutionized our understanding of protein dynamics. Nuclear magnetic resonance spectroscopy provides atomic-resolution insights, while molecular dynamics simulations offer detailed trajectories of protein motions. Computational methods applied to X-ray crystallography and cryo-electron microscopy (cryo-EM) have enabled the exploration of protein dynamics, capturing conformational ensembles that were previously unattainable. The integration of machine learning, exemplified by AlphaFold2, has accelerated structure prediction and dynamics analysis. These approaches have revealed the importance of protein dynamics in allosteric regulation, enzyme catalysis, and intrinsically disordered proteins. The shift towards ensemble representations of protein structures and the application of single-molecule techniques have further enhanced our ability to capture the dynamic nature of proteins. Understanding protein dynamics is essential for elucidating biological mechanisms, designing drugs, and developing novel biocatalysts, marking a significant paradigm shift in structural biology and drug discovery.

Quantifying Induced Dipole Effects in Small Molecule Permeation in a Model Phospholipid Bilayer

The cell membrane functions as a semipermeable barrier that governs the transport of materials into and out of cells. The bilayer features a distinct dielectric gradient due to the amphiphilic nature of its lipid components. This gradient influences various aspects of small molecule permeation and the folding and functioning of membrane proteins. Here, we employ polarizable molecular dynamics simulations to elucidate the impact of the electronic environment on the permeation process. We simulated eight distinct amino-acid side chain analogs within a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer using the Drude polarizable force field (FF). Our approach includes both unbiased and umbrella sampling simulations. By using a polarizable FF, we sought to investigate explicit dipole responses in relation to local electric fields along the membrane normal. We evaluate molecular dipole moments, which exhibit variation based on their localization within the membrane, and compare the outcomes with analogous simulations using the nonpolarizable CHARMM36 FF. This comparative analysis aims to discern characteristic differences in the free energy surfaces of permeation for the various amino-acid analogs. Our results provide the first systematic quantification of the impact of employing an explicitly polarizable FF in this context compared to the fixed-charge convention inherent to nonpolarizable FFs, which may not fully capture the influence of the membrane dielectric gradient.

Curvature Footprints of Transmembrane Proteins in Simulations with the Martini Force Field

Membranes play essential roles in biological systems and are tremendously diverse in the topologies and chemical and elastic properties that define their functions. In many cases, a given membrane may display considerable heterogeneity, with localized clusters of lipids and proteins exhibiting distinct characteristics compared to adjoining regions. These lipid-protein assemblies can span nanometers to micrometers and are associated with cellular processes such as transport and signaling. While lipid-protein assemblages are dynamic, they can be stabilized by coupling between local membrane composition and shape. Due to the inherent difficulty in resolving atomistic details of membrane proteins in their native lipid environments, these complexes are notoriously challenging to study experimentally; however, molecular dynamics (MD) simulations might be a viable alternative. Here, we aim to assess the utility of coarse-grained (CG) MD simulations with the Martini force field for studying membrane curvature induced by transmembrane (TM) proteins that are reported to generate local curvature. The direction and magnitude of curvature induced by five different TM proteins, as well as certain lipid-protein and protein-protein interactions, were found to be in good agreement with available reference data.

Analyzing lipid distributions and curvature in molecular dynamics simulations of complex membranes

We describe methods to analyze lipid distributions and curvature in membranes with complex lipid mixtures and embedded membrane proteins. We discuss issues involved in these analyses, available tools to calculate curvature preferences of lipids and proteins, and focus on tools developed in our group for visual analysis of lipid-protein interactions and the analysis of membrane curvature.

Shocker─A Molecular Dynamics Protocol and Tool for Accelerating and Analyzing the Effects of Osmotic Shocks

The process of osmosis, a fundamental phenomenon in life, drives water through a semipermeable membrane in response to a solute concentration gradient across this membrane. In vitro, osmotic shocks are often used to drive shape changes in lipid vesicles, for instance, to study fission events in the context of artificial cells. While experimental techniques provide a macroscopic picture of large-scale membrane remodeling processes, molecular dynamics (MD) simulations are a powerful tool to study membrane deformations at the molecular level. However, simulating an osmotic shock is a time-consuming process due to slow water diffusion across the membrane, making it practically impossible to examine its effects in classic MD simulations. In this article, we present Shocker, a Python-based MD tool for simulating the effects of an osmotic shock by selecting and relocating water particles across a membrane over the course of several pumping cycles. Although this method is primarily aimed at efficiently simulating volume changes in vesicles, it can also handle membrane tubes and double bilayer systems. Additionally, Shocker is force field-independent and compatible with both coarse-grained and all-atom systems. We demonstrate that our tool is applicable to simulate both hypertonic and hypotonic osmotic shocks for a range of vesicular and bilamellar setups, including complex multicomponent systems containing membrane proteins or crowded internal solutions.

Electrostatic switch mechanisms of membrane protein trafficking and regulation

Lipid-protein interactions are normally classified as either specific or general. Specific interactions refer to lipid binding to specific binding sites within a membrane protein, thereby modulating the protein's thermal stability or kinetics. General interactions refer to indirect effects whereby lipids affect membrane proteins by modulating the membrane's physical properties, e.g., its fluidity, thickness, or dipole potential. It is not widely recognized that there is a third distinct type of lipid-protein interaction. Intrinsically disordered N- or C-termini of membrane proteins can interact directly but nonspecifically with the surrounding membrane. Many peripheral membrane proteins are held to the cytoplasmic surface of the plasma membrane via a cooperative combination of two forces: hydrophobic anchoring and electrostatic attraction. An acyl chain, e.g., myristoyl, added post-translationally to one of the protein's termini inserts itself into the lipid matrix and helps hold peripheral membrane proteins onto the membrane. Electrostatic attraction occurs between positively charged basic amino acid residues (lysine and arginine) on one of the protein's terminal tails and negatively charged phospholipid head groups, such as phosphatidylserine. Phosphorylation of either serine or tyrosine residues on the terminal tails via regulatory protein kinases allows for an electrostatic switch mechanism to control trafficking of the protein. Kinase action reduces the positive charge on the protein's tail, weakening the electrostatic attraction and releasing the protein from the membrane. A similar mechanism regulates many integral membrane proteins, but here only electrostatic interactions are involved, and the electrostatic switch modulates protein activity by altering the stabilities of different protein conformational states.

A high-throughput structural dynamics approach for identification of potential agonists of FFAR4 for type 2 diabetes mellitus therapy

Diabetes mellitus is a metabolic disorder that persists as a global threat to the world. A G-protein coupled receptor (GPCR), free fatty acid receptor 4 (FFAR4), has emerged as a potential target for type 2 diabetes mellitus (T2DM) and obesity-related disorders. The current study has investigated the FFAR4, deploying 3-dimensional structure modeling, molecular docking, machine learning, and high-throughput virtual screening methods to unravel the receptor's crucial and non-crucial binding site residues. We screened four lakh compounds and shortlisted them based on binding energy, stereochemical considerations, non-bonded interactions, and pharmacokinetic profiling. Out of the screened compounds, four compounds were selected for ligand-bound simulations. The molecular dynamic simulations were carried out for 1µs for native FFAR4 and 500 ns each for complexes of FFAR4 with compound 1, compound 2, compound 3, and compound 4. Our findings showed that in addition to reported binding site residues ARG99, ARG183, and VAL98 in known agonists like TUG-891, the amino acids ARG22, ARG24, THR23, TRP305, and GLU43 were also critical binding site residues. These amino acids impart stability to the FFAR4 complexes and contribute to the stronger binding affinity of the compounds. The study also indicated that aromatic residues like PHE211 are crucial for recognizing the active site's pi-pi and C-C double bonds. Since FFAR4 is a membrane protein, the simulation studies give an insight into the mechanisms of the crucial protein-lipid and lipid-water interactions. The analysis of the molecular dynamics trajectories showed all four compounds as potential hit molecules that can be developed further into potential agonists for T2DM therapy. Amongst the four compounds, compound 4 showed relatively better binding affinity, stronger non-bonded interactions, and a stable complex.Communicated by Ramaswamy H. Sarma.

Exploring the World of Membrane Proteins: Techniques and Methods for Understanding Structure, Function, and Dynamics

In eukaryotic cells, membrane proteins play a crucial role. They fall into three categories: intrinsic proteins, extrinsic proteins, and proteins that are essential to the human genome (30% of which is devoted to encoding them). Hydrophobic interactions inside the membrane serve to stabilize integral proteins, which span the lipid bilayer. This review investigates a number of computational and experimental methods used to study membrane proteins. It encompasses a variety of technologies, including electrophoresis, X-ray crystallography, cryogenic electron microscopy (cryo-EM), nuclear magnetic resonance spectroscopy (NMR), biophysical methods, computational methods, and artificial intelligence. The link between structure and function of membrane proteins has been better understood thanks to these approaches, which also hold great promise for future study in the field. The significance of fusing artificial intelligence with experimental data to improve our comprehension of membrane protein biology is also covered in this paper. This effort aims to shed light on the complexity of membrane protein biology by investigating a variety of experimental and computational methods. Overall, the goal of this review is to emphasize how crucial it is to understand the functions of membrane proteins in eukaryotic cells. It gives a general review of the numerous methods used to look into these crucial elements and highlights the demand for multidisciplinary approaches to advance our understanding.

Lipid-mediated prestin organization in outer hair cell membranes and its implications in sound amplification

Prestin is a high-density motor protein in the outer hair cells (OHCs), whose conformational response to acoustic signals alters the shape of the cell, thereby playing a major role in sound amplification by the cochlea. Despite recent structures, prestin's intimate interactions with the membrane, which are central to its function remained unresolved. Here, employing a large set (collectively, more than 0.5 ms) of coarse-grained molecular dynamics simulations, we demonstrate the impact of prestin's lipid-protein interactions on its organization at densities relevant to the OHCs and its effectiveness in reshaping OHCs. Prestin causes anisotropic membrane deformation, which mediates a preferential membrane organization of prestin where deformation patterns by neighboring copies are aligned constructively. The resulting reduced membrane rigidity is hypothesized to maximize the impact of prestin on OHC reshaping. These results demonstrate a clear case of protein-protein cooperative communication in membrane, purely mediated by interactions with lipids.

COGRIMEN: Coarse-Grained Method for Modeling of Membrane Proteins in Implicit Environments

The presented methodology is based on coarse-grained representation of biomolecules in implicit environments and is designed for the molecular dynamics simulations of membrane proteins and their complexes. The membrane proteins are not only found in the cell membrane but also in all membranous compartments of the cell: Golgi apparatus, mitochondria, endosomes and lysosomes, and they usually form large complexes. To investigate such systems the methodology is proposed based on two independent approaches combining the coarse-grained MARTINI model for proteins and the effective energy function to mimic the water/membrane environments. The latter is based on the implicit environment developed for all-atom simulations in the IMM1 method. The force field solvation parameters for COGRIMEN were initially calculated from IMM1 all-atom parameters and then optimized using Genetic Algorithms. The new methodology was tested on membrane proteins, their complexes and oligomers. COGRIMEN method is implemented as a patch for NAMD program and can be useful for fast and brief studies of large membrane protein complexes.

Studying KcsA Channel Clustering Using Single Channel Voltage-Clamp Fluorescence Imaging*

Oligomerization and complex formation play a key role for many membrane proteins and has been described to influence ion channel function in both neurons and the heart. In this study, we observed clustering of single KcsA channels in planar lipid bilayer using single molecule fluorescence, while simultaneously measuring single channel currents. Clustering coincided with cooperative opening of KcsA. We demonstrate that clustering was not caused by direct protein-protein interactions or hydrophobic mismatch with the lipid environment, as suggested earlier, but was mediated via microdomains induced by the channel in the lipid matrix. We found that single channel activity of KcsA requires conically-shaped lipids in the lamellar liquid-crystalline (L_α) phase, and the need for a negative spontaneous curvature seem to lead to the deformations in the membrane that cause the clustering. The method introduced here will be applicable to follow oligomerization of a wide range of membrane proteins.

Curvature-based sorting of eight lipid types in asymmetric buckled plasma membrane models

Curvature is a fundamental property of biological membranes and has essential roles in cellular function. Bending of membranes can be induced by their lipid and protein compositions, as well as peripheral proteins, such as those that make up the cytoskeleton. An important aspect of membrane function is the grouping of lipid species into microdomains, or rafts, which serve as platforms for specific biochemical processes. The fluid mosaic model of membranes has evolved to recognize the importance of curvature and leaflet asymmetry, and there are efforts towards evaluating their functional roles. This work investigates the effect of curvature on the sorting of lipids in buckled asymmetric bilayers containing eight lipid types, approximating an average mammalian plasma membrane, through coarse-grained (CG) molecular dynamics (MD) simulations with the Martini force field. The simulations reveal that i) leaflet compositional asymmetry can induce curvature asymmetry, ii) lipids are sorted by curvature to different extents, and iii) curvature-based partitioning trends show moderate to strong correlations with lipid molecular volumes and head to tail bead ratios, respectively. The findings provide unique insights into the role of curvature in membrane organization, and the curvature-based sorting trends should be useful references for later investigations, and potentially interpreting the functional roles of specific lipids.

Spontaneous local membrane curvature induced by transmembrane proteins

The (local) curvature of cellular membranes acts as a driving force for the targeting of membrane-associated proteins to specific membrane domains, as well as a sorting mechanism for transmembrane proteins, e.g., by accumulation in regions of matching spontaneous curvature. The latter measure was previously experimentally employed to study the curvature induced by the potassium channel KvAP and by aquaporin AQP0. However, the direction of the reported spontaneous curvature levels as well as the molecular driving forces governing the membrane curvature induced by these integral transmembrane proteins could not be addressed experimentally. Here, using both coarse-grained and atomistic molecular dynamics (MD) simulations, we report induced spontaneous curvature values for the homologous potassium channel Kv 1.2/2.1 Chimera (KvChim) and AQP0 embedded in unrestrained lipid bicelles that are in very good agreement with experiment. Importantly, the direction of curvature could be directly assessed from our simulations: KvChim induces a strong positive membrane curvature (≈0.036 nm^-1) whereas AQP0 causes a comparably small negative curvature (≈-0.019 nm^-1). Analyses of protein-lipid interactions within the bicelle revealed that the potassium channel shapes the surrounding membrane via structural determinants. Differences in shape of the protein-lipid interface of the voltage-gating domains between the extracellular and cytosolic membrane leaflets induce membrane stress and thereby promote a protein-proximal membrane curvature. In contrast, the water pore AQP0 displayed a high structural stability and an only faint effect on the surrounding membrane environment that is connected to its wedge-like shape.

Biophysical Reviews- 2021, the year that was

The current issue (volume 13 issue 6, 2021) is a Special Issue jointly dedicated to scientific content presented at the 20th triennial IUPAB Congress that was held in conjunction with both the 45th Annual Meeting of the Brazilian Biophysical Society (Sociedade Brasileira de Biofísica - SBBf) and the 50th Annual Meeting of the Brazilian Society for Biochemistry and Molecular Biology (Sociedade Brasileira de Bioquímica e Biologia Molecular - SBBq). In addition to describing the scientific and nonscientific content arising from the meeting this sub-editorial also provides a look back at some of the high points for Biophysical Reviews in the year 2021 before going on to describe a number of matters of interest to readers of the journal in relation to the coming year of 2022.