Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

Effect of a low fat versus a low carbohydrate weight loss dietary intervention on biomarkers of long term survival in breast cancer patients ('CHOICE'): study protocol

<p>Abstract</p> <p>Background</p> <p>Weight loss in overweight or obese breast cancer patients is associated with an improved prognosis for long term survival. However, it is not clear whether the macronutrient composition of the chosen weight loss dietary plan imparts further prognostic benefit. A study protocol is presented for a dietary intervention to investigate the effects of weight loss dietary patterns that vary markedly in fat and carbohydrate contents on biomarkers of exposure to metabolic processes that may promote tumorigenesis and that are predictive of long term survival. The study will also determine how much weight must be lost for biomarkers to change in a favorable direction.</p> <p>Methods/Design</p> <p>Approximately 370 overweight or obese postmenopausal breast cancer survivors (body mass index: 25.0 to 34.9 kg/m²) will be accrued and assigned to one of two weight loss intervention programs or a non-intervention control group. The dietary intervention is implemented in a free living population to test the two extremes of popular weight loss dietary patterns: a high carbohydrate, low fat diet versus a low carbohydrate, high fat diet. The effects of these dietary patterns on biomarkers for glucose homeostasis, chronic inflammation, cellular oxidation, and steroid sex hormone metabolism will be measured. Participants will attend 3 screening and dietary education visits, and 7 monthly one-on-one dietary counseling and clinical data measurement visits in addition to 5 group visits in the intervention arms. Participants in the control arm will attend two clinical data measurement visits at baseline and 6 months. The primary outcome is high sensitivity C-reactive protein. Secondary outcomes include interleukin-6, tumor necrosis factor-α, insulin-like growth factor-1 (IGF), IGF binding protein-3, 8-isoprostane-F2-alpha, estrone, estradiol, progesterone, sex hormone binding globulin, adiponectin, and leptin.</p> <p>Discussion</p> <p>While clinical data indicate that excess weight for height is associated with poor prognosis for long term survival, little attention is paid to weight control in the clinical management of breast cancer. This study will provide information that can be used to answer important patient questions about the effects of dietary pattern and magnitude of weight loss on long term survival following breast cancer treatment.</p> <p>Clinical Trial Registration</p> <p>CA125243</p

Membrane remodeling by the M2 amphipathic helix drives influenza virus membrane scission

Membrane scission is a crucial step in all budding processes, from endocytosis to viral budding. Many proteins have been associated with scission, though the underlying molecular details of how scission is accomplished often remain unknown. Here, we investigate the process of M2-mediated membrane scission during the budding of influenza viruses. Residues 50–61 of the viral M2 protein bind membrane and form an amphipathic α-helix (AH). Membrane binding requires hydrophobic interactions with the lipid tails but not charged interactions with the lipid headgroups. Upon binding, the M2AH induces membrane curvature and lipid ordering, constricting and destabilizing the membrane neck, causing scission. We further show that AHs in the cellular proteins Arf1 and Epsin1 behave in a similar manner. Together, they represent a class of membrane-induced AH domains that alter membrane curvature and fluidity, mediating the scission of constricted membrane necks in multiple biological pathways

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta