Structure, molecular dynamics, and stress in a linear polymer

We present a study correlating uniaxial stress in a polymer with its underlying structure when it is strained. The uniaxial stress is significantly influenced by the mean-square bond length and mean bond angle. In contrast, the size and shape of the polymer, typically represented by the end-to-end length, mass ratio, and radius of gyration, contribute negligibly. Among externally set control variables, density and polymer chain length play a critical role in influencing the anisotropic uniaxial stress. Short chain polymers more or less behave like rigid molecules. Temperature and rate of loading, in the range considered, have a very mild effect on the uniaxial stress.

Structural characterization and molecular order of rodlike mesogens with three- and four-ring core by XRD and C-13 NMR spectroscopy

Structural characterizations using XRD and C-13 NMR spectroscopy of two rodlike mesogens consisting of (i) three phenyl ring core with a polar cyano terminal and (ii) four phenyl ring core with flexible dodecyl terminal chain are presented. The three-ring-core mesogen with cyano terminal exhibits enantiotropic smectic A phase while the four-ring mesogen reveals polymesomorphism and shows enantiotropic nematic, smectic C, and tilted hexatic phases. The molecular organization in the three-ring mesogen is found to be partial bilayer smectic Ad type, and the interdigitation of the molecules in the neighboring layers is attributed to the presence of the polar terminal group. For the four-ring mesogen, the XRD results confirm the existence of the smectic C and the tilted hexatic mesophases. A thermal variation of the layer spacing across the smectic C phase followed by a discrete jump at the transition to the tilted hexatic phase is also observed. The tilt angles have been estimated to be about 45 degrees in the smectic C phase and about 40 degrees in tilted hexatic phase. C-13 NMR results indicate that in the mesophase the molecules are aligned parallel to the magnetic field. From the C-13-H-1 dipolar couplings determined from the 2D experiments, the overall order parameter for the three-ring mesogen in its smectic A phase has been estimated to be 0.72 while values ranging from 0.88 to 0.44 have been obtained for the four-ring mesogen as it passes from the tilted hexatic to the nematic phase. The orientations of the different rings of the core unit with respect to each other and also with respect to the long axis of the molecule have also been obtained.

Multiple molecular alterations in phosphodegrons 1-3 within AAV2 capsid demonstrates higher hepatic gene transfer efficiency

The success of AAV2 mediated hepatic gene transfer in human trials for diseases such as hemophilia has been hampered by a combination of low transduction efficiency and a robust immune response directed against these vectors. We have previously shown that AAV2 is targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal degradation machinery and modification of the serine(S)/threonine(T) kinase and lysine(K) targets on AAV capsid is beneficial. Thus targeted single mutations of S/T>A(S489A, S498A, T251A) and K>R (K532R) improved the efficiency of gene transfer in vivo as compared to wild type (WT)-AAV2 vectors (∼6-14 fold). In the present study, we evaluated if combined alteration of the phosphodegrons (PD), which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, improves further the gene transfer efficiency. Thus, we generated four multiple mutant vectors (PD: 1+3, S489A+K532R, PD: 1+3, S489A+K532R together with T251 residue which did not lie in any of the phosphodegrons but had shown increased transduction efficiency compared to the WT-AAV2 vector (∼6 fold) and was also conserved in 9 out of 10 AAV serotypes (AAV 1 to 10), PD: 1+3, S489A+K532R+S498A and a fourth combination of PD: 3, K532R+T251. We then evaluated them in vitro and in vivo and compared their gene transfer efficiency with either the WT-AAV2 or the best single mutant S489A-AAV2 vector. The novel multiple mutations on the AAV2 capsid did not affect the overall vector packaging efficiency. All the multiple AAV2 mutants showed superior transduction efficiency in HeLa cells in vitro when compared to either the WT (62-72% Vs 21%) or the single mutant S489A (62-72% Vs 50%) AAV2 vectors as demonstrated by FACS analysis (Fig. 1A). On hepatic gene transfer with 5x10^10 vgs per animal in C57BL/6 mice, all the multiple mutants showed increased transgene expression compared to either the WT-AAV2 (∼15-23 fold) or the S489A single mutant vector (∼2-3 fold) (Fig.1B and C). These novel multiple mutant AAV2 vectors also showed higher vector copy number in murine hepatocytes 4 weeks post transduction, as compared to the WT-AAV2 (∼5-6 Vs 1.4 vector copies/diploid genome) and further higher when compared to the single mutant S489A(∼5-6 fold Vs 3.8 fold) (Fig.1D). Further ongoing studies will demonstrate the therapeutic benefit of one or more of the multiple mutants vectors in preclinical models of hemophilia.

Crystal structures and binding studies of atovaquone and its derivatives with cytochrome bc(1): a molecular basis for drug design

Crystal structure of trans-atovaquone (antimalarial drug), its polymorph and its stereoisomer (cis) along with five other derivatives with different functional groups have been analyzed. Based on the conformational features of these compounds and the characteristics of the nature of intermolecular interactions, valuable insights into the atomistic details of protein-inhibitor interactions have been derived by docking studies. Atovaquone and its derivatives pack in the crystal lattice using intermolecular O-H center dot center dot center dot O hydrogen bond dimer motifs supported by surrogate weak interactions including C-H center dot center dot center dot O and C-H center dot center dot center dot Cl hydrogen bonds. The docking results of these molecules with cytochrome bc(1) show preferences to form N-H center dot center dot center dot O, O-H center dot center dot center dot O and O-H center dot center dot center dot Cl hydrogen bonds. The involvement of halogen atoms in the binding pocket appears to be significant and is contrary to the theoretically predicted mechanism of protein-ligand docking reported earlier based on mimicking experimental binding results of stigmatellin with cytochrome bc(1). The significance of subtle energy factors controlled by weak intermolecular interactions appears to play a major role in drug binding.

Analysis of designed beta-hairpin peptides: molecular conformation and packing in crystals

The crystal structures of several designed peptide hairpins have been determined in order to establish features of molecular conformations and modes of aggregation in the crystals. Hairpin formation has been induced using a centrally positioned (D)Pro-Xxx segment (Xxx = (L)Pro, Aib, Ac(6)c, Ala; Aib = alpha-aminoisobutyric acid; Ac(6)c = 1-aminocyclohexane-1-carboxylic acid). Structures of the peptides Boc-Leu-Phe-Val-(D)Pro-(L)Pro-Leu-Phe-Val-OMe (1), Boc-Leu-Tyr-Val-(D)Pro-(L)Pro-Leu-Phe-Val-OMe (2, polymorphic forms labeled as 2a and 2b), Boc-Leu-Val-Val-(D)Pro-(L)Pro-Leu-Val-Val-OMe (3), Boc-Leu-Phe-Val-(D)Pro-Aib-Leu-Phe-Val-OMe (4, polymorphic forms labeled as 4a and 4b), Boc-Leu-Phe-Val-(D)Pro-Ac(6)c-Leu-Phe-Val-OMe (5) and Boc-Leu-Phe-Val-(D)Pro-Ala-Leu-Phe-Val-OMe (6) are described. All the octapeptides adopt type II' beta-turn nucleated hairpins, stabilized by three or four cross-strand intramolecular hydrogen bonds. The angle of twist between the two antiparallel strands lies in the range of -9.8 degrees to -26.7 degrees. A detailed analysis of packing motifs in peptide hairpin crystals is presented, revealing three broad modes of association: parallel packing, antiparallel packing and orthogonal packing. An attempt to correlate aggregation modes in solution with observed packing motifs in crystals has been made by indexing of crystal faces in the case of three of the peptide hairpins. The observed modes of hairpin aggregation may be of relevance in modeling multiple modes of association, which may provide insights into the structure of insoluble polypeptide aggregates.

Phylogenetic analysis and molecular dating suggest that hemidactylus anamallensis is not a member of the hemidactylus radiation and has an ancient late cretaceous origin

Background of the Work: The phylogenetic position and evolution of Hemidactylus anamallensis (family Gekkonidae) has been much debated in recent times. In the past it has been variously assigned to genus Hoplodactylus (Diplodactylidae) as well as a monotypic genus `Dravidogecko' (Gekkonidae). Since 1995, this species has been assigned to Hemidactylus, but there is much disagreement between authors regarding its phylogenetic position within this genus. In a recent molecular study H. anamallensis was sister to Hemidactylus but appeared distinct from it in both mitochondrial and nuclear markers. However, this study did not include genera closely allied to Hemidactylus, thus a robust evaluation of this hypothesis was not undertaken. Methods: The objective of this study was to investigate the phylogenetic position of H. anamallensis within the gekkonid radiation. To this end, several nuclear and mitochondrial markers were sequenced from H. anamallensis, selected members of the Hemidactylus radiation and genera closely allied to Hemidactylus. These sequences in conjunction with published sequences were subjected to multiple phylogenetic analyses. Furthermore the nuclear dataset was also subjected to molecular dating analysis to ascertain the divergence between H. anamallensis and related genera. Results and Conclusion: Results showed that H. anamallensis lineage was indeed sister to Hemidactylus group but was separated from the rest of the Hemidactylus by a long branch. The divergence estimates supported a scenario wherein H. anamallensis dispersed across a marine barrier to the drifting peninsular Indian plate in the late Cretaceous whereas Hemidactylus arrived on the peninsular India after the Indian plate collided with the Eurasian plate. Based on these molecular evidence and biogeographical scenario we suggest that the genus Dravidogecko should be resurrected.

Calibration of etched fiber bragg grating sensor arrays for measurement of molecular surface adsorption

Etched Fiber Bragg Grating (EFBG) sensors are attractive from the point of the inherently high multiplexing ability of fiber based sensors. However, the strong dependence of the sensitivity of EFBG sensors on the fiber diameter requires robust methods for calibration when used for distributed sensing in a large array format. Using experimental data and numerical modelling, we show that knowledge of the wavelength shift during the etch process is necessary for high-fidelity calibration of EFBG arrays. However as this approach requires the monitoring of every element of the sensor array during etching, we also proposed and demonstrated a calibration scheme using data from bulk refractometry measurements conducted post-fabrication without needing any information about the etching process. Although this approach is not as precise as the first one, it may be more practical as there is no requirement to monitor each element of the sensor array. We were able to calibrate the response of the sensors to within 3% with the approach using information acquired during etching and to within 5% using the post-fabrication bulk refractometry approach in spite of the sensitivities of the array element differing by more than a factor of 4. These two approaches present a tradeoff between accuracy and practicality.

Molecular dynamics simulations of PPI dendrimer-drug complexes

Dendrimeric nanoparticles are potential drug delivery devices which can enhance the solubility of hydrophobic drugs, thus increasing their bioavailability and sustained release action. A quantitative understanding of the dendrimer-drug interactions can give valuable insight into the solubility and release profile of hydrophobic drug molecules in various solvent conditions. Fully atomistic molecular dynamics (MD) simulations have been performed to study the interactions of G5 PPIEDA (G5 ethylenediamine cored poly(propylene imine)) dendrimer and two well known drugs (Famotidine and Indomethacin) at different pH conditions. The study suggested that at low pH the dendrimer-drug complexes are thermodynamically unstable as compared to neutral and high pH conditions. Calculated Potential of Mean Force (PMF) by umbrella sampling showed that the release of drugs from the dendrimer at low pH is spontaneous, median release at neutral pH and slow release at high pH. In addition, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding free energy calculations were also performed at each umbrella sampling window to identify the various energy contributions. To understand the effect of dendrimer chemistry and topology on the solubility and release profile of drugs, this study is extended to explore the solubility and release profile of phenylbutazone drug complexed with G3 poly(amidoamine) and G4 diaminobutane cored PPI dendrimers. The results indicate that the pH-induced conformational changes in dendrimer, ionization states, dendrimer type and pK(a) of the guest molecules influence the free energy barrier and stability of complexation, and thus regulate drug loading, solubility and release.

H/F isosteric substitution to attest different equi-energetic molecular conformations in crystals

The sequential replacement of aromatic H-atoms by F-atoms in 1,6-bis(phenylcarbonate) hexa-2,4-diyne allows access to its possible iso-energetic ``syn'', ``gauche'' and ``anti'' conformations.

Solvation dynamics of tryptophan in water-dimethyl sulfoxide binary mixture: In search of molecular origin of composition dependent multiple anomalies

Experimental and simulation studies have uncovered at least two anomalous concentration regimes in water-dimethyl sulfoxide (DMSO) binary mixture whose precise origin has remained a subject of debate. In order to facilitate time domain experimental investigation of the dynamics of such binary mixtures, we explore strength or extent of influence of these anomalies in dipolar solvation dynamics by carrying out long molecular dynamics simulations over a wide range of DMSO concentration. The solvation time correlation function so calculated indeed displays strong composition dependent anomalies, reflected in pronounced non-exponential kinetics and non-monotonous composition dependence of the average solvation time constant. In particular, we find remarkable slow-down in the solvation dynamics around 10%-20% and 35%-50% mole percentage. We investigate microscopic origin of these two anomalies. The population distribution analyses of different structural morphology elucidate that these two slowing down are reflections of intriguing structural transformations in water-DMSO mixture. The structural transformations themselves can be explained in terms of a change in the relative coordination number of DMSO and water molecules, from 1DMSO:2H(2)O to 1H(2)O:1DMSO and 1H(2)O:2DMSO complex formation. Thus, while the emergence of first slow down (at 15% DMSO mole percentage) is due to the percolation among DMSO molecules supported by the water molecules (whose percolating network remains largely unaffected), the 2nd anomaly (centered on 40%-50%) is due to the formation of the network structure where the unit of 1DMSO:1H(2)O and 2DMSO:1H(2)O dominates to give rise to rich dynamical features. Through an analysis of partial solvation dynamics an interesting negative cross-correlation between water and DMSO is observed that makes an important contribution to relaxation at intermediate to longer times.

Molecular flexibility tuned emission in ``V'' shaped naphthalimides: Hg(II) detection and aggregation-induced emission enhancement (AIEE)

Four ``V'' shaped 1,8-naphthalimides (1-4) have been synthesized and their fluorescence quantum-yields correlated to their molecular flexibility. The correlation was used for detection of Hg(II) via a chemodosimetric approach. 4 was found to be an AIE active molecule with the formation of fluorescent nanoaggregates.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.

Impact of substrate nitridation on the photoluminescence and photovoltaic characteristics of GaN grown on p-Si (100) by molecular beam epitaxy

This report focuses on the structural and optical properties of the GaN films grown on p-Si (100) substrates along with photovoltaic characteristics of GaN/p-Si heterojunctions fabricated with substrate nitridation and in absence of substrate nitridation. The high resolution X-ray diffraction (HRXRD), atomic force microscopy (AFM), Raman and photoluminescence (PL) spectroscopic studies reveal that the significant enhancement in the structural as well as in the optical properties of GaN epifilms grown with silicon nitride buffer layer when compared with the sample grown without silicon nitride buffer layer. The low temperature PL shows a free excitonic (FX) emission peak at 3.51 eV at the temperature of 5 K with a very narrow line width of 35 meV. Temperature dependent PL spectra follow the Varshni equation well and peak energy blue shifts by similar to 63 meV from 300 to 5 K. Raman data confirms the strain free nature and reasonably good crystallinity of the films. The GaN/p-Si heterojunctions fabricated without substrate nitridation show a superior photovoltaic performance compared to the devices fabricated in presence of substrate nitridation. The discussions have been carried out on the junction properties. Such single junction devices exhibit a promising fill factor and conversion efficiency of 23.36 and 0.12 %, respectively, under concentrated AM1.5 illumination.

Gap renormalization of molecular crystals from density-functional theory

Fundamental gap renormalization due to electronic polarization is a basic phenomenon in molecular crystals. Despite its ubiquity and importance, all conventional approaches within density-functional theory completely fail to capture it, even qualitatively. Here, we present a new screened range-separated hybrid functional, which, through judicious introduction of the scalar dielectric constant, quantitatively captures polarization-induced gap renormalization, as demonstrated on the prototypical organic molecular crystals of benzene, pentacene, and C-60. This functional is predictive, as it contains system-specific adjustable parameters that are determined from first principles, rather than from empirical considerations.

Kinetic proofreading at single molecular level: aminoacylation of tRNA(Ile) and the role of water as an editor

Proofreading/editing in protein synthesis is essential for accurate translation of information from the genetic code. In this article we present a theoretical investigation of efficiency of a kinetic proofreading mechanism that employs hydrolysis of the wrong substrate as the discriminatory step in enzyme catalytic reactions. We consider aminoacylation of tRNA(Ile) which is a crucial step in protein synthesis and for which experimental results are now available. We present an augmented kinetic scheme and then employ methods of stochastic simulation algorithm to obtain time dependent concentrations of different substances involved in the reaction and their rates of formation. We obtain the rates of product formation and ATP hydrolysis for both correct and wrong substrates (isoleucine and valine in our case, respectively), in single molecular enzyme as well as ensemble enzyme kinetics. The present theoretical scheme correctly reproduces (i) the amplitude of the discrimination factor in the overall rates between isoleucine and valine which is obtained as (1.8x10(2)).(4.33x10(2)) = 7.8x10(4), (ii) the rates of ATP hydrolysis for both Ile and Val at different substrate concentrations in the aminoacylation of tRNA(Ile). The present study shows a non-michaelis type dependence of rate of reaction on tRNA(Ile) concentration in case of valine. The overall editing in steady state is found to be independent of amino acid concentration. Interestingly, the computed ATP hydrolysis rate for valine at high substrate concentration is same as the rate of formation of Ile-tRNA(Ile) whereas at intermediate substrate concentration the ATP hydrolysis rate is relatively low. We find that the presence of additional editing domain in class I editing enzyme makes the kinetic proofreading more efficient through enhanced hydrolysis of wrong product at the editing CP1 domain.