Technical note: GAMMORA,a free,open-source,and validated GATE-based model for Monte-Carlo simulations of the Varian TrueBeam

PurposeMonte Carlo (MC) is the reference computation method for medical physics. In radiotherapy, MC computations are necessary for some issues (such as assessing figures of merit, double checks, and dose conversions). A tool based on GATE is proposed to easily create full MC simulations of the Varian TrueBeam STx.MethodsGAMMORA is a package that contains photon phase spaces as a pre-trained generative adversarial network (GAN) and the TrueBeam’s full geometry. It allows users to easily create MC simulations for simple or complex radiotherapy plans such as VMAT. To validate the model, the characteristics of generated photons are first compared to those provided by Varian (IAEA format). Simulated data are also compared to measurements in water and heterogeneous media. Simulations of 8 SBRT plans are compared to measurements (in a phantom). Two examples of applications (a second check and interplay effect assessment) are presented.ResultsThe simulated photons generated by the GAN have the same characteristics (energy, position, and direction) as the IAEA data. Computed dose distributions of simple cases (in water) and complex plans delivered in a phantom are compared to measurements, and the Gamma index (3%/3mm) was always superior to 98%. The feasibility of both clinical applications is shown.ConclusionsThis model is now shared as a free and open-source tool that generates radiotherapy MC simulations. It has been validated and used for five years. Several applications can be envisaged for research and clinical purposes.     

Some ways of looking at compensatory kosmotropes and different water environments

Hydration of macromolecular structures determines biological activity. Stabilizing solutes are kosmotropic (increase order of water) rather than chaotropic (decrease order). Preferential hydration of surfaces is a thermodynamic consequence of the solution behavior of kosmotropic solutes, but inconsistencies imply interactions such as the hydration of specific sites within macromolecules. Thermodynamic measures require bulk pure solutes; here simpler measures of the effects on bulk water, water at surfaces and hydration water of probes have been applied to solutes including natural stabilizers, analogues and example chaotropes. Changes in the near-infrared spectra, water proton NMR chemical shifts and relaxation times measure changes in the bulk liquid; HPLC-column retention of solutes indicate interactions with hydration water at different surfaces, and fluorescence probes detect effects on functional group hydration water. Ab initio calculations and Monte-Carlo simulations of the solutes in water measure the energetics of the solute-water interactions, the dipole moments of these molecules, their charge distributions and the effect of the solute molecules on the structure of water. The rankings of the test solutes by these measures are not consistent. Thus, stabilizing solutes are not interchangeable in biological systems and the intracellular replacement of one by another could affect the integration of cell metabolism.     

PKQuest: volatile solutes - application to enflurane,nitrous oxide,halothane, methoxyflurane and toluene pharmacokinetics

BACKGROUND: The application of physiologically based pharmacokinetic models (PBPK) to human studies has been limited by the lack of the detailed organ information that is required for this analysis. PKQuest is a new generic PBPK that is designed to avoid this problem by using a set of "standard human" default parameters that are applicable to most solutes. RESULTS: PKQuest is used to model the human pharmacokinetics of the volatile solutes. A "standard human" value for the lipid content of the blood and each organ (klip) was chosen. This set of klip and the oil/water partition coefficient then specifies the organ/blood partition for each organ. Using this approach, the pharmacokinetics of inert volatile solute is completely specified by just 2 parameters: the water/air and oil/water partition coefficients. The model predictions of PKQuest were in good agreement with the experimental data for the inert solutes enflurane and nitrous oxide and the metabolized solutes halothane and toluene. METHODS: The experimental data that was modeled was taken from previous publications. CONCLUSIONS: This approach greatly increases the predictive power of the PBPK. For inert volatile solutes the pharmacokinetics are determined just from the water/air and oil/water partition coefficient. Methoxyflurane cannot be modeled by this PBPK because the arterial and end tidal partial pressures are not equal (as assumed in the PBPK). This inequality results from the "washin-washout" artifact in the large airways that is established for solutes with large water/air partition coefficients.PKQuest and the worked examples are available on the web www.pkquest.com.     

A cohesion/tension mechanism explains the gating of water channels (aquaporins) in Chara internodes by high concentration

Isolated internodes of Chara corallina have been used to study the gating of aquaporins (water channels) in the presence of high concentrations of osmotic solutes of different size (molecular weight). Osmolytes were acetone and three glycol ethers: ethylene glycol monomethyl ether (EGMME), diethylene glycol monomethyl ether (DEGMME), and triethylene glycol monoethyl ether (TEGMEE). The 'osmotic efficiency' of osmolytes was quite different. Their reflection coefficients ranged between 0.15 (acetone), 0.59 (EGMME), 0.78 (DEGMME), and 0.80 (TEGMEE). Bulk water permeability (Lp) and diffusive permeabilities (Ps) of heavy water (HDO), hydrogen peroxide (H2O2), acetone, and glycol ethers (EGMME, DEGMME, and TEGMEE) were measured using a cell pressure probe. Cells were treated with different concentrations of osmotic solutes of up to 800 mM ( approximately 2.0 MPa of osmotic pressure). Inhibition of aquaporin activity increased with both increasing concentration and size of solutes (reflection coefficients). As cell Lp decreased, Ps increased, indicating that water and solutes used different passages across the plasma membrane. Similar to earlier findings of an osmotic gating of ion channels, a cohesion/tension model of the gating of water channels in Chara internodes by high concentration is proposed. According to the model, tensions (negative pressures) within water channels affected the open/closed state by changing the free energy between states and favoured a distorted/collapsed rather than the open state. They should have differed depending on the concentration and size of solutes that are more or less excluded from aquaporins. The bigger the solute, the lower was the concentration required to induce a reversible closure of aquaporins, as predicted by the model.     

Hydration and conformational equilibria of simple hydrophobic and amphiphilic solutes.

We consider whether the continuum model of hydration optimized to reproduce vacuum-to-water transfer free energies simultaneously describes the hydration free energy contributions to conformational equilibria of the same solutes in water. To this end, transfer and conformational free energies of idealized hydrophobic and amphiphilic solutes in water are calculated from explicit water simulations and compared to continuum model predictions. As benchmark hydrophobic solutes, we examine the hydration of linear alkanes from methane through hexane. Amphiphilic solutes were created by adding a charge of +/-1e to a terminal methyl group of butane. We find that phenomenological continuum parameters fit to transfer free energies are significantly different from those fit to conformational free energies of our model solutes. This difference is attributed to continuum model parameters that depend on solute conformation in water, and leads to effective values for the free energy/surface area coefficient and Born radii that best describe conformational equilibrium. In light of these results, we believe that continuum models of hydration optimized to fit transfer free energies do not accurately capture the balance between hydrophobic and electrostatic contributions that determines the solute conformational state in aqueous solution.     

Hydration heat capacity of nucleic acid constituents determined from the random network model

The heat capacities of hydration (dCp) of the five nucleic acid bases A, G, C, T, and U, the sugars ribose and deoxyribose, and the phosphate backbone were determined using Monte Carlo simulations and the random network model. Solute-induced changes in the mean length and root mean square angle of hydrogen bonds between hydration shell waters were used to compute dCp for these solutes. For all solutes the dCp is significantly more positive than predicted from accessible surface area (ASA) models of heat capacity. In ASA models, nitrogen, oxygen, and phosphorus atoms are considered as uniformly polar, therefore making a negative contribution to dCp. However, the simulations show that many of these polar atoms are hydrated by water whose hydrogen bonds are less distorted than in bulk, leading to a positive dCp. This is in contrast to the effect of polar groups seen previously in small molecules and amino acids, which increase the water H-bond distortion, giving negative dCp contributions. Our results imply that dCp accompanying DNA dehydration in DNA-ligand and DNA-protein binding reactions may be significantly more negative than previously believed and that dehydration is a significant contributor to the large decrease in heat capacity seen in experiments.     

Nitrogen Forms Influence Microcystin Concentration and Composition via Changes in Cyanobacterial Community Structure

The eutrophication of freshwaters is a global health concern as lakes with excess nutrients are often subject to toxic cyanobacterial blooms. Although phosphorus is considered the main element regulating cyanobacterial biomass, nitrogen (N) concentration and more specifically the availability of different N forms may influence the overall toxicity of blooms. In this study of three eutrophic lakes prone to cyanobacterial blooms, we examined the effects of nitrogen species and concentrations and other environmental factors in influencing cyanobacterial community structure, microcystin (MC) concentrations and MC congener composition. The identification of specific MC congeners was of particular interest as they vary widely in toxicity. Different nitrogen forms appeared to influence cyanobacterial community structure leading to corresponding effects on MC concentrations and composition. Total MC concentrations across the lakes were largely explained by a combination of abiotic factors: dissolved organic nitrogen, water temperature and ammonium, but Microcystis spp. biomass was overall the best predictor of MC concentrations. Environmental factors did not appear to affect MC congener composition directly but there were significant associations between specific MC congeners and particular species. Based on redundancy analyses (RDA), the relative biomass of Microcystis aeruginosa was associated with MC-RR, M. wesenbergii with MC-LA and Aphanizomenon flos-aquae with MC-YR. The latter two species are not generally considered capable of MC production. Total nitrogen, water temperature, ammonium and dissolved organic nitrogen influenced the cyanobacterial community structure, which in turn resulted in differences in the dominant MC congener and the overall toxicity.     

Why do water molecules around small hydrophobic solutes form stronger hydrogen bonds than in the bulk?

Molecular solutes are known to have a strong effect on the structural and dynamical properties of the surrounding water. In our recent study (PNAS, 114, 322 (2017)) we have identified the presence of strengthened water hydrogen bonds near hydrophobic solutes by using both IR spectroscopy and ab-initio molecular dynamics simulations. The water molecules involved in the enhanced hydrogen bonding have been shown to display extensive structural ordering and restricted mobility. We observed that an individual pair of water molecules can make stronger hydrogen bond to each other if it is not surrounded by intercalating water molecules. Here we present compelling simulation results which unravel a simple mechanistic picture of the emergence of the hydrogen bond (HB) strengthening around solvated methane. We show explicitly that actual absence of water molecules within the excluded volume due to the hydrophobic molecule assures smaller residual torque on neighboring water molecules enabling the formation of stronger HBs between them.     

Local Partition Coefficients Govern Solute Permeability of Cholesterol-Containing Membranes

The permeability of lipid membranes for metabolic molecules or drugs is routinely estimated from the solute’s oil/water partition coefficient. However, the molecular determinants that modulate the permeability in different lipid compositions have remained unclear. Here, we combine scanning electrochemical microscopy and molecular-dynamics simulations to study the effect of cholesterol on membrane permeability, because cholesterol is abundant in all animal membranes. The permeability of membranes from natural lipid mixtures to both hydrophilic and hydrophobic solutes monotonously decreases with cholesterol concentration [Chol]. The same is true for hydrophilic solutes and planar bilayers composed of dioleoyl-phosphatidylcholine or dioleoyl-phosphatidyl-ethanolamine. However, these synthetic lipids give rise to a bell-shaped dependence of membrane permeability on [Chol] for very hydrophobic solutes. The simulations indicate that cholesterol does not affect the diffusion constant inside the membrane. Instead, local partition coefficients at the lipid headgroups and at the lipid tails are modulated oppositely by cholesterol, explaining the experimental findings. Structurally, these modulations are induced by looser packing at the lipid headgroups and tighter packing at the tails upon the addition of cholesterol.     

System-forming functions of bound water in the mechanism of topochemical reactions of formation of ultrathin layers on water surface

The system-forming role of the surface layer of water in the process of formation of ultrathin layers from the epoxide oligomer and triethylenetetramine or phosphotungstic acid on the surface of the aqueous phase is shown. Previously, this process was studied experimentally. The surface layer of the aqueous phase plays the role of a matrix, on which an epoxide oligomer monolayer and triethylenetetramine or phosphotungstic acid molecules are immobilized by hydrogen-bonding with water molecules from above or from below this matrix, respectively. Thus, topochemical reactions between the epoxide oligomer and triethylenetetramine or between the epoxide oligomer molecules in the presence of phosphotungstic acid become possible, so that ultrathin network epoxide-triethylenetetramine or epoxide layers are formed on the surface of water.     

Local Partition Coefficients Govern Solute Permeability of Cholesterol-Containing Membranes

The permeability of lipid membranes for metabolic molecules or drugs is routinely estimated from the solute’s oil/water partition coefficient. However, the molecular determinants that modulate the permeability in different lipid compositions have remained unclear. Here, we combine scanning electrochemical microscopy and molecular-dynamics simulations to study the effect of cholesterol on membrane permeability, because cholesterol is abundant in all animal membranes. The permeability of membranes from natural lipid mixtures to both hydrophilic and hydrophobic solutes monotonously decreases with cholesterol concentration [Chol]. The same is true for hydrophilic solutes and planar bilayers composed of dioleoyl-phosphatidylcholine or dioleoyl-phosphatidyl-ethanolamine. However, these synthetic lipids give rise to a bell-shaped dependence of membrane permeability on [Chol] for very hydrophobic solutes. The simulations indicate that cholesterol does not affect the diffusion constant inside the membrane. Instead, local partition coefficients at the lipid headgroups and at the lipid tails are modulated oppositely by cholesterol, explaining the experimental findings. Structurally, these modulations are induced by looser packing at the lipid headgroups and tighter packing at the tails upon the addition of cholesterol.     

Molecular dynamics simulation of amyloid beta dimer formation

Recent experiments with amyloid beta (Abeta) peptide indicate that formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Abeta oligomers depends on their structure, which is governed by assembly dynamics. Due to limitations of current experimental techniques, a detailed knowledge of oligomer structure at the atomic level is missing. We introduce a molecular dynamics approach to study Abeta dimer formation. 1), We use discrete molecular dynamics simulations of a coarse-grained model to identify a variety of dimer conformations; and 2), we employ all-atom molecular mechanics simulations to estimate thermodynamic stability of all dimer conformations. Our simulations of a coarse-grained Abeta peptide model predicts 10 different planar beta-strand dimer conformations. We then estimate the free energies of all dimer conformations in all-atom molecular mechanics simulations with explicit water. We compare the free energies of Abeta(1-42) and Abeta(1-40) dimers. We find that 1), dimer conformations have higher free energies compared to their corresponding monomeric states; and 2), the free-energy difference between the Abeta(1-42) and the corresponding Abeta(1-40) dimer conformation is not significant. Our results suggest that Abeta oligomerization is not accompanied by the formation of thermodynamically stable planar beta-strand dimers.     

Origins of protein denatured state compactness and hydrophobic clustering in aqueous urea: inferences from nonpolar potentials of mean force

Free energies of pairwise hydrophobic association are simulated in aqueous solutions of urea at concentrations ranging from 0-8 M. Consistent with the expectation that hydrophobic interactions are weakened by urea, the association of relatively large nonpolar solutes is destabilized by urea. However, the association of two small methane-sized nonpolar solutes in water has the opposite tendency of being slightly strengthened by the addition of urea. Such size effects and the dependence of urea-induced stability changes on the configuration of nonpolar solutes are not predicted by solvent accessible surface area approaches based on energetic parameters derived from bulk-phase solubilities of model compounds. Thus, to understand hydrophobic interactions in proteins, it is not sufficient to rely solely on transfer experiment data that effectively characterize a single nonpolar solute in an aqueous environment but not the solvent-mediated interactions among two or more nonpolar solutes. We find that the m-values for the rate of change of two-methane association free energy with respect to urea concentration is a dramatically nonmonotonic function of the spatial separation between the two methanes, with a distance-dependent profile similar to the corresponding two-methane heat capacity of association in pure water. Our results rationalize the persistence of residual hydrophobic contacts in some proteins at high urea concentrations and explain why the heat capacity signature (DeltaC(P)) of a compact denatured state can be similar to DeltaC(P) values calculated by assuming an open random-coil-like unfolded state.     

Cellular-automata models of solid-liquid interfaces

A series of cellular-automata (CA) models have been created, simulating relationships between water (or aqueous solutions) and solid surfaces of differing hydropathic (i.e., hydrophilic or hydrophobic) nature. Both equilibrium- and dynamic-flow models were examined, employing simple breaking and joining rules to simulate the hydropathic interactions. The CA simulations show that water accumulates near hydrophilic surfaces and avoids hydrophobic surfaces, forming concave-up and concave-down meniscuses, resp., under equilibrium conditions. In the dynamic-flow simulations, the flow rate of water was found to increase past a wall surface as the surface became less hydrophilic, reaching a maximum rate when the solid surface was of intermediate hydropathic state, and then declining with further increase in the hydrophobicity of the surface. Solution simulations show that non-polar solutes tend to achieve higher concentrations near hydrophobic-wall surfaces, whereas other hydrophobic/hydrophilic combinations of solutes and surfaces do not show such accumulations. Physical interpretations of the results are presented, as are some possible biological consequences.     

32S Tubulin oligomer. The resistance to factors suppressing the microtubule formation]

The resistance of brain 32S tubulin oligomer to factors suppressing the microtubules formation: colchicine, CaCl2, cooling and the absence of GTP is studied. The content of oligomer in the preparation and the polymerization degree were estimated by means of analytical centrifugation. Colchicine at 25 degrees and at a concentration of 10 muM does not change and at a concentration of 100 muM only slightly decreases the content of oligomer. The oligomer dissociated (but not completely) in 1 mM colchicine. Tubulin polymerization was partly suppressed by 10 muM and completely--by 100 muM colchicine. CaCl2 at 1 and 10 mM concentrations did not destroy the oligomer but inhibited its polymerization even in lesser of these concentrations. The cooling of the incubation medium to 14 degrees C or 4 degrees C completely inhibited the polymerization and did not affect the content of oligomer in the preparations. Tbulin preparations with low amount of exogenous GTP (less than or equal 3.10(-6) M) had a usual oligomer content, whereas GTP is necessary for polymerization at concentrations exceeding 10(-4) M. Thus, the reaction of tubulin oligomerization is relatively resistant to factors preventing the microtubules assembly. This probably means that there are at least two types of intereaction between tubulin molecules: 1) bonds in microtubules which are sensitive to colchicine, Ca2+ and cold, and which are formed only in the presence of nucleosidetriphosphates; 2) bonds in 32S tubulin oligomer which are more stable and do not need in exogenous nucleotides.     

Identification of a high-affinity glycine betaine transport system in Staphylococcus aureus.

Staphylococcus aureus accumulates proline and glycine betaine when cells are grown at low water activity. In the present study, we have identified a high-affinity glycine betaine transport system in this bacterium. Optimal activity for this transport system was measured in the presence of high NaCl concentrations, but transport activity was not stimulated by high concentrations of other solutes.     

The membrane action of antidiuretic hormone (ADH) on toad urinary bladder.

Radioactive tracer and electrical techniques were used to study the transport of nonelectrolytes and sodium, respectively, across toad urinary bladders in the presence and absence of ADH. The permeability of lipophilic molecules was roughly proportional to bulk phase oil/water partition coefficients both in the presence and absence of hormone; i.e., ADH elicited a general nonselective increase in the permeation of all nine solutes tested. The branched nonelectrolyte, isobutyramide, was less permeable than its straight-chain isomer, n-butyramide, in control tissues. ADH reduced the discrimination between these structural isomers. Hydrophilic solutes permeated more rapidly than expected. In the presence of hormone, there was no change in the permeation of large hydrophilic solutes considered to move via an extracellular pathway, but there was a marked increase in the permeability of water and other small hydrophilic solutes. Collectively, these results suggest that ADH acts to increase the motional freedom or fluidity of lipids in the cell membrane which is considered to be the preferred pathway for the permeation of lipophilic and small hydrophilic molecules. At concentrations of cAMP and ADH which elicit equivalent increments in the shortcircuit current, the effects of these agents on nonelectrolyte transport and membrane electrical conductance are divergent. Such observations suggest that some membrane effects of ADH may not be directly dependent upon cAMP. ADH in the mucosal solution increased the permeability of the toad bladder when the surface charge on the outer surface of the apical membrane was screened with the polyvalent cation, La-3+. These experiments emphasize that interaction of ADH with membranes of toad urinary bladder may account for at least some effects of this hormone.     

The stability of cylindrin β‐barrel amyloid oligomer models—A molecular dynamics study

Small‐soluble amyloid oligomers are believed to play a significant role in the pathology of amyloid diseases. Recently, the atomic structure of a toxic oligomer formed by an 11 residue and its tandem repeat was found to have an out‐off register antiparallel β‐strands in the shape of a β‐barrel. In the present article we investigate the effect of mutations in the hydrophobic cores on the structure and dynamic of the β‐barrels using all atom multiple molecular dynamics simulations with an explicit solvent. Extending previous experiments with molecular dynamics simulations we systematically test how stability and formation of cylindrin depends on the interplay between hydrophobicity and steric effects of the core residues. We find that strong hydrophobic interactions between geometrically fitting residues keep the strands (both in register and out‐off‐register interface) in close proximity, which in turn stabilizes the side‐chain and main‐chain hydrogen bonds, and the salt bridges on the outer surface along the weak out‐of‐register interface. Our simulations also indicate presence of water molecules in the hydrophobic interior of the cylindrin β‐barrel.Proteins 2013. © 2013 Wiley Periodicals, Inc.     

Further quantification of the role of internal unstirred layers during the measurement of transport coefficients in giant internodes of Chara by a new stop-flow technique

A new stop-flow technique was employed to quantify the impact of internal unstirred layers on the measurement of the solute permeability coefficient (P(s)) across the plasma membrane of internodes of the giant-celled alga Chara corallina using a cell pressure probe. During permeation experiments with rapidly permeating solutes (acetone, 2-propanol, and dimethylformamide), the solute concentration inside the cell was estimated and the external medium was adjusted to stop solute transport across the membrane, after which responses in turgor were measured. This allowed estimation of the solute concentration right at the membrane. Stop-flow experiments were also simulated with a computer. Both the stop-flow experiments and simulations provided quantitative data about internal concentration gradients and the contribution of unstirred layers to overall measured values of P(meas)(s) for the three solutes. The stop-flow experimental results agreed with stop-flow simulations assuming that solutes diffused into a completely stagnant cell interior. The effects of internal unstirred layers on the underestimation of membrane P(s) declined with decreasing P(s). They were no bigger than 37% in the presence of the most rapidly permeating solute, acetone (P(meas)(s) =4.2 x 10(-6) m s(-1)), and 14% for the less rapidly permeating dimethylformamide (P(meas)(s) =1.6x10(-6) m s(-1)). It is concluded that, even in the case of rapidly permeating solutes such as isotopic water and, even when making pessimistic assumptions about the internal mixing of solutes, an upper limit for the underestimation of P(s) due to internal unstirred layers was 37%. The data are discussed in terms of recent theoretical estimates of the effect of internal unstirred layers and in terms of some recent criticism of cell pressure probe measurements of water and solute transport coefficients. The current stop-flow data are in line with earlier estimations of the role of unstirred layers in the literature on cell water relations.