Modeling the molecular weight distribution of block copolymer formation in a reversible addition-fragmentation chain transfer mediated living radical polymerization

Block copolymers have become an integral part of the preparation of complex architectures through self-assembly. The use of reversible addition-fragmentation chain transfer (RAFT) allows blocks ranging from functional to nonfunctional polymers to be made with predictable molecular weight distributions. This article models block formation by varying many of the kinetic parameters. The simulations provide insight into the overall polydispersities (PDIs) that will be obtained when the chain-transfer constants in the main equilibrium steps are varied from 100 to 0.5. When the first dormant block [polymer-S-C(Z)=S] has a PDI of 1 and the second propagating radical has a low reactivity to the RAFT moiety, the overall PDI will be greater than 1 and dependent on the weight fraction of each block. When the first block has a PDI of 2 and the second propagating radical has a low reactivity to the RAFT moiety, the PDI will decrease to around 1.5 because of random coupling of two broad distributions. It is also shown how we can in principle use only one RAFT agent to obtain block copolymers with any desired molecular weight distribution. We can accomplish this by maintaining the monomer concentration at a constant level in the reactor over the course of the reaction. (c) 2005 Wiley Periodicals, Inc.

Determination of thermal and optical parameters of melanins by photopyroelectric spectroscopy

Photopyroelectric spectroscopy (PPE) was used to study the thermal and optical properties of melanins. The photopyroelectric intensity signal and its phase were independently measured as a function of wavelength and chopping frequency for a given wavelength in the saturation part of the PPE spectrum. Equations for both the intensity and the phase of the PPE signal were used to fit the experimental results. From these fits we obtained for the first time, with great accuracy, the thermal diffusivity coefficient, the thermal conductivity, and the specific heat of the samples, as well as a value for the condensed phase optical gap, which we found to be 1.70 eV. (c) 2005 American Institute of Physics.

Enhancing the ratio of molecular ions to non-covalent compounds in the electrospray interface of LC-MS in quantitative analysis

A common problem encountered during the development of MS methods for the quantitation of small organic molecules by LGMS is the formation of non-covalently bound species or adducts in the electrospray interface. Often the population of the molecular ion is insignificant compared to those of all other forms of the analyte produced in the electrospray, making it difficult to obtain the sensitivity required for accurate quantitation. We have investigated the effects of the following variables: orifice potential, nebulizer gas flow, temperature, solvent composition and the sample pH on the relative distributions of ions of the types MH+, MNa+, MNH+, and 2MNa(+), where M represents a 4 small organic molecule: BAY 11-7082 ((E)-3-[4-methylphenylsulfonyl]-2-propenenitrile). Orifice potential, solvent composition and the sample pH had the greatest influence on the relative distributions of these ions, making these parameters the most useful for optimizing methods for the quantitation of small molecules.

Design strategies for controlling the molecular weight and rate using reversible addition-fragmentation chain transfer mediated living radical polymerization

Living radical polymerization has allowed complex polymer architectures to be synthesized in bulk, solution, and water. The most versatile of these techniques is reversible addition-fragmentation chain transfer (RAFT), which allows a wide range of functional and nonfunctional polymers to be made with predictable molecular weight distributions (MWDs), ranging from very narrow to quite broad. The great complexity of the RAFT mechanism and how the kinetic parameters affect the rate of polymerization and MWD are not obvious. Therefore, the aim of this article is to provide useful insights into the important kinetic parameters that control the rate of polymerization and the evolution of the MWD with conversion. We discuss how a change in the chain-transfer constant can affect the evolution of the MWD. It is shown how we can, in principle, use only one RAFT agent to obtain a poly-mer with any MWD. Retardation and inhibition are discussed in terms of (1) the leaving R group reactivity and (2) the intermediate radical termination model versus the slow fragmentation model. (c) 2005 Wiley Periodicals, Inc.

Adaptation is not required to explain the long-term response of axons to molecular gradients

It has been suggested that growth cones navigating through the developing nervous system might display adaptation, so that their response to gradient signals is conserved over wide variations in ligand concentration. Recently however, a new chemotaxis assay that allows the effect of gradient parameters on axonal trajectories to be finely varied has revealed a decline in gradient sensitivity on either side of an optimal concentration. We show that this behavior can be quantitatively reproduced with a computational model of axonal chemotaxis that does not employ explicit adaptation. Two crucial components of this model required to reproduce the observed sensitivity are spatial and temporal averaging. These can be interpreted as corresponding, respectively, to the spatial spread of signaling effects downstream from receptor binding, and to the finite time over which these signaling effects decay. For spatial averaging, the model predicts that an effective range of roughly one-third of the extent of the growth cone is optimal for detecting small gradient signals. For temporal decay, a timescale of about 3 minutes is required for the model to reproduce the experimentally observed sensitivity.

Heterozygote effects in mice with partial truncations in the growth hormone receptor cytoplasmic domain: Assessment of growth parameters and phenotype

The GH receptor (GHR) is essential for normal postnatal growth and development, and the molecular basis of GHR action has been studied intensively. Clinical case studies and more recently mouse models have revealed the extensive phenotype of impaired GH action. We recently reported two new mouse models, possessing cytoplasmic truncations at position 569 (plus Y539/545-F) and 391, which were created to identify functional subdomains within the cytoplasmic signaling domain. In the homozygous state, these animals show progressively impaired postnatal growth coupled with complex changes in gene expression. We describe here an extended phenotype analysis encompassing the heterozygote state to identify whether single copies of these mutant receptors bring about partial or dominant-negative phenotypes. It appears that the retention of the ubiquitin-dependent endocytosis motif the N-terminal cytoplasmic domain permits turnover of these mutant receptors because no dominant-negative phenotype is seen. Nonetheless, we do observe partial impairment of postnatal growth in heterozygotes supporting limited haploinsufficiency. Reproductive function is impaired in these models in a progressive manner, in parallel with loss of signal transducer and activator of transcription-5 activation ability. In summary, we describe a more comprehensive phenotypic analysis of these mouse models, encompassing overall and longitudinal body growth, reproductive function, and hormonal status in both the heterozygote and homozygote state. Our results suggest that patients expressing single copies of similarly mutated GHRs would not display an obvious clinical phenotype.

Negative second virial coefficients as predictors of protein crystal growth: Evidence from sedimentation equilibrium studies that refutes the designation of those light scattering parameters as osmotic virial coefficients

The effects of ammonium sulphate concentration on the osmotic second virial coefficient (B-AA/M-A) for equine serum albumin (pH 5.6, 20 degrees C) have been examined by sedimentation equilibrium. After an initial steep decrease with increasing ammonium sulphate concentration, B-AA/M-A assumes an essentially concentration-independent magnitude of 8-9 ml/g. Such behaviour conforms with the statistical-mechanical prediction that a sufficient increase in ionic strength should effectively eliminate the contributions of charge interactions to B-AA/M-A but have no effect on the covolume contribution (8.4 ml/g for serum albumin). A similar situation is shown to apply to published sedimentation equilibrium data for lysozyme (pH 4.5). Although termed osmotic second virial coefficients and designated as such (B-22), the negative values obtained in published light scattering studies of both systems have been described incorrectly because of the concomitant inclusion of the protein-salt contribution to thermodynamic nonideality of the protein. Those negative values are still valid predictors of conditions conducive to crystal growth inasmuch as they do reflect situations in which there is net attraction between protein molecules. However, the source of attraction responsible for the negative virial coefficient stems from the protein-salt rather than the protein-protein contribution, which is necessarily positive. (c) 2005 Elsevier B.V. All rights reserved.

A Bethe ansatz solvable model for superpositions of Cooper pairs and condensed molecular bosons

We introduce a general Hamiltonian describing coherent superpositions of Cooper pairs and condensed molecular bosons. For particular choices of the coupling parameters, the model is integrable. One integrable manifold, as well as the Bethe ansatz solution, was found by Dukelsky et al. [J. Dukelsky, G.G. Dussel, C. Esebbag, S. Pittel, Phys. Rev. Lett. 93 (2004) 050403]. Here we show that there is a second integrable manifold, established using the boundary quantum inverse scattering method. In this manner we obtain the exact solution by means of the algebraic Bethe ansatz. In the case where the Cooper pair energies are degenerate we examine the relationship between the spectrum of these integrable Hamiltonians and the quasi-exactly solvable spectrum of particular Schrodinger operators. For the solution we derive here the potential of the Schrodinger operator is given in terms of hyperbolic functions. For the solution derived by Dukelsky et al., loc. cit. the potential is sextic and the wavefunctions obey PT-symmetric boundary conditions. This latter case provides a novel example of an integrable Hermitian Hamiltonian acting on a Fock space whose states map into a Hilbert space of PE-symmetric wavefunctions defined on a contour in the complex plane. (c) 2006 Elsevier B.V. All rights reserved.

Flexibility of short-strand RNA in aqueous solution as revealed by molecular dynamics simulation:are A′-RNA and A-RNA distinct conformational structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A'-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA namely major groove width, inclination and the number of base pairs in a helical twist over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A'-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors). © CSIRO 2009.

Controlling protein molecular dynamics:how to accelerate folding while preserving the native state

The dynamics of peptides and proteins generated by classical molecular dynamics (MD) is described by using a Markov model. The model is built by clustering the trajectory into conformational states and estimating transition probabilities between the states. Assuming that it is possible to influence the dynamics of the system by varying simulation parameters, we show how to use the Markov model to determine the parameter values that preserve the folded state of the protein and at the same time, reduce the folding time in the simulation. We investigate this by applying the method to two systems. The first system is an imaginary peptide described by given transition probabilities with a total folding time of 1 micros. We find that only small changes in the transition probabilities are needed to accelerate (or decelerate) the folding. This implies that folding times for slowly folding peptides and proteins calculated using MD cannot be meaningfully compared to experimental results. The second system is a four residue peptide valine-proline-alanine-leucine in water. We control the dynamics of the transitions by varying the temperature and the atom masses. The simulation results show that it is possible to find the combinations of parameter values that accelerate the dynamics and at the same time preserve the native state of the peptide. A method for accelerating larger systems without performing simulations for the whole folding process is outlined.

Solvent induced NMR chemical shifts that arise from molecular encounters

Recently Homer and Percival have postulated that intermolecular van der Waals dispersion forces can be characterized by three mechanisms. The first arises via the mean square reaction field < R1; 2> due to the transient dipole of a particular solute molecule that is considered situated in a cavity surrounded by solvent molecules; this was characterized by an extended Onsager approach. The second stems from the extra cavity mean square reaction field < R2; 2> of the near neighbour solvent molecules. The third originates from square field electric fields E2BI due to a newly characterized effect in which solute atoms are `buffeted' by the peripheral atoms of adjacent solvent molecules. The present work concerns more detailed studies of the buffeting screening, which is governed by sterically controlled parameter (2T - T)2, where and are geometric structural parameters. The original approach is used to characterise the buffeting shifts induced by large solvent molecules and the approach is found to be inadequate. Consequently, improved methods of calculating and are reported. Using the improved approach it is shown that buffeting is dependent on the nature of the solvent as well as the nature of the solute molecule. Detailed investigation of the buffeting component of the van der Waals chemical shifts of selected solutes in a range of solvents containing either H or Cl as peripheral atoms have enabled the determination of a theoretical acceptable value for the classical screening coefficient B for protons. 1H and 13C resonance studies of tetraethylmethane and 1H, 13C and 29Si resonance studies of TMS have been used to support the original contention that three (< R1; 2> , < R2; 2> and E2BI) components of intermolecular van der Waals dispersion fields are required to characterise vdW chemical shifts.

Studies of molecular interactions using physical techniques

The investigations described in this thesis concern the molecular interactions between polar solute molecules and various aromatic compounds in solution. Three different physical methods were employed. Nuclear magnetic resonance (n.m.r.) spectroscopy was used to determine the nature and strength of the interactions and the geometry of the transient complexes formed. Cryoscopic studies were used to provide information on the stoichiometry of the complexes. Dielectric constant studies were conducted in an attempt to confirm and supplement the spectroscopic investigations. The systems studied were those between nitromethane, chloroform, acetonitrile (solutes) and various methyl substituted benzenes. In the n.m.r. work the dependence of the solute chemical shift upon the compositions of the solutions was determined. From this the equilibrium quotients (K) for the formation of each complex and the shift induced in the solute proton by the aromatic in the complex were evaluated. The thermodynamic parameters for the interactions were obtained from the determination of K at several temperatures. The stoichiometries of the complexes obtained from cryoscopic studies were found to agree with those deduced from spectroscopic investigations. For most systems it is suggested that only one type of complex, of 1:1 stiochiometry, predominates except that for the acetonitrile-benzene system a 1:2 complex is formed. Two sets of dielectric studies were conducted, the first to show that the nature of the interaction is dipole-induced dipole and the second to calculate K. The equilibrium quotients obtained from spectroscopic and dielectric studies are compared. Time-averaged geometries of the complexes are proposed. The orientation of solute, with respect to the aromatic for the 1:1 complexes, appears to be the one in which the solute lies symmetrically about the aromatic six-fold axis whereas for the 1:2 complex, a sandwich structure is proposed. It is suggested that the complexes are formed through a dipole-induced dipole interaction and steric factors play some part in the complex formation.

T-cell epitope prediction and immune complex simulation using molecular dynamics:state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

Adsorption of dimethyl ether (DME) on zeolite molecular sieves

In recent years there has been growing interest in the use of dimethyl ether (DME) as an alternative fuel. In this study, the adsorption of DME on molecular sieves 4Å (Mol4A) and 5Å (Mol5A) has been experimentally investigated using the volumetric adsorption method. Data on the adsorption isotherms, heats of adsorption, and adsorption kinetic have been obtained and used to draw conclusions and compare the performance of the two adsorbents. Within the conditions considered, the adsorption capacity of Mol5A was found to be around eight times higher than the capacity of Mol4A. Low temperature adsorption and thermal pre-treatment of the adsorbents in vacuum were observed to be favourable for increased adsorption capacity. The adsorption isotherms for both adsorbent were fitted to the Freundlich model and the corresponding model parameters are proposed. The adsorption kinetic analysis suggest that the DME adsorption on Mol5A is controlled by intracrystalline diffusion resistance, while on Mol4A it is mainly controlled by surface layering resistance with the diffusion only taking place at the start of adsorption and for a very limited short time. The heats of adsorption were calculated by a calorimetric method based on direct temperature measurements inside the adsorption cell. Isosteric heats, calculated by the thermodynamic approach (Clasius-Clapeyron equation), have consistently shown lower values. The maximum heat of adsorption was found to be 25.9kJmol-1 and 20.1kJmol-1 on Mol4A and Mol5A, respectively; thus indicating a physisorption type of interactions. © 2014 Elsevier B.V.

Bridging large and small scales of water models using hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework

This thesis presents a two-dimensional water model investigation and development of a multiscale method for the modelling of large systems, such as virus in water or peptide immersed in the solvent. We have implemented a two-dimensional ‘Mercedes Benz’ (MB) or BN2D water model using Molecular Dynamics. We have studied its dynamical and structural properties dependence on the model’s parameters. For the first time we derived formulas to calculate thermodynamic properties of the MB model in the microcanonical (NVE) ensemble. We also derived equations of motion in the isothermal–isobaric (NPT) ensemble. We have analysed the rotational degree of freedom of the model in both ensembles. We have developed and implemented a self-consistent multiscale method, which is able to communicate micro- and macro- scales. This multiscale method assumes, that matter consists of the two phases. One phase is related to micro- and the other to macroscale. We simulate the macro scale using Landau Lifshitz-Fluctuating Hydrodynamics, while we describe the microscale using Molecular Dynamics. We have demonstrated that the communication between the disparate scales is possible without introduction of fictitious interface or approximations which reduce the accuracy of the information exchange between the scales. We have investigated control parameters, which were introduced to control the contribution of each phases to the matter behaviour. We have shown, that microscales inherit dynamical properties of the macroscales and vice versa, depending on the concentration of each phase. We have shown, that Radial Distribution Function is not altered and velocity autocorrelation functions are gradually transformed, from Molecular Dynamics to Fluctuating Hydrodynamics description, when phase balance is changed. In this work we test our multiscale method for the liquid argon, BN2D and SPC/E water models. For the SPC/E water model we investigate microscale fluctuations which are computed using advanced mapping technique of the small scales to the large scales, which was developed by Voulgarakisand et. al.