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.     

Quantum chemical study of the structure and properties of citrinin

Detailed structures and electronic properties of three tautomeric forms of the toxin citrinin were investigated using several quantum calculation methods. Energetic preference of the predominant p- and o-quinone methide tautomeric forms is dependent on the method of calculation. A previously unstudied carboxylic acid enol tautomer was calculated to be surprisingly stable in vacuo, being within 2.5 kcal mol^? 1 at the B3LYP/6-311++G(2d,2p) level of theory. Despite differences in bond nature and connectivity of tautomers, the natural bond orbital analysis revealed that tautomeric forms share similar natural charges and natural electron configurations. Calculated bond lengths corresponded with experimentally observed values and assignments for the calculated infrared vibrational frequencies are reported.     

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.     

Structure and hydrogen bonds of cyclohexapeptide RA-VII by molecular dynamics simulations and quantum chemical calculations

We computationally examined the structure of anti-tumour bicyclic hexapeptide RA-VII. This peptide adopts three conformations (confs.), A, B and C, in dimethyl sulfoxide (DMSO). Although it was experimentally reported that the structure of conf. A is important for anti-tumour activity, the dynamics of confs. A, B and C are not well known. We performed quantum chemical calculations and molecular dynamics (MD) simulations of RA-VII in DMSO. The MD simulations indicated two different local stable structures for conf. C: a structure containing a bent 18-membered ring and another structure containing a rotated peptide bond between Tyr⁶ and d-Ala¹. The root-mean-square fluctuation of the 14-membered ring for conf. A was larger than that for confs. B and C. Ala⁴ formed intramolecular hydrogen bonds more often in conf. A than in the other conformations. A large number of hydrogen bonds and large structural fluctuations are important for the anti-tumour activity of RA-VII. Our results for the structural change of conf. C and the analysis of the dynamics for confs. A, B and C may contribute to the design of new analogues of cyclic peptides.     

Molecular dynamics studies of a DNA-binding protein: 2. An evaluation of implicit and explicit solvent models for the molecular dynamics simulation of the Escherichia coli trp repressor.

Although aqueous simulations with periodic boundary conditions more accurately describe protein dynamics than in vacuo simulations, these are computationally intensive for most proteins. Trp repressor dynamic simulations with a small water shell surrounding the starting model yield protein trajectories that are markedly improved over gas phase, yet computationally efficient. Explicit water in molecular dynamics simulations maintains surface exposure of protein hydrophilic atoms and burial of hydrophobic atoms by opposing the otherwise asymmetric protein-protein forces. This properly orients protein surface side chains, reduces protein fluctuations, and lowers the overall root mean square deviation from the crystal structure. For simulations with crystallographic waters only, a linear or sigmoidal distance-dependent dielectric yields a much better trajectory than does a constant dielectric model. As more water is added to the starting model, the differences between using distance-dependent and constant dielectric models becomes smaller, although the linear distance-dependent dielectric yields an average structure closer to the crystal structure than does a constant dielectric model. Multiplicative constants greater than one, for the linear distance-dependent dielectric simulations, produced trajectories that are progressively worse in describing trp repressor dynamics. Simulations of bovine pancreatic trypsin were used to ensure that the trp repressor results were not protein dependent and to explore the effect of the nonbonded cutoff on the distance-dependent and constant dielectric simulation models. The nonbonded cutoff markedly affected the constant but not distance-dependent dielectric bovine pancreatic trypsin inhibitor simulations. As with trp repressor, the distance-dependent dielectric model with a shell of water surrounding the protein produced a trajectory in better agreement with the crystal structure than a constant dielectric model, and the physical properties of the trajectory average structure, both with and without a nonbonded cutoff, were comparable.     

Vibrational (FT-IR and FT-Raman) spectral investigations of 7-aminoflavone with density functional theoretical simulations

FT-Raman and FT-IR spectra of the 7-aminoflavone have been recorded and analysed. The detailed interpretation of the vibrational spectra has been carried out with the aid of normal coordinate analysis following the scaled quantum mechanical force field methodology. The various intramolecular interactions that are responsible for stabilisation of the molecule were revealed by natural bond orbital analysis. The obtained vibrational wavenumbers and optimised geometric parameters were observed to be in good agreement with the experimental data. The carbonyl stretching vibrations have been lowered due to conjugation and hydrogen bonding in the molecules.     

Characterization of the active site of DNA polymerase beta by molecular dynamics and quantum chemical calculation

It is well established that the fully formed polymerase active site of the DNA repair enzyme, polymerase beta (pol beta), including two bound Mg2+ cations and the nucleoside triphosphate (dNTP) substrate, exists at only one point in the catalytic cycle just prior to the chemical nucleotidyl transfer step. The structure of the active conformation has been the subject of much interest as it relates to the mechanism of the chemical step and also to the question of fidelity assurance. Although crystal structures of ternary pol beta-(primer-template) DNA-dNTP complexes have provided the main structural features of the active site, they are necessarily incomplete due to intentional alterations (e.g., removal of the 3'OH groups from primer and substrate) needed to obtain a structure from midcycle. Working from the crystal structure closest to the fully formed active site [Protein Data Bank (PDB) code: 1bpy], two molecular dynamics (MD) simulations of the solvated ternary complex were performed: one with the missing 3'OHs restored, via modeling, to the primer and substrate, and the other without restoration of the 3'OHs. The results of the simulations, together with ab initio optimizations on simplified active-site models, indicate that the missing primer 3'OH in the crystal structure is responsible for a significant perturbation in the coordination sphere of the catalytic cation and allow us to suggest several corrections and additions to the active-site structure as observed by crystallography. In addition, the calculations help to resolve questions raised regarding the protonation states of coordinating ligands.     

Molecular dynamics simulations of trehalose as a 'dynamic reducer' for solvent water molecules in the hydration shell

Systematic computational work for a series of 13 disaccharides was performed to provide an atomic-level insight of unique biochemical role of the alpha,alpha-(1-->1)-linked glucopyranoside dimer over the other glycosidically linked sugars. Superior osmotic and cryoprotective abilities of trehalose were explained on the basis of conformational and hydration characteristics of the trehalose molecule. Analyses of the hydration number and radial distribution function of solvent water molecules showed that there was very little hydration adjacent to the glycosidic oxygen of trehalose and that the dynamic conformation of trehalose was less flexible than any of the other sugars due to this anisotropic hydration. The remarkable conformational rigidity that allowed trehalose to act as a sugar template was required for stable interactions with hydrogen-bonded water molecules. Trehalose made an average of 2.8 long-lived hydrogen bonds per each MD step, which was much larger than the average of 2.1 for the other sugars. The stable hydrogen-bond network is derived from the formation of long-lived water bridges at the expense of decreasing the dynamics of the water molecules. Evidence for this dynamic reduction of water by trehalose was also established based on each of the lowest translational diffusion coefficients and the lowest intermolecular coulombic energy of the water molecules around trehalose. Overall results indicate that trehalose functions as a 'dynamic reducer' for solvent water molecules based on its anisotropic hydration and conformational rigidity, suggesting that macroscopic solvent properties could be modulated by changes in the type of glycosidic linkages in sugar molecules.     

Molecular dynamics studies in nanoscale liquid structures: geometry and thermal effects on nanojet development

Conventional macroscopic jet theory relies heavily on experimental correlations which cannot be easily extended to the nanoscale regime. Moreover, the fluid dynamic effects at small length scales and their contribution to the development of nanoscale liquid structures are fundamentally different from their macroscopic counterparts. This coupled with the high spatial and temporal resolution requirements at nanoscale domains make molecular dynamics (MD) an excellent tool for studying such structures. In this study, the formation and breakup of nanojets (NJs) developing from high pressure into vacuum is investigated using MD based on non-Hamiltonian formulations. By ejecting the equilibrated argon atoms through various nozzle geometries and diameters, nanoscale jet flows were generated. The dependence of the jet structure on nozzle geometry and diameter is studied. The influence of geometry on NJ formation is also studied along with issues involved in the equilibration and thermostat coupling parameter. Various thermostats are compared to understand the role they play in MD simulations of liquid nanostructures. Tuning of the thermostat coupling parameter has also been discussed. The jet breakup phenomenon is analysed and a comparative study, vis-à-vis, well-established continuum and stochastic models, is attempted.     

Identification of InhA inhibitors: A combination of virtual screening,molecular dynamics simulations and quantum chemical studies

In the present work, multiple pharmacophore-based virtual screening of the SPECS natural product database was carried out to identify novel inhibitors of the validated biological target, InhA. The pharmacophore models were built from the five different groups of the co-crystallized ligands present within the active site. The generated models with the same features from each group were pooled and subjected to the test set validation, receiver–operator characteristic analysis and Güner–Henry studies. A set of five hypotheses with sensitivity > 0.5, specificity > 0.5, area under curve (AUC) > 0.7, and goodness of hit score > 0.7 were retrieved and exploited for the virtual screening. The common hits (87 molecules) obtained from these hypotheses were processed via drug-likeness filters. The filtered molecules (27 molecules) were compared for the binding modes and the top scored molecules (12 molecules) along with the reference (triclosan (TCL), docking score = ?11.65 kcal/mol) were rescored and reprioritized via molecular mechanics-generalized Born surface area approach. Eventually, the stability of reprioritized (10 molecules) docked complexes was scrutinized via molecular dynamics simulations. Moreover, the quantum chemical studies of the dynamically stable compounds (9 molecules) were performed to understand structural features essential for the activity. Overall, the protocol resulted in the recognition of nine lead compounds that can be targeted against InhA.     

An overview about the impact of hinge region towards the anticancer binding affinity of the Ck2 ligands: a quantum chemical analysis

Abstract

Casein kinase 2 (CK2) is extremely preserved and universally uttered serine/threonine kinase, vital for cellular feasibility. The present study aimed to analyse the binding strength of CK2 ligands specifically in the hinge region, as it is aware that most of the existing drugs are targeted to bind the hinge of the corresponding protein. The analysis will give a clear picture about the role of hinge region with ligand, which will be useful for scientist community in drug designing. To predict the binding strength of CK2 ligands, the role of halogen bond, hydrogen bond interaction at the hinge region was depicted in detail through interaction energy calculations at M062Z/def2-QZVP level of theory. Highest occupied molecular orbital (HOMO) map plotted for CK2 ligands gives a clear pictorial representation of orbitals, which induce for interaction. Ligand properties discussed in detail through Lipinski’s five rules predict that almost all the ligands satisfy the rule, except 3KXG, which violates Lipinski’s two rules, i.e. molecular mass exceeds 500?Da, i.e. 512.61?Da, and Log P value is high of 5.09. The natural bond orbital analysis deliberates that the hydrogen/halogen bonds figuring out within the complexes are observed to have moderate stabilization energy, but those hydrogen/halogen bonds that exist with close contacts have high stabilization energy. Overall, this computational work will give an understandable depiction for modelling anticancer ligands along the hinge region in CK2 protein; also, it will give a new path for the choice of side chains on the ligand.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics studies on structure and dynamics of phospholamban monomer and pentamer in membranes

Phospholamban (PLB) is an integral membrane protein of 52 residues that regulates the activity of the sarcoplasmic reticulum calcium pump in cardiac muscle cells through reversible phosphorylation of Ser16. To explore its possible conformations and dynamics in a monomeric state, we have performed comparative molecular dynamics simulations of unphosphorylated and phosphorylated PLB (pPLB) with various orientations in POPC membranes. The simulations indicate that dynamics of the cytoplasmic domain is highly dependent on its interactions with membranes, that is, large conformational changes in the absence of membrane interactions, but very restricted dynamics in their presence. pPLB shows more structural flexibility in its cytoplasmic domain, which is consistent with experimental observations. We have also performed a simulation of a PLB pentameric structure (the so‐called bellflower model), recently determined in micelles, to investigate its behaviors in a POPC membrane. The cytoplasmic domain in each monomer shows uncorrelated dynamics and undergoes large conformational changes toward the membrane surface during the simulation, which supports the so‐called pinwheel model of the PLB pentamer structure. The hydrophobic nature of the pentameric pore excludes water molecules in the pore region, which illustrates that the pore appears to be an energetic barrier for ion and water translocation. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Molecular dynamics simulation of sucrose- and trehalose-coated carboxy-myoglobin

We performed a room temperature molecular dynamics (MD) simulation on a system containing 1 carboxy-myoglobin (MbCO) molecule in a sucrose-water matrix of identical composition (89% [sucrose/(sucrose + water)] w/w) as for a previous trehalose-water-MbCO simulation (Cottone et al., Biophys J 2001;80:931-938). Results show that, as for trehalose, the amplitude of protein atomic mean-square fluctuations, on the nanosecond timescale, is reduced with respect to aqueous solutions also in sucrose. A detailed comparison as a function of residue number evidences mobility differences along the protein backbone, which can be related to a different efficacy in bioprotection. Different heme pocket structures are observed in the 2 systems. The joint distribution of the magnitude of the electric field at the CO oxygen atom and of the angle between the field and the CO unit vector shows a secondary maximum in sucrose, absent in trehalose. This can explain the CO stretching band profile (A substates distribution) differences evidenced by infrared spectroscopy in sucrose- and trehalose-coated MbCO (Giuffrida et al., J Phys Chem B 2004;108:15415-15421), and in particular the appearance of a further substate in sucrose. Analysis of hydrogen bonds at the protein-solvent interface shows that the fraction of water molecules shared between the protein and the sugar is lower in sucrose than in trehalose, in spite of a larger number of water molecules bound to the protein in the former system, thus indicating a lower protein-matrix coupling, as recently observed by Fourier transform infrared (FTIR) experiments (Giuffrida et al., J Phys Chem B 2004;108:15415-15421).     

Molecular dynamics studies on the NMR structures of rabbit prion protein wild type and mutants: surface electrostatic charge distributions

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of mammalian species such as sheep and goats, cattle, deer and elk, and humans. But for rabbits, studies have shown that they have a low susceptibility to be infected by prion diseases. This paper does molecular dynamics (MD) studies of rabbit NMR structures (of the wild type and its two mutants of two surface residues), in order to understand the specific mechanism of rabbit prion proteins (RaPrP^C). Protein surface electrostatic charge distributions are specially focused to analyze the MD trajectories. This paper can conclude that surface electrostatic charge distributions indeed contribute to the structural stability of wild-type RaPrP^C; this may be useful for the medicinal treatment of prion diseases.     

Molecular dynamics insights on temperature and pressure effects on electroporation

Electroporation is a cell-level phenomenon caused by an ionic imbalance in the membrane, being of great relevance in various fields of knowledge. A dependence of the pore formation kinetics on the environmental conditions (temperature and pressure) of the cell membrane has already been reported, but further clarification regarding how these variables affect the pore formation/resealing dynamics and the transport of molecules through the membrane is still lacking. The objective of the present study was to investigate the temperature (288–348 K) and pressure (1–5000 atm) effects on the electroporation kinetics using coarse-grained molecular dynamics simulations. Results shown that the time for pore formation and resealing increased with pressure and decreased with temperature, whereas the maximum pore radius increased with temperature and decreased with pressure. This behavior influenced the ion migration through the bilayer, and the higher ionic mobility was obtained in the 288 K/1000 atm simulations, i.e., a combination of low temperature and (not excessively) high pressure. These results were used to discuss some experimental observations regarding the extraction of intracellular compounds applying this technique. This study contributes to a better understanding of electroporation under different thermodynamic conditions and to an optimal selection of processing parameters in practical applications which exploit this phenomenon.     

Conformational dynamics of recoverin's Ca2+-myristoyl switch probed by 15N NMR relaxation dispersion and chemical shift analysis

Recoverin, a member of the neuronal calcium sensor (NCS) branch of the calmodulin superfamily, serves as a calcium sensor in retinal rod cells. Ca²⁺‐induced conformational changes in recoverin promote extrusion of its covalently attached myristate, known as the Ca²⁺‐myristoyl switch. Here, we present nuclear magnetic resonance (NMR) relaxation dispersion and chemical shift analysis on ¹⁵N‐labeled recoverin to probe main chain conformational dynamics. ¹⁵N NMR relaxation data suggest that Ca²⁺‐free recoverin undergoes millisecond conformational dynamics at particular amide sites throughout the protein. The addition of trace Ca²⁺ levels (0.05 equivalents) increases the number of residues that show detectable relaxation dispersion. The Ca²⁺‐dependent chemical shifts and relaxation dispersion suggest that recoverin has an intermediate conformational state (I) between the sequestered apo state (T) and Ca²⁺ saturated extruded state (R): T ? I ? R. The first step is a fast conformational equilibrium ([T]/[I] < 100) on the millisecond time scale (τ_exδω < 1). The final step (I ? R) is much slower (τ_exδω > 1). The main chain structure of I is similar in part to the structure of half‐saturated E85Q recoverin with a sequestered myristoyl group. We propose that millisecond dynamics during T ? I may transiently increase the exposure of Ca²⁺‐binding sites to initiate Ca²⁺ binding that drives extrusion of the myristoyl group during I ? R. Proteins 2011; © 2011 Wiley‐Liss, Inc.     

13C and 1H solid state NMR investigation of hydration effects on gluten dynamics

The effect of hydration on the molecular dynamics of soft wheat gluten was investigated by solid state NMR. For this purpose, we recorded static and MAS ¹H spectra and SPE, CP, and other selective ¹³C spectra under MAS and dipolar decoupling conditions on samples of dry and H₂O and D₂O hydrated gluten. Measurements of carbon-proton CP times and several relaxation times (proton T₁, T_1ρ and T₂, and carbon T₁) were also performed. The combination of these techniques allowed both site-specific and domain-averaged motional information to be obtained in different characteristic frequency ranges. Domains with different structural and dynamic behaviour were identified and the changes induced by hydration on the dynamics of different domains could be monitored. The proton spin diffusion process was exploited to get information on the degree of mixing among different gluten domains. The results are consistent with the “loop and train” model proposed for hydrated gluten.     

Structure and quantum chemical analysis of NAD+-dependent isocitrate dehydrogenase: hydride transfer and co-factor specificity

Mutations in human coagulation factor IX cause an X-linked bleeding disorder Hemophilia B, which can be classified as severe, moderately severe and mild based on the plasma levels of factor IX among affected individuals with respect to normal factor IX activity assayed in the patients' plasma (<1%, 2-5%, 6-30%, respectively). Recently, we identified hemophilia B to be a disease with mutations showing clinical variation and speculated that this phenotypic heterogeneity might be a replacement-specific property. Here, we have analyzed the differences in sequence and structural properties among identical mutations with varying phenotypes (IMVPs) by comparing with mutations with uniform phenotypes (MUPs), with recurring reports in Haemophilia B mutation database. Classification of mutations into IMVPs and MUPs has been done based on rigorous systematic evaluation of the clotting activity each mutation is associated with. IMVPs (n = 51) occur in less conserved mutant sites with more tolerated substitutions compared to MUPs (n = 100). A preponderance of CpG site mutations and Arg as the mutated residue in IMVPs compared to Cys in MUPs was observed. Hence, a CpG site substitution at less conserved Arg site might have an increased propensity of expressing variable phenotypes. The changes in intrinsic properties associated with the mutation are less drastic for IMVPs than for MUPs, though no significant differences were observed in structural properties. Based on this study and available literature we speculate that modifier genes at other loci, epigenetic interactions and environment may serve individually or cumulatively to bring about the clinical variation implicating hemophilia B to be deviation from classical Mendelian disorder with complete penetrance. We demonstrate that phenotypic heterogeneity appears to be site-specific also owing to the lesser conservation of the mutant site.