Andrew H. Beaven

How Tolerant are Membrane Simulations with Mismatch in Area per Lipid between Leaflets

Difficulties in estimating the correct number of lipids in each leaflet of complex bilayer membrane simulation systems make it inevitable to introduce a mismatch in lipid packing (i.e., area per lipid) and thus alter the lateral pressure of each leaflet. To investigate potential impacts of such mismatch on simulation results, we performed molecular dynamics simulations of saturated and monounsaturated lipid bilayers with and without gramicidin A or WALP23 at various mismatches by adjusting the number of lipids in the lower leaflet from no mismatch to a 25% reduction compared to that in the upper leaflet. All simulations were stable under the constant pressure barostat, but the mismatch induces asymmetric lipid packing between the leaflets, so that the upper leaflet becomes more ordered, and the lower leaflet becomes less ordered. The mismatch impacts on various bilayer properties are mild up to 5–10% mismatch, and bilayers with fully saturated chains appear to be more prone to these impacts than those with unsaturated tails. The nonvanishing leaflet surface tensions and the free energy derivatives with respect to the bilayer curvature indicate that the bilayer would be energetically unstable in the presence of mismatch. We propose a quantitative criterion for allowable mismatch based on the energetics derived from a continuum elastic model, which grows as a square root of the number of the lipids in the system. On the basis of this criterion, we infer that the area per lipid mismatch up to 5% would be tolerable in various membrane simulations of reasonable all-atom system sizes (40–160 lipids per leaflet).

CHARMM-GUI 10 years for biomolecular modeling and simulation.

CHARMM‐GUI, http://www.charmm-gui.org, is a web‐based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM‐GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM‐GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all‐atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse‐grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ‐Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM‐GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram‐negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM‐GUI development project.

An extensive simulation study of lipid bilayer properties with different head groups, acyl chain lengths, and chain saturations

Previous MD simulations of six phosphocholine (PC) lipid bilayers demonstrated the accuracy of the CHARMM36 force field (C36FF) for PC bilayer simulation at varied temperatures (BBA-Biomembranes, 1838 (2014): 2520-2529). In this work, we further examine the accuracy of C36FF over a wide temperature range for a broader range of lipid types such as various head groups (phosphatidic acid (PA), PC, phosphoethanolamine (PE), phosphoglycerol (PG), and phosphoserine (PS)), and tails (saturated, mono-, mixed- and poly-unsaturated acyl chains with varied chain lengths). The structural properties (surface area per lipid (SA/lip), overall bilayer thickness, hydrophobic thickness, headgroup-to-headgroup thickness, deuterium order parameter (SCD), and spin-lattice relaxation time (T1)) obtained from simulations agree well with nearly all available experimental data. Our analyses indicate that PS lipids have the most inter-lipid hydrogen bonds, while PG lipids have the most intra-lipid hydrogen bonds, which play the main role in their low SA/lip in PS lipids and low thicknesses in PG lipids, respectively. PS, PE, and PA lipids have the largest contact clusters with on average 5-8 lipids per cluster, while PC and PG have clusters of 4 lipids based on a cutoff distance of 6.5Å. PS lipids have much slower lipid wobble (i.e., higher correlation time) than other head groups at a given temperature as the hydrogen bonded network significantly reduces a lipids mobility, and the rate of lipid wobble increases dramatically as temperature increases. These in-depth analyses facilitate further understanding of lipid bilayers at the atomic level.

Characterizing Residue-Bilayer Interactions Using Gramicidin A as a Scaffold and Tryptophan Substitutions as Probes

Previous experiments have shown that the lifetime of a gramicidin A dimer channel (which forms from two nonconducting monomers) in a lipid bilayer is modulated by mutations of the tryptophan (Trp) residues at the bilayer-water interface. We explore this further using extensive molecular dynamics simulations of various gA dimer and monomer mutants at the Trp positions in phosphatidylcholine bilayers with different tail lengths. gA interactions with the surrounding bilayer are strongly modulated by mutating these Trp residues. There are three principal effects: eliminating residue hydrogen bonding ability (i.e., reducing the channel-monolayer coupling strength) reduces the extent of the bilayer deformation caused by the assembled dimeric channel; a residues size and geometry affects its orientation, leading to different hydrogen bonding partners; and increasing a residues hydrophobicity increases the depth of gA monomer insertion relative to the bilayer center, thereby increasing the lipid bending frustration.

Biophysical and functional characterization of Norrin signaling through Frizzled4

Significance Wnt signaling has a broad spectrum of effects on cellular physiology and diseases like cancer. Much of the specificity and regulation mechanism of the Frizzled (Fzd) family, the main receptor of the Wnt pathway, remains elusive, but exploiting the activation mechanism of Fzd is pivotal to understanding the Wnt signaling pathway. Here, by using biophysical and biochemical techniques, we identify the key conformational changes that occur as the Norrin ligand binds to Fzd4, and we show the functional implications of the involved regions. We also reveal that the linker region plays an important role in communicating ligand binding with cytoplasmic signaling. Our observations offer insight into the activation mechanism of the Fzd family and the regulation of the Wnt signaling pathway. Wnt signaling is initiated by Wnt ligand binding to the extracellular ligand binding domain, called the cysteine-rich domain (CRD), of a Frizzled (Fzd) receptor. Norrin, an atypical Fzd ligand, specifically interacts with Fzd4 to activate β-catenin–dependent canonical Wnt signaling. Much of the molecular basis that confers Norrin selectivity in binding to Fzd4 was revealed through the structural study of the Fzd4CRD–Norrin complex. However, how the ligand interaction, seemingly localized at the CRD, is transmitted across full-length Fzd4 to the cytoplasm remains largely unknown. Here, we show that a flexible linker domain, which connects the CRD to the transmembrane domain, plays an important role in Norrin signaling. The linker domain directly contributes to the high-affinity interaction between Fzd4 and Norrin as shown by ∼10-fold higher binding affinity of Fzd4CRD to Norrin in the presence of the linker. Swapping the Fzd4 linker with the Fzd5 linker resulted in the loss of Norrin signaling, suggesting the importance of the linker in ligand-specific cellular response. In addition, structural dynamics of Fzd4 associated with Norrin binding investigated by hydrogen/deuterium exchange MS revealed Norrin-induced conformational changes on the linker domain and the intracellular loop 3 (ICL3) region of Fzd4. Cell-based functional assays showed that linker deletion, L430A and L433A mutations at ICL3, and C-terminal tail truncation displayed reduced β-catenin–dependent signaling activity, indicating the functional significance of these sites. Together, our results provide functional and biochemical dissection of Fzd4 in Norrin signaling.