Polarizable atomic multipole-based force field for DOPC and POPE membrane lipids

Evaluating Polarizable Biomembrane Simulations against Experiments

Owing to the increase of available computational capabilities and the potential for providing a more accurate description, polarizable molecular dynamics force fields are gaining popularity in modeling biomolecular systems. It is, however, crucial to evaluate how much precision is truly gained with increasing cost and complexity of the simulation. Here, we leverage the NMRlipids open collaboration and Databank to assess the performance of available polarizable lipid models─the CHARMM-Drude and the AMOEBA-based parameters─against high-fidelity experimental data and compare them to the top-performing nonpolarizable models. While some improvement in the description of ion binding to membranes is observed in the most recent CHARMM-Drude parameters, and the conformational dynamics of AMOEBA-based parameters are excellent, the best nonpolarizable models tend to outperform their polarizable counterparts for each property we explored. The identified shortcomings range from inaccuracies in describing the conformational space of lipids to excessively slow conformational dynamics. Our results provide valuable insights for the further refinement of polarizable lipid force fields and for selecting the best simulation parameters for specific applications.

Water Nanoconfined in a Hydrophobic Pore: Molecular Dynamics Simulations of Transmembrane Protein 175 and the Influence of Water Models

Water molecules within biological ion channels are in a nanoconfined environment and therefore exhibit behaviors which differ from that of bulk water. Here, we investigate the phenomenon of hydrophobic gating, the process by which a nanopore may spontaneously dewet to form a "vapor lock" if the pore is sufficiently hydrophobic and/or narrow. This occurs without steric occlusion of the pore. Using molecular dynamics simulations with both rigid fixed-charge and polarizable (AMOEBA) force fields, we investigate this wetting/dewetting behavior in the transmembrane protein 175 ion channel. We examine how a range of rigid fixed-charge and polarizable water models affect wetting/dewetting in both the wild-type structure and in mutants chosen to cover a range of nanopore radii and pore-lining hydrophobicities. Crucially, we find that the rigid fixed-charge water models lead to similar wetting/dewetting behaviors, but that the polarizable water model resulted in an increased wettability of the hydrophobic gating region of the pore. This has significant implications for molecular simulations of nanoconfined water, as it implies that polarizability may need to be included if we are to gain detailed mechanistic insights into wetting/dewetting processes. These findings are of importance for the design of functionalized biomimetic nanopores (e.g., sensing or desalination) as well as for furthering our understanding of the mechanistic processes underlying biological ion channel function.

Multiscale simulations of large complexes in conjunction with cryo-EM analysis

The cellular environment is highly crowded with most proteins and RNA/DNA forming homomeric and heteromeric complexes. Essential questions regarding how these complexes switch between functional, rest, and abnormal states with regulators or modifications remain challenging and complicated. Here, we review the recent progress integrating cryoelectron microscopy and multiscale molecular modeling to understand the dynamics and function-related mechanism in protein-RNA/DNA complexes, protein-protein complexes/assemblies, and membrane protein complexes. One future direction of multiscale simulations will be to interpret the large complex multibody regulation in assembly-induced function enhancement in conjunction with advanced atomic resolution structural-biology techniques and specialized computing architectures.

Thermodynamics of ion binding and occupancy in potassium channels

Potassium channels modulate various cellular functions through efficient and selective conduction of K⁺ ions. The mechanism of ion conduction in potassium channels has recently emerged as a topic of debate. Crystal structures of potassium channels show four K⁺ ions bound to adjacent binding sites in the selectivity filter, while chemical intuition and molecular modeling suggest that the direct ion contacts are unstable. Molecular dynamics (MD) simulations have been instrumental in the study of conduction and gating mechanisms of ion channels. Based on MD simulations, two hypotheses have been proposed, in which the four-ion configuration is an artifact due to either averaged structures or low temperature in crystallographic experiments. The two hypotheses have been supported or challenged by different experiments. Here, MD simulations with polarizable force fields validated by ab initio calculations were used to investigate the ion binding thermodynamics. Contrary to previous beliefs, the four-ion configuration was predicted to be thermodynamically stable after accounting for the complex electrostatic interactions and dielectric screening. Polarization plays a critical role in the thermodynamic stabilities. As a result, the ion conduction likely operates through a simple single-vacancy and water-free mechanism. The simulations explained crystal structures, ion binding experiments and recent controversial mutagenesis experiments. This work provides a clear view of the mechanism underlying the efficient ion conduction and demonstrates the importance of polarization in ion channel simulations.

Polarization shapes the energy landscape of ion conduction in potassium channels.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.

Induced Polarization in Molecular Dynamics Simulations of the 5-HT3 Receptor Channel

Ion channel proteins form water-filled nanoscale pores within lipid bilayers, and their properties are dependent on the complex behavior of water in a nanoconfined environment. Using a simplified model of the pore of the 5-HT3 receptor (5HT3R) which restrains the backbone structure to that of the parent channel protein from which it is derived, we compare additive with polarizable models in describing the behavior of water in nanopores. Molecular dynamics simulations were performed with four conformations of the channel: two closed state structures, an intermediate state, and an open state, each embedded in a phosphatidylcholine bilayer. Water density profiles revealed that for all water models, the closed and intermediate states exhibited strong dewetting within the central hydrophobic gate region of the pore. However, the open state conformation exhibited varying degrees of hydration, ranging from partial wetting for the TIP4P/2005 water model to complete wetting for the polarizable AMOEBA14 model. Water dipole moments calculated using polarizable force fields also revealed that water molecules remaining within dewetted sections of the pore resemble gas phase water. Free energy profiles for Na+ and for Cl- ions within the open state pore revealed more rugged energy landscapes using polarizable force fields, and the hydration number profiles of these ions were also sensitive to induced polarization resulting in a substantive reduction of the number of waters within the first hydration shell of Cl- while it permeates the pore. These results demonstrate that induced polarization can influence the complex behavior of water and ions within nanoscale pores and provides important new insights into their chemical properties.

Intramembrane ionic protein-lipid interaction regulates integrin structure and function

Protein transmembrane domains (TMDs) are generally hydrophobic, but our bioinformatics analysis shows that many TMDs contain basic residues at terminal regions. Physiological functions of these membrane-snorkeling basic residues are largely unclear. Here, we show that a membrane-snorkeling Lys residue in integrin αLβ2 (also known as lymphocyte function-associated antigen 1 [LFA-1]) regulates transmembrane heterodimer formation and integrin adhesion through ionic interplay with acidic phospholipids and calcium ions (Ca²⁺) in T cells. The amino group of the conserved Lys ionically interacts with the phosphate group of acidic phospholipids to stabilize αLβ2 transmembrane association, thus keeping the integrin at low-affinity conformation. Intracellular Ca²⁺ uses its charge to directly disrupt this ionic interaction, leading to the transmembrane separation and the subsequent extracellular domain extension to increase adhesion activity. This Ca²⁺-mediated regulation is independent on the canonical Ca²⁺ signaling or integrin inside-out signaling. Our work therefore showcases the importance of intramembrane ionic protein–lipid interaction and provides a new mechanism of integrin activation.

Integrin αLβ2 is the major integrin in T cells and plays a vital role in regulating T-cell activation, adhesion, and migration. The transmembrane association of αL and β2 is crucial for maintaining the integrin at low-affinity conformation. Here, we find that the conserved basic residue (K702) in the transmembrane domain of β2 contributes to transmembrane association through ternary ionic interaction with acidic phospholipid and αL cytoplasmic residue. Upon T-cell activation, influxed calcium ions (Ca²⁺) can directly disrupt the ionic K702–lipid interaction through its positive charges, which leads to transmembrane separation and subsequent extracellular domain extension to switch αLβ2 to high-affinity conformation. This Ca²⁺-mediated regulation is through the modulation of the ionic Lys–lipid interaction but not through the canonical Ca²⁺ signaling or integrin inside-out signaling. Our study thus reports a new regulatory mechanism of integrin activation and showcases the importance of intramembrane ionic protein–lipid interaction. This finding might have general relevance, as bioinformatics analysis shows the presence of membrane-snorkeling basic residue is a common feature of transmembrane proteins.