Potentials for Molecular Dynamics Simulation of Silicate Glasses

A new potential model has been developed for the simulation of amorphous silica based on the ab initio potential model of Pyper. This model promises to be of value in the simulation of silica at high pressures.     

Theoretical mechanism of the formation of cholesteryl chloride from cholesterol and thionyl chloride

This work describes a theoretical approach to the substitution reaction mechanism involving the conversion of cholesterol to cholesteryl chloride. Two chlorosulfite ester molecules were formed as intermediates. An iso-steroid was found as the transition state. The final product was cholesteryl chloride and the side products were HCl and SO₂. Calculations were carried out at high level Hartree–Fock theory, using the 6–31G* basis set. From the electronic structure of the reactants, the most important physicochemical properties involved in the reaction pathway were used. Thus, to determine the participation of each molecule and to explain the mechanism of reaction; the total energy, HOMO and LUMO, atomic orbital contribution to frontier orbitals formation, electrostatic potentials, atomic charges, hardness and dipole moment were used. Characterization of intermediates and transition state was supported by their respective energy minima, fundamental frequencies and equilibrium geometry.Figure Synopsis of the reaction pathway. The reaction starts when the lone pair of the Ch oxygen interacts with the sulfur atom, releasing a chloride ion. As a result, the first intermediate is formed. Next, in the first intermediate the nucleophilic chloride ion bonds the electrophilic hydrogen atom, releasing HCl and yielding the second intermediate. In the second intermediate, the electrophilic H-atom from HCl bonds with the lone pair of the Cl atom adjacent to the sulfur atom, restoring HCl. Concurrently, SO₂ is liberated and causes the formation of the C3-C5 partial bond and breaking of the C5-C6 -bond leading to the transition state. In the transition state, the electrophilic H from HCl bonds with the Cl lone pair at C6-Cl, forming HCl again and leaving the C6 atom electron-deficient, which restores the C5-C6 -bond and breaks the C3-C5 partial bond. Finally, the electrophilic C3 atom and the nucleophilic Cl atom form a bond, yielding cholesteryl chloride. HCl and SO₂ are also formed as side products. The arrows show the rearrangement of electrons.     

The Implementation of the Iterative Diagonalization Scheme and ab initio Molecular Dynamics Simulation with the LAPW Method

Abstract

A new ab initio molecular dynamics method based on the full-potential linearized-augmented-plane-wave (LAPW) basis set has been implemented. The LAPW basis set has been successfully employed for systems containing localized electrons such as first row atoms and transition metals. In our implementation of the LAPW-MD scheme, iterative residual minimization algorithm is used to solve the electronic states problem. The atoms are moved according to forces derived from the Hellman–Feynman theorem and incomplete basis set correction terms. The performance of the program is further enhanced by parallelization. We will discuss technical details of the program implementation and present results obtained from this code to the equilibrium structures and vibrational properties of simple diatomic molecules.     

Density functional theory and ab initio studies of vibrational spectra of 2-bis (2-chloroethyl) aminoperhydro-1,3,2-oxazaphosphorinane-2-oxide

The Fourier transform Raman (FTR) and Fourier transform infrared (FTIR) spectra of 2-bis (2-chloroethyl) aminoperhydro-1,3,2-oxazaphosphorinane-2-oxide were recorded in the regions 4000–100 cm^? 1 and 4000–400 cm¹, respectively, in the solid phase. Molecular electronic energy, geometrical structure, harmonic vibrational spectra, infrared intensities and Raman scattering activities, highest occupied molecular orbital, lowest unoccupied molecular orbital energy, energy gaps and thermodynamical properties such as zero-point vibrational energies, rotational constants, entropies and dipole moment were computed at the Hartree–Fock/6-31G(d,p) and three parameter hybrid functional Lee–Yang–Parr/6-31G(d,p) levels of theory. The vibrational studies were interpreted in terms of potential energy distribution. The results were compared with experimental values with the help of scaling procedures. The observed wave number in FTIR and FTR spectra was analysed and assigned to different normal modes of the molecule. Most of the modes have wave numbers in the expected range and are in good agreement with computed values.     

Molecular dynamics study of the binding of phenylalanine stereoisomers to thermolysin

The stereospecificity in binding of phenylalanine as inhibitor in the active site of the thermolysin, has been investigated by means of molecular dynamics simulations using free energy integration techniques. The difference in the free energy of binding was found to be 2.0 ± 1.8 kJ/mol in favour of the D-form. This agrees with the experimental value, 2.8 kJ/mol. The result was obtained using a standard empirical force field (that of GROMOS). A different force field with 30% bigger charges (more like ab initio charges) was also tried. This resulted in less fluctuations and a more precise binding, but in a free energy. difference that was clearly larger than the experimental one. The phenylalanine backbone is located close to the zinc atom and the ring stays in the hydrophobic pocket in both the cases. The two stereoisomers differ mainly in the orientation of the backbone plane with respect to the active site and the rotational state of the dihedral around the C---CO^β bond.     

Distance geometry generates native-like folds for small helical proteins using the consensus distances of predicted protein structures.

For successful ab initio protein structure prediction, a method is needed to identify native-like structures from a set containing both native and non-native protein-like conformations. In this regard, the use of distance geometry has shown promise when accurate inter-residue distances are available. We describe a method by which distance geometry restraints are culled from sets of 500 protein-like conformations for four small helical proteins generated by the method of Simons et al. (1997). A consensus-based approach was applied in which every inter-Calpha distance was measured, and the most frequently occurring distances were used as input restraints for distance geometry. For each protein, a structure with lower coordinate root-mean-square (RMS) error than the mean of the original set was constructed; in three cases the topology of the fold resembled that of the native protein. When the fold sets were filtered for the best scoring conformations with respect to an all-atom knowledge-based scoring function, the remaining subset of 50 structures yielded restraints of higher accuracy. A second round of distance geometry using these restraints resulted in an average coordinate RMS error of 4.38 A.     

Molecular mechanic study of nerve agent O-ethyl S-[2-(diisopropylamino)ethyl]methylphosphonothioate (VX) bound to the active site of Torpedo californica acetylcholinesterase

Herein a molecular mechanic study of the interaction of a lethal chemical warfare agent, O-ethyl S-[2-(diisopropylamino)ethyl]methylphosphonothioate (also called VX), with Torpedo californica acetylcholinesterase (TcAChE) is discussed. This compound inhibits the enzyme by phosphonylating the active site serine. The chirality of the phosphorus atom induces an enantiomeric inhibitory effect resulting in an enhanced anticholinesterasic activity of the S_P isomer (VX^S) versus its R_P counterpart (VX^R). As formation of the enzyme-inhibitor Michaelis complex is known to be a crucial step in the inhibitory pathway, this complex was addressed by stochastic boundary molecular dynamics and quantum mechanical calculations. For this purpose two models of interaction were analyzed: in the first, the leaving group of VX was oriented toward the anionic subsite of TcAChE, in a similar way as it has been suggested for the natural substrate acetylcholine; in the second, it was oriented toward the gorge entrance, placing the active site serine in a suitable position for a backside attack on the phosphorus atom. This last model was consistent with experimental data related to the high inhibitory effect of this compound and the difference in activity observed for the two enantiomers. Proteins 28:543–555, 1997. © 1997 Wiley-Liss, Inc.     

Molecular dynamics - potential of mean force calculations as a tool for understanding ion permeation and selectivity in narrow channels

Ion channels catalyze the permeation of charged molecules across cell membranes and are essential for many vital physiological functions, including nerve and muscle activity. To understand better the mechanisms underlying ion conduction and valence selectivity of narrow ion channels, we have employed free energy techniques to calculate the potential of mean force (PMF) for ion movement through the prototypical gramicidin A channel. Employing modern all-atom molecular dynamics (MD) force fields with umbrella sampling methods that incorporate one hundred 1-2 ns trajectories, we find that it is possible to achieve semi-quantitative agreement with experimental binding and conductance measurements. We also examine the sensitivity of the MD-PMF results to the choice of MD force field and compare PMFs for potassium, calcium and chloride ions to explore the basis for the valence selectivity of this narrow and uncharged ion channel. A large central barrier is observed for both anions and divalent ions, consistent with lack of experimental conductance. Neither anion or divalent cation is seen to be stabilized inside the channel relative to the bulk electrolyte and each leads to large disruptions to the protein and membrane structure when held deep inside the channel. Weak binding of calcium ions outside the channel corresponds to a free energy well that is too shallow to demonstrate channel blocking. Our findings emphasize the success of the MD-PMF approach and the sensitivity of ion energetics to the choice of biomolecular force field.     

Kinetics and thermodynamics of the P870+Q−A → P870+Q−B reaction in isolated reaction centers from the photosynthetic bacterium Rhodopseudomonas sphaeroides

The temperature dependences of the P870⁺Q^?_A → P870Q_A and P870⁺Q^?_B → P870Q_B recombination reactions were measured in reaction centers from Rhodopseudomonas sphaeroides. The data indicate that the P870⁺Q^?_B state decays by thermal repopulation of the P870⁺Q^?_A state, followed by recombination. ΔG° for the P870⁺Q^?_A → P870⁺Q^?_B reaction is ?6.89 kJ · mol^?1, while ΔH° = ?14.45 kJ · mol^?1 and ?TΔS° = + 7.53 kJ · mol^?1. The activation ethalpy, $\text{H}{}^3$, for the P870⁺Q^?_A Δ P870⁺Q^?_B reaction is +56.9 kJ · mol^?1, while the activation entropy is near zero. The results permit an estimate of the shape of the potential energy curve for the P870⁺Q^?_A → P870⁺Q^?_B electron transfer reaction.     

Ab initio investigations on the geometric and electronic structures of a diblock molecular diode under the influence of an external bias

Theoretical investigations on the diblock molecular diode, thiophene–thiazole compound, have been carried out at the Hartree–Fock (HF) level by considering the interaction under the external bias. They demonstrate that the electronic structures of this kind of diode molecule are essentially different from those based on the Aviram and Ratner model, in which donor and acceptor π-conjugated segments are separated by an insulating σ-bonded segment, in terms of the energy levels of the frontier molecular orbitals as well as their spatial distributions. The introduction of the external bias modifies both the geometric and electronic structures. In particular, the spatial distributions of the frontier molecular orbitals are also shifted under the external bias. Moreover, all these features show a strong dependence on the polarity of the applied bias due to the build in intrinsic molecular asymmetric structures, which could be used to intuitively interpret the asymmetrical current–voltage behaviours of molecules.     

Electrostatic interactions in proteins: Calculations of the electrostatic term of free energy and the electrostatic potential field

The pH-dependence of the electrostatic energy of interactions between titratable groups is calculated for some well studied globular proteins: basic pancreatic trypsin inhibitor, sperm whale myoglobin and tuna cytochrome c. The calculations are carried out using a semi-empirical appraach in terms of the macroscopic model based on the Kirkwood-Tanford theory. The results are discussed in the light of their physicochemical and biological properties. It was found that the pH-dependence of the electrostatic energy correlates with the III–IV transition of cytochrome c. The electrostatic field of the cysteine proteinase inhibitor, cystatin, was calculated in two ways. In the first one, the electrostatic field created by the pH dependent charges of the ionizable groups and peptide dipoles was calculated using the approach proposed. In the second one, the finite-difference method was used. The results obtained by the two methods are in overall agreement. The calculated field was discussed in terms of the binding of cystatin to papain.     

DFT study of a substrate and inhibitors of 1-deoxy-2-xylulose-5-phosphate reductoisomerase – the potential novel target molecule for anti-malaria drug development

1-Deoxy-2-xylulose-5-phosphate (DOXP) reductoisomerase is a novel target for developing anti-malaria drugs. The determination of structural and electronic properties of the inhibitor molecules is of crucial importance for analyzing the interactions between DOXP-reductoisomerase and its inhibitors. Geometry-optimizations and single point calculations at the B3LYP/3-21G*//B3LYP/3-21G** and B3LYP/3-21G*//MP2/3-21G** levels were performed to determine the structures and charge distributions of an enzyme substrate (1-deoxy-D-xylulose 5-phosphate) and the two inhibitors (fosmidomycin and FR-900098). The theoretically derived bond lengths are in excellent agreement with the corresponding experimental values reported for similar structures. Partial charges and dipole moments are assigned using the Mulliken and natural population analyses. The calculated structures and partial charge distributions can readily be used for the further development of biologically active inhibitors of DOXP-reductoisomerase as well as parameters for docking experiments.Electronic Supplementary Material available.     

Effect of methylation on the stability and solvation free energy of amylose and cellulose fragments: a molecular dynamics study

Molecular dynamics (MD) simulations were used to study the stability and solvation of amylose and cellulose fragments. The recently developed gromos carbohydrate force field was further tested by simulating maltose, cellobiose, and maltoheptaose. The MD simulations reproduced fairly well the favorable conformations of disaccharides defined by the torsional angles related with the glycosidic bond and the radius gyration of maltoheptaose. The effects of methylation at different hydroxyl groups on the stability of amylose and cellulose fragments were investigated. The methylations of O-2 and O-3 reduce the stability of a single helix more than methylation at O-6, while the latter reduces the stability of a double helix more. Solvation free-energy differences between the unsubstituted amylose and cellulose fragments and the methylated species were studied using the single-step perturbation method. It was found that methylation at O-2 has the biggest effect, in agreement with experiment.     

Molecular dynamics study of the dissociation of an antigen-antibody complex in solution

Preliminary results are reported of a molecular dynamics calculation of free energy variations during the dissociation of an antigen-antibody complex, hen egg-white lysozyme — Fab D1.3, using atomic coordinates determined by the group of R. J. Poljak, and explicit handling of solvent molecules. After equilibration of the complex in solution at 300 K, a dissociation path was generated by a directed dynamics protocol. Then the thermodynamic perturbation method was used for computing the derivative of the free energy of the system with respect to dissociation coordinate, both for the undissociated complex and in two points along the path. 200-ps molecular dynamics simulations were carried out at each of these points. The results obtained are discussed, with special emphasis on the role of interstitial water in the appearance of a hydrophobic activation free energy.Université Paul Sabatier, and URA 505 of C.N.R.S.     

The importance of the protein in controlling the electrochemistry of heme metalloproteins: methods of calculation and analysis

 The importance of electrostatic effects in determining the free energy of redox reactions in proteins such as cytochromes and iron-sulfur complexes is well established. Several theoretical techniques have been used to analyze how the protein and its environment combine to produce the observed electrochemical midpoints. The free energy of changing the cofactor charge is influenced by the distribution of charges and dipoles in the protein, solvent and ions surrounding the protein, and by the redistribution of these charges and dipoles coupled to the reaction. An outline of a consistent view for calculating these effects will be presented and compared with other theoretical models. Heme redox potentials in yeast cytochrome c and the cytochrome subunit of photosynthetic reaction centers will be calculated to show how these protein structures produce the observed electrochemistry. Received, accepted: 26 November 1996     

Pyramidal inversion mechanism of simple chiral and achiral sulfoxides: a theoretical study

The pyramidal inversion mechanism of simple sulfoxides was studied, employing ab initio and DFT methods. The convergence of the geometrical and energetic parameters of H2SO and DMSO with respect to the Hamiltonian and basis set was analyzed in order to determine a computational level suitable for methyl phenyl sulfoxide (3), methyl 4-cyanophenyl sulfoxide (4), diphenyl sulfoxide (5), 4,4'-dicyanodiphenyl sulfoxide (6), benzyl methyl sulfoxide (7) and benzyl phenyl sulfoxide (8). The DFT B3LYP/6-311G(d,p) level was chosen for further calculations of larger sulfoxides. The barriers DeltaE calculated for the pyramidal inversion mechanism of sulfoxides 3-8 are in the range of 38.7-47.1 kcal/mol. These values are in good agreement with the experimental barriers for racemization via the pyramidal inversion mechanism. A resonance effect of a phenyl ring selectively stabilizes the transition state conformations, decreasing the energy barrier for pyramidal inversion by about 3 kcal/mol, compared to a similar molecule without a phenyl substituent. Introducing electron withdrawing groups (cyano) at the para positions of the phenyl ring(s) causes a further decrease of the energy barrier.     

Complete basis set,hybrid-DFT study and NBO interpretation of conformational analysis of 2-methoxytetrahydropyran and its thiopyran and selenopyran analogues in relation to the anomeric effect

2-Methoxytetrahydropyran (1), -thiopyran (2) and -selenopyran (3) have been chosen as model compounds to investigate the origin of the anomeric effect (AE). The impacts of the hyperconjugation, electrostatic and steric interactions on the conformational preferences of compounds 1–3 have been analysed by means of complete basis set-4, hybrid-density functional theory (B3LYP/6-311+G**) based methods and natural bond orbital (NBO) interpretation. Both levels of theory showed that the axial conformations of compounds 1–3 are more stable than their equatorial conformations. The Gibbs free energy difference (G _eq–G _ax) values (i.e. ΔG _eq–ax) between the axial and equatorial conformations increase from compound 1 to compound 2 but decrease from compound 2 to compound 3. Based on the NBO results obtained, the AE associated with the electron delocalisation [i.e. Σ(endo-AE_eq + exo-AE_eq) ? Σ(endo-AE_ax + exo-AE_ax)] increase slightly from compound 1 to compound 2 but decrease from compound 2 to compound 3. Similar trend is also observed for the differences between the calculated total steric exchange energy values [i.e. Δ(TSEE)_eq–ax]. On the other hand, the calculated differences between the dipole moment values of the axial and equatorial conformations [i.e. Δ(μ_eq–μ_ax)] decrease from compound 1 to compound 3. These findings led to the proposal that the AE associated with the electron delocalisation (the hyperconjugation effect) is more significant for the explanation of the conformational preferences of compounds 1–3 than the electrostatic model. The correlations between the AE associated with the electron delocalisation, bond orders, TSEE, ΔG _eq–ax, dipole–dipole interactions, structural parameters and conformational behaviours of compounds 1–3 have been investigated.     

Molecular modeling and in-silico engineering of Cardamom mosaic virus coat protein for the presentation of immunogenic epitopes of Leptospira LipL32

Expression of Cardamom mosaic virus (CdMV) coat protein (CP) in E. coli forms virus-like particles. In this study, the structure of CdMV CP was predicted and used as a platform to display epitopes of the most abundant surface-associated protein, LipL32 of Leptospira at C, N, and both the termini of CdMV CP. In silico, we have mapped sequential and conformational B-cell epitopes from the crystal structure of LipL32 of Leptospira interrogans serovar Copenhageni str. Fiocruz L1-130 using IEDB Elipro, ABCpred, BCPRED, and VaxiJen servers. Our results show that the epitopes displayed at the N-terminus of CdMV CP are promising vaccine candidates as compared to those displayed at the C-terminus or at both the termini. LipL32 epitopes, EP2, EP3, EP4, and EP6 are found to be promising B-cell epitopes for vaccine development. Based on the type of amino acids, length, surface accessibility, and docking energy with CdMV CP model, the order of antigenicity of the LipL32 epitopes was found to be EP4 > EP3 > EP2 > EP6.     

Theoretical modeling of the O-intermediate structure of bacteriorhodopsin

Bacteriorhodopsin is a prototype of efficient molecular machinery functioning as a light-activated proton pump. Among the five distinct intermediates (K, L, M, N, and O) of the photocycle, there is less structural information on the later stages compared with the early intermediates. Here, we report the structural modeling of the O-intermediate for which the determination of experimental structure remains difficult. Hypothetical conformational change of the molecule from the light-adapted state to the O-intermediate state was simulated by gradually changing the protonation state of two residues. To achieve accurate molecular modeling, we carefully constructed a realistic system of the native purple membrane. The modeled structure of the O-intermediate has some implications about proton transfer in the later stages of the photocycle and the structural response of bacteriorhodopsin to the inner charge distribution.