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

Impact of carbohydrate restriction with and without fatty acid loading on myocardial 18F-FDG uptake during PET: A randomized controlled trial

Low-carbohydrate (LC) and high-fat, low-carbohydrate (HFLC) dietary preparations may enhance 18F-FDG-PET-based imaging of small, inflamed structures near the heart by suppressing myocardial FDG signal. We compared myocardial 18F-FDG uptake in patients randomized to LC, HFLC, and unrestricted (UR) preparations prior to 18F-FDG-PET. We randomized 63 outpatients referred for oncologic 18F-FDG-PET to LC, HFLC, or UR dietary preparations (1:1:1 allocation) starting the evening before PET. After eating dinner according to instructions, UR and LC patients fasted until FDG injection (mean time 745 minutes for UR, 899 minutes for LC), and HFLC patients drank a fatty drink 60-70 minutes prior to FDG injection. Attenuation-corrected PET imaging was performed 60 minutes after FDG administration. Maximal myocardial standard uptake values (MyoSUVmax) were systematically measured in axial view and compared between the three groups. Using UR patients as reference, mean MyoSUVmax was lower in LC patients (3.3 ± 2.7 vs 6.2 ± 5.2, P = .03) but not in HFLC patients (5.5 ± 4.2, P = .63). Ratios of MyoSUVmax to liver SUVmax, calculated to control for background uptake, were not significantly different amongst the groups (1.9 ± 2.1 LC, 2.6 ± 2.3 HFLC, 3.6 ± 3.5 UR). In this small randomized controlled trial using UR diet as reference, LC dietary preparation followed by extended fasting resulted in significant myocardial uptake suppression

Complete mitochondrial DNA sequences provide new insights into the Polynesian motif and the peopling of Madagascar

More than a decade of mitochondrial DNA (mtDNA) studies have given the 'Polynesian motif' renowned status as a marker for tracing the late-Holocene expansion of Austronesian speaking populations. Despite considerable research on the Polynesian motif in Oceania, there has been little equivalent work on the western edge of its expansion - leaving major issues unresolved regarding the motif's evolutionary history. This has also led to considerable uncertainty regarding the settlement of Madagascar. In this study, we assess mtDNA variation in 266 individuals from three Malagasy ethnic groups: the Mikea, Vezo, and Merina. Complete mtDNA genome sequencing reveals a new variant of the Polynesian motif in Madagascar; two coding region mutations define a Malagasy-specific sub-branch. This newly defined 'Malagasy motif' occurs at high frequency in all three ethnic groups (13-50%), and its phylogenetic position, geographic distribution, and estimated age all support a recent origin, but without conclusively identifying a specific source region. Nevertheless, the haplotype's limited diversity, similar to those of other mtDNA haplogroups found in our Malagasy groups, best supports a small number of initial settlers arriving to Madagascar through the same migratory process. Finally, the discovery of this lineage provides a set of new polymorphic positions to help localize the Austronesian ancestors of the Malagasy, as well as uncover the origin and evolution of the Polynesian motif itself

Quantifying similarity of pore-geometry in nanoporous materials

In most applications of nanoporous materials the pore structure is as important as the chemical composition as a determinant of performance. For example, one can alter performance in applications like carbon capture or methane storage by orders of magnitude by only modifying the pore structure. For these applications it is therefore important to identify the optimal pore geometry and use this information to find similar materials. However, the mathematical language and tools to identify materials with similar pore structures, but different composition, has been lacking. We develop a pore recognition approach to quantify similarity of pore structures and classify them using topological data analysis. This allows us to identify materials with similar pore geometries, and to screen for materials that are similar to given top-performing structures. Using methane storage as a case study, we also show that materials can be divided into topologically distinct classes requiring different optimization strategies

Complete Mitochondrial Genome Sequencing Reveals Novel Haplotypes in a Polynesian Population

The high risk of metabolic disease traits in Polynesians may be partly explained by elevated prevalence of genetic variants involved in energy metabolism. The genetics of Polynesian populations has been shaped by island hoping migration events which have possibly favoured thrifty genes. The aim of this study was to sequence the mitochondrial genome in a group of Maoris in an effort to characterise genome variation in this Polynesian population for use in future disease association studies. We sequenced the complete mitochondrial genomes of 20 non-admixed Maori subjects using Affymetrix technology. DNA diversity analyses showed the Maori group exhibited reduced mitochondrial genome diversity compared to other worldwide populations, which is consistent with historical bottleneck and founder effects. Global phylogenetic analysis positioned these Maori subjects specifically within mitochondrial haplogroup - B4a1a1. Interestingly, we identified several novel variants that collectively form new and unique Maori motifs – B4a1a1c, B4a1a1a3 and B4a1a1a5. Compared to ancestral populations we observed an increased frequency of non-synonymous coding variants of several mitochondrial genes in the Maori group, which may be a result of positive selection and/or genetic drift effects. In conclusion, this study reports the first complete mitochondrial genome sequence data for a Maori population. Overall, these new data reveal novel mitochondrial genome signatures in this Polynesian population and enhance the phylogenetic picture of maternal ancestry in Oceania. The increased frequency of several mitochondrial coding variants makes them good candidates for future studies aimed at assessment of metabolic disease risk in Polynesian populations

Chemogenetic fingerprinting by analysis of cellular growth dynamics

<p>Abstract</p> <p>Background</p> <p>A fundamental goal in chemical biology is the elucidation of on- and off-target effects of drugs and biocides. To this aim chemogenetic screens that quantify drug induced changes in cellular fitness, typically taken as changes in composite growth, is commonly applied.</p> <p>Results</p> <p>Using the model organism <it>Saccharomyces cerevisiae </it>we here report that resolving cellular growth dynamics into its individual components, growth lag, growth rate and growth efficiency, increases the predictive power of chemogenetic screens. Both in terms of drug-drug and gene-drug interactions did the individual growth variables capture distinct and only partially overlapping aspects of cell physiology. In fact, the impact on cellular growth dynamics represented functionally distinct chemical fingerprints.</p> <p>Discussion</p> <p>Our findings suggest that the resolution and quantification of all facets of growth increases the informational and interpretational output of chemogenetic screening. Hence, by facilitating a physiologically more complete analysis of gene-drug and drug-drug interactions the here reported results may simplify the assignment of mode-of-action to orphan bioactive compounds.</p

Novel, synergistic antifungal combinations that target translation fidelity

There is an unmet need for new antifungal or fungicide treatments, as resistance to existing treatments grows. Combination treatments help to combat resistance. Here we develop a novel, effective target for combination antifungal therapy. Different aminoglycoside antibiotics combined with different sulphate-transport inhibitors produced strong, synergistic growth-inhibition of several fungi. Combinations decreased the respective MICs by ≥8 fold. Synergy was suppressed in yeast mutants resistant to effects of sulphate-mimetics (like chromate or molybdate) on sulphate transport. By different mechanisms, aminoglycosides and inhibition of sulphate transport cause errors in mRNA translation. The mistranslation rate was stimulated up to 10-fold when the agents were used in combination, consistent with this being the mode of synergistic action. A range of undesirable fungi were susceptible to synergistic inhibition by the combinations, including the human pathogens Candida albicans, C. glabrata and Cryptococcus neoformans, the food spoilage organism Zygosaccharomyces bailii and the phytopathogens Rhizoctonia solani and Zymoseptoria tritici. There was some specificity as certain fungi were unaffected. There was no synergy against bacterial or mammalian cells. The results indicate that translation fidelity is a promising new target for combinatorial treatment of undesirable fungi, the combinations requiring substantially decreased doses of active components compared to each agent alone

Noisy-threshold control of cell death

<p>Abstract</p> <p>Background</p> <p>Cellular responses to death-promoting stimuli typically proceed through a differentiated multistage process, involving a lag phase, extensive death, and potential adaptation. Deregulation of this chain of events is at the root of many diseases. Improper adaptation is particularly important because it allows cell sub-populations to survive even in the continuous presence of death conditions, which results, among others, in the eventual failure of many targeted anticancer therapies.</p> <p>Results</p> <p>Here, I show that these typical responses arise naturally from the interplay of intracellular variability with a threshold-based control mechanism that detects cellular changes in addition to just the cellular state itself. Implementation of this mechanism in a quantitative model for T-cell apoptosis, a prototypical example of programmed cell death, captures with exceptional accuracy experimental observations for different expression levels of the oncogene Bcl-x_Land directly links adaptation with noise in an ATP threshold below which cells die.</p> <p>Conclusions</p> <p>These results indicate that oncogenes like Bcl-x_L, besides regulating absolute death values, can have a novel role as active controllers of cell-cell variability and the extent of adaptation.</p

In vivo efficacy of XR9051, a potent modulator of P-glycoprotein mediated multidrug resistance

Overexpression of P-glycoprotein (P-gp) is a potential cause of multidrug resistance (MDR) in tumours. We have previously reported that XR9051 (N-(4-(2-(6,7-dimethoxy-1,2,3,4-tetrahydro-2-isoquinolyl)ethyl)phenyl)-3-((3Z,6Z)-6-benzylidene-1-methyl-2,5-dioxo-3-piperazinylidene)methylbenzamide) is a potent and specific inhibitor of P-gp, which reverses drug resistance in several murine and human MDR cell lines. In this study we have evaluated the in vivo efficacy of this novel modulator in a panel of murine and human tumour models and examined its pharmacokinetic profile. Efficacy studies in mice bearing MDR syngeneic tumours (P388/DX Johnson, MC26) or human tumour xenografts (A2780AD, CH1/DOXr, H69/LX) demonstrated that co-administration of XR9051 significantly potentiated the anti-tumour activity of a range of cytotoxic drugs. This modulatory activity was observed following parenteral and oral co-administration of XR9051. In addition, the combination schedules were well-tolerated. Following intravenous administration in mice, XR9051 is rapidly distributed and accumulates in tumours and other tissues. In addition, the compound is well-absorbed after oral administration. These data suggest that XR9051 has the potential for reversing clinical MDR mediated by P-glycoprotien. © 1999 Cancer Research Campaig

Melanesian mtDNA Complexity

Melanesian populations are known for their diversity, but it has been hard to grasp the pattern of the variation or its underlying dynamic. Using 1,223 mitochondrial DNA (mtDNA) sequences from hypervariable regions 1 and 2 (HVR1 and HVR2) from 32 populations, we found the among-group variation is structured by island, island size, and also by language affiliation. The more isolated inland Papuan-speaking groups on the largest islands have the greatest distinctions, while shore dwelling populations are considerably less diverse (at the same time, within-group haplotype diversity is less in the most isolated groups). Persistent differences between shore and inland groups in effective population sizes and marital migration rates probably cause these differences. We also add 16 whole sequences to the Melanesian mtDNA phylogenies. We identify the likely origins of a number of the haplogroups and ancient branches in specific islands, point to some ancient mtDNA connections between Near Oceania and Australia, and show additional Holocene connections between Island Southeast Asia/Taiwan and Island Melanesia with branches of haplogroup E. Coalescence estimates based on synonymous transitions in the coding region suggest an initial settlement and expansion in the region at ∼30–50,000 years before present (YBP), and a second important expansion from Island Southeast Asia/Taiwan during the interval ∼3,500–8,000 YBP. However, there are some important variance components in molecular dating that have been overlooked, and the specific nature of ancestral (maternal) Austronesian influence in this region remains unresolved