Combined theoretical studies on solvation and hydrogen bond interactions in glycine tripeptide

The theoretical conformational analysis of glycine tripeptide (GT) has been carried out by molecular dynamics (MD) method in order to find minimum energy conformations. The MD studies on GT with water have been carried out for over 10 ns with a time step of 2 fs using fixed charge force field (AMBER ff03). By adding the solvation effect using water as a solvent, the GT conformers identified in this study exhibit α-helical conformation. Compared with the earlier reports, this MD study is able to identify the energetically favourable GT conformations. The obtained geometry of the five most stable GT conformations was optimised using the density functional theory method at B3LYP/6-311G** level of theory. Subsequently, the effects of solvation on the conformational characteristics of five most stable GT conformers with four water molecules (the number of water molecules in the first solvation shell of GT obtained from MD study) were investigated using the same method and the same level of theory. The effect of microsolvation on the fifth GT conformer has been studied with a cluster of 11 water molecules as the first hydration shell which generates folded structure. The interaction energies of all the complexes are calculated by correcting the basis set superposition error. The strong hydrogen bond mainly contributes to the interaction energies. The atoms in molecules theory and natural bond orbital analysis were used to study the origin of H-bonds. A good correlation between the structural parameters and the properties of charge density is found. NMR calculations show that the C = O carbons of the amine groups of the first and middle glycine fragments have maximum chemical shifts.     

Molecular dynamics simulations and density functional theory studies of NALMA and NAGMA dipeptides

Classical molecular dynamics (MD) simulations using fixed charged force field (AMBER ff03) and density functional theory method using the M05-2X/6-31G^?? level of theory have been used to investigate the plasticity of the hydrogen bond formed between dipeptides of N-Acetyl-Leucine-MethylAmide (NALMA), N-Acetyl-Glycine-MethylAmide (NAGMA), and vicinity of water molecules at temperature of 300?K. We have noticed that 2–3 water molecules contribute to change in the conformations of dipeptides NAGMA and NALMA. The self-assembly of 11 water molecules leads to the formation of water bridge at vicinity of the dipeptides and it constrain the conformations of dipeptides. We have found that the energy balance between breaking of the C?=?O…H–N H bonds and the formation of the C?=?O…H–O (wat) H bonds may be one of the determining factors to control the dynamics of the folding process of protein molecules.     

Free energy surface for rotamers of cis-enol malonaldehyde in aqueous solution studied by molecular dynamics calculations

Two-dimensional free energy surfaces for four rotamers of cis-enol malonaldehyde in water have been investigated by umbrella sampling molecular dynamics (MD) calculations. Biasing potential used in the umbrella sampling calculation was adopted to be the minus of conformational free energy preliminary obtained by the thermodynamic integration MD calculations for the rigid malonaldehyde whose stretching and bending were all fixed. The calculated free energy surface shows that, in water, a rotamer that has an intramolecular hydrogen bond is most stable among the rotamers. This is the same as that in vacuum, while order of relative stability of the other three rotamers is different in water and in vacuum. Inclusion of intramolecular vibrations changed the free energy surface little, i.e. at most 2.6 kJ/mol, which is much smaller than the solvation free energy. Free energy barriers from the most stable intramolecular hydrogen bonded rotamer to the others are lowered by hydration but they are still very high, >50 kJ/mol, such that the malonaldehyde molecule spends most of its time in water taking this conformation. Thus, reaction coordinate for intramolecular proton transfer reaction in water may be constructed assuming this rotamer.     

Theoretical calculation on the reaction of alkene molecules on H-terminated Si(100)-3×1 surface

A study on mechanisms of radical initiated surface chain reaction of ethylene molecule on H-terminated Si(100)-3 × 1 has been carried out in a supercell approach by using density functional theory and ab initio molecular dynamic method. On the H-terminated Si(100)-3 × 1 surface, one of the crucial steps of the surface chain reaction, namely, the abstraction of a H atom from a nearby surface hydride unit, is found to have a somewhat smaller activation energy from the nearest silicon site than from the next-nearest silicon site. From the intermediate state to the final state, the transition state has bigger activation energy. Ab initio molecular dynamics (MD) shows that the H-abstraction on Si(100)-3 × 1 surface bound organic group with a carbon-centered radical is very easy to be obtained from the transition state, and it also shows that the C…H bond at methyl group is formed in a very short MD time, and the Si…C bond between the Si surface and the alkyl chain oscillates with time evolution on Si(100)-3 × 1 surface.     

A molecular simulation of the compatibility of chitosan and poly(vinyl pyrrolidone)

The compatibility of chitosan (CS) and poly(vinyl pyrrolidone) was investigated by molecular dynamic (MD) simulations using the Flory–Huggins theory. The specific interactions in blends were studied by the radial distribution function (RDF). The Flory–Huggins interaction parameter, χ, was calculated at 298 K to assess the blend compatibility at different component ratios in the polymers. Miscibility was observed for blends with more than 50% of CS in the molar fraction, while immiscibility was prevalent at the molar fraction of CS between 10 and 50% of CS. Miscibility between poly(N-vinyl-2-pyrrolidone) (PVP) and CS polymers is attributed to the hydrogen bond formation of the –C = O group of PVP and the –CH₂OH groups of CS. This was further confirmed by MD simulations of RDFs for groups or atoms that are involved in interactions. These results are correlated well to obtain more realistic information on interactions involved as a function of blend composition.     

The solvation of bromide anion in acetonitrile: a structural study based on the combination of theoretical calculations and X-ray absorption spectroscopy

This work studies the solvation of bromide in acetonitrile by combining quantum mechanics, computer simulations and X-ray absorption near edge structure (XANES) spectroscopy. Three different sets of interaction potentials were tested, one of them derived from literature and the other two are simple modifications of the previous one to include specificities of the bromide–acetonitrile interactions. Results for microsolvation of bromide were obtained by quantum mechanical optimization and classical minimization of small clusters [Br(ACN) _n ]^?  (n = 9, 20). Analysis of molecular dynamics (MD) simulations has provided structural, dynamic and energetic aspects of the solvation phenomenon. The theoretical computation of Br K-edge XANES spectrum in solution using the structural information obtained from the different simulations has allowed the comparison among the three different potentials, as well as the examination of the main structural and dynamic factors determining the shape of the experimental spectrum.     

Investigation of rubidium(I) ion solvation in liquid ammonia using QMCF-MD simulation and NBO analysis of first solvation shell structure

Rb(I) ion solvation in liquid ammonia has been studied by an ab initio quantum mechanical charge field molecular dynamics simulation, and the first solvation shell structure has been analyzed using natural bond orbital. The simulation was performed for an ion and 593 ammonia molecules in a box with a length of 29.03 Å corresponding to a liquid ammonia density of 0.69 g/mL at 235.16 K. The quantum mechanical calculation was carried out for atomic interactions in the radius of 6.4 Å from the ion using LANL2DZ ECP and DZP (Dunning) basis sets for Rb(I) ion and ammonia respectively. The trajectories of the simulation were analyzed in terms of radial, angular, and coordination number distribution functions, vibration, and mean residence time (MRT). Two solvation shell regions are observed for the Rb(I)-N as well as the Rb(I)-H. The maximum distance of Rb(I)-N in the first solvation shell is in accordance with experimental data where a coordination number of 8 is favorable. A non-single coordination number of the first and second shell indicates dynamic solvation structure. It is confirmed by frequent exchange ligand processes observed within a simulation time of 15 ps. The low stabilization energy of donor acceptor ion-ligand interaction with a small Wiberg bond index affirms that the Rb(I)-NH₃ interaction is weak electrostatically.     

Exploring the structure and dynamics of nano-confined water molecules using molecular dynamics simulations

Confinement effects can lead to drastic changes in the structural and dynamical properties of water molecules. In this work, we have performed classical molecular dynamics simulations of endohedral fullerenes of type (H₂O)_n@C_m (n = 1, 12, 21, 62, 108 and m = 60, 180, 240, 500 and 720) to explore the effects of spherical confinement on water properties. It is shown that these confined water molecules can form distinct solvation pattern depending upon the available space inside the fullerene cavity. For the systems with smaller diameter, cage-like structure is predominant whereas bulk-like structure is observed for larger fullerenes. The orientational relaxation of these confined water molecules showed slower relaxation as the cavity diameter increases except for the (H₂O)₂₁@C₂₄₀. In this case, stable cage-like structure hinders the overall dynamics of the trapped water molecules. Finally, we have calculated the hydrogen bond lifetimes from the hydrogen bond time correlation functions and compared with that of bulk water.     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Solvation free energies of nucleic acid bases in ionic liquids

The solvation free energies of five nucleic acid bases in [C_nbim]Br (where n = 2, 4, 6) ionic liquids (ILs) were computed using the Bennett acceptance ratio (BAR) method employing molecular dynamics simulations. The computed free energies using BAR were in agreement with other methods. The large and negative predicted free energies of the bases in ILs indicated that the bases were better solvated in the ILs rather than in water. Hydrogen bonding interactions between polar sites of the bases and ILs’ ions significantly contributed to the solvation mechanism.     

Structure,stability and water permeation of ([D-Leu-L-Lys-(D-Gln-L-Ala)3]) cyclic peptide nanotube: a molecular dynamics study

The structural stability of 8 × ([D-Leu-L-Lys-(D-Gln-L-Ala)₃]) cyclic peptide nanotube (CPN) in water and different phospholipid bilayers were explored by 100 ns independent molecular dynamics (MD) simulations. The role of non-bonded interaction energy between the side and main chains of cyclic peptide rings in different membrane environments assessed, wherein the repulsive electrostatic interaction energy between neighbouring cyclic peptide rings was found adequate to break hydrogen bond energy thereby to crumple CPN. Further, the water permeation across the CPN channel was studied in four types of phospholipid bilayers- DMPG (1,2-Dimyristoyl-sn-glycero-3-phosphorylglycerol), DMPS (1,2-Dimyristoyl-sn-glycero-3-phosphoserine), POPC (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and POPE (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) from MD simulations. DMPS membrane shows higher non-bonded interaction energies (?1913.06 kJ/mol of electrostatic interaction energy and ?994.13 kJ/mol of van der Waals interaction energy) with CPN due to the presence of polar molecules in lipid structure. Thusly, the non-bonded interaction energies were essential towards the stability of CPN than hydrogen bonds between the nearby cyclic peptides. The result also reveals the role of side chains, hydrogen bonds and non-bonded interaction energies in an aqueous environment. The diffusion coefficient of water obtained from means square deviation calculation shows similar coefficients irrespective of the lipid surroundings. However, the permeation coefficients demonstrate water flow in the channel relies upon the environment.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Peptides and Peptoids - A Systematic Structure Comparison

A systematic analysis of the conformational space of the basic structure unit of peptoids in comparison to the corresponding peptide unit was performed based on ab initio MO theory and complemented by molecular mechanics (MM) and molecular dynamics (MD) calculations both in the gas phase and in aqueous solution.The calculations show three minimum conformations denoted as C_7ß, a_D and a that do not correspond to conformers on the gas phase peptide potential energy hypersurface. The influence of aqueous solvation was estimated by means of continuum models. The MD simulations indicate the a_D form as the preferred conformation in solution both in cis and trans peptide bond orientations.     

Extending the host range of the bioherbicidal fungus Colletotrichum gloeosporioides f. sp. aeschynomene

Spore formulations of the bioherbicidal fungus Colletotrichum gloeosporioides f. sp. aeschynomene (ATCC No. 20358) (CGA) were evaluated for control of three weed species: northern jointvetch (Aeschynomene virginica), Indian jointvetch (A. indica), and hemp sesbania (Sesbania exaltata) in greenhouse experiments. Mortality, dry weight reduction, and plant height reduction of A. virginica seedlings ranged from 98% to 100%, 15 days after inoculation with CGA in water, in an invert emulsion, or in Silwet L-77 surfactant. However, CGA in water caused no effects of these parameters on S. exaltata, and only slight effects on A. indica. A. indica and S. exaltata were also severely injured (mortality, dry weight and plant height reduction, 98–100%) by CGA in the invert emulsion or in Silwet L-77. The CGA in Silwet formulation incited severe disease development more rapidly than the invert emulsion or water formulations of CGA in all species. These results suggest that the host range of CGA can be expanded though formulation modification to enable this bioherbicide to control multiple weeds, thus making this product more economically acceptable.     

Structural analysis,spectroscopic characterization and molecular docking studies of the cyclic heptapeptide

The theoretically possible stable conformer of the cyclic heptapeptide, that has significant anti-metastatic activity, was examined by conformational analysis followed by DFT calculations. Experimental infrared and Raman spectroscopy, together with theoretical DFT (6-31G (d,p) basis set)-based quantum chemical calculations, have been used to understand the structural and spectral characteristics of cyclo(Gly-Arg-Gly-Asp-Ser-Pro-Ala) {cyclo(GRGDSPA)}. A complete analysis of the vibrational spectrum has been reported on the basis of potential energy distribution (PED%) data of the vibrational modes. Finally, the calculation results were applied to simulate infrared and Raman spectra of the title compound. The simulated spectra satisfactorily coincide with the experimental spectra. In addition, molecular electrostatic potential and frontier molecular orbital analysis were investigated using theoretical calculations. The stability of the molecule, arising from hyperconjugative interaction and charge delocalization, has been analyzed using natural bond orbital analysis and a high E(2) value reveals the presence of strong interaction between donors and acceptors. Molecular docking studies with fibronectin were performed on cyclo(GRGDSPA) in order to understand its inhibitory nature. The results indicate that the docked ligand {cyclo(GRGDSPA)} forms a stable complex with human fibronectin and gives a binding affinity value of ?7.7 kcal/mol, which points out that cyclo(GRGDSPA) might exhibit inhibitory activity against the attachment of melanoma cells to human fibronectin.     

Structural and dynamical studies of all-trans and all-cis cyclo[(1R,3S)-γ-Acc-Gly]<Subscript>3</Subscript> peptides

The quantum chemical and molecular dynamics studies have been performed to infer the structural changes of all-trans and all-cis forms of cyclo[(1R,3S)-3-aminocyclohexanecarboxylicacid(γ-Acc)-α-Glycine(Gly)]₃ hexapeptide. The backbone conformations of the above peptide have been analyzed using the valence and peptide deformation angles applying B3LYP/6–311G** level of theory. The conformational preference of the backbone of all-trans and all-cis cyclo[(1R,3S)-γ-Acc-Gly]₃ hexapeptides is found to depend on the puckering of cyclohexane rings. The non-uniform distribution of water inside the cavity is observed, where sometimes water molecules formed a chain like conformation through hydrogen bond networks while traversing the pore of all-cis cyclo[(1R,3S)-γ-Acc-Gly]₃ peptide. Larger relaxation times of the order of a hundred to two hundred pico seconds for active site…water hydrogen bond interactions were noticed. The hydrophobic nature of the cavity of all-trans cyclo[(1R,3S)-γ-Acc-Gly]₃ due to the presence of cyclohexane moiety has been analyzed. Further this investigation emphasized on the non-transport of molecules through the pore of all-trans cyclo[(1R,3S)-γ-Acc-Gly]₃ peptide due to the obstruction produced by cyclohexane groups.     

Drug targeted virtual screening and molecular dynamics of LipU protein of Mycobacterium tuberculosis and Mycobacterium leprae

The lipolytic protein LipU was conserved in mycobacterium sp. including M. tuberculosis (MTB LipU) and M. leprae (MLP LipU). The MTB LipU was identified in extracellular fraction and was reported to be essential for the survival of mycobacterium. Therefore to address the problem of drug resistance in pathogen, LipU was selected as a drug target and the viability of finding out some FDA approved drugs as LipU inhibitors in both the cases was explored. Three-dimensional (3D) model structures of MTB LipU and MLP LipU were generated and stabilized through molecular dynamics (MD). FDA approved drugs were screened against these proteins. The result showed that the top-scoring compounds for MTB LipU were Diosmin, Acarbose and Ouabain with the Glide XP score of ?12.8, ?11.9 and ?11.7 kcal/mol, respectively, whereas for MLP LipU protein, Digoxin (?9.2 kcal/mol), Indinavir (?8.2 kcal/mol) and Travoprost (?8.2 kcal/mol) showed highest affinity. These drugs remained bound in the active site pocket of MTB LipU and MLP LipU structure and interaction grew stronger after dynamics. RMSD, RMSF and Rg were found to be persistent throughout the simulation period. Hydrogen bonds along with large number of hydrophobic interactions stabilized the complex structures. Binding free energies obtained through Prime/MM-GBSA were found in the significant range from ?63.85 kcal/mol to ?34.57 kcal/mol for MTB LipU and ?71.33 kcal/mol to ?23.91 kcal/mol for MLP LipU. The report suggested high probability of these drugs to demolish the LipU activity and could be probable drug candidates to combat TB and leprosy disease.     

Structures and stability of salt-bridge in aqueous solution

Structures and stability of salt-bridges in aqueous solutions were investigated using a complex formed from the guanidinium (Gdm+) and formate (FmO-) ions as a model system. The Test-particle model (T-model) potentials to describe the interactions in the Gdm+-H2O, FmO(-)-H2O and Gdm+-FmO- complexes were constructed, tested and applied in molecular dynamics (MD) simulations of the aqueous solutions at 298 K. The three-dimensional structures and energetic of the hydrogen bond (H-bond) networks of water in the first hydration shells of the Gdm+ and FmO- ions, as well as the Gdm+-FmO- complex, were visualized and analyzed using various probability distribution (PD) maps. The structures of the average potential energy landscapes at the H-bond networks were employed to characterize the stability and dynamic behavior of water molecules in the first hydration shells of the solutes. It was observed that water molecules in the first hydration shell of the close-contact Gdm+-FmO- complex form associated H-bond networks, which introduce a net stabilization effect to the ion-pair, whereas those in the interstitial H-bond network destabilize and break the solvent-separated Gdm+-FmO- complex. The present results showed that, in order to provide complete insights into the structures and stability of ion-pairs in aqueous solutions, explicit water molecules have to be included in the model calculations.     

In-silico analysis of Sirt2 from Schistosoma mansoni: structures,conformations, and interactions with inhibitors

Sirtuins are NAD+-dependent lysine deacetylases member of the class III HDAC family. These are demonstrated to be therapeutic targets in parasitic diseases like schistosomiasis. Observations suggested that sirtuin enzyme is necessary for the functionality of fe/male reproductive system, due to which SmSirt2 is treated as a potential therapeutic target. There are no structural and molecular features of SmSirt2 have been reported yet. In this study, homology modeling has been used to determine the three-dimensional features of the SmSITRT2. Further, structure validation has been performed by energy minimization and Ramachandran plot. Validated structures are further subjected to molecular docking and virtual screening to find the best lead molecules for downstream analysis. Ten lead molecules were selected while comparing virtual screening of hSirt2 and SmSirt2 both. These leads are further compared with AKG2 which is known inhibitor of hSirt2 (?8.8 kcal/mol). Out of selected 10 leads, four of them (ZINC23995485 (?9.5 kcal/mol), ZINC53298162 (?9.4 kcal/mol), ZINC70927268 (?10.0 kcal/mol), ZINC89878705 (?11.2 kcal/mol)) have shown better interaction with SmSirt2, in which ZINC89878705 (?11.2 kcal/mol) shows a more compact packing as compared to AKG2 and rest of ligands. These molecules could be further subject to in vitro study and model of SmSirt2 has been proposed for further structure-based drug design projects concerning sirtuins from Schistosoma mansoni.     

Insight into the nature of the interactions of furan and thiophene with hydrogen halides and lithium halides: ab initio and QTAIM studies

The nature of the interactions of furan and thiophene with hydrogen halides and lithium halides has been investigated using ab initio calculations and QTAIM analysis. The concept of molecule formation density difference (MFDD) is introduced to study weak hydrogen bond (HB) and lithium bond (LB) interactions. The results have shown the molecular electrostatic potentials of furan and thiophene, as well as of the hydrogen halides and lithium halides, determine the geometries of the complexes. Both the studied HB and LB interactions can be classified as “closed-shell” weak interactions. The topological properties and energy properties at the bond critical points of HB and LB have been shown to be exponentially dependent on intermolecular distances d(H-bond) and d(Li-bond), which enables interpretation of the strength of the HB and LB interactions in terms of these ρ(r) properties. Electron transfer plays a more important role in the formation of HB than in that of LB, while electrostatic interaction in LB is more dominant than that in HB.