Polarizable Water Model for the Coarse-Grained MARTINI Force Field

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models

Tonsillar metastasis of gastric cancer

Metastasis from a malignant tumor to the palatine tonsils is rare, with only 100 cases reported in the English-language literature. Tonsillar metastasis from a gastric cancer is very rare. We report here a case of palatine tonsillar metastasis after gastric cancer surgery. The patient was an 88-year-old woman who had gastric cancer with abdominal wall invasion. She had undergone a distal gastrectomy with abdominal wall resection and D2 lymph node dissection. Histologically, the tumor was primarily a moderately differentiated adenocarcinoma. It was stage IV (T4, N1, M0) using TNM clinical classification. The patient developed pharyngeal discomfort and abdominal pain and was hospitalized during the follow-up period, 1 year 9 months post-operatively. Multiple lung metastases, Virchow’s lymph node metastasis, and adrenal metastasis were observed. A mass of 2.5 cm was also observed in the right palatine tonsil. It was diagnosed as a moderately differentiated adenocarcinoma, a metastasis from gastric cancer. There was a concern of asphyxiation due to hemorrhage of the tumor; however, the tumor dislodged approximately 10 days after biopsy and tonsillar recurrence was not observed. The patient died 1 year 10 months post-operatively. In the literature there are cases with tonsillar metastases where surgical treatment, radiotherapy, and chemotherapy were performed and extension of survival was seen. Tonsillar metastasis is a form of systemic metastasis of a malignant tumor, and there is a high risk for asphyxiation from tumor dislodgement or hemorrhage. Thus, it is important to recognize tonsillar metastasis as an oncologic emergency

Conformational Changes and Slow Dynamics through Microsecond Polarized Atomistic Molecular Simulation of an Integral Kv1.2 Ion Channel

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins—not to mention numerous applications in drug design. Here, we present a full 1 µs atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120° rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation (∼35°) of the extracellular end of all S4 segments is present also in a reference 0.5 µs simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 310 helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4–lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5–1 µs). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations

Calcium Triggered Lα-H2 Phase Transition Monitored by Combined Rapid Mixing and Time-Resolved Synchrotron SAXS

BACKGROUND: Awad et al. reported on the Ca(2+)-induced transitions of dioleoyl-phosphatidylglycerol (DOPG)/monoolein (MO) vesicles to bicontinuous cubic phases at equilibrium conditions. In the present study, the combination of rapid mixing and time-resolved synchrotron small-angle X-ray scattering (SAXS) was applied for the in-situ investigations of fast structural transitions of diluted DOPG/MO vesicles into well-ordered nanostructures by the addition of low concentrated Ca(2+) solutions. METHODOLOGY/PRINCIPAL FINDINGS: Under static conditions and the in absence of the divalent cations, the DOPG/MO system forms large vesicles composed of weakly correlated bilayers with a d-spacing of approximately 140 A (L(alpha)-phase). The utilization of a stopped-flow apparatus allowed mixing these DOPG/MO vesicles with a solution of Ca(2+) ions within 10 milliseconds (ms). In such a way the dynamics of negatively charged PG to divalent cation interactions, and the kinetics of the induced structural transitions were studied. Ca(2+) ions have a very strong impact on the lipidic nanostructures. Intriguingly, already at low salt concentrations (DOPG/Ca(2+)>2), Ca(2+) ions trigger the transformation from bilayers to monolayer nanotubes (inverted hexagonal phase, H(2)). Our results reveal that a binding ratio of 1 Ca(2+) per 8 DOPG is sufficient for the formation of the H(2) phase. At 50 degrees C a direct transition from the vesicles to the H(2) phase was observed, whereas at ambient temperature (20 degrees C) a short lived intermediate phase (possibly the cubic Pn3m phase) coexisting with the H(2) phase was detected. CONCLUSIONS/SIGNIFICANCE: The strong binding of the divalent cations to the negatively charged DOPG molecules enhances the negative spontaneous curvature of the monolayers and causes a rapid collapsing of the vesicles. The rapid loss of the bilayer stability and the reorganization of the lipid molecules within ms support the argument that the transition mechanism is based on a leaky fusion of the vesicles

Membrainy: a ‘smart’, unified membrane analysis tool

BACKGROUND: The study of biological membranes using Molecular Dynamics has become an increasingly popular means by which to investigate the interactions of proteins, peptides and potentials with lipid bilayers. These interactions often result in changes to the properties of the lipids which can modify the behaviour of the membrane. Membrainy is a unified membrane analysis tool that contains a broad spectrum of analytical techniques to enable: measurement of acyl chain order parameters; presentation of 2D surface and thickness maps; determination of lateral and axial headgroup orientations; measurement of bilayer and leaflet thickness; analysis of the annular shell surrounding membrane-embedded objects; quantification of gel percentage; time evolution of the transmembrane voltage; area per lipid calculations; and quantification of lipid mixing/demixing entropy. RESULTS: Each analytical component within Membrainy has been tested on a variety of lipid bilayer systems and was found to be either comparable to or an improvement upon existing software. For the analytical techniques that have no direct comparable software, our results were confirmed with experimental data. CONCLUSIONS: Membrainy is a user-friendly, intelligent membrane analysis tool that automatically interprets a variety of input formats and force fields, is compatible with both single and double bilayers, and capable of handling asymmetric bilayers and lipid flip-flopping. Membrainy has been designed for ease of use, requiring no installation or configuration and minimal user-input to operate