Development of a Cell Size Relaxed Scheme for the Direct Simulation Monte Carlo Method

The conventional direct simulation Monte Carlo (DSMC) method has a strong restriction on the cell size because simulated particles are selected randomly within the cell for collisions. Cells with size larger than the molecular mean free path are generally not allowed in correct DSMC simulations. However, the cell-size induced numerical error can be controlled if the gradients of flow properties are properly involved during collisions. In this study, a large cell DSMC scheme is proposed to relax the cell size restriction. The scheme is applied to simulate several test problems and promising results are obtained even when the cell size is greater than 10 mean free paths of gas molecules. However, it is still necessary, of course, that the cell size be small with respect to the flow field structures that must be resolved.

Identification of larval sea basses (Centropristis spp.) using ribosomal DNA-specific molecular assays

The identification of sea bass (Centropristis) larvae to species is difficult because of similar morphological characters, spawning times, and overlapping species ranges. Black sea bass (Centropristis striata) is an important fishery species and is currently considered to be overfished south of Cape Hatteras, North Carolina. We describe methods for identifying three species of sea bass larvae using polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) assays based on species-specific amplification of rDNA internal transcribed spacer regions. The assays were tested against DNA of ten other co-occurring reef fish species to ensure the assay's specificity. Centropristis larvae were collected on three cruises during cross-shelf transects and were used to validate the assays. Seventy-six Centropristis larva were assayed and 69 (91%) were identified successfully. DNA was not amplified from 5% of the larvae and identification was inconclusive for 3% of the larvae. Those assays can be used to identify sea bass eggs and larvae and will help to assess spawning locations, spawning times, and larval dispersal.

The impacts of coastal protection structures in California’s Monterey Bay National Marine Sanctuary

This report outlines the potential impacts of coastal protection structures on the resources of the Monterey Bay National Marine Sanctuary. At least 15 miles of the Sanctuary’s 300-mile shoreline are currently armored with seawalls and riprap revetments. Most of these coastal protection structures are placed above the mean high tide line, the official boundary of the Sanctuary, yet some influences of armoring impinge on the marine realm and on recreational use. In addition, continued sea level rise and accompanying coastal retreat will force many of these structures below the high tide line over time. The Monterey Bay National Marine Sanctuary staff has recognized the significance of coastal armoring, identifying it as a critical issue in the Coastal Armoring Action Plan of the draft Joint Management Plan. This summary is intended to provide general background information for Sanctuary policies on coastal armoring. The impacts discussed include: aesthetic depreciation, beach loss due to placement, access restriction, loss of sand supply from eroding cliffs, passive erosion, and active erosion. In addition, the potential biological impacts are explored. Finally, an appraisal of how differing armor types compare in relation to impacts, expense and engineering is presented. While the literature cited in this report focus predominantly on the California coast, the framework for this discussion could have implications for other actively eroding coastlines. (PDF contains 26 pages.)

A Note on Risk Acceptance, Bankruptcy Avoidance and Riskiness Measures

In this work we clarify the relationships between riskiness, risk acceptance and bankruptcy avoidance. We distinguish between the restriction on the current wealth required to make a gamble acceptable to the decision maker and the restriction on the current wealth required to guarantee no bankruptcy if a gamble is accepted. We focus on the measure of riskiness proposed by Foster and Hart.

Online exercise for the design and simulation of PCR and PCR-RFLP experiments

[EN] Background: Polymerase Chain Reaction (PCR) and Restriction Fragment Length Polymorphism of PCR products (PCR-RFLP) are extensively used molecular biology techniques. An exercise for the design and simulation of PCR and PCR-RFLP experiments will be a useful educational tool. Findings: An online PCR and PCR-RFLP exercise has been create that requires users to find the target genes,compare them, design primers, search for restriction endonucleases, and finally to simulate the experiment. Each user of the service is randomly assigned a gene from Escherichia coli; to complete the exercise, users must design an experiment capable of distinguishing among E. coli strains. By applying the experimental procedure to all completely sequenced E. coli, a basic understanding of strain comparison and clustering can also be acquired. Comparison of results obtained in different experiments is also very instructive. Conclusions: The exercise is freely available at http://insilico.ehu.es/edu.

Emergence of Cooperation in Heterogeneous Population: A Discrete-Time Replicator Dynamics Analysis

The emergence of cooperation is analyzed in heterogeneous populations where individuals can be classified in two groups according to their phenotypic appearance. Phenotype recognition is assumed for all individuals: individuals are able to identify the type of every other individual, but fail to recognize their own type, and thus behave under partial information conditions. The interactions between individuals are described by 2 × 2 symmetric games where individuals can either cooperate or defect. The evolution of such populations is studied in the framework of evolutionary game by means of the replicator dynamics. Overlapping generations are considered, so the replicator equations are formulated in discrete-time form. The well-posedness conditions of the system are derived. Depending on the parameters of the game, a restriction may exist for the generation length. The stability analysis of the dynamical system is carried out and a detailed description of the behavior of trajectories starting from the interior of the state-space is given. We find that, provided the conditions of well-posedness are verified, the linear stability of monomorphic states in the discrete-time replicator coincides with the one of the continuous case. Specific from the discrete-time case, a relaxed restriction for the generation length is derived, for which larger time-steps can be used without compromising the well-posedness of the replicator system.

Is General Relativity a simplified theory?

IARD 8th Biennial Conference on Classical and Quantum Relativistic Dynamics of Particles and Fields - Galileo Galilei Inst Theoret Phys (GGI), Florence, ITALY - MAY 29-JUN 01, 2012. Edited by:Horowitz, LP

From transistor to trapped-ion computers for quantum chemistry

Over the last few decades, quantum chemistry has progressed through the development of computational methods based on modern digital computers. However, these methods can hardly fulfill the exponentially-growing resource requirements when applied to large quantum systems. As pointed out by Feynman, this restriction is intrinsic to all computational models based on classical physics. Recently, the rapid advancement of trapped-ion technologies has opened new possibilities for quantum control and quantum simulations. Here, we present an efficient toolkit that exploits both the internal and motional degrees of freedom of trapped ions for solving problems in quantum chemistry, including molecular electronic structure, molecular dynamics, and vibronic coupling. We focus on applications that go beyond the capacity of classical computers, but may be realizable on state-of-the-art trapped-ion systems. These results allow us to envision a new paradigm of quantum chemistry that shifts from the current transistor to a near-future trapped-ion-based technology.

The COMT Val158 Met polymorphism as an associated risk factor for Alzheimer disease and mild cognitive impairment in APOE 4 carriers

Background: The aim of this study is to examine the influence of the catechol-O-methyltranferase (COMT) gene (polymorphism Val158 Met) as a risk factor for Alzheimer's disease (AD) and mild cognitive impairment of amnesic type (MCI), and its synergistic effect with the apolipoprotein E gene (APOE). A total of 223 MCI patients, 345 AD and 253 healthy controls were analyzed. Clinical criteria and neuropsychological tests were used to establish diagnostic groups. The DNA Bank of the University of the Basque Country (UPV-EHU) (Spain) determined COMT Val158 Met and APOE genotypes using real time polymerase chain reaction (rtPCR) and polymerase chain reaction (PCR), and restriction fragment length polymorphism (RFLPs), respectively. Multinomial logistic regression models were used to determine the risk of AD and MCI. Results: Neither COMT alleles nor genotypes were independent risk factors for AD or MCI. The high activity genotypes (GG and AG) showed a synergistic effect with APOE epsilon 4 allele, increasing the risk of AD (OR = 5.96, 95% CI 2.74-12.94, p < 0.001 and OR = 6.71, 95% CI 3.36-13.41, p < 0.001 respectivily). In AD patients this effect was greater in women. In MCI patients such as synergistic effect was only found between AG and APOE epsilon 4 allele (OR = 3.21 95% CI 1.56-6.63, p = 0.02) and was greater in men (OR = 5.88 95% CI 1.69-20.42, p < 0.01). Conclusion: COMT (Val158 Met) polymorphism is not an independent risk factor for AD or MCI, but shows a synergistic effect with APOE epsilon 4 allele that proves greater in women with AD.

Structure-function studies of serine hydrolases: synthesis of a gene for ∝-lytic protease and the purification and characterization of a mutant β-lactamase

The author has constructed a synthetic gene for ∝-lytic protease. Since the DNA sequence of the protein is not known, the gene was designed by using the reverse translation of ∝-lytic protease's amino acid sequence. Unique restriction sites are carefully sought in the degenerate DNA sequence to aid in future mutagenesis studies. The unique restriction sites are designed approximately 50 base pairs apart and their appropriate codons used in the DNA sequence. The codons used to construct the DNA sequence of ∝-lytic protease are preferred codons in E-coli or used in the production of β-lactamase. Codon usage is also distributed evenly to ensure that one particular codon is not heavily used. The gene is essentially constructed from the outside in. The gene is built in a stepwise fashion using plasmids as the vehicles for the ∝-lytic oligomers. The use of plasmids allows the replication and isolation of large quantities of the intermediates during gene synthesis. The ∝-lytic DNA is a double-stranded oligomer that has sufficient overhang and sticky ends to anneal correctly in the vector. After six steps of incorporating ∝-lytic DNA, the gene is completed and sequenced to ensure that the correct DNA sequence is present and that no mutations occurred in the structural gene.

β-lactamase is the other serine hydrolase studied in this thesis. The author used the class A RTEM-1 β- lactamase encoded on the plasmid pBR322 to investigate the roll of the conserved threonine residue at position 71. Cassette mutagenesis was previously used to generate all possible amino acid substitutions at position 71. The work presented here describes the purification and kinetic characterization of a T71H mutant previously constructed by S. Schultz. The mutated gene was transferred into plasmid pJN for expression and induced with IPTG. The enzyme is purified by column chromatography and FPLC to homogeneity. Kinetic studies reveal that the mutant has lower k_(cat) values on benzylpenicillin, cephalothin and 6-aminopenicillanic acid but no changes in k_m except for cephalothin which is approximately 4 times higher. The mutant did not change siginificantly in its pH profile compared to the wild-type enzyme. Also, the mutant is more sensitive to thermal denaturation as compared to the wild-type enzyme. However, experimental evidence indicates that the probable generation of a positive charge at position 71 thermally stabilized the mutant.

Oscillatory integral operators related to pointwise convergence of Schrödinger operators

In this thesis we consider smooth analogues of operators studied in connection with the pointwise convergence of the solution, u(x,t), (x,t) ∈ ℝ^n x ℝ, of the free Schrodinger equation to the given initial data. Such operators are interesting examples of oscillatory integral operators with degenerate phase functions, and we develop strategies to capture the oscillations and obtain sharp L^2 → L^2 bounds. We then consider, for fixed smooth t(x), the restriction of u to the surface (x,t(x)). We find that u(x,t(x)) ∈ L^2(D^n) when the initial data is in a suitable L^2-Sobolev space H^8 (ℝ^n), where s depends on conditions on t.

Distinct intron DNA structures in simian virus 40 T-antigen and adenovirus 2 E1A genes

Distinct structures delineating the introns of Simian Virus 40 T-antigen and Adenovirus 2 E1A genes have been discovered. The structures, which are centered around the branch points of the genes inserted in supercoiled double-stranded plasmids, are specifically targeted through photoactivated strand cleavage by the metal complex tris(4,7-diphenyl-1,10-phenanthroline)rhodium(III). The DNA sites that are recognized lack sequence homology but are similar in demarcating functionally important sites on the RNA level. The single-stranded DNA fragments corresponding to the coding strands of the genes were also found to fold into a structure apparently identical to that in the supercoiled genes based on the recognition by the metal complex. Further investigation of different single-stranded DNA fragments with other structural probes, such as another metal complex bis(1,10-phenanthroline)(phenanthrenequinone diimine)rhodium(III), AMT (4'aminomethyl-4,5',8 trimethylpsoralen), restriction enzyme Mse I, and mung bean nuclease, showed that the structures require the sequ ences at both ends of the intron plus the flanking sequences but not the middle of the intron. The two ends form independent helices which interact with each other to form the global tertiary structures. Both of the intron structures share similarities to the structure of the Holliday junction, which is also known to be specifically targeted by the former metal complex. These structures may have arisen from early RNA intron structures and may have been used to facilitate the evolution of genes through exon shuffling by acting as target sites for recombinase enzymes.

Characterizing distribution rules for cost sharing games

In noncooperative cost sharing games, individually strategic agents choose resources based on how the welfare (cost or revenue) generated at each resource (which depends on the set of agents that choose the resource) is distributed. The focus is on finding distribution rules that lead to stable allocations, which is formalized by the concept of Nash equilibrium, e.g., Shapley value (budget-balanced) and marginal contribution (not budget-balanced) rules.

Recent work that seeks to characterize the space of all such rules shows that the only budget-balanced distribution rules that guarantee equilibrium existence in all welfare sharing games are generalized weighted Shapley values (GWSVs), by exhibiting a specific 'worst-case' welfare function which requires that GWSV rules be used. Our work provides an exact characterization of the space of distribution rules (not necessarily budget-balanced) for any specific local welfare functions remains, for a general class of scalable and separable games with well-known applications, e.g., facility location, routing, network formation, and coverage games.

We show that all games conditioned on any fixed local welfare functions possess an equilibrium if and only if the distribution rules are equivalent to GWSV rules on some 'ground' welfare functions. Therefore, it is neither the existence of some worst-case welfare function, nor the restriction of budget-balance, which limits the design to GWSVs. Also, in order to guarantee equilibrium existence, it is necessary to work within the class of potential games, since GWSVs result in (weighted) potential games.

We also provide an alternative characterization—all games conditioned on any fixed local welfare functions possess an equilibrium if and only if the distribution rules are equivalent to generalized weighted marginal contribution (GWMC) rules on some 'ground' welfare functions. This result is due to a deeper fundamental connection between Shapley values and marginal contributions that our proofs expose—they are equivalent given a transformation connecting their ground welfare functions. (This connection leads to novel closed-form expressions for the GWSV potential function.) Since GWMCs are more tractable than GWSVs, a designer can tradeoff budget-balance with computational tractability in deciding which rule to implement.

Isolation and characterization of proteins that bind to yeast origins of DNA replication

Yeast chromosomes contain sequences called ARSs which function as origins of replication in vitro and in vivo. We have carried out a systematic deletion analysis of ARS1, allowing us to define three functionally distinct domains, designated A, B, and C. Domain A is a sequence of 11 to 19bp, containing the core consensus element that is required for replication. The core consensus sequence, A/TTTTATPuTTTA/T, is conserved at all ARSs sequenced to date. A fragment containing only element A and 8 flanking nucleotides enables autonomous replication of centromeric plasmids. These plasmids replicate very inefficiently, suggesting that flanking sequences must be important for ARS function. Domain B also provides important sequences needed for efficient replication. Deletion of domain B drastically increases the doubling times of transformants and reduces plasmid stability. Domain B contains a potential consensus sequence conserved at some ARSs which overlaps a region of bent DNA. Mutational analysis suggests this bent DNA may be important for ARS function. Deletion of domain C has only a slight effect on replication of plasmids carrying those deletions.

We have identified a protein called ARS binding factor I (ABF-I) that binds to the HMR-E ARS and ARS1. We have purified this protein to homogeneity using conventional and oligonucleotide affinity chromatography. The protein has an apparent molecular weight of 135kDa and is present at about 700 molecules per diploid cell, based on the yield of purified protein and in situ antibody staining. DNaseI footprinting reveals that ABF-I binds sequence-specifically to an approximately 24bp sequence that overlaps element Bat ARSl. This same protein binds to and protects a similar size region at the HMR-E ARS.

We also find evidence for another ARS binding protein, ABF-III, based on DN asei footprint analysis and gel retardation assays. The protein protects approximately 22bp adjacent to the ABF-I site. There appears to be no interaction between ABF-I and ABF-III despite the proximity of their binding sites.

To address the function of ABF-I in DNA replication, we have cloned the ABF-I gene using rabbit polyclonal anti-sera and murine monoclonal antibodies against ABF-I to screen a λgt11 expression library. Four EcoRI restriction fragments were isolated which encoded proteins that were recognized by both polyclonal and monoclonal antibodies. A gene disruption can now be constructed to determine the in vivo function of ABF-I.

Structural and mechanistic motifs in membrane proteins: the three-dimensional modelling of rhodopsin, band 3, and the nicotinic acetylcholine receptor

Because so little is known about the structure of membrane proteins, an attempt has been made in this work to develop techniques by which to model them in three dimensions. The procedures devised rely heavily upon the availability of several sequences of a given protein. The modelling procedure is composed of two parts. The first identifies transmembrane regions within the protein sequence on the basis of hydrophobicity, β-turn potential, and the presence of certain amino acid types, specifically, proline and basic residues. The second part of the procedure arranges these transmembrane helices within the bilayer based upon the evolutionary conservation of their residues. Conserved residues are oriented toward other helices and variable residues are positioned to face the surrounding lipids. Available structural information concerning the protein's helical arrangement, including the lengths of interhelical loops, is also taken into account. Rhodopsin, band 3, and the nicotinic acetylcholine receptor have all been modelled using this methodology, and mechanisms of action could be proposed based upon the resulting structures.

Specific residues in the rhodopsin and iodopsin sequences were identified, which may regulate the proteins' wavelength selectivities. A hinge-like motion of helices M3, M4, and M5 with respect to the rest of the protein was proposed to result in the activation of transducin, the G-protein associated with rhodopsin. A similar mechanism is also proposed for signal transduction by the muscarinic acetylcholine and β-adrenergic receptors.

The nicotinic acetylcholine receptor was modelled with four trans-membrane helices per subunit and with the five homologous M2 helices forming the cation channel. Putative channel-lining residues were identified and a mechanism of channel-opening based upon the concerted, tangential rotation of the M2 helices was proposed.

Band 3, the anion exchange protein found in the erythrocyte membrane, was modelled with 14 transmembrane helices. In general the pathway of anion transport can be viewed as a channel composed of six helices that contains a single hydrophobic restriction. This hydrophobic region will not allow the passage of charged species, unless they are part of an ion-pair. An arginine residue located near this restriction is proposed to be responsible for anion transport. When ion-paired with a transportable anion it rotates across the barrier and releases the anion on the other side of the membrane. A similar process returns it to its original position. This proposed mechanism, based on the three-dimensional model, can account for the passive, electroneutral, anion exchange observed for band 3. Dianions can be transported through a similar mechanism with the additional participation of a histidine residue. Both residues are located on M10.