Molecular modeling of Bt Cry1Ac (DI–DII)–ASAL (Allium sativum lectin)–fusion protein and its interaction with aminopeptidase N (APN) receptor of Manduca sexta

Genetic engineering of Bacillus thuringiensis (Bt) Cry proteins has resulted in the synthesis of various novel toxin proteins with enhanced insecticidal activity and specificity towards different insect pests. In this study, a fusion protein consisting of the DI–DII domains of Cry1Ac and garlic lectin (ASAL) has been designed in silico by replacing the DIII domain of Cry1Ac with ASAL. The binding interface between the DI–DII domains of Cry1Ac and lectin has been identified using protein–protein docking studies. Free energy of binding calculations and interaction profiles between the Cry1Ac and lectin domains confirmed the stability of fusion protein. A total of 18 hydrogen bonds was observed in the DI–DII–lectin fusion protein compared to 11 hydrogen bonds in the Cry1Ac (DI–DII–DIII) protein. Molecular mechanics/Poisson–Boltzmann (generalized-Born) surface area [MM/PB (GB) SA] methods were used for predicting free energy of interactions of the fusion proteins. Protein–protein docking studies based on the number of hydrogen bonds, hydrophobic interactions, aromatic–aromatic, aromatic–sulphur, cation–pi interactions and binding energy of Cry1Ac/fusion proteins with the aminopeptidase N (APN) of Manduca sexta rationalised the higher binding affinity of the fusion protein with the APN receptor compared to that of the Cry1Ac–APN complex, as predicted by ZDOCK, Rosetta and ClusPro analysis. The molecular binding interface between the fusion protein and the APN receptor is well packed, analogously to that of the Cry1Ac–APN complex. These findings offer scope for the design and development of customized fusion molecules for improved pest management in crop plants.

Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.

Sequence Independent Recombination Triple Helices: A Molecular Dynamics Study

Recent experimental studies have shown that the Rec-A mediated homologous recombination reaction involves a triple helical intermediate, in which the third strand base forms hydrogen bonds with both the bases in the major groove of the Watson-Crick duplex. Such 'mixed' hydrogen bonds allow formation of sequence independent triplexes. DNA triple helices involving 'mixed' hydrogen bonds have been studied, using model building, molecular mechanics (MM) and molecular dynamics (MD). Models were built for a tripler comprising all four possible triplets viz., G.C*C, C.G*G, A.T*T and T.A*A. To check the stability of all the 'mixed' hydrogen bonds in such triplexes and the conformational preferences of such tripler structures, MD studies were carried out starting from two structures with 30 degrees and 36 degrees twist between the basepairs. It was observed that though the two triplexes converged towards a similar structure, the various hydrogen bonds between the WC duplex and the third strand showed differential stabilities. An MD simulation with restrained hydrogen bonds showed that the resulting structure was stable and remained close to the starting structure. These studies help us in defining stable hydrogen bond geometries involving the third strand and the WC duplex. It was observed that in the C.G*G triplets the N7 atom of the second strand is always involved in hydrogen bonding. In the G.C*C triplets, either N3 or O2 in the third strand cytosine can interchangeably act as a hydrogen bond acceptor.

Polyelectrolyte microcapsules for sustained delivery of water-soluble drugs

Polyelectrolyte capsules composed of weak polyelectrolytes are introduced as a simple and efficient system for spontaneous encapsulation of low molecular weight water-soluble drugs. Polyelectrolyte capsules were prepared by layer-by-layer (LbL) assembling of weak polyelectrolytes, poly(allylamine hydrochloride) (PAH) and poly (methacrylic acid) (PMA) on polystyrene sulfonate (PSS) doped CaCO3 particles followed by core removal with ethylene-diaminetetraacetic add (EDTA). The loading process was observed by confocal laser scanning microscopy (CLSM) using tetramethylrhodamineisothiocyanate labeled dextran (TRITC-dextran) as a fluorescent probe. The intensity of fluorescent probe inside the capsule decreased with increase in cross-linking time. Ciprofloxacin hydrochloride (a model water-soluble drug) was spontaneously deposited into PAH/PMA capsules and their morphological changes were investigated by scanning electron microscopy (SEM) and atomic force microscopy (AFM). The quantitative study of drug loading was also elucidated which showed that drug loading increased with initial drug concentration, but decreased with increase in pH. The loaded drug was released in a sustained manner for 6 h, which could be further extended by cross-linking the capsule wall. The released drug showed significant antibacterial activity against E. coli. These findings indicate that such capsules can be potential carriers for water-soluble drugs in sustained/controlled drug delivery applications. (C) 2010 Elsevier B.V. All rights reserved.

Chain folding and A:T pairing in human telomeric DNA: a model-building and molecular dynamics study

The various types of chain folding and possible intraloop as well as interloop base pairing in human telomeric DNA containing d(TTAG(3)) repeats have been investigated by model-building, molecular mechanics, and molecular dynamics techniques. Model-building and molecular mechanics studies indicate that it is possible to build a variety of energetically favorable folded-back structures with the two TTA loops on same side and the 5' end thymines in the two loops forming TATA tetrads involving a number of different intraloop as well as interloop A:T pairing schemes. In these folded-back structures, although both intraloop and interloop Watson-Crick pairing is feasible, no structure is possible with interloop Hoogsteen pairing. MD studies of representative structures indicate that the guanine-tetraplex stem is very rigid and, while the loop regions are relatively much more flexible, most of the hydrogen bonds remain intact throughout the 350-ps in vacuo simulation. The various possible TTA loop structures, although they are energetically similar, have characteristic inter proton distances, which could give rise to unique cross-peaks in two-dimensional nuclear Overhauser effect spectroscopy (NOESY) experiments. These folded-back structures with A:T pairings in the loop region help in rationalizing the data from chemical probing and other biochemical studies on human telomeric DNA.

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.

Molecular dynamics research of G249 and S249 substitutions of p53 protein.

,The molecular dynamics research of the core domain of p53 protein crystal structure shows that besides the stability in biochemistry this domain also shows a high stability in molecular mechanics. Based on that work, the residue R249 was substituted with amino acids Gly and Ser respectively, and molecular dynamics researches were performed separately. The results show that these substitutions cause a relax tendency between loop2 and 3 domains, leading to an alteration of the whole conformation of p53 core domain and ruining its stability. The results visually explains the mechanism of p53 changes in immunological and biochemical reactions, which are caused by 249 residue substitutions from 3-D structure variations.

Molecular modeling on some human interleukins sharing gamma c of interleukin-2 receptor, and structure-function relationships

Five models for human interleukin-7 (HIL-7), HIL-9, HIL-13, HIL-15 and HIL-17 have been generated by SYBYL software package. The primary models were optimized using molecular dynamics and molecular mechanics methods. The final models were optimized using a steepest descent algorithm and a subsequent conjugate gradient method. The complexes with these interleukins and the common gamma chain of interleukin-2 receptor (IL-2R) were constructed and subjected to energy minimization. We found residues, such as Gln127 and Tyr103, of the common gamma chain of IL-2R are very important. Other residues, e.g. Lys70, Asn128 and Glu162, are also significant. Four hydrophobic grooves and two hydrophilic sites converge at the active site triad of the gamma chain. The binding sites of these interleukins interaction with the common gamma chain exist in the first helical and/or the fourth helical domains. Therefore, we conclude that these interleukins binds to the common gamma chain of IL-2R by the first and the fourth helix domain. Especially at the binding sites of some residues (lysine, arginine, asparagine, glutamic acid and aspartic acid), with a discontinuous region of the common gamma chain of IL-2R, termed the interleukins binding sites (103-210). The study of these sites can be important for the development of new drugs. (C) 2000 Elsevier Science B.V. All rights reserved.

Stability and molecular modeling of triplex DNA inhibiting DNA binding protein binding to the core promoter of hepatitis B virus.

Two three-dimensional structure models of the 21nt oligodeoxyribonucleotides, CPI (G3TG-2TGT2G5TG2TGT) and CP3 (TGTG2TGST2GTG2TG3), were constructed by InsightII (MSI) software in IRIS Indigo2 (SGI) workstation using the crystal structure of TAT tripler formation as the template. The initial structures subsequently were minimized by molecular mechanics. The final structures were believed as the dominant conformation. The results showed that the energy of CP1 is lower than that of CP3, and the former is more stable than the latter. Moreover, the results further proved that the 21nt oligodeoxyribo-nucleotide CP1 stably combines with the core promoter (Cp) fragment of hepatitis B virus (HBV) to form a tripler DNA, and CP1 specifically inhibits a specific cellular factor (DNA binding protein) binding to Cp fragment. These results indicated that specific repression of gene transcription of HBV DNA might be possible by tripler-formation DNA.

Novel inclusion complex of ibuprofen tromethamine with cyclodextrins: Physico-chemical characterization

Guest-host interactions of ibuprofen tromethamine salt (Ibu.T) with native and modified cyclodextrins (CyDs) have been investigated using several techniques, namely phase solubility diagrams (PSDs), proton nuclear magnetic resonance (H-1 NMR), differential scanning calorimetry (DSC), Fourier transform infrared spectroscopy (FT-IR), X-ray powder diffractometry (XRPD). scanning-electron microscopy (SEM) and molecular mechanics (MM). From the analysis of PSD data (A(L)-type) it is concluded that the anionic tromethamine salt of ibuprofen (pK(a) = 4.55) forms 1: 1 soluble complexes with all CyDs investigated in buffered water at pH 7.0, while the neutral form of Ibu forms an insoluble complex with beta-CyD (B-S-type) in buffered water at pH 2.0. Ibu.T has a lower tendency to complex with beta-CyD (K-11 = 58 M-1 at pH 7.0) compared with the neutral Ibu (K-11 = 4200 M (1)) in water. Complex formation of Ibu.T with beta-CyD (Delta G degrees = -20.4 kJ/mol) is enthalpy driven (Delta H degrees = -22.9 kJ/mol) and is accompanied by a small unfavorable entropy (Delta S degrees = -8.4 J/mol K) change. H-1 NMR studies and MM computations revealed that, on complexation, the hydrophobic central benzene ring of lbu.T and part of the isobutyl group reside within the beta-CyD cavity leaving the peripheral groups (carboxylate, tromethamine and methyl groups) located near the hydroxyl group networks at either rim of beta-CyD. PSD, H-1 NMR, DSC, FT-IR, XRPD, SEM and MM studies confirmed the formation of Ibu.T/beta-CyD inclusion complex in solution and the solid state. (C) 2009 Elsevier B.V. All rights reserved.

Cation complexation by chemically modified calixarenes. 11. Complexation and extraction of alkali cations by calix[5]- and -[6]arene ketones. Crystal and molecular structures of calix[5]arene ketones and Na+ and Rb+ complexes

A series of four calix[5]arenes and three calix[6]arenes (R-calixarene-OCH2COR1) (R = H or Bu-t) with alkyl ketone residues (R-1 = Me or Bu-t) on the lower rim have been synthesized, and their affinity for complexation of alkali cations has been assessed through phase-transfer experiments and stability constant measurements. The conformations of these ketones have been probed by H-1 NMR and X-ray diffraction analysis, and by molecular mechanics calculations. Pentamer 3 (R R-1 = Bu-t) possesses a symmetrical cone conformation in solution and a very distorted cone conformation in the solid state. Pentamer 5 (R = H, R-1 = Bu-t) exists in a distorted 1,2-alternate conformation in the solid state, but in solution two slowly interconverting conformations, one a cone and the other presumed to be 1,2-alternate, can be detected. X-ray structure analysis of the sodium and rubidium perchlorate complexes of 3 reveal the cations deeply encapsulated by the ethereal and carbonyl oxygen atoms in distorted cone conformations which can be accurately reproduced by molecular mechanics calculations. The phase-transfer and stability constant data reveal that the extent of complexation depends on calixarene size and the nature of the alkyl residues adjacent to the ketonic carbonyls with tert-butyl much more efficacious than methyl.

Molecular modeling of calixarenes with group I metal ions

Molecular mechanics calculations have been used to model the geometries of the complexes of Group I metal ions with calix[n]arenes (n = 4,5). A simple procedure in which the calixarene atoms are assigned partial charges on the basis of AM1 calculations and the metal ions are allowed to bind electrostatically to the calixarenes produces surprising good results when the resulting structures are compared to known crystallographic data on the complexes. Encapsulated solvent molecules and/or counterions can be included in the calculations and, indeed, are necessary to reproduce the X-ray data.

Volumetric Properties, Viscosities, and Isobaric Heat Capacities of Imidazolium Octanoate Protic Ionic Liquid in Molecular Solvents

Densities (F), viscosities (?), and isobaric heat molar capacities (Cp) of binary mixtures containing imidazolium octanoate, [Im][C7CO2], a protic ionic liquid (PIL), with four molecular solvents, water, acetonitrile, ethanol, and 1-octanol, are determined as a function of temperature from (298.15 to 323.15) K and within the whole composition range at atmospheric pressure. Excess molar volumes, VE, excess molar heat capacities, Cp E, and the deviation from additivity rules of viscosities, ??, of imidazolium octanoate solutions were then deduced from the experimental results, as well as apparent molar volumes, Vfi, and partial molar volumes, V j m,i. Results are discussed according to the nature of the interaction between the PIL and the molecules and the effect of temperature. The excess Gibbs energies of activation of viscous flow (?G*E) for these systems were then calculated at 298.15 K. The excess isobaric heat capacities, Cp E, of binary ([Im][C7CO2] + solvent) systems, depend also of the nature of the molecular solvent in mixture. The excess properties were then correlated, at each temperature, as a function of composition by a Redlich-Kister-type equation. Finally results have been discussed in terms of molecular interactions and molecular structures in these binary mixtures, and thermodynamic properties of investigated binary mixtures were then compared to literature values together to investigate the impact of the nature of the solvent on these reported properties.

Volumetric properties, viscosity and refractive index of the protic ionic liquid, pyrrolidinium octanoate, in molecular solvents

Densities ([rho]) and viscosities ([eta]) of binary mixtures containing the Protic Ionic Liquid (PIL), pyrrolidinium octanoate with five molecular solvents: water, methanol, ethanol, n-butanol, and acetonitrile are determined at the atmospheric pressure as a function of the temperature and within the whole composition range. The refractive index of all mixtures (nD) is measured at 298.15†K. The excess molar volumes VE and deviation from additivity rules of viscosities [eta]E and refractive index [Delta][phi]n, of pyrrolidinium octanoate solutions were then deduced from the experimental results as well as apparent molar volumes V[phi]i, partial molar volumes and thermal expansion coefficients [alpha]p. The excess molar volumes VE are negative over the entire mole fraction range for mixture with water, acetonitrile, and methanol indicating strong hydrogen-bonding interaction for the entire mole fraction. In the case of longest carbon chain alcohols (such as ethanol and n-butanol)†+†pyrrolidinium octanoate solutions, the VE variation as a function of the composition describes an S shape. The deviation from additivity rules of viscosities is negative over the entire composition range for the acetonitrile, methanol, ethanol, and butanol, and becomes less negative with increasing temperature. Whereas, [eta]E of the {[Pyrr][C7CO2]†+†water} binary mixtures is positive in the whole mole fraction range and decreases with increasing temperature. the excess Gibbs free energies of activation of viscous flow ([Delta]G*E) for these systems were calculated. The deviation from additivity rules of refractive index [Delta][phi]n are positive over the whole composition range and approach a maximum of 0.25 in PIL mole fraction for all systems. The magnitude of deviation for [Delta][phi]n describes the following order: water†>†methanol†>†acetonitrile†>†ethanol. Results have been discussed in terms of molecular interactions and molecular structures in these binary mixtures.

Rotation-electronic interaction in the 3p-complex Rydberg state of water

An energy theory is formulated for the rotational energy levels in a p-complex Rydberg state of an asymmetric top molecule of symmetry C2v. The effective Hamiltonian used consists of the usual rigid rotor Hamiltonian augmented with terms representing electronic spin and orbital angular momentum effects. Criteria for assigning symmetry species to the rotational energy levels, following Houganfs scheme that uses the full molecular group,are established and given in the form of a table. This is particularly suitable when eigenvectors are calculated on a digital computer. Also, an intensity theory for transitions to the Rydberg p-complex singlet states is presented and selection rules in terms of symmetry species of energy states are established. Finally, applications to HpO and DpO are given.