Trophic microRNA: Post-transcriptional regulation of target genes and larval development impairment in Plutella xylostella upon precursor and mature microRNA ingestion.

MicroRNAs (miRNAs) are post-transcriptional gene regulators. In the miRNA pathway's cytoplasmic part, the miRNA is processed from a hairpin-structured precursor to a double-stranded (ds) mature RNA and ultimately to a single-stranded mature miRNA. In insects, ingesting these two ds forms can regulate the target gene expression; this inspired the trophic miRNA's use as a functional genomics and pest management tool. However, systematic studies enabling comparisons of pre- and mature forms, dosages, administration times and instar-wise effects on target transcripts and phenotypes, which can help develop a miRNA administration method, are unavailable due to the different focuses of the previous investigations. We investigated the impact of trophically delivered Px-let-7 miRNA on the lepidopteran pest Plutella xylostella, to compare the efficacies of its pre- and ds-mature forms. Continuous feeding on the miRNA-supplemented diet suppressed expressions of FTZ-F1 and E74, the target ecdysone pathway genes. Both the pre-let-7 and mature let-7 miRNA forms similarly downregulated the target transcripts in all four larval instars. Pre-let-7 and let-7 ingestions decreased larval mass and instar duration and increased mortality in all instars, exhibiting adverse effects on larval growth and development. miRNA processing Dicer-1 and AGO-1's upregulations upon miRNA ingestion denoted the systemic miRNA spread in larval tissues. The scrambled sequence controls did not affect the target transcripts, suggesting the sequence-specific targeting by the mature miRNA and hairpin cassette's non-involvement in the target downregulation. This work provides a framework for miRNA and target gene function analyses and potentiates the trophic miRNA's utility in pest management.

Chemical ecology of Himalayan eggplant variety's antixenosis: identification of geraniol as an oviposition deterrent against the eggplant shoot and fruit borer.

Eggplant (Solanum melongena) suffers severe losses due to a multi-insecticide-resistant lepidopteran pest, shoot and fruit borer (SFB, Leucinodes orbonalis). Heavy and combinatorial application of pesticides for SFB control renders eggplant risky for human consumption. We observed that gravid SFB females do not oviposit on Himalayan eggplant variety RC-RL-22 (RL22). We hypothesized that RL22 contained an antixenosis factor. Females' behavior indicated that the RL22 cue they perceived was olfactory. To identify it, leaf volatile blends of seven eggplant varieties were profiled using solid phase microextraction and gas chromatography mass spectrometry. Seven RL22-specific compounds were detected in the plant headspace. In choice assays, oviposition deterrence efficacies of these candidate compounds were independently tested by their foliar application on SFB-susceptible varieties. Complementation of geraniol, which was exclusively found in RL22, reduced oviposition (> 90%). To validate geraniol's role in RL22's SFB-deterrence, we characterized RL22's geraniol synthase and silenced its gene in planta, using virus-induced gene silencing. Geraniol biosynthesis suppression rendered RL22 SFB-susceptible; foliar geraniol application on the geraniol synthase-silenced plants restored oviposition deterrence. We infer that geraniol is RL22's SFB oviposition deterrent. The use of natural compounds like geraniol, which influence the chemical ecology of oviposition, can reduce the load of hazardous synthetic larvicides.

High-Dimensional Parameter Search Method to Determine Force Field Mixing Terms in Molecular Simulations.

Molecular dynamics (MD) force fields for lipids and ions are typically developed independently of one another. In simulations consisting of both lipids and ions, lipid-ion interaction energies are estimated using a predefined set of mixing rules for Lennard-Jones (LJ) interactions. This, however, does not guarantee their reliability. In fact, compared to the quantum mechanical reference data, Lorentz-Berthelot mixing rules substantially underestimate the binding energies of Na+ ions with small-molecule analogues of lipid headgroups, yielding errors on the order of 80 and 130 kJ/mol, respectively, for methyl acetate and diethyl phosphate. Previously, errors associated with mixing force fields have been reduced using approaches such as "NB-fix" in which LJ interactions are computed using explicit cross terms rather than those from mixing rules. Building on this idea, we derive explicit lipid-ion cross terms that also may implicitly include many-body cooperativity effects. Additionally, to account for the interdependency between cross terms, we optimize all cross terms simultaneously by performing high-dimensional searches using our ParOpt software. The cross terms we obtain reduce the errors due to mixing rules to below 10 kJ/mol. MD simulation of the lipid bilayer conducted using these optimized cross terms resolves the structural discrepancies between our previous simulations and small-angle X-ray and neutron scattering experiments. These results demonstrate that simulations of lipid bilayers with ions that are accurate up to structural data from scattering experiments can be performed without explicit polarization terms. However, it is worth noting that such NB-fix cross terms are not based on any physical principle; a polarizable lipid model would be more realistic and is still desired. Our approach is generic and can be applied to improve the accuracies of simulations employing mixed force fields.

Inclusion of High-Field Target Data in AMOEBA's Calibration Improves Predictions of Protein-Ion Interactions.

The reliability of molecular mechanics simulations to predict effects of ion binding to proteins depends on their ability to simultaneously describe ion-protein, ion-water, and protein-water interactions. Force fields (FFs) to describe protein-water and ion-water interactions have been constructed carefully and have also been refined routinely to improve accuracy. Descriptions for ion-protein interactions have also been refined, although in an a posteriori manner through the use of "nonbonded-fix (NB-fix)" approaches in which parameters from default Lennard-Jones mixing rules are replaced with those optimized against some reference data. However, even after NB-fix corrections, there remains a significant need for improvement. This is also true for polarizable FFs that include self-consistent inducible moments. Our recent studies on the polarizable AMOEBA FF suggested that the problem associated with modeling ion-protein interactions could be alleviated by recalibrating polarization models of cation-coordinating functional groups so that they respond better to the high electric fields present near ions. Here, we present such a recalibration of carbonyls, carboxylates, and hydroxyls in the AMOEBA protein FF and report that it does improve predictions substantially─mean absolute errors in Na+-protein and K+-protein interaction energies decrease from 8.7 to 5.3 and 9.6 to 6.3 kcal/mol, respectively. Errors are computed with respect to estimates from van der Waals-inclusive density functional theory benchmarked against high-level quantum mechanical calculations and experiments. While recalibration does improve ion-protein interaction energies, they still remain underestimated, suggesting that further improvements can be made in a systematic manner through modifications in classical formalism. Nevertheless, we show that by applying our many-body NB-fix correction to Lennard-Jones components, these errors are further reduced to 2.7 and 2.6 kcal/mol, respectively, for Na+ and K+ ions. Finally, we show that the recalibrated AMOEBA protein FF retains its intrinsic reliability in predicting protein structure and dynamics in the condensed phase.

Molecular basis for higher affinity of SARS-CoV-2 spike RBD for human ACE2 receptor.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused substantially more infections, deaths, and economic disruptions than the 2002-2003 SARS-CoV. The key to understanding SARS-CoV-2's higher infectivity lies partly in its host receptor recognition mechanism. Experiments show that the human angiotensin converting enzyme 2 (ACE2) protein, which serves as the primary receptor for both CoVs, binds to the receptor binding domain (RBD) of CoV-2's spike protein stronger than SARS-CoV's spike RBD. The molecular basis for this difference in binding affinity, however, remains unexplained from X-ray structures. To go beyond insights gained from X-ray structures and investigate the role of thermal fluctuations in structure, we employ all-atom molecular dynamics simulations. Microseconds-long simulations reveal that while CoV and CoV-2 spike-ACE2 interfaces have similar conformational binding modes, CoV-2 spike interacts with ACE2 via a larger combinatorics of polar contacts, and on average, makes 45% more polar contacts. Correlation analysis and thermodynamic calculations indicate that these differences in the density and dynamics of polar contacts arise from differences in spatial arrangements of interfacial residues, and dynamical coupling between interfacial and non-interfacial residues. These results recommend that ongoing efforts to design spike-ACE2 peptide blockers will benefit from incorporating dynamical information as well as allosteric coupling effects.

Detoxification of plant defensive glucosinolates by an herbivorous caterpillar is beneficial to its endoparasitic wasp.

Plant chemical defences impact not only herbivores, but also organisms in higher trophic levels that prey on or parasitize herbivores. While herbivorous insects can often detoxify plant chemicals ingested from suitable host plants, how such detoxification affects endoparasitoids that use these herbivores as hosts is largely unknown. Here, we used transformed plants to experimentally manipulate the major detoxification reaction used by Plutella xylostella (diamondback moth) to deactivate the glucosinolate defences of its Brassicaceae host plants. We then assessed the developmental, metabolic, immune, and reproductive consequences of this genetic manipulation on the herbivore as well as its hymenopteran endoparasitoid Diadegma semiclausum. Inhibition of P. xylostella glucosinolate metabolism by plant-mediated RNA interference increased the accumulation of the principal glucosinolate activation products, the toxic isothiocyanates, in the herbivore, with negative effects on its growth. Although the endoparasitoid manipulated the excretion of toxins by its insect host to its own advantage, the inhibition of herbivore glucosinolate detoxification slowed endoparasitoid development, impaired its reproduction, and suppressed the expression of genes of a parasitoid-symbiotic polydnavirus that aids parasitism. Therefore, the detoxification of plant glucosinolates by an herbivore lowers its toxicity as a host and benefits the parasitoid D. semiclausum at multiple levels.

Transferable interactions of Li+ and Mg2+ ions in polarizable models.

Therapeutic implications of Li+, in many cases, stem from its ability to inhibit certain Mg2+-dependent enzymes, where it interacts with or substitutes for Mg2+. The underlying details of its action are, however, unknown. Molecular simulations can provide insights, but their reliability depends on how well they describe relative interactions of Li+ and Mg2+ with water and other biochemical groups. Here, we explore, benchmark, and recommend improvements to two simulation approaches: the one that employs an all-atom polarizable molecular mechanics (MM) model and the other that uses a hybrid quantum and MM implementation of the quasi-chemical theory (QCT). The strength of the former is that it describes thermal motions explicitly and that of the latter is that it derives local contributions from electron densities. Reference data are taken from the experiment, and also obtained systematically from CCSD(T) theory, followed by a benchmarked vdW-inclusive density functional theory. We find that the QCT model predicts relative hydration energies and structures in agreement with the experiment and without the need for additional parameterization. This implies that accurate descriptions of local interactions are essential. Consistent with this observation, recalibration of local interactions in the MM model, which reduces errors from 10.0 kcal/mol to 1.4 kcal/mol, also fixes aqueous phase properties. Finally, we show that ion-ligand transferability errors in the MM model can be reduced significantly from 10.3 kcal/mol to 1.2 kcal/mol by correcting the ligand's polarization term and by introducing Lennard-Jones cross-terms. In general, this work sets up systematic approaches to evaluate and improve molecular models of ions binding to proteins.

Interactions of Monovalent and Divalent Cations at Palmitoyl-Oleoyl-Phosphatidylcholine Interface.

Li+ is a biologically active and medically important cation. Experiments show that Li+ modulates some phospholipid bilayer properties in a manner similar to divalent cations, rather than other monovalent cations. We previously performed a comparative simulation study of the interaction of several monovalent cations with palmitoyl-oleoyl-phosphatidylcholine bilayers and reported that Li+ exhibited the highest association with lipids and formed a unique tetrahedral coordinated structure with lipid head groups. Here we extend these studies to two biologically important divalent cations, Mg2+ and Ca2+, and observe that, just like monovalent cations, Mg2+ and Ca2+ reduce bilayer areas and increase chain order. Bilayer area changes induced by cations are strongly correlated with the amount of charge inside the headgroup region; however, Mg2+ and Li+ are clear outliers. At the same time though, Mg2+ adsorption in the bilayer is the smallest among all cations, which is in contrast to Li+ that binds strongly to lipids. In fact, in contrast to all other cations, Mg2+ remains fully hydrated in the lipid headgroup region. However, Li+ and Mg2+ share high overlap between their inner-shell coordination topologies. This suggests that Li+ can structurally replace Mg2+, which is bound to other biomolecules with up to fourfold coordination, provided such replacement is energetically feasible. We compute structural topologies and compare them quantitatively using a new weighted-graphs-based method. Finally, we find that the specificity of cation interaction with lipid head groups exhibit consistent trend with the solvation shell energetics of ions in lipid headgroup and bulk water regions.

Molecular dynamics simulations of ether- and ester-linked phospholipids.

Dissimilarities in the bulk structure of bilayers composed of ether- vs ester-linked lipids are well-established; however, the atomistic interactions responsible for these differences are not well known. These differences are important in understanding of why archaea have a different bilayer composition than the other domains of life and why humans have larger concentrations of plasmalogens in specialized membranes? In this paper, we simulate two lipid bilayers, the ester linked dipalmitoylphosphatidylcholine (DPPC) and the ether lined dihexadecylphosphatidylcholine (DHPC), to study these variations. The structural analysis of the bilayers reveals that DPPC is more compressible than DHPC. A closer examination of dipole potential shows DHPC, despite having a smaller dipole potential of the bilayer, has a higher potential barrier than DPPC at the surface. Analysis of water order and dynamics suggests DHPC has a more ordered, less mobile layer of water in the headgroup. These results seem to resolve the issue as to whether the decrease in permeability of DHPC is due to of differences in minimum area per lipid (A0) or diffusion coefficient of water in the headgroup region (Dhead) (Guler et al., 2009) since we have shown significant changes in the order and mobility of water in that region.

Towards a Unified Understanding of Lithium Action in Basic Biology and its Significance for Applied Biology.

Lithium has literally been everywhere forever, since it is one of the three elements created in the Big Bang. Lithium concentration in rocks, soil, and fresh water is highly variable from place to place, and has varied widely in specific regions over evolutionary and geologic time. The biological effects of lithium are many and varied. Based on experiments in which animals are deprived of lithium, lithium is an essential nutrient. At the other extreme, at lithium ingestion sufficient to raise blood concentration significantly over 1 mM/, lithium is acutely toxic. There is no consensus regarding optimum levels of lithium intake for populations or individuals-with the single exception that lithium is a generally accepted first-line therapy for bipolar disorder, and specific dosage guidelines for sufferers of that condition are generally agreed on. Epidemiological evidence correlating various markers of social dysfunction and disease vs. lithium level in drinking water suggest benefits of moderately elevated lithium compared to average levels of lithium intake. In contrast to other biologically significant ions, lithium is unusual in not having its concentration in fluids of multicellular animals closely regulated. For hydrogen ions, sodium ions, potassium ions, calcium ions, chloride ions, and magnesium ions, blood and extracellular fluid concentrations are closely and necessarily regulated by systems of highly selective channels, and primary and secondary active transporters. Lithium, while having strong biological activity, is tolerated over body fluid concentrations ranging over many orders of magnitude. The lack of biological regulation of lithium appears due to lack of lithium-specific binding sites and selectivity filters. Rather lithium exerts its myriad physiological and biochemical effects by competing for macromolecular sites that are relatively specific for other cations, most especially for sodium and magnesium. This review will consider what is known about the nature of this competition and suggest using and extending this knowledge towards the goal of a unified understanding of lithium in biology and the application of that understanding in medicine and nutrition.

Effects of Lithium and Other Monovalent Ions on Palmitoyl Oleoyl Phosphatidylcholine Bilayer.

Interactions of monovalent salts with lipid membranes are explored with molecular dynamics (MD) simulations. The simulations included the monovalent ions Na+ and K+, for their importance in physiology, Li+ for its small size and importance in several medical conditions including bipolar disorder, and Rb+ for its large size. All simulations included Cl- as counterions. One bilayer was simulated without salt as a control. Palmitoyl oleoyl phosphatidylcholine (POPC) bilayers experienced reductions in area per lipid with the addition of salt; the smaller the ion the smaller the area, with the exception of Li+. Li+ exhibited unique binding affinities between phosphates and sn-2 carbonyls that lowered the order of the top part of sn-2 chain, which increased the area per lipid, compared to other ionic simulations. Further, we observe that monovalent salts alter bilayer properties through structural changes and not so much through the changes in surface potential.

Natural history-driven, plant-mediated RNAi-based study reveals CYP6B46's role in a nicotine-mediated antipredator herbivore defense.

Manduca sexta (Ms) larvae are known to efficiently excrete ingested nicotine when feeding on their nicotine-producing native hostplant, Nicotiana attenuata. Here we describe how ingested nicotine is co-opted for larval defense by a unique mechanism. Plant-mediated RNAi was used to silence a midgut-expressed, nicotine-induced cytochrome P450 6B46 (CYP6B46) in larvae consuming transgenic N. attenuata plants producing MsCYP6B46 dsRNA. These and transgenic nicotine-deficient plants were planted into native habitats to study the phenotypes of larvae feeding on these plants and the behavior of their predators. The attack-behavior of a native wolf spider (Camptocosa parallela), a major nocturnal predator, provided the key to understanding MsCYP6B46's function: spiders clearly preferred CYP6B46-silenced larvae, just as they had preferred larvae fed nicotine-deficient plants. MsCYP6B46 redirects a small amount (0.65%) of ingested nicotine from the midgut into hemolymph, from which nicotine is exhaled through the spiracles as an antispider signal. CYP6B46-silenced larvae were more susceptible to spider-attack because they exhaled less nicotine because of lower hemolymph nicotine concentrations. CYP6B46-silenced larvae were impaired in distributing ingested nicotine from midgut to hemolymph, but not in the clearing of hemolymph nicotine or in the exhalation of nicotine from hemolymph. MsCYP6B46 could be a component of a previously hypothesized pump that converts nicotine to a short-lived, transportable, metabolite. Other predators, big-eyed bugs, and antlion larvae were insensitive to this defense. Thus, chemical defenses, too toxic to sequester, can be repurposed for defensive functions through respiration as a form of defensive halitosis, and predators can assist the functional elucidation of herbivore genes.

Atomically detailed lipid bilayer models for the interpretation of small angle neutron and X-ray scattering data.

We present a new atom density profile (ADP) model and a statistical approach for extracting structural characteristics of lipid bilayers from X-ray and neutron scattering data. Models for five lipids with varying head and tail chemical composition in the fluid phase, 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG), are optimized using a simplex based method to simultaneously reproduce both neutron and X-ray scattering data. Structural properties are determined using statistical analysis of multiple optimal model structures. The method and models presented make minimal assumptions regarding the atomic configuration, while taking into account the underlying physical properties of the system. The more general model and statistical approach yield data with well defined uncertainties, indicating the precision in determining density profiles, atomic locations, and bilayer structural characteristics. Resulting bilayer structures include regions exhibiting large conformational variation. Due to the increased detail in the model, the results demonstrate the possibility of a distinct hydration layer within the interfacial (backbone) region.

Why do chimpanzee males attack the females of neighboring communities?

Our closest nonhuman primate relatives, chimpanzees, engage in potentially lethal between-group conflict; this collective aggressive behavior shows parallels with human warfare. In some communities, chimpanzee males also severely attack and even kill females of the neighboring groups. This is surprising given their system of resource defense polygyny, where males are expected to acquire potential mates. We develop a simple mathematical model based on reproductive skew among primate males to solve this puzzle. The model predicts that it is advantageous for high-ranking males but not for low-ranking males to attack females. It also predicts that more males gain a benefit from attacking females as the community's reproductive skew decreases, i.e., as mating success is more evenly distributed. Thus, fatal attacks on females should be concentrated in communities with low reproductive skew. These attacks should also concur with between-community infanticide. A review of the chimpanzee literature provides enough preliminary support for this prediction to warrant more detailed testing.

Effectiveness of interventions for emergency care of hypoglycaemia and diabetic ketoacidosis: A systematic review.

AIM: This systematic review aims to provide evidence on effectiveness of interventions used in emergency care of hypoglycaemia and diabetic ketoacidosis (DKA). METHODOLOGY: This is a systematic review of randomized controlled trials and analytical studies. We selected studies based on eligibility criteria. The databases Medline, Cochrane library and Embase were searched from their inception till November 2, 2022, using search strategy. We used the term such as "diabetes mellitus", "treatment", "hypoglycaemia", "diabetic ketoacidosis", "low blood sugar", "high blood sugar" and Mesh terms like "disease management", "hypoglycaemia", "diabetic ketoacidosis", and "diabetes mellitus" to form search strategy. RESULTS: Hypoglycemia: Both 10 % dextrose (D10) and 50 % dextrose (D50) are effective options with similar hospital mortality D10 (4.7 %) and D50 (6.2 %). DKA: Low dose insulin is non-inferior to standard dose with time till resolution of DKA 16.5 (7.2) hours and 17.2 (7.7) hours (p value = 0.73) respectively. In children, subcutaneous insulin was associated with reduced ICU admissions and hospital readmissions (67.8 % to 27.9 %). Plasmalyte (PL) is noninferior to sodium chloride (SC), with ICU length of stay 49 h (IQR 23-72) and 55 h (IQR 41-80) respectively, hyperchloremia was associated with longer in-hospital length of stay and longer time to resolution of DKA. And potassium replacement at < 10 mmol/L was associated with higher mortality (n = 72). CONCLUSION: We conclude either of the 10 % or 50 % dextrose is effective for management of hypoglycaemia. For DKA subcutaneous insulin and intravenous insulin, chloride levels ≤ 109 mEq/L, potassium above 10 mmol/l, IV fluids like Plasmalyte and normal saline are effective.

Approximate Algorithms for Computing Spatial Distance Histograms with Accuracy Guarantees.

Particle simulation has become an important research tool in many scientific and engineering fields. Data generated by such simulations impose great challenges to database storage and query processing. One of the queries against particle simulation data, the spatial distance histogram (SDH) query, is the building block of many high-level analytics, and requires quadratic time to compute using a straightforward algorithm. Previous work has developed efficient algorithms that compute exact SDHs. While beating the naive solution, such algorithms are still not practical in processing SDH queries against large-scale simulation data. In this paper, we take a different path to tackle this problem by focusing on approximate algorithms with provable error bounds. We first present a solution derived from the aforementioned exact SDH algorithm, and this solution has running time that is unrelated to the system size N. We also develop a mathematical model to analyze the mechanism that leads to errors in the basic approximate algorithm. Our model provides insights on how the algorithm can be improved to achieve higher accuracy and efficiency. Such insights give rise to a new approximate algorithm with improved time/accuracy tradeoff. Experimental results confirm our analysis.

Mating strategies in primates: a game theoretical approach to infanticide.

Infanticide by newly immigrated or newly dominant males is reported among a variety of taxa, such as birds, rodents, carnivores and primates. Here we present a game theoretical model to explain the presence and prevalence of infanticide in primate groups. We have formulated a three-player game involving two males and one female and show that the strategies of infanticide on the males' part and polyandrous mating on the females' part emerge as Nash equilibria that are stable under certain conditions. Moreover, we have identified all the Nash equilibria of the game and arranged them in a novel hierarchical scheme. Only in the subspace spanned by the males are the Nash equilibria found to be strict, and hence evolutionarily stable. We have therefore proposed a selection mechanism informed by adaptive dynamics to permit the females to transition to, and remain in, optimal equilibria after successive generations. Our model concludes that polyandrous mating by females is an optimal strategy for the females that minimizes infanticide and that infanticide confers advantage to the males only in certain regions of parameter space. We have shown that infanticide occurs during turbulent changes accompanying male immigration into the group. For changes in the dominance hierarchy within the group, we have shown that infanticide occurs only in primate groups where the chance for the killer to sire the next infant is high. These conclusions are confirmed by observations in the wild. This model thus has enabled us to pinpoint the fundamental processes behind the reproductive decisions of the players involved, which was not possible using earlier theoretical studies.

Multiscale simulations of heterogeneous model membranes.

This review will focus on computer modeling aimed at providing insights into the existence, structure, size, and thermodynamic stability of localized domains in membranes of heterogeneous composition. Modeling the lateral organization within a membrane is problematic due to the relatively slow lateral diffusion rate for lipid molecules so that microsecond or longer time scales are needed to fully model the formation and stability of a raft in a membrane. Although atomistic simulations currently are not able to reach this scale, they can provide data on the intermolecular forces and correlations that are involved in lateral organization. These data can be used to define coarse grained models that are capable of predictions of lateral organization in membranes. In this paper, we review modeling efforts that use interaction data from MD simulations to construct coarse grained models for heterogeneous bilayers. In this review we will discuss MD simulations done with the aim of gaining the information needed to build accurate coarse-grained models. We will then review some of the coarse-graining work, emphasizing modeling that has resulted from or has a basis in atomistic simulations.

Lateral organization of complex lipid mixtures from multiscale modeling.

The organizational properties of complex lipid mixtures can give rise to functionally important structures in cell membranes. In model membranes, ternary lipid-cholesterol (CHOL) mixtures are often used as representative systems to investigate the formation and stabilization of localized structural domains ("rafts"). In this work, we describe a self-consistent mean-field model that builds on molecular dynamics simulations to incorporate multiple lipid components and to investigate the lateral organization of such mixtures. The model predictions reveal regions of bimodal order on ternary plots that are in good agreement with experiment. Specifically, we have applied the model to ternary mixtures composed of dioleoylphosphatidylcholine:18:0 sphingomyelin:CHOL. This work provides insight into the specific intermolecular interactions that drive the formation of localized domains in these mixtures. The model makes use of molecular dynamics simulations to extract interaction parameters and to provide chain configuration order parameter libraries.

A reactive molecular dynamics simulation of the silica-water interface.

We report our study of a silica-water interface using reactive molecular dynamics. This first-of-its-kind simulation achieves length and time scales required to investigate the detailed chemistry of the system. Our molecular dynamics approach is based on the ReaxFF force field of van Duin et al. [J. Phys. Chem. A 107, 3803 (2003)]. The specific ReaxFF implementation (SERIALREAX) and force fields are first validated on structural properties of pure silica and water systems. Chemical reactions between reactive water and dangling bonds on a freshly cut silica surface are analyzed by studying changing chemical composition at the interface. In our simulations, reactions involving silanol groups reach chemical equilibrium in approximately 250 ps. It is observed that water molecules penetrate a silica film through a proton-transfer process we call "hydrogen hopping," which is similar to the Grotthuss mechanism. In this process, hydrogen atoms pass through the film by associating and dissociating with oxygen atoms within bulk silica, as opposed to diffusion of intact water molecules. The effective diffusion constant for this process, taken to be that of hydrogen atoms within silica, is calculated to be 1.68 x 10(-6) cm(2)/s. Polarization of water molecules in proximity of the silica surface is also observed. The subsequent alignment of dipoles leads to an electric potential difference of approximately 10.5 V between the silica slab and water.