Computer simulation of polypeptide translocation through a nanopore

A simplified model of polypeptide chains was designed and studied by means of computer simulations. Chains were represented by a sequence of united atoms located at the positions of the -carbons. A further assumption was the lattice approximation for the chains. We used a (310) lattice, which was found useful for studying properties of proteins. The force field used consisted of a long-range contact potential between amino-acid residues and a local preference for forming -helical states. The chain consisted of two kinds of residues: hydrophilic (P) and hydrophobic (H) ones forming model helical septets –HHPPHPP– in a sequence. The chains were placed near an impenetrable surface with a square hole in it. The size of the hole was comparable or smaller than the size of a chain. The properties of these model chains were determined using the Monte-Carlo simulation method. During the simulations, translocation of the chain through the hole in the wall was observed. The influence of the chain length, the temperature differences on both sides of the wall and the force field on the chain properties were investigated. It was shown that the translocation time scales as N^2.2 and it was found that the presence of the local helical potential significantly slows down the process of translocation.Figure: The snapshots of typical chains conformation obtained during the simulation for chain consisted of N = 60. The values of the local potential _loc = -8.     

Theoretical potential scans and vibrational spectra of vinyl selenonyl halides CH2=CH–SeO2X (X is F, Cl and Br)

The structure and conformational stability of vinyl selenonyl fluoride, chloride and bromide CH₂=CH–SeO₂X (X is F, Cl and Br) were investigated using density functional B3LYP/6-311+G** and ab initio MP2/6-311+G** calculations. From the calculations the molecules were predicted to exist only in the non-planar gauche conformation with the vinyl C=C group almost eclipsing one of the selenonyl Se=O bonds as a result of conjugation between the two moieties. Single-minimum potential scans were calculated at the DFT level for the molecules. The vibrational frequencies were computed using B3LYP/6-311+G**. Normal coordinate calculations were then carried out and potential energy distributions were calculated for the three molecules in the gauche conformation.Figure Potential function for the asymmetric torsion in vinyl selenonyl fluoride (dotted line), chloride (dashed line) and bromide (solid line) as determined at the DFT-B3LYP/6-311+G** level     

An unusual feature of end-substituted model carbon (6,0) nanotubes

We have examined the effects of substituents on the computed electrostatic potentials V_S(r) and average local ionization energies on the surfaces of model carbon nanotubes of the types (5,5), (6,1) and (6,0). For the (5,5) and the (6,1), the effects upon both V_S(r) and of substituting a hydroxyl group at one end are primarily localized to that part of the system. For the (6,0) tube, however, a remarkable change is observed over its entire length, with V_S(r) showing a marked gradation from strongly positive at the substituted end to strongly negative at the other; correspondingly goes from higher to lower values. Replacing OH by another resonance- donor, NH₂, produces similar results in the (6,0) system, while the resonance withdrawing NO₂ does the opposite, but in equally striking fashion. We explain these observations by noting that the arrangement of the C–C bonds in the (6,0) tube facilitates charge delocalization over the full length and entire surface of the tube. Substituting NH₂ and NO₂ at opposite ends of the (6,0) tube greatly strengthens the gradations in both V_S(r) and The first hyperpolarizability of this system was found to be nine times that of para-nitroaniline, suggesting possible nonlinear optical applications.Figure HF/STO-5G electrostatic potential on outer surface of open (6,0) C₇₂H₁₀NH₂NO₂. The nitro group is at the right end of the tube, the amino group at the left. In eV: purple is less than 14, blue is between 14 and 15, green is between 15 and 16.5, yellow is between 16.5 and 17.5, and red is more than 17.5.      

Structure of glutathione S-transferase of the filarial parasite Wuchereria bancrofti: a target for drug development against adult worm

A three dimensional structural model of Glutathione-S-transferase (GST) of the lymphatic filarial parasite Wuchereria bancrofti (wb) was constructed by homology modeling. The three dimensional X-ray crystal structure of porcine -class GST with PDB ID: 2gsr-A chain protein with 42% sequential and functional homology was used as the template. The model of wbGST built by MODELLER6v2 was analyzed by the PROCHECK programs. Ramachandran plot analysis showed that 93.5% of the residues are in the core region followed by 5.4 and 1.1% residues in the allowed and generously allowed regions, respectively. None of the non-glycine residues is in disallowed regions. The PROSA II z-score and the energy graph for the final model further confirmed the quality of the modeled structure. The computationally modeled three-dimensional (3D) structure of wbGST has been submitted to the Protein Data Bank (PDB) (PDB ID: 1SFM and RCSB ID: RCSB021668). 1SFM was used for docking with GST inhibitors by Hex4.2 macromolecular docking using spherical polar Fourier correlations.Figure: A three-dimensional (3D) structure of Glutathione-S-transferase (GST) of the lymphatic filarial parasite Wuchereria bancrofti (wb) was constructed by homology modeling. This modeled 3D structure of wbGST has been submitted to the Protein Data Bank (PDB) (PDB ID: 1SFM and RCSB ID: RCSB021668).     

Organization of rhodopsin molecules in native membranes of rod cells–an old theoretical model compared to new experimental data

It has been shown that rhodopsin forms an oligomer in the shape of long double rows of monomers. Because of the importance of rhodopsin as a template for all G protein-coupled receptors, its dimeric, tetrameric and higher-oligomeric structures also provide a useful pattern for similar structures in GPCRs. New experimental data published recently are discussed in the context of a proposed model of the rhodopsin oligomer 1N3M deposited in the protein data bank. The new rhodopsin structure at 2.2 Å resolution with all residues resolved as well as an electron cryomicroscopy structure from 2D crystals of rhodopsin are in agreement with the 1N3M model. Accommodation of movement of transmembrane helix VI, regarded as a major event during the activation of rhodopsin, in a steady structure of the oligomer is also discussed.Figure Superimposition of the 1U19 (red wire), 1GZM (purple wire) and 1N3M (blue wire) rhodopsin structures. Size of the wires is proportional to thermal factors of backbone C atoms, view parallel to the membrane.      

Monte Carlo studies of two-dimensional polymer–solvent systems

The static properties of two-dimensional athermal polymer solutions were studied by performing Monte Carlo lattice simulations using the cooperative motion algorithm (CMA) and taking into account the presence of explicit solvent molecules. The simulations were performed for a wide range of polymer chain lengths N (16–1024) and concentrations φ (0.0156–1). The results obtained for short chains (N?<?256) were in good agreement with those given by previous simulations. For the longest chains (512 or 1024 beads), some unexpected behavior was observed in the dilute and semidilute regimes. A pronounced change in the concentration dependence of chain size and shape was observed below a certain critical concentration (0.6 for the longest chains under consideration, consisting of 1024 beads). Longer chains became more extended below this concentration. The behavior of the single-chain structure factor confirmed these changes in the fractal dimension of the chain as a function of the concentration. The observed phenomena are related to the excluded volume of solvent molecules, which causes the chain statistics to be modified in the vicinity of other chains; this effect is important in strictly 2D systems.

Open image in new window

Graphical abstract Extended long chains at moderate density with solvent molecules inside coils.

     

Theoretical studies of the conformational behavior of chain molecules containing polar groups: simulations of a poly(vinylidene fluoride) model

Atomistic Monte Carlo simulations have been conducted to elucidate the conformational behavior of a single chain molecule containing polar functional groups. Here, we resort to an atomistic poly(vinylidene fluoride) (PVDF) chain model as a representative example. The model is modified in such a way that bond lengths and bond angles are fixed, aiming to manifest the role of dipolar interactions. For a given chain length, chain conformation is sensitive to two environmental parameters, temperature and dielectric constant. The mean chain size increases when temperature and/or dielectric constant are increased. The conformational behavior is further characterized by chain size distribution function, and our findings show that temperature induced conformational transition for a chain molecule can be discrete or continuous, depending on its chain length. Also, the dipolar interactions in PVDF are effectively attractive, and enhance chain contraction. As a result, when the strength of dipolar interactions is increased, the discrete conformational transition shifts toward longer chains; and for a given chain length, such a transition occurs at higher temperatures.Figure Variation of R² with temperature for different dielectric constants =1 and 8, denoted by dotted and solid lines, respectively, and for different chain lengths M=8 and 12, as marked. Lines are meant for eye guidance     

The structure of polymer chains in confinement. A Monte Carlo study

A coarse-grained model of polymer star chains confined in two parallel impenetrable surfaces, which were attractive for polymer beads was studied. The flexible homopolymer chains were built of united atoms whose positions in space were restricted to vertices of a simple cubic lattice. The chains were modeled in good solvent conditions and, thus, there were no long-range specific interactions between polymer beads—only the excluded volume was present. The influence of the polymer density and the distances between the confining surfaces on the properties of star-branched polymers was studied. It is shown that the chains adsorbed on one surface could change their position so that they swap between both surfaces with frequency depending on the size of the slit and on the density of the system only. The increase of the polymer density diminished the frequency of jumps and caused that chains became only partially adsorbed. The analysis of structural elements of chains showed that the increase of the density of the system leads to increase of the number of bridges connecting the two adsorbing surfaces, thus, the frequency of jumps between them decreases.     

Analysis of the effect of electrostatic energy truncation in molecular dynamics simulations of immunoglobulin G light chain dimer

Molecular dynamics (MD) simulations of immunoglobulin G (IgG) light chain dimer using particle mesh Ewald (PME) and cutoff methods of treating electrostatic interactions were performed. The results indicate that structural parameters (RMSD, radius of gyration, solvent accessible surface) are very similar for both schemes; however, PME simulation shows increased mobility of side chains. This leads to larger fluctuations in the distance between the monomers in the dimer molecule, and, as a consequence, results in decreased number of interactions across the dimer interface. The wall clock time of the simulations was also compared. It was shown that the PME method is approximately 30% faster than the cutoff method for the system studied on a single processor.Figure Backbone order parameters for PME (red) and cutoff (green) calculations. Thick, horizontal lines show stable secondary structures     

Characterization of molecular structures and properties of polyurethanes using molecular dynamics simulations

Thermotropic polyurethanes with mesogenic groups in side chains were prepared from two diisocyanates and four diols with stoichiometric ratios of reactive isocyanate (NCO) and hydroxy (OH) groups. Their thermal behavior was determined by differential scanning calorimetry. The effect of structure modifications of the diisocyanates and diols, in particular changes in the mesogen, were investigated. Introduction of mesogenic segments into the polymers suppresses the ordering. Stiff end substituents (phenyl and alkoxy groups) of the mesogens stabilize the mesophases to such an extent that the negative influence of long polymer chains is compensated and the liquid-crystalline properties are recovered. All-atom molecular dynamics simulations in the Cerius² modeling environment were carried out to characterize the structures of the polymers. Analysis of the dynamic trajectories at 20, 100, 120 and 170 °C revealed changes in conformation of macromolecules, which correlate with DSC measurements.Figure Example of structure relaxation of D4/TDI molecule at indicated simulation times (temperature 20 °C): a complete structure; b backbone structure; c top view of molecule     

Local molecular properties and their use in predicting reactivity

Expressions for the local electron affinity, electronegativity and hardness are derived in analogy to the local ionization energy introduced by Sjoberg, Murray and Politzer. The local polarizability is also defined based on an additive atomic orbital polarizability model that uses Rivail's variational technique. The characteristics of these local properties at molecular surfaces and their relevance to electrophilic aromatic substitution, to S_N2 reactivity and to the nucleophilicity of enolate ions are discussed.Figure The local ionization energy at the SES surfaces of methyl benzoate. The color scale ranges from 375 (blue) to 550 kcal mol^–1 (red). The blue areas are those for which interaction with an acceptor is most favorable.     

The design of organic catalysis for epoxidation by hydrogen peroxide

The potential of various organic species to catalyze epoxidation of ethene by hydrogen peroxide is explored with B3LYP/6-31G* DFT calculations. Electronic Supplementary Materials Supplementary material is available for this article at http://dx.doi.org/10.1007/s00894-005-0044-4.Dedicated to Professor Dr. Paul von Ragué Schleyer on the occasion of his 75th birthday     

Reinvestigation of molecular structure and barrier to internal rotation of pyridinium N-phenolate betaine dye

In the present paper, the results of a systematic theoretical study of the molecular structure of 4-(1-pyridinium-1-yl)phenolate betaine are reported. The ground-state molecular structure and the barrier to internal rotation of the betaine dye molecule were calculated ab inito (with Hartree–Fock theory and the second-order of Möller–Plesset method) and with density functional theory (DFT). In order to estimate the complete basis set limit, the calculations of barriers to internal rotations were performed using correlation–consistent basis sets with a maximal cardinal number of four. It was determined that electron correlation is crucial in order to obtain reliable geometries and rotational barriers of the molecule investigated. For the sake of comparison, the results of calculations using the AM1 Hamiltonian are also presented.Figure Investigated betaine dye.     

On the nature of bonding in HCOOH...Ar and HCOOH...Kr complexes

The chemical interaction in HCOOH...Ng (Ng=Ar, Kr) complex was analyzed by topological analysis of the electron density based on Atoms-In-Molecules theory. For all computationally stable equilibrium structures of 1:1 HCOOH...Ng complexes, an intermolecular bond path with a bond critical point was found and perturbation of formic acid (FA) atomic basins and electron density was observed. The intermolecular interaction between the two complex subunits can be classified, according to its nature, as a closed-shell van der Waals type of interaction. However, one of the computed structures (complex II), pictures a noble gas atom attached linearly to the acidic O–H tail of FA. In this particular case, the electron density at the intermolecular bond critical point was found to resemble a hydrogen-bonded system and thus, may be termed a hydrogen-bond-like interaction. This change in the nature of the interaction is also shown by large perturbations of the FA properties found for this complex structure. The structural and vibrational perturbations are larger than for the other three structures and they increase for the Kr complexes compared to the Ar complex.Figure Electron density analysis of HCOOH...Ng (Ng=Ar,Kr) complex.     

A study on the influence of molecular properties in the psychoactivity of cannabinoid compounds

Several molecular properties are calculated for a set of 26 cannabinoid compounds with the goal of connecting the psychoactivity of the compounds with an appropriate set of calculated properties. For this purpose we used quantum chemical (the AM1 semi-empirical method) and chemometric methods. The AM1 method was employed to calculate the set of quantum chemical molecular properties and the chemometric methods were employed with the aim of selecting the most relevant properties to be correlated with psychoactivity. The chemometric methods used were Principal Component Analysis (PCA), Hierarchical Cluster Analysis (HCA) and the K-Nearest Neighbor (KNN) method. The chemometric analysis showed that an electronic property (energy of LUMO), a hydrophobic property (log P), a steric property (volume of the substituent at the C4 position) and a topological property (Lovasz–Pelikan index) were the most important variables for the separation between the psychoactive and psychoinactive compounds. In order to validate our PCA, HCA and KNN results, eight new cannabinoid compounds (with known psychoactivity) were used in a prediction study and were classified correctly by the methods used in this work, indicating that our PCA, HCA and KNN models are able to predict reliable psychoactivity of cannabinoid compounds. Figure: ⁹-THC This revised version was published online in June 2005 with corrections to Table 1.     

Quantitative structure–activity relationship study on some benzodiazepine derivatives as anti-Alzheimer agents

A QSAR study was performed in an attempt to explore the pharmacophore of some benzodiazepine derivatives as anti-Alzheimer agents for the inhibition of -secretase. The study, which used the electrotopological state atom (ETSA) index, which encodes electronic and topological information, reveals the importance of two phenyl rings—one substituted and another unsubstituted, for the inhibition of the enzyme. Fluorine substitution on the substituted phenyl ring has an important contribution to the activity. R-configurations of the aliphatic chain substituents provide the exact conformation of the molecules to enter into the binding pockets of the receptor(s).Figure General structure of benzodiazepine containing -secretase inhibitors     

Effect of Pt and Sn on the adsorption of n-heptane in γ-Al2O3 catalyst models

The Grand Canonical Monte-Carlo (GCMC) method has been used to carry out simulations of the adsorption of n-heptane in models of naphtha-reforming catalysts. Models used in the study differed in the number and distribution of metal atoms—Pt and Sn. The number of adsorbed n-heptane molecules grows linearly with increasing number of metal atoms. The effect of Pt content on the adsorption of n-heptane molecules is most distinct at approximately 100 kPa and within the lower range of the temperatures investigated. In the models of bimetallic catalysts, the effect of the two metals is additive.Figure Effect of Pt and Sn on number of n-heptane molecules adsorbed in Al₂O₃ catalyst in 773 K and 1000 kPa.      

Immature erythroid cells with novel morphology and cytoskeletal organization in adultXenopus

Summary Nucleated erythrocytes of non-mammalian vertebrates are a useful model system for studying the correlation between changes in cell shape and cytoskeletal organization during cellular morphogenesis. They are believed to transform from spheres to flattened discs to ellipsoids. Our previous work on developing erythroblasts suggested that pointed cells containing incomplete, pointed marginal bands (MBs) of microtubules might be intermediate stages in the larval axolotl. To test whether the occurrence of such pointed cells was characteristic of amphibian erythrogenesis, we have utilized phenylhydrazine (PHZ)-induced anemia in adultXenopus. In this system, circulating erythrocytes are destroyed and replaced by erythroblasts that differentiate in the blood, making them experimentally accessible. Thus, we followed the time-course of morphological and cytoskeletal changes in the new erythroid population during recovery. During days 7–9 post-PHZ, pointed cells did indeed begin to appear, as did spherical and discoidal cells. The percentage of pointed cells peaked at days 11–13 in different animals, subsequently declining as the percentage of elliptical cells increased. Since degenerating old erythrocytes were still present when pointed cells appeared, we tested directly whether pointed ones were old or new cells. Blood was removed via the dorsal tarsus vein, and the erythrocytes washed, fluorescently tagged, and re-injected. In different animals, 2–8% of circulating erythrocytes were labeled. Subsequent to induction of anemia in these frogs, time-course sampling showed that no pointed cells were labeled, identifying them as new cells. Use of propidium iodide revealed large nuclei and cytoplasmic staining indicative of immaturity, and video-enhanced phase contrast and anti-tubulin immunofluorescence showed that the pointed cells contained pointed MBs. The results show that pointed cells, containing incomplete, pointed MBs are a consistent feature of amphibian erythrogenesis. These cells may represent intermediate stages in the formation of elliptical erythrocytes.Abbreviations MB marginal band - MS membrane skeleton - PHZ phenylhydrazine     

EMPIRE: a highly parallel semiempirical molecular orbital program: 2: periodic boundary conditions

ᅟ

ᅟ

     

Relationships between aqueous acidities and computed surface-electrostatic potentials and local ionization energies of substituted phenols and benzoic acids

Electrostatic potentials and average local ionization energies on the molecular surfaces of 19 phenols, 17 benzoic acids and their respective conjugate bases were computed at the HF/STO-5G(d)//B3LYP/6-311G(d,p) level. Good correlations were found between pK_as and the V_S,max values of the neutral acids and the V_S,min and of the conjugate bases for both sets of molecules. V_S,max is the most positive value of the electrostatic potential on the molecular surface and is an indicator of the ease with which the phenols and benzoic acids lose their acidic hydrogens. V_S,min and are the minimum values of the electrostatic potential and the local ionization energy computed on the molecular surface; the V_S,min and of the conjugate bases of the phenoxides and benzoates are indicative, respectively, of the tendencies of electrophiles to approach the anions (V_S,min) and to react with the anions ( ) to reform the original acids. The correlations observed between the computed molecular surface quantities, taken as single parameters, and the experimental pK_a values ranged from R=0.938 to R=0.970 for the two classes of compounds.Figure Correlation between pK_a and the computed V_S,max of the benzoic acids listed in Table 2. The linear correlation coefficient is 0.970