Independent experimental and computational approaches show agreement concerning arginine/membrane interactions when a single arginine is introduced at selected positions within the membrane-spanning region of acetyl-GGALW5LALALAL12AL14ALALW19LAGA-ethanolamide, designated GWALP23. Peptide sequence isomers having Arg in position 12 or position 14 display markedly different behaviors, as deduced by both solid-state NMR experiments and coarse-grained molecular dynamics (CG-MD) simulations. With respect to the membrane normal of DOPC or DPPC lipid bilayer membranes, GWALP23-R14 shows one major state whose apparent average tilt is ~10° greater than that of GWALP23. The presence of R14 furthermore induces bilayer thinning and peptide displacement to “lift” the charged guanidinium toward the bilayer surface. By contrast, GWALP23-R12 exhibits multiple states that are in slow exchange on the NMR time scale, with CG-MD simulations indicating two distinct positions with different screw rotation angles in the membrane, along with an increased tendency to exit the lipid bilayer.

The use of molecular simulation to estimate the strength of macromolecular binding free energies is becoming increasingly widespread, with goals ranging from lead optimization and enrichment in drug discovery to personalizing or stratifying treatment regimes. In order to realize the potential of such approaches to predict new results, not merely to explain previous experimental findings, it is necessary that the methods used are reliable and accurate, and that their limitations are thoroughly understood. However, the computational cost of atomistic simulation techniques such as molecular dynamics (MD) has meant that until recently little work has focused on validating and verifying the available free energy methodologies, with the consequence that many of the results published in the literature are not reproducible. Here, we present a detailed analysis of two of the most popular approximate methods for calculating binding free energies from molecular simulations, molecular mechanics Poisson–Boltzmann surface area (MMPBSA) and molecular mechanics generalized Born surface area (MMGBSA), applied to the nine FDA-approved HIV-1 protease inhibitors. Our results show that the values obtained from replica simulations of the same protease–drug complex, differing only in initially assigned atom velocities, can vary by as much as 10 kcal mol–1, which is greater than the difference between the best and worst binding inhibitors under investigation. Despite this, analysis of ensembles of simulations producing 50 trajectories of 4 ns duration leads to well converged free energy estimates. For seven inhibitors, we find that with correctly converged normal mode estimates of the configurational entropy, we can correctly distinguish inhibitors in agreement with experimental data for both the MMPBSA and MMGBSA methods and thus have the ability to rank the efficacy of binding of this selection of drugs to the protease (no account is made for free energy penalties associated with protein distortion leading to the over estimation of the binding strength of the two largest inhibitors ritonavir and atazanavir). We obtain improved rankings and estimates of the relative binding strengths of the drugs by using a novel combination of MMPBSA/MMGBSA with normal mode entropy estimates and the free energy of association calculated directly from simulation trajectories. Our work provides a thorough assessment of what is required to produce converged and hence reliable free energies for protein–ligand binding.

Lymph node (LN) stromal cells, particularly fibroblastic reticular cells (FRCs), provide critical structural support and regulate immunity, tolerance and transport properties of LNs. In many tumors, LN metastasis is predictive of poor prognosis. However, the stromal contribution to the evolving microenvironment of tumor draining LNs (TDLNs) remains poorly understood. Here we show that FRCs specifically of TDLNs proliferate in response to tumor-derived cues and that the network they form is remodeled. Comparative transcriptional analysis of non-draining and TDLN FRCs demonstrated reprogramming of key pathways including matrix remodeling, chemokine/cytokine signaling and immune functions including leukocyte recruitment, migration and activation. In particular, downregulation of FRC-derived CCL21 and IL-7 were accompanied by altered immune composition and aberrant localization. These data imply that stroma of TDLNs adapt on multiple levels, following exposure to tumor-derived factors, to exhibit features typically associated with immune suppression.

Voltage-gated sodium channels have essential roles in electrical signalling. Prokaryotic sodium channels are tetramers consisting of transmembrane (TM) voltage-sensing and pore domains, and a cytoplasmic carboxy-terminal domain. Previous crystal structures of bacterial sodium channels revealed the nature of their TM domains but not their C-terminal domains (CTDs). Here, using electron paramagnetic resonance (EPR) spectroscopy combined with molecular dynamics, we show that the CTD of the NavMs channel from Magnetococcus marinus includes a flexible region linking the TM domains to a four-helix coiled-coil bundle. A 2.9 Å resolution crystal structure of the NavMs pore indicates the position of the CTD, which is consistent with the EPR-derived structure. Functional analyses demonstrate that the coiled-coil domain couples inactivation with channel opening, and is enabled by negatively charged residues in the linker region. A mechanism for gating is proposed based on the structure, whereby splaying of the bottom of the pore is possible without requiring unravelling of the coiled-coil.