The dynamics of the calcium‐induced chain–chain association in the polyuronate systems

The calcium‐induced formation of strong, hydrophilic gels is the important feature of polyuronates, connected with most of their practical applications. The insight into the molecular details of gelling process dynamics is hardly feasible for both experimental and theoretical methods. Here, the application of the transition path sampling method for studying this problem is reported; the focus was on the poly(α‐L ‐guluronate) systems, treated as the representative for all polyuronate‐containing systems. The results allowed for identifying several distinct local minima of the free energy lying on the transition paths and visited by the system during the process of chain–chain association. These minima usually correspond to the intermediate structures in which the water molecules bridge calcium ion and carboxyl groups. This work emphasizes the importance of water and provides more complete understanding of the calcium binding by the polyuronate chains. © 2012 Wiley Periodicals, Inc.     

诺氟沙星-DNA复合物的分子动力学模拟

采用分子模建的方法构建了诺氟沙星-DNA复合物的初始结构, 通过2 ns的分子动力学(MD)模拟研究表明: 诺氟沙星能够和双螺旋d[ATATCGATAT]₂形成稳定的复合物, 药物分子可紧密结合在DNA的小沟区域, 并且能够与DNA的鸟嘌呤碱基形成两个稳定的氢键. 在分子水平上提供了诺氟沙星直接与双螺旋DNA相互作用的结构及复合物的动态变化情况.     

Molecular mobility of nitroxide spin probes in glassy polymers. Quasi‐libration model

ESR spectra of three spin probes with different molecular volumes: 2,2,6,6‐tetramethyl‐4‐oxopiperidine‐1‐oxyl, di‐p‐anisylnitroxide, and nitroxide derivative of fullerene in glassy polystyrene, polyvinyl trimethylsilane, and Teflon AF‐2400 were calculated numerically within the model of quasi‐libration motions. Temperature ranges, where the model is capable to reproduce spectra within experimental errors, were defined. It was found that simulation of X‐band ESR spectra allows to determine quasi‐libration amplitudes around molecular axes X and Y with accuracy ～ 3° and around Z axis with accuracy ～ 15–20°. A shape of distribution of quasi‐libration amplitudes was also determined qualitatively by ESR spectra simulations. It was established that the average amplitude of quasi‐libration motion depends on the free volume of each polymer and geometrical molecular volume of a spin probe. Quasi‐libration amplitudes increase as the temperature increases, and reach the value of 40 degrees. We found that upon further temperature increase, quasi‐libration model becomes inapplicable for quantitative numerical spectra simulation. © 2008 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 47: 107–120, 2009     

Molecular Modeling of the Three‐Dimensional Structure of Human Sphingomyelin Synthase

Sphingomyelin synthase (SMS) produces sphingomyelin and diacylglycerol from ceramide and phosphatidylcholine. It plays an important role in cell survival and apoptosis, inflammation, and lipid homeostasis, and therefore has been noticed in recent years as a novel potential drug target. In this study, we combined homology modeling, molecular docking, molecular dynamics simulation, and normal mode analysis to derive a three‐dimensional structure of human sphingomyelin synthase (hSMS1) in complex with sphingomyelin. Our model provides a reasonable explanation on the catalytic mechanism of hSMS1. It can also explain the high selectivity of hSMS1 towards phosphocholine and sphingomyelin as well as some other known experimental results about hSMS1. Moreover, we also derived a complex model of D609, the only known small‐molecule inhibitor of hSMS1 so far. Our hSMS1 model may serve as a reasonable structural basis for the discovery of more effective small‐molecule inhibitors of hSMS1.     

Comparison of thermodynamic integration and Bennett acceptance ratio for calculating relative protein‐ligand binding free energies

The performances of Bennett's acceptance ratio method and thermodynamic integration (TI) for the calculation of free energy differences in protein simulations are compared. For the latter, the standard trapezoidal rule, Simpson's rule, and Clenshaw‐Curtis integration are used as numerical integration methods. We evaluate the influence of the number and definition of intermediate states on the precision, accuracy, and efficiency of the free energy calculations. Our results show that non‐equidistantly spaced intermediate states are in some cases beneficial for the TI methods. Using several combinations of softness parameters and the λ power dependence, it is shown that these benefits are strongly dependent on the shape of the integrand. Although TI is more user‐friendly due to its simplicity, it was found that Bennett's acceptance ratio method is the more efficient method. It is also the least dependent on the choice of the intermediate states, making it more robust than TI. © 2013 Wiley Periodicals, Inc.     

Exploration of the Ca2+ interaction modes of the nifedipine calcium channel antagonist.

A comprehensive study is carried out using quantum chemical computation and molecular dynamics (MD) simulations to gain insight into the interaction between Ca(2+) ions and the most important class of calcium channel antagonists--nifedipine. First, the chelating structures and energetic characters of nifedipine-Ca(2+) in the gas phase are explored, and 25 isomers are found. The most favorable chelating mode is a tridentate one, that is, Ca(2+) binds to two carbonyl O atoms and one nitryl O atom, where Ca(2+) is above the plane of the three O atoms to form a pyramidal structure. Accurate geometric structures, relative stabilities, vertical and adiabatic binding energies, and charge distributions are discussed. The differences in the geometries and energies among these isomers are analyzed from the contributions of chelating sites, electrostatics and polarizations, steric repulsions, and charge distributions. The interconversions among isomers with similar geometries and energies are also investigated because of the importance of the geometric transformation in the biological system. Furthermore, certain numbers of water molecules are added to the nifedipine-Ca(2+) system to probe the effect of water. A detailed study is performed on the hydrated geometries on the basis of the most stable isomer 1. Stepwise hydration can weaken the nifedipine-Ca(2+) interaction, and the chelating sites of nifedipine are gradually replaced by the added water molecules. Hexacoordination is found to be the most favorable geometry no matter how many water molecules were added, which can be verified by the MD simulations. The transfer of water molecules from the inner shell to the outer shell is also supported by MD simulations of the hexahydrated complexes.     

Molecular dynamics simulations and free energy calculations on the enzyme 4‐hydroxyphenylpyruvate dioxygenase

4‐Hydroxyphenylpyruvate dioxygenase is a relevant target in both pharmaceutical and agricultural research. We report on molecular dynamics simulations and free energy calculations on this enzyme, in complex with 12 inhibitors for which experimental affinities were determined. We applied the thermodynamic integration approach and the more efficient one‐step perturbation. Even though simulations seem well converged and both methods show excellent agreement between them, the correlation with the experimental values remains poor. We investigate the effect of slight modifications on the charge distribution of these highly conjugated systems and find that accurate models can be obtained when using improved force field parameters. This study gives insight into the applicability of free energy methods and current limitations in force field parameterization. © 2011 Wiley Periodicals, Inc. J Comput Chem 2011     

Modeling of Structural,Energetic, and Dynamic Properties of Few‐Atom Silver Clusters Embedded in Polynucleotide Strands by Using Molecular Dynamics

This work concerns the study of the structural, energetic, and dynamic properties of fluorescent systems composed of silver clusters stabilized by polynucleotide strands. To do so, classical interaction potentials relative to silver, neutral and cationic, were introduced in the AMBER force field. Molecular dynamics simulations allowed analysis of the nature and force of the interactions between the various parts of the nucleic oligomers and the silver clusters. Conformational analyses were necessary to explore the flexibility of the supramolecular assemblies, specifically by radial distribution functions and Ramachandran‐type maps.     

Postgel properties in the statistical crosslinking of heterochains. I. Systems with N types of chains

The heterochain crosslinking model describes nonrandom crosslinking of polymer chains and is an extension of the classical Flory/Stockmayer gelation theory. We consider the postgelation relationship for the system consisting of N types of polymer chains, in which the probability that a crosslink point on an i‐type chain is connected to a j‐type chain is explicitly given by p_ij. The analytical solutions for the weight fraction of the sol, the number‐average and weight‐average molecular weights within the sol fraction, and the crosslinking density within the sol and gel fractions are derived for the systems, with each type of chain conforming to the Schulz–Zimm distribution. Illustrative calculations are shown for the systems consisting of two and three types of chains, and the obtained results agree with those from the Monte Carlo method. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 2333–2341, 2000     

On the influence of charged side chains on the folding-unfolding equilibrium of beta-peptides: a molecular dynamics simulation study

The influence of charged side chains on the folding-unfolding equilibrium of beta-peptides was investigated by means of molecular dynamics simulations. Four different peptides containing only negatively charged side chains, positively charged side chains, both types of charged side chains (with the ability to form stabilizing salt bridges) or no charged side chains were studied under various conditions (different simulation temperatures, starting structures and solvent environment). The NMR solution structure in methanol of one of the peptides (A) has already been published; the synthesis and NMR analysis of another peptide (B) is described here. The other peptides (C and D) studied herein have hitherto not been synthesized. All four peptides A-D are expected to adopt a left-handed 3(14)-helix in solution as well as in the simulations. The resulting ensembles of structures were analyzed in terms of conformational space sampled by the peptides, folding behavior, structural properties such as hydrogen bonding, side chain-side chain and side chain-backbone interactions and in terms of the level of agreement with the NMR data available for two of the peptides. It was found that the presence of charged side chains significantly slows down the folding process in methanol solution due to the stabilization of intermediate conformers with side chain-backbone interactions. In water, where the solvent competes with the solute-solute polar interactions, the folding process to the 3(14)-helix is faster in the simulations.     

A Central Partition of Molecular Conformational Space. III. Combinatorial Determination of the Volume Spanned by a Molecular System

In the first work of this series (Gabarro-Arpa, Comp. Biol. Chem. 27 (2003) 153–159) it was shown that the conformational space of a molecule could be described to a fair degree of accuracy by means of a central hyperplane arrangement. The hyperplanes divide the space into a hierarchical set of cells that can be encoded by the face lattice poset of the arrangement. The model however, lacked explicit rotational symmetry, which made impossible to distinguish rotated structures in conformational space. This problem was solved in a second work (Gabarro-Arpa, Proc. 26th Ann. Int. Conf. of the IEEE EMBS (San Franciso, 2004) 3007–3010) by sorting the elementary 3-dimensional components of the molecular system into a set of morphological classes that can be properly oriented in a standard 3-D reference frame. This also made possible to find a solution to the problem that is being addressed in the present work: for a molecular system immersed in a heat bath we want to enumerate the subset of cells in conformational space that are visited by the molecule in its thermal wandering. If each visited cell is a vertex on a graph with edges to the adjacent cells, here it is explained how such graph can be built.     

A computational model of the nicotinic acetylcholine binding site

We have derived a model of the nicotinic acetylcholine binding site. This was accomplished by using three known agonists (acetylcholine, nicotine and epibatidine) as templates around which polypeptide side chains, found to be part of the receptor cavity from published molecular biology studies, are allowed to flow freely in molecular dynamics simulations and mold themselves around these templates. The resulting supramolecular complex should thus be a complement, both in terms of steric effects as well as electronic effects, to the agonists and it should be a good estimation of the true receptor cavity structure. The shapes of those minireceptor cavities equilibrated rapidly on the simulation time scale and their structural congruence is very high, implying that a satisfactory model of the nicotinic acetylcholine binding site has been achieved. The computational methodology was internally tested against two rigid and specific antagonists (dihydro--erytroidine and erysoidine), that are expected to give rise to a somewhat differently shaped binding site compared to that derived from the agonists. Using these antagonists as templates there were structural reorganizations of the initial receptor cavities leading to distinctly different cavities compared to agonists. This indicates that adequate times and temperatures were used in our computational protocols to achieve equilibrium structures for the agonists. Overall, both minireceptor geometries for agonists and antagonists are similar with the exception of one amino acid (ARG209).     

Validation by Molecular Dynamics of the Major Components of Sugarcane Vinasse,On a Surface of Calcium Carbonate (Calcite)

There is ongoing interest in the alcohol industry to significantly reduce and/or add value to the liquid residue, vinasse, produced after the distillation and rectification of ethanol from sugar cane. Vinasse contains potassium, glycerol, and a protein component that can cause environmental issues if improperly disposed of. Currently, some industries have optimized their processes to reduce waste, and a significant proportion of vinasse is being considered for use as an additive in other industrial processes. In the manufacture of cement and asphalt, vinasse has been used in the mixtures at low concentrations, albeit with some physical and mechanical problems. This work is the first molecular approximation of the components of the sugar cane vinasse in an industrial context, and it provides atomic details of complex molecular events. In the current study, the major components of sugar cane vinasse, alone or complexed on the surface of calcium carbonate, were modeled and simulated using molecular dynamics. The results showed that the protein component, represented by the mannoprotein Mp1p, has a high affinity for forming hydrogen bonds with potassium and glycerol in the vinasse. Additionally, it provides atomic stability to the calcium carbonate surface, preserving the calcite crystalline structure in the same way potassium ions interact with the carbonate group through ion–dipole interactions to improve the cohesion of the modeled surface. On the contrary, when the glycerol molecule interacts with calcium carbonate using more than two hydrogen bonds, it triggers the breakdown of the crystalline structure of calcite expanding the ionic pair.     

Molecular dynamics simulation exploration of cooperative migration mechanism of calcium ions in sarcoplasmic reticulum Ca2+‐ATPase

Calcium ATPase is a member of the P‐type ATPase, and it pumps calcium ions from the cytoplasm into the reticulum against a concentration gradient. Several X‐ray structures of different conformations have been solved in recent years, providing basis for elucidating the active transport mechanism of Ca²⁺ ions. In this work, molecular dynamics (MD) simulations were performed at atomic level to investigate the dynamical process of calcium ions moving from the outer mouth of the protein to their binding sites. Five initial locations of Ca²⁺ ions were considered, and the simulations lasted for 2 or 6 ns, respectively. Specific pathways leading to the binding sites and large structural rearrangements around binding sites caused by uptake of calcium ions were identified. A cooperative binding mechanism was observed from our simulation. Firstly, the first Ca²⁺ ion binds to site I , and then, the second Ca²⁺ ion approaches. The interactions between the second Ca²⁺ and the residues around site I disturb the binding state of site I and weaken its binding ability for the first bound Ca²⁺. Because of the electrostatic repulsion of the second Ca²⁺ and the electrostatic attraction of site II , the first bound Ca²⁺ shifts from site I to site II . Concertedly, the second Ca²⁺ binds to site I , forming a binding state with two Ca²⁺ ions, one at site I and the other at site II . Both of Glu908 and Asp800 coordinate with the two Ca²⁺ ions simultaneously during the concerted binding process, which is believed to be the hinge to achieve the concerted binding. In our simulations, four amino acid residues that serve as the channel to link the outer mouth and the binding sites during the binding process were recognized, namely Tyr837, Tyr763, Asn911, and Ser767. The analyses regarding the activity of the proteins via mutations of some key residues also supported our cooperative mechanism. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

Effects of Moisture on the Mechanical Properties of Microcrystalline Cellulose and the Mobility of the Water Molecules as Studied by the Hybrid Molecular Mechanics–Molecular Dynamics Simulation Method

A hybrid molecular mechanics–molecular dynamics simulation method has been performed to study the effects of moisture content on the mechanical properties of microcrystalline cellulose (MCC) and the mobility of the water molecules. The specific volume and diffusion coefficient of the water increase with increasing moisture content in the range studied of 1.8–25.5 w/w%, while the Young's modulus decreases. The simulation results are in close agreement with the published experimental data. Both the bound scission and free‐volume mechanisms contribute to the plasticization of MCC by water. The Voronoi volume increases with increasing moisture content. It is related to the free volume and the increase enhances the mobility of the water molecules and thus increases the coefficient of diffusion of the water. Moreover, with increasing moisture content, the hydrogen bonding per water molecule between MCC–water molecules decreases, thus increasing the water mobility and number of free water molecules. © 2019 The Authors. Journal of Polymer Science Part B: Polymer Physics published by Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 454–464     

Modeling of saccharide utilization in primary beer fermentation with yeasts immobilized in calcium alginate

Immobilized beer fermentation was studied using an industrial bottom-fermenting yeast strain Saccharomyces cerevisiae. The yeast cells were immobilized in 2.5% calcium alginate gel and used for brewing in a five-vessel cascade reactor. The fermentation was performed at 15°C at various flow rates. A nonstructured mathematical model was developed to simulate the performance of continuous primary fermentation of lager beer. The model was based on the following variables: maltose, maltotriose, glucose, fructose, ethanol, and cell concentration. Experimental values of these variables were determined in samples taken at regular intervals. For experimental data fitting a nonlinear regression was used. Substrate consumption was characterized by specific substrate consumption rate and saturation constant. The values of these two parameters were optimized for all four substrates. Inhibition effects of substrates and product were analyzed using various inhibition patterns. Only the inhibition effect of maltose on maltose consumption was clearly identified. A good-fitting relationship for maltose inhibition was found, and inhibition constants were calculated.     

Estimating the temperature dependence of peptide folding entropies and free enthalpies from total energies in molecular dynamics simulations

The temperature dependence of thermodynamic quantities, such as heat capacity, entropy and free enthalpy, may be obtained by using well-known equations that relate these quantities to the enthalpy of the molecular system of interest at a range of temperatures. In turn, the enthalpy of a molecular system can be estimated from molecular dynamics simulations of an appropriate model. To demonstrate this, we have investigated the temperature dependence of the enthalpy, heat capacity, entropy and free enthalpy of a system that consists of a beta-heptapeptide in methanol and have used the statistical mechanics relationships to describe the thermodynamics of the folding/unfolding equilibrium of the peptide. The results illustrate the power of current molecular simulation force fields and techniques in establishing the link between thermodynamic quantities and conformational distributions.