Many important physiological processes operate at time and space scales far beyond those accessible to atom-realistic simulations, and yet discrete stochastic rather than continuum methods may best represent finite numbers of molecules interacting in complex cellular spaces. We describe and validate new tools and algorithms developed for a new version of the MCell simulation program (MCell3), which supports generalized Monte Carlo modeling of diffusion and chemical reaction in solution, on surfaces representing membranes, and combinations thereof. A new syntax for describing the spatial directionality of surface reactions is introduced, along with optimizations and algorithms that can substantially reduce computational costs (e.g., event scheduling, variable time and space steps). Examples for simple reactions in simple spaces are validated by comparison to analytic solutions. Thus we show how spatially realistic Monte Carlo simulations of biological systems can be far more cost-effective than often is assumed, and provide a level of accuracy and insight beyond that of continuum methods.

Most of the cellular ATP in living organisms is synthesized by the enzyme F(1)F(o)-ATP synthase. The water soluble F(1) part of the enzyme can also work in reverse and utilize the chemical energy released during ATP hydrolysis to generate mechanical motion. Despite the availability of a large amount of biochemical data and several x-ray crystallographic structures of F(1), there still remains a considerable lack of understanding as to how this protein efficiently converts the chemical energy released during the reaction ATP + H(2)O --> ADP + P(i) into mechanical motion of the stalk. We report here an ab initio QM/MM study of ATP hydrolysis in the beta(TP) catalytic site of F(1). Our simulations provide an atomic level description of the reaction path, its energetics, and the interaction of the nucleotide with the protein environment during catalysis. The simulations suggest that the reaction path with the lowest potential energy barrier proceeds via nucleophilic attack on the gamma-phosphate involving two water molecules. Furthermore, the ATP hydrolysis reaction in beta(TP) is found to be endothermic, demonstrating that the catalytic site is able to support the synthesis of ATP and does not promote ATP hydrolysis in the particular conformation studied.

The enzyme F1-adenosine triphosphatase (ATPase) is a molecular motor that converts the chemical energy stored in the molecule adenosine triphosphate (ATP) into mechanical rotation of its gamma-subunit. During steady-state catalysis, the three catalytic sites of F1 operate in a cooperative fashion such that at every instant each site is in a different conformation corresponding to a different stage along the catalytic cycle. Notwithstanding a large amount of biochemical and, recently, structural data, we still lack an understanding of how ATP hydrolysis in F1 is coupled to mechanical motion and how the catalytic sites achieve cooperativity during rotatory catalysis. In this publication, we report combined quantum mechanical/molecular mechanical simulations of ATP hydrolysis in the betaTP and betaDP catalytic sites of F1-ATPase. Our simulations reveal a dramatic change in the reaction energetics from strongly endothermic in betaTP to approximately equienergetic in betaDP. The simulations identify the responsible protein residues, the arginine finger alphaR373 being the most important one. Similar to our earlier study of betaTP, we find a multicenter proton relay mechanism to be the energetically most favorable hydrolysis pathway. The results elucidate how cooperativity between catalytic sites might be achieved by this remarkable molecular motor.

LOV domains are the light-sensitive portion of plant phototropins. They absorb light through a flavin cofactor, photochemically form a covalent bond between the chromophore and a cysteine residue in the protein, and proceed to mediate activation of an attached kinase domain. Although the photoreaction itself is now well-characterized experimentally and computationally, it is still unclear how the formation of the adduct leads to kinase activation. We have performed molecular dynamics simulations on the LOV1 domain of Chlamydomonas reinhardtii and the LOV2 domain of Avena sativa, both before and after the photoreaction, to answer this question. The extensive simulations, over 240 ns in duration, reveal significant differences in how the LOV1 and LOV2 domains respond to photoactivation. The simulations indicate that LOV1 activation is likely caused by a change in hydrogen bonding between protein and ligand that destabilizes a highly conserved salt bridge, whereas LOV2 activation seems to result from a change in the flexibility of a set of protein loops. Results of electrostatics calculations, principal component analysis, sequence alignments, and root mean-square deviation analysis corroborate the above findings.

Plants use sophisticated photosensing mechanisms to maximize their utilization of the available sunlight and to control developmental processes. The plant blue-light receptors of the phot family mediate plant phototropism and contain two light, oxygen, and voltage (LOV) sensitive domains as photoactive elements. Here, we report combined quantum mechanical/molecular mechanical simulations of the photocycle of a complete Phot-LOV1 domain from C. reinhardtii. We have investigated the electronic properties and structural changes that follow blue-light absorption. This permitted us to characterize the pathway for flavin-cysteinyl adduct formation, which was found to proceed via a neutral radical state generated by hydrogen atom transfer from the reactive cysteine residue, Cys57, to the chromophore flavin mononucleotide. Interestingly, we find that adduct formation does not cause any larger scale conformational changes in Phot-LOV1 which suggests that dynamical effects mediate signal transmission following the initial photoexcitation event.

We used high-resolution fluorescence imaging and single-pixel optical fluctuation analysis to estimate the opening probability of individual voltage-gated calcium (Ca2+) channels during an action potential and the number of such Ca2+ channels within active zones of frog neuromuscular junctions. Analysis revealed ~36 Ca2+ channels within each active zone, similar to the number of docked synaptic vesicles but far less than the total number of transmembrane particles reported based on freeze-fracture analysis (~200–250). The probability that each channel opened during an action potential was only ~0.2. These results suggest why each active zone averages only one quantal release event during every other action potential, despite a substantial number of docked vesicles. With sparse Ca2+ channels and low opening probability, triggering of fusion for each vesicle is primarily controlled by Ca2+ influx through individual Ca2+ channels. In contrast, the entire synapse is highly reliable because it contains hundreds of active zones.

Despite decades of intense experimental studies, we still lack a detailed understanding of synaptic function. Fortunately, using computational approaches, we can obtain important new insights into the inner workings of these important neural systems. Here, we report the development of a spatially realistic computational model of an entire frog active zone in which we constrained model parameters with experimental data, and then used Monte Carlo simulation methods to predict the Ca(2+)-binding stoichiometry and dynamics that underlie neurotransmitter release. Our model reveals that 20-40 independent Ca(2+)-binding sites on synaptic vesicles, only a fraction of which need to bind Ca(2+) to trigger fusion, are sufficient to predict physiological release. Our excess-Ca(2+)-binding-site model has many functional advantages, agrees with recent data on synaptotagmin copy number, and is the first (to our knowledge) to link detailed physiological observations with the molecular machinery of Ca(2+)-triggered exocytosis. In addition, our model provides detailed microscopic insight into the underlying Ca(2+) dynamics during synapse activation.

The objective of the study was to investigate nerve ultrasound (US) in comparison to nerve conduction studies (NCS) for differential diagnosis of amyotrophic lateral sclerosis with predominant lower motoneuron disease(ALS/LMND) and multifocal motor neuropathy(MMN). A single-center, prospective, examiner-blinded cross-sectional diagnostic study in two cohorts was carried out. Cohort I: convenience sample of subjects diagnosed with ALS/LMND or MMN (minimal diagnostic criteria:possible ALS (revised EL-Escorial criteria), possible MMN (European Federation of Neurosciences guidelines).Cohort II: consecutive subjects with suspected diagnosis of either ALS/LMND or MMN. Diagnostic US and NCS models were developed based on ROC analysis of 28 different US and 32 different NCS values measured in cohort I. Main outcome criterion was sensitivity/specificity of these models between ALS/LMND and MMN in cohort II.Cohort I consisted of 16 patients with ALS/LMND and 8 patients with MMN. For cohort II, 30 patients were recruited, 8 with ALS/LMND, 5 with MMN, and 17 with other diseases. In cohort I, the three best US measures showed higher mean ± SD areas under the curve than the respective NCS measures (0.99 ± 0.01 vs. 0.79 ± 0.03, p<0.001; two-sided t test). The US model with highest measurement efficacy (8 values) and diagnostic quality reached 100 % sensitivity and 92 % specificity for MMN in cohort II, while the respective NCS model (6 values, including presence of conduction blocks) reached 100 and 52 %. Nerve US is of high diagnostic accuracy for differential diagnosis of ALS/LMND and MMN. It might be superior to NCS in the diagnosis of MMN in hospital-admitted patients with this differential diagnosis.