Coarse-Grained Molecular Dynamics Simulation of the Aggregation Properties of Multiheaded Cationic Surfactants in Water

The aggregation property of multiheaded surfactants has been investigated by constant pressure molecular dynamics (MD) simulation in aqueous medium. The model multiheaded surfactants contain more than one headgroup (x = 2, 3, and 4) for a single tail group. This increases the hydrophilic charge progressively over the hydrophobic tail which has dramatic consequences in the aggregation behavior. In particular, we have looked at the change in the aggregation property such as critical micellar concentration (cmc), aggregation number, and size of the micelles for the multiheaded surfactants in water. We find with increasing number of headgroups of the Multiheaded surfactants that the cmc values increase and the aggregation numbers as well as the size of the micelles decrease. These trends are in agreement with the experimental findings as reported earlier with x = 1, 2, and 3. We also predict the aggregation properties of multiheaded surfactant With four headgroups (x = 4) for which no experimental studies exist yet.

X-ray crystal structure of a ternary copper(II) vitamin B6 complex, hydroxo(2,2'-bipyridyl)(pyridoxinato) copper(II) monohydrate. A rare example of monodentate coordination of copper(II) by the hydroxyl ion

A ternary metal complex involving Vitamin B6 with the formula [Cu(bipy)(pn) (OH)]H2O (bipy = 2,2'²-bipyridine, PN = anionic pyridoxine) has been synthesized and studied in the solid state by means of spectroscopy and X-ray crystallography. The geometry around copper(II) is distorted square pyramidal, two oxygens from phenolic and 4-(hydroxymethyl) groups of pn, two nitrogens from bipy and an axial OH- ion forming the coordination sphere. In this structure pn exists in a new anionic form with deprotonation of the phenolic group. The structure also provides a rare example of monodentate hydroxyl coordination to copper.

A study of the nature of coordination in metal monothiocarbamates based on infrared spectroscopy and the HSAB concept

The nature of coordination in metal monothiocarbamates is shown to depend on the hardness or softness of the metal ton. Thus, the monothiocarbamate ion acts as a monodentate ligand with metal-sulphur bending when the metal ion is a soft acid while it acts as a bidentate ligand when the metal ion is a hard acid; it can exhibit either behaviour when the metal ion is a borderline acid. In dialkyltin and dialkylmonocholorotin complexes, the monothiocarbamate ion acts as a bidentate ligand with strong Sn-S bonding while in trialkyl-or triaryl-tin complexes it acts essentially as a monodentate ligand. Thus, R3Sn(I) seems to be a soft or borderline acid while R2Sn(II) is a hard acid.

Synthesis, characterization, thermal degradation, and comparative chain dynamics studies of weak-link polysulfide polymers

Synthesis, spectroscopic and thermal characterization of two new classes of polysulfide polymers: poly[1(phenoxymethyl) ethylene polysulfide] (PPMEP), and poly [1-(phenoxy) ethylene polysulfide] (PPEP) is presented. The direct pyrolysis mass spectrometry (DP-MS) technique, used to study the thermal degradation behavior of these polysulfide polymers, indicated that the polymers underwent degradation through the weak-links scission. The thermal stability of the polysulfide polymers decreased as the ``rank'' (number of sulfur atoms in the polysulfide linkage; n=1, 2, 4) increased. The main-chain flexibility of these polysulfide polymers in terms of their C-13 NMR spinlattice relaxation time (T-1) measurements on the backbone methine (-CH-) and methylene (-CH2-) carbons are reported here for the first time. A comparative study of the solution chain dynamics indicated that it increased as ``rank'' of the polysulfide linkages decreased as well as by introducing side chain spacers such as, ether (-O-) or methyleneoxy (-CH2O-) groups.

Metal-Nucleotide Interactions: Crystal Structures of Alkali (Li+, Na+, K+) and Alkaline Earth (Ca2+, Mg2+) Metal Complexes of Adenosine 2'-Monophosphate

Crystal structures of lithium, sodium, potassium, calcium and magnesium salts of adenosine 2'-monophosphate (2'-AMP) have been obtained at atomic resolution by X-ray crystallographic methods. 2'-AMP.Li belongs to the monoclinic space group P21 with a = 7.472(3)Å, b = 26.853(6) Å, c = 9.184(1)Å, b = 113.36(1)Å and Z= 4. 2'-AMP.Na and 2'-AMP.K crystallize in the trigonal space groups P31 and P3121 with a = 8.762(1)Å, c = 34.630(5)Å, Z= 6 and a = 8.931(4), Åc = 34.852(9)Å and Z= 6 respectively while 2'-AMP.Ca and 2'-AMP.Mg belong to space groups P6522 and P21 with cell parameters a = 9.487(2), c = 74.622(13), Z = 12 and a = 4.973(1), b = 10.023(2), c = 16.506(2), beta = 91.1(0) and Z = 2 respectively. All the structures were solved by direct methods and refined by full matrix least-squares to final R factors of 0.033, 0.028, 0.075, 0.069 and 0.030 for 2'-AMP.Li, 2'-AMP.Na, 2'- AMP.K, 2'-AMP.Ca and 2'-AMP.Mg, respectively. The neutral adenine bases in all the structures are in syn conformation stabilized by the O5'-N3 intramolecular hydrogen bond as in free acid and ammonium complex reported earlier. In striking contrast, the adenine base is in the anti geometry (cCN = -156.4(2)°) in 2'-AMP.Mg. Ribose moieties adopt C2'-endo puckering in 2'-AMP.Li and 2'-AMP.Ca, C2'-endo-C3'-exo twist puckering in 2'-AMP.Na and 2'-AMP.K and a C3'-endo-C2'-exo twist puckering in 2'-AMP.Mg structure. The conformation about the exocyclic C4'-C5' bond is the commonly observed gauche-gauche (g+) in all the structures except the gauche- trans (g-) conformation observed in 2'-AMP.Mg structure. Lithium ions coordinate with water, ribose and phosphate oxygens at distances 1.88 to 1.99Å. Na+ ions and K+ ions interact with phosphate and ribose oxygens directly and with N7 indirectly through a water oxygen. A distinct feature of 2'-AMP.Na and 2'-AMP.K structures is the involvement of ribose O4' in metal coordination. The calcium ion situated on a two-fold axis coordinates directly with three oxygens OW1, OW2 and O2 and their symmetry mates at distances 2.18 to 2.42Å forming an octahedron. A classic example of an exception to the existence of the O5'-N3 intramolecular hydorgen bond is the 2'-AMP.Mg strucure. Magnesium ion forms an octahedral coordination with three water and three phosphate oxygens at distances ranging from 2.02 to 2.11Å. A noteworthy feature of its coordination is the indirect link with N3 through OW3 oxygen resulting in macrochelation between the base and the phosphate group. Greater affnity of metal clays towards 5' compared to 2' and 3' nucleotides (J. Lawless, E. Edelson, and L. Manring, Am. Chem. Soc. Northwest Region Meeting, Seattle. 1978) due to macrochelation infered from solution studies (S. S. Massoud, H. Sigel, Eur. J. Biochem. 179, 451-458 (1989)) and interligand hydrogen bonding induced by metals postulated from metal-nucleotide structures in solid state (V. Swaminathan and M. Sundaralingam, CRC. Crit. Rev. Biochem. 6, 245-336 (1979)) are borne out by our structures also. The stacking patterns of adenine bases of both 2'-AMP.Na and 2'-AMP.K structures resemble the 2'-AMP.NH4 structure reported in the previous article. 2'-AMP.Li, 2'-AMP.Ca and 2'-AMP.Mg structures display base-ribose O4' stacking. An overview of interaction of monovalent and divalent cations with 2' and 5'-nucleotides has been presented.

Ab initio molecular dynamics simulations of excited hydrogen halides and methyl halides

Using excited-state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field approach, we study the mechanism of photodissociation in terms of time evolution of structure, kinetic energy, charges and potential energy for the first excited state of hydrogen halides and methyl halides. Although the hydrogen halides and methyl halides are similar in the photodissociation mechanism, their dynamics are slightly different. The presence of the methyl group causes delay in photodissociation as compared to hydrogen halides.

Structure, dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X-ray studies

Molecular dynamics simulations have been carried out on all the jacalin-carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X-ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin-carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate-binding residues are consistent with the known thermodynamic parameters of jacalin-carbohydrate interactions. The simulations, along with X-ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin-carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit. Proteins 2009; 77:760-777.

Deciphering the Distinct Role for the Metal Coordination Motif in the Catalytic Activity of Mycobacterium smegmatis Topoisomerase I

Mycobacterium smegmatis topoisomerase I (Mstopol) is distinct from typical type IA topoisomerases. The enzyme binds to both single- and double-stranded DNA with high affinity, making specific contacts. The enzyme comprises conserved regions similar to type IA topoisomerases from Escherichia coli and other eubacteria but lacks the typically found zinc fingers in the carboxy-terminal domain. The enzyme can perform DNA cleavage m the absence of Mg2+ but religation needs exogenously added Mg2+. One molecule of Mg2+ tightly bound to the enzyme has no role in DNA cleavage but is needed only for the religation reaction. The toprim. (topoisomerase-primase) domain in MstopoI comprising the Mg2+ binding pocket, conserved in both type IA and type II topoisomerases, was subjected to mutagenesis to understand the role of Mg2+, in different steps of the reaction. The residues D108, D110, and E112 of the enzyme, which form the acidic triad in the DXDXE motif, were changed to alanines. D108A mutation resulted in an enzyme that is Mg2+ dependent for DNA cleavage unlike Mstopol and exhibited enhanced DNA cleavage property and reduced religation activity. The mutant was toxic for cell growth, most likely due to the imbalance in cleavage-religation equilibrium. In contrast, the E112A mutant behaved like wild-type enzyme, cleaving DNA in a Mg2+-independent fashion, albeit to a reduced extent. Intra- and intermolecular religation assays indicated specific roles for D108 and E112 residues during the reaction. Together, these results indicate that the D108 residue has a major role during cleavage and religation, while E112 is important for enhancing the efficiency of cleavage. Thus, although architecturally and mechanistically similar to topoisomerase I from E. coli, the metal coordination pattern of the mycobacterial enzyme is distinct, opening up avenues to exploit the enzyme to develop inhibitors.

Ligand dependent intra and inter subunit communication in human tryptophanyl tRNA synthetase as deduced from the dynamics of structure networks

Homodimeric protein tryptophanyl tRNA synthetase (TrpRS) has a Rossmann fold domain and belongs to the 1c subclass of aminoacyl tRNA synthetases. This enzyme performs the function of acylating the cognate tRNA. This process involves a number of molecules (2 protein subunits, 2 tRNAs and 2 activated Trps) and thus it is difficult to follow the complex steps in this process. Structures of human TrpRS complexed with certain ligands are available. Based on structural and biochemical data, mechanism of activation of Trp has been speculated. However, no structure has yet been solved in the presence of both the tRNA(Trp) and the activated Trp (TrpAMP). In this study, we have modeled the structure of human TrpRS bound to the activated ligand and the cognate tRNA. In addition, we have performed molecular dynamics (MD) simulations on these models as well as other complexes to capture the dynamical process of ligand induced conformational changes. We have analyzed both the local and global changes in the protein conformation from the protein structure network (PSN) of MD snapshots, by a method which was recently developed in our laboratory in the context of the functionally monomeric protein, methionyl tRNA synthetase. From these investigations, we obtain important information such as the ligand induced correlation between different residues of this protein, asymmetric binding of the ligands to the two subunits of the protein as seen in the crystal structure analysis, and the path of communication between the anticodon region and the aminoacylation site. Here we are able to elucidate the role of dimer interface at a level of detail, which has not been captured so far.

A comparative molecular dynamics study of diffusion of n-decane and 3-methyl pentane in Y zeolite

Molecular dynamics simulations are reported on the structure and dynamics of n-decane and 3-methylpentane in zeolite NaY. We have calculated several properties such as the center of mass-center of mass rdf, the end-end distance distribution, bond angle distribution and dihedral angle distribution. We have also analysed trajectory to obtain diffusivity and velocity autocorrelation function (VACF). Surprisingly, the diffusivity of 3-methylpentane which is having larger cross-section perpendicular to the long molecular axis is higher than n-decane at 300 K. Activation energies have been obtained from simulations performed at 200 K, 300 K, 350 K, 400 K and 450 K in the NVE ensemble. These results can be understood in terms of the previously known levitation effect. Arrhenious plot has higher value of slope for n-decane (5 center dot 9 kJ/mol) than 3-methylpentane (3 center dot 7 kJ/mol) in agreement with the prediction of levitation effect.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

Conformational study of valinomycin: a molecular dynamics approach

Valinomycin is a highly flexible cyclic dodecadepsipeptide that transports ions across membranes. Such a flexibility in the conformation is required for its biological function since it has to encounter a variety of environments and liganding state. Exploration of conformational space of this molecule is therefore important and is one of the objectives of the present study that has been carried out by means of high temperature Molecular Dynamics. Further, the stability of the known bracelet-like structure of the uncomplexed valinomycin and the inherent flexibility around this structure has been investigated. The uncomplexed form of valinomycin has been simulated at 75–100 K for 1 ns in order to elucidate the average conformational properties. An alanine-analog of valinomycin has been simulated under identical conditions in order to evaluate the effect of sidechain on the conformational properties, The studies confirm the effect of sidechain on conformational equilibrium.

Buckling, stiffening, and negative dissipation in the dynamics of a biopolymer in an active medium

We present a generic theory for the dynamics of a stiff filament under tension, in an active medium with orientational correlations, such as a microtubule in contractile actin. In sharp contrast to the case of a passive medium, we find the filament can stiffen, and possibly oscillate or buckle, depending on both the contractile or tensile nature of the activity and the filament-medium anchoring interaction. We also demonstrate a strong violation of the fluctuation-dissipation (FD) relation in the effective dynamics of the filament, including a negative FD ratio. Our approach is also of relevance to the dynamics of axons, and our model equations bear a remarkable formal similarity to those in recent work [Martin P, Hudspeth AJ, Juelicher F (2001) Proc Natl Acad Sci USA 98: 14380-14385] on auditory hair cells. Detailed tests of our predictions can be made by using a single filament in actomyosin extracts or bacterial suspensions.

Novel 3-dimensional sixfold interpenetrating diamondoid networks of copper(I) coordination polymers of polypyridyl ligands - Syntheses,characterization and crystal structures

Two new coordination polymers [Cu(L-1)(2)](n)(ClO4)(n)center dot 2nH(2)O (1), [Cu(L-2)(2)](n)(ClO4)(n)center dot 2nH(2)O (2) of polydentate imine/pyridyl ligands, L-1 and L-2 with Cu(I) ion have been synthesized and characterized by single crystal X-ray diffraction studies, elemental analyses, IR' UV-vis and NMR spectroscopy. They represent 3-dimensional, sixfold interpenetrating diamondoid network structures having large pores of dimension, 35 x 21 angstrom(2) in 1 and 38 x 19 angstrom(2) in 2, respectively.

Femtosecond carrier dynamics and saturable absorption in graphene suspensions

Nonlinear optical properties and carrier relaxation dynamics in graphene, suspended in three different solvents, are investigated using femtosecond (80 fs pulses) Z-scan and degenerate pump-probe spectroscopy at 790 nm. The results demonstrate saturable absorption property of graphene with a nonlinear absorption coefficient, beta of (similar to 2-9) x 10(-8) cm/W. Two distinct time scales associated with the relaxation of photoexcited carriers, a fast one in the range of 130-330 fs (related to carrier-carrier scattering) followed by it slower one in 3.5-4.9 ps range (associated with carrier-phonon scattering) are observed. (C) 2009 American Institute of Physics.