Restoration of native folding of single-stranded DNA sequences through reverse mutations: an indication of a new epigenetic mechanism

We used in vivo (biological), in silico (computational structure prediction), and in vitro (model sequence folding) analyses of single-stranded DNA sequences to show that nucleic acid folding conservation is the selective principle behind a high-frequency single-nucleotide reversion observed in a three-nucleotide mutated motif of the Maize streak virus replication associated protein (Rep) gene. In silico and in vitro studies showed that the three-nucleotide mutation adversely affected Rep nucleic acid folding, and that the single-nucleotide reversion [C(601)A] restored wild-type-like folding. In vivo support came from infecting maize with mutant viruses: those with Rep genes containing nucleotide changes predicted to restore a wild-type-like fold [A(601)/G(601)] preferentially accumulated over those predicted to fold differently [C(601)/T(601)], which frequently reverted to A(601) and displaced the original population. We propose that the selection of native nucleic acid folding is an epigenetic effect, which might have broad implications in the evolution of plants and their viruses.

A multi-scaled approach for simulating chemical reaction systems

In this paper we give an overview of some very recent work, as well as presenting a new approach, on the stochastic simulation of multi-scaled systems involving chemical reactions. In many biological systems (such as genetic regulation and cellular dynamics) there is a mix between small numbers of key regulatory proteins, and medium and large numbers of molecules. In addition, it is important to be able to follow the trajectories of individual molecules by taking proper account of the randomness inherent in such a system. We describe different types of simulation techniques (including the stochastic simulation algorithm, Poisson Runge-Kutta methods and the balanced Euler method) for treating simulations in the three different reaction regimes: slow, medium and fast. We then review some recent techniques on the treatment of coupled slow and fast reactions for stochastic chemical kinetics and present a new approach which couples the three regimes mentioned above. We then apply this approach to a biologically inspired problem involving the expression and activity of LacZ and LacY proteins in E coli, and conclude with a discussion on the significance of this work. (C) 2004 Elsevier Ltd. All rights reserved.

Contextual models and the non-Newtonian paradigm

Biological systems exhibit a wide range of contextual effects, and this often makes it difficult to construct valid mathematical models of their behaviour. In particular, mathematical paradigms built upon the successes of Newtonian physics make assumptions about the nature of biological systems that are unlikely to hold true. After discussing two of the key assumptions underlying the Newtonian paradigm, we discuss two key aspects of the formalism that extended it, Quantum Theory (QT). We draw attention to the similarities between biological and quantum systems, motivating the development of a similar formalism that can be applied to the modelling of biological processes.

A simplified approach for calculating the moments of action for linear reaction-diffusion equations

The mean action time is the mean of a probability density function that can be interpreted as a critical time, which is a finite estimate of the time taken for the transient solution of a reaction-diffusion equation to effectively reach steady state. For high-variance distributions, the mean action time under-approximates the critical time since it neglects to account for the spread about the mean. We can improve our estimate of the critical time by calculating the higher moments of the probability density function, called the moments of action, which provide additional information regarding the spread about the mean. Existing methods for calculating the nth moment of action require the solution of n nonhomogeneous boundary value problems which can be difficult and tedious to solve exactly. Here we present a simplified approach using Laplace transforms which allows us to calculate the nth moment of action without solving this family of boundary value problems and also without solving for the transient solution of the underlying reaction-diffusion problem. We demonstrate the generality of our method by calculating exact expressions for the moments of action for three problems from the biophysics literature. While the first problem we consider can be solved using existing methods, the second problem, which is readily solved using our approach, is intractable using previous techniques. The third problem illustrates how the Laplace transform approach can be used to study coupled linear reaction-diffusion equations.

Multiscale investigation of microfilament networks disruption in living cells

Actin is the most abundantly distributed protein in living cells which plays critical roles in the cell interior force generation and transmission. The fracture mechanism of microfilament networks, whose principle component is actin, would provide insights which can contribute to the understandings of self-protective characters of cytoskeleton. In this study, molecular simulations are conducted to investigate the molecular mechanisms of disruption of microfilament networks from the viewpoint of biophysics. By employing a coarse-grained (CG) model of actin filament networks, we focused on the ultimate strength and crack growth mode of microfilament networks that have dependency on the crack length. It can be found that, the fracture mechanism of microfilament network has dependency on the structural properties of microfilament networks. The structure flaws marginally change the strength of microfilament networks which would explain the self-protective characters of cytoskeleton.

Prospects for studying how high-intensity compression waves cause damage in human blast injuries

Since World War I, explosions have accounted for over 70% of all injuries in conflict. With the development of improved personnel protection of the torso, improved medical care and faster aeromedical evacuation, casualties are surviving with more severe injuries to the extremities. Understanding the processes involved in the transfer of blast-induced shock waves through biological tissues is essential for supporting efforts aimed at mitigating and treating blast injury. Given the inherent heterogeneities in the human body, we argue that studying these processes demands a highly integrated approach requiring expertise in shock physics, biomechanics and fundamental biological processes. This multidisciplinary systems approach enables one to develop the experimental framework for investigating the material properties of human tissues that are subjected to high compression waves in blast conditions and the fundamental cellular processes altered by this type of stimuli. Ultimately, we hope to use the information gained from these studies in translational research aimed at developing improved protection for those at risk and improved clinical outcomes for those who have been injured from a blast wave.

Molecular sliding filament model for muscular contraction based on multiscale investigation

A multiscale approach that bridges the biophysics of the actin molecules at nanoscale and the biomechanics of actin filament at microscale level is developed and used to evaluate the mechanical performances of actin filament bundles. In order to investigate the contractile properties of skeletal muscle which is induced by the protein motor of myosin, a molecular model is proposed in the prediction of the dynamic behaviors of skeletal muscle based on classic sliding filament model. Randomly distributed myosin motors are applied on a 2.2 μm long sarcomere, whose principal components include actin and myosin filaments. It can be found that, the more myosin motors on the sarcomere, the faster the sarcomere contracts. The result demonstrates that the sarcomere shortening speed cannot increase infinitely by the modulation of myosin, thus providing insight into the self-protective properties of skeletal muscles. This molecular filament sliding model provides a theoretical way to evaluate the properties of skeletal muscles, and contributes to the understandings of the molecular mechanisms in the physiological phenomenon of muscular contraction.

MAS NMR measurements of intact articular bovine cartilage

Articular cartilage (AC), an avascular connective tissue lining articulating surfaces of the long bones, comprises extracellular biopolymers. In functionally compromised states such as osteoarthritis, thinned or lost AC causes reduced mobility and increased health-care costs. Understanding of the characteristics responsible for the load bearing efficiency of AC and the factors leading to its degradation are incomplete. DTI shows the structural alignment of collagen in AC [1] and T2 relaxation measurements suggest that the average director of reorientational motion of water molecules depends on the degree of alignment of collagen in AC [2]. Information on the nature of the chemical interactions involved in functional AC is lacking. The need for AC structural integrity makes solid state NMR an ideal tool to study this tissue. We examined the contribution of water in different functional ‘compartments’ using 1H-MAS, 13C-MAS and 13C-CPMAS NMR of bovine patellar cartilage incubated in D2O. 1H-MAS spectra signal intensity was reduced due to H/D exchange without a measureable redistribution of relative signal intensity. Chemical shift anisotropy was estimated by lineshape analysis of multiple peaks in the 1H-MAS spinning sidebands. These asymmetrical sidebands suggested the presence of multiple water species in AC. Therefore, water was added in small aliquots to D2O saturated AC and the influence of H2O and D2O on organic components was studied with 13C-MAS-NMR and 13C-CPMAS-NMR. Signal intensity in 13C-MAS spectra showed no change in relative signal intensity throughout the spectrum. In 13C-CPMAS spectra, displacement of water by D2O resulted in a loss of signal in the aliphatic region due to a reduction in proton availability for cross-polarization. These results complement dehydration studies of cartilage using osmotic manipulation [3] and demonstrate components of cartilage that are in contact with mobile water.

Production of biodegradable nano and micro particles via ultrasonic atomization for biopharmaceutical delivery

Biopharmaceuticals have been shown to have low delivery and transformation efficiencies. To over come this, larger doses are administered in order to obtain the desired response which may lead to toxicity and drug resistance. This paper reports upon a continuous particle production method utilizing surface acoustic wave atomization to reliably produce micro and nanoparticles with physical characteristics to facilitate the cellular uptake of biopharmaceuticals. By producing particles of an optimal size for cellular uptake, the efficacy and specificity of drug loaded nanoparticles will be increased. Better delivery methods will result in dosage reduction (hence lower costs per dose), reduced toxicity, and reduced problems associated with multidrug resistance due to over dosing.

SecB binds only to a late native-like intermediate in the folding pathway of barstar and not to the unfolded state

SecB is a cytosolic, tetrameric chaperone of Escherichia coli which maintains precursor proteins in a translocation competent state. We have investigated the effect of SecB on the refolding kinetics of the small protein barstar in I M guanidine hydrochloride at pH 7.0 and 25 degrees C using fluorescence spectroscopy. We show that SecB does not bind either the native or the unfolded states of barstar but binds to a late near-native intermediate along the folding pathway. For barstar, polypeptide collapse and formation of a hydrophobic surface are required for binding to SecB. SecB does not change the apparent rate constant of barstar refolding. The kinetic data for SecB binding to barstar are not consistent with simple kinetic partitioning models.

Chiral molecular self-assembly of phospholipid tubules: A circular dichroism study

We report on spectroscopic studies of the chiral structure in phospholipid tubules formed in mixtures of alcohol and water. Synthetic phospholipids containing diacetylenic moieties in the acyl chains self-assemble into hollow, cylindrical tubules in appropriate conditions. Circular dichroism provides a direct measure of chirality of the molecular structure. We find that the CD spectra of tubules formed in mixtures of alcohol and water depends strongly on the alcohol used and the lipid concentration. The relative spectral intensity of different circular dichroism bands correlates with the number of bilayers observed using microscopy. The results provide experimental evidence that tubule formation is based on chiral packing of the lipid molecules and that interbilayer interactions are important to the tubule structure

Beta-Turn Conformations In Crystal-Structures Of Model Peptides Containing Alpha,Alpha-Di-N-Propylglycine And Alpha,Alpha-Di-N-Butylglycine

The crystal state conformations of three peptides containing the alpha, alpha-dialkylated residues, alpha,alpha-di-n-propylglycine (Dpg) and alpha,alpha-di-n-butylglycine (Dbg), have been established by x-ray diffraction. Boc-Ala-Dpg-Ala-OMe (I) and Boc-Ala-Dbg-Ala-OMe (III) adopt distorted type II beta-turn conformations with Ala (1) and Dpg/Dbg (2) as the corner residues. In both peptides the conformational angles at the Dxg residue (I: phi = 66.2 degrees, psi = 19.3 degrees; III: phi = 66.5 degrees, psi = 21.1 degrees) deviate appreciably from ideal values for the i + 2 residue in a type II beta-turn. In both peptides the observed (N...O) distances between the Boc CO and Ala(3) NH groups are far too long (I: 3.44 Angstrom; III: 3.63 Angstrom) for an intramolecular 4 --> 1 hydrogen bond. Boc-Ala-Dpg-Ala-NHMe (II) crystallizes with two independent molecules in the asymmetric unit. Both molecules IIA and IIB adopt consecutive beta-turn (type III-III in IIA and type III-I in IIB) or incipient 3(10)-helical structures, stabilized by two intramolecular 4 --> 1 hydrogen bonds. In all four molecules the bond angle N-C-alpha-C' (tau) at the Dxg residues are greater than or equal to 110 degrees. The observation of conformational angles in the helical region of phi,psi space at these residues is consistent with theoretical predictions.

X-ray Structure of Gelonin at 1.8 Å Resolution

Gelonin is a single chain ribosome inactivating protein (RIP) with potential application in the treatment of cancer and AIDS. Diffraction quality crystals grown using PEG3350, belong to the space group P2(1), with it a = 49.4 Angstrom b = 44.9 Angstrom, c = 137.4 Angstrom and beta = 98.4 degrees, and contain two molecules in the asymmetric unit. Diffraction data collected to 1.8 Angstrom resolution has a R(m) value of 7.3%. Structure of gelonin has been solved by the molecular replacement method, using ricin A chain as the search model. Crystallographic refinement using X-PLOR resulted in a model for which the r.m.s deviations from ideal bond lengths and bond angles are 0.012 Angstrom and 2.7 degrees, respectively The final R-factor is 18.4% for 39,806 reflections for which I > 1.0 sigma(I).The C-alpha atoms of the two molecules in the asymmetric unit superpose to within 0.38 Angstrom for 247 atom pairs. The overall fold of gelonin is similar to that of other RIPs such as ricin A chain and alpha-momorcharin, the r.m.s.d. for C-alpha superpositions being 1.3 and 1.4 Angstrom, respectively The-catalytic residues (Glu166, Arg169 and Tyr113) in the active site form a hydrogen bond scheme similar to that observed in other RIPs. The conformation of Tyr74 in the active site, however, is significantly different from that in alpha-momorcharin. Three well defined water molecules are located in the active site cavity and one of them, X319, superposes to within 0.2 Angstrom of a corresponding water molecule in the structure of alpha-momorcharin. Any of the three could be the substrate water molecule in the hydrolysis reaction catalysed by gelonin.Difference electron density for a N-linked sugar moiety has been observed near only one of the two potential glycosylation sites in the sequence. The amino acid at position 239 has been established as Lys by calculation of omit electron density maps.The two cysteine residues in the sequence, Cys44 and Cys50, form a disulphide bond, and are therefore not available for disulphide conjugation with antibodies. Based on the structure, the region of the molecule that is involved in intradimer interactions is suggested to be suitable for introducing a Cys residue for purposes of conjugation with an antibody to produce useful immunotoxins.