Molecular dynamics study of Pseudomonas aeruginosa lectin-II complexed with monosaccharides

We present the results of a series of 10-ns molecular dynamics simulations on Pseudomonas aeruginosa lectin-II (PA-IIL) and its complexes with four different monosaccharides. We compare the saccharide-free, saccharide-occupied, and saccharide- and ion-free forms of the lectin. The results are coupled with analysis of the water density map and calcium coordination. The water density pattern around the binding site in the free lectin molecular dynamics was fitted with that in the X-ray and with the hydroxyl groups of the monosaccharide within the lectin/monosaccharide complexes and the best ligand was predicted based on the best fit. Interestingly, the water density pattern around the binding site in the uncomplexed lectin exactly fitted the O2, O3, and O4 hydroxyl groups of the fucose complex with the lectin. This observation could lead to a hypothesis that the replacement of these three water molecules from the binding site by the monosaccharide decreases the entropy of the complex and increases the entropy of the water molecules, which favors the binding. It suggests that the high density peaks of the solvent around the binding site in the free protein could be the tool to predict hydroxyl group orientation of the sugar in the protein/sugar complexes. The high affinity of PA-IIL binding site is also attributed to the presence of two calcium ions, each of them making five to six coordinations with the protein part and two coordinations with either water or the monosaccharide. When the calcium ions are removed from the simulated system, they are replaced by sodium ions from the solvent. These observations rationalize the high binding affinity of PA-IIL towards fucose.     

Enhanced crystal packing due to solvent reorganization through reductive methylation of lysine residues in oxidoreductase from <Emphasis Type="Italic">Streptococcus pneumoniae</Emphasis>

Protein crystallization is in part driven by the changes in the entropy of the system, but opinions differ as to whether the solute (protein) or solvent (water) molecules make more of a contribution to the overall entropic change. Methylation of lysine residues in proteins has been used to enhance protein crystallization. We investigated using molecular dynamics simulations with explicit solvent molecules, the behavior of several native proteins and their methylated counterparts chosen from an earlier large-scale study. Methylated lysines are capable of making a variety of interactions including H-bonds with protein residues and solvent. We demonstrate that methylation on the lysine slightly increases its side chain conformational entropy by about 3.5 J mol⁻¹ K⁻¹. Analysis of the radial and spatial distributions of the water molecules around the methylated lysine surface in oxidoreductase from Streptococcus pneumoniae revealed a larger sphere of water molecules with low entropy, as compared with solvent associated with unmethylated lysine. If methylated lysine were to make interactions at the protein–protein interface, the low-entropy water molecules associated with methylated lysines would be released, resulting in a gain of entropy. We show that this gain more than compensates for the loss of protein entropy. Therefore, we propose that lysine methylation favors the formation of crystals through solvent entropic gain.     

Reconstructing the protein–water interface

Using molecular dynamics simulations of fully hydrated proteins and analysis of crystal structures contained in the Protein Data Bank, we develop a transferable set of perpendicular radial distribution functions for water molecules around globular proteins. These universal functions may be used to reconstruct the unique three-dimensional solvent density distribution around every individual protein with a modest error. We discuss potential applications of this solvent treatment in protein x-ray crystallographic refinements and in theoretical modeling. We also present a fast, grid-based algorithm for construction of the perpendicular solvent density distributions. © 1998 John Wiley & Sons, Inc. Biopoly 45: 469–478, 1998     

一步法制备羧甲基茯苓多糖的工艺研究

本实验对在有机溶剂中一步法半合成羧甲基茯苓多糖的合成条件进行了研究。结果表明,乙醇是作为羧甲基化反应的合适介质。反应温度提高能加快反应速度;反应时间延长能提高取代度。茯苓多糖葡萄糖当量与氢氧化钠和一氯乙酸的摩尔比调配适当,能减少副产物羟乙酸钠的产生。     

Hydration of cis and trans N-methylformamide as revealed by the use of 17O-NMR,molecular mechanics,and ab initio calculations

The solvation of cis and trans N-methylformamide (NMF) by water was investigated using a combination of ¹⁷O-nmr spectroscopy, classical molecular mechanics [MM2(77) and MM2(87)] force field, and ab initio 4-31G* gradient optimization calculations. In dilute aqueous solution, the ¹⁷O-nmr spectra of NMF indicate strong shielding by 66.9 and 66.1 ppm for the cis and trans amide oxygens, respectively, compared to those values obtained in dilute toluene solution. This demonstrates that both isomers are equally solvated by molecules of water, which are further hydrogen bonded to molecules of water of the bulk solvent. Molecular mechanics simulations were carried out for cis and trans NMF in a cluster of water molecules. Radial distribution functions show structural contacts by several water molecules at the amide CO and NH group, which are significantly more pronounced with MM2 (87) calculations. Ab initio 4-31G* gradient optimization calculations on the supermolecule trans NMF-(H₂O)₃ indicates the presence of more than two hydrogen-bond contacts at the carbonyl oxygen. This is in agreement with MM2 calculations and provides further evidence for multiple acceptor properties of the amide oxygen and an out of the amide plane arrangement of the bound molecules of water. Comparison of the integration data to the first radial distribution function (rdf) minima shows that the local solvation of the CO and NH groups is very similar for both cis and trans isomers. The intermolecular geometric parameters of the supermolecule trans NMA–(H₂O)₃ and the first rdf maxima resulting from MM2 (87) and MM2 (77) calculations are compared with distribution of water molecules around the CO and NH groups of peptides and proteins resulting from x-ray and neutron diffraction experiments. The rdfs involving the methyl group of NMF demonstrate the nonrandom distribution of solvent sites with first maxima in reasonable agreement with distribution of water molecules around the apolar side chain of amino acid residues in proteins. © 1995 John Wiley & Sons, Inc.     

Modeling high-resolution hydration patterns in correlation with DNA sequence and conformation

Hydration around the DNA fragment d(C5T5).(A5G5) is presented from two molecular dynamics simulations of 10 and 12 ns total simulation time. The DNA has been simulated as a flexible molecule with both the CHARMM and AMBER force fields in explicit solvent including counterions and 0.8 M additional NaCl salt. From the previous analysis of the DNA structure B-DNA conformations were found with the AMBER force-field and A-DNA conformations with CHARMM parameters. High-resolution hydration patterns are compared between the two conformations and between C.G and T.A base-pairs from the homopolymeric parts of the simulated sequence. Crystallographic results from a statistical analysis of hydration sites around DNA crystal structures compare very well with the simulation results. Differences between the crystal sites and our data are explained by variations in conformation, sequence, and limitations in the resolution of water sites by crystal diffraction. Hydration layers are defined from radial distribution functions and compared with experimental results. Excellent agreement is found when the measured experimental quantities are compared with the equivalent distribution of water molecules in the first hydration shell. The number of water molecules bound to DNA was found smaller around T.A base-pairs and around A-DNA as compared to B-DNA. This is partially offset by a larger number of water molecules in hydrophobic contact with DNA around T.A base-pairs and around A-DNA. The numbers of water molecules in minor and major grooves have been correlated with helical roll, twist, and inclination angles. The data more fully explain the observed B-->A transition at low humidity.     

Contribution of solvent water to the solution X-ray scattering profile of proteins

A theoretical framework is presented to analyze how solvent water contributes to the X-ray scattering profile of protein solution. Molecular dynamics simulations were carried out on pure water and an aqueous solution of myoglobin to determine the spatial distribution of water molecules in each of them. Their solution X-ray scattering (SXS) profiles were numerically evaluated with obtained atomic-coordinate data. It is shown that two kinds of contributions from solvent water must be considered to predict the SXS profile of a solution accurately. One is the excluded solvent scattering originating in exclusion of water molecules from the space occupied by solutes. The other is the hydration effect resulting from formation of a specific distribution of water around solutes. Explicit consideration of only two molecular layers of water is practically enough to incorporate the hydration effect. Care should be given to using an approximation in which an averaged electron density distribution is assumed for the structure factor because it may predict profiles considerably deviating from the correct profile at large K.     

Hydration of oligonucleotides in crystals

The solvent molecules found around crystallized oligonucleotides after X-ray refinement are analysed in terms of interaction sites to bases, phosphates and sugars in the three main forms of nucleic acid structures, the A-form, the B-form and the Z-form. The average numbers of contacts to nucleic acid atoms made by solvent molecules are identical in the three forms, but it appears that the average number of contacts solvent molecules make with each other depends on the resolution of the structure. The phosphate anionic oxygen atoms are the most hydrated, while the O(3′) and O(5′) backbone atoms and the ring oxygen atom O(4′) are the least hydrated. Among the hydrophilic atoms of the bases, there is a modulation of the relative water affinities with the nucleic acid form. Numerous hydration sites are such that water molecules can bridge hydrophilic atoms of the same residue, of adjacent residues on the same strand, of distant residues on the two strands, or belonging to symmetry-related residues. Through the helical periodicity of the nucleic acid structure, those bridges can lead to regular and striking hydration networks involving several water molecules and characteristic of the nucleic acid form. Solvent dynamics, as seen by temperature factor versus occupancy plots, seems intimately related to nucleic acid structure and dynamics, since they depend on hydration sites around the nucleic acids.     

Inspecting fluctuation and coordination around chromophore inside green fluorescent protein from water to nonpolar solvent

Green fluorescent protein (GFP) is a widely used biomarker that demands systematical rational approaches to its structure function redesign. In this work, we mainly utilized atomistic molecular dynamics simulations to inspect and visualize internal fluctuation and coordination around chromophore inside GFP, from water to nonpolar octane solvent. We found that GFP not only maintains its β-barrel structure well into the octane, but also sustains internal residue and water coordination to position the chromophore stably while suppress dihedral fluctuations of the chromophore, so that functional robustness of GFP is achieved. Our accompanied fluorescence microscope measurements accordingly confirmed the GFP functioning into the octane. Furthermore, we identified that crucial water sites inside GFP along with permeable pores on the β-barrel of the protein are largely preserved from the water to the octane solvent, which allows sufficiently fast exchanges of internal water with the bulk or with the water layer kept on the surface of the protein. By additionally pulling GFP from bulk water to octane, we suggest that the GFP function can be well maintained into the nonpolar solvent as long as, first, the protein does not denature in the nonpolar solvent nor across the polar-nonpolar solvent interface; second, a minimal set of water molecules are in accompany with the protein; third, the nonpolar solvent molecules may need to be large enough to be nonpermeable via the water pores on the β-barrel.     

Molecular dynamics study of onset of water gelation around the collagen triple helix

The onset of water gelation around a collagen-like triple helix peptide was studied at ambient temperature and pressure by performing Molecular Dynamics simulations. The radial distribution functions of the oxygen and hydrogen atoms of water are distorted below 4 A from the peptide. The distortion is accompanied by the breakdown of the tetrahedral coordination of the hydrogen-bonded network of water molecules. The water shell around the peptide consists of alternating regions of higher and lower density. In agreement with experiments we find that the first hydration shell is kinetically labile, with a residence time in the order of picoseconds for a water molecule. From the computed diffusion coefficient, a key measure of the collective dynamics, we estimate the average diffusion speed decreases by a factor of 1.5 close to the peptide compared to the liquid. Our results give new insight in gel formation and structure on a molecular level.     

Myosin-induced volume increase of the hyper-mobile water surrounding actin filaments

Microwave dielectric spectroscopy can measure the rotational mobility of water molecules that hydrate proteins and the hydration-shell volume. Using this technique, we have recently shown that apart from typical hydrating water molecules with lowered mobility there are other water molecules around the actin filaments (F-actin) which have a much higher mobility than that of bulk water [Biophys. J. 85 (2003) 3154]. We report here that the volume of this water component (hyper-mobile water) markedly increases without significant change of the volume of the ordinary hydration shell when the myosin motor-domain (S1, myosin subfragment-1) binds to F-actin. No hyper-mobile component was found in the hydration shell of S1 itself. The present results strongly suggest that the solvent space around S1 bound to F-actin is diffusionally asymmetric, which supports our model of force generation by actomyosin proposed previously [op. cit.].     

Molecular Dynamics Study on Protein and it's Water Structure at High Pressure

Abstract

We have performed NPT molecular dynamics simulations (Langevin Piston Method) on two types of solvated proteins-‘denaturation-unfavorable’ protein (insulin) and ‘denaturation-favorable protein’ (ribonuclease A) at high pressure (from 1 bar up to 20 kbar). The method is based on the extended system formalism introduced by Andersen, where the deterministic equations of motion for the piston degree of freedom are replaced by Langevin equation. We report the structural changes of proteins (ribonuclease A and insulin) and water molecules through radius of gyration, solvent accessible surface area, hydrogen bond pattern, and the topology of water clusters connected by the hydrogen bonded circular network. The solvent accessibility of ribonuclease A is mainly decreased by hydrophilic residues rather than hydrophobic residues under high pressure. From the results of hydrogen bond analysis, we have found that α-helix is more stable than β-sheet under high pressure. In addition, from the analysis of the water cluster, we have observed that for ribonuclease A, 5-membered ring structure is more favorable than 6-membered ring at higher pressure. However, for insulin, the ratio of 5 to 6-ring is constant over the pressure ranges for which we have performed MD simulation. This indicates that the water structure around insulin does not change under high pressure.     

Structural and dynamics analysis of matrix metalloproteinases MMP-2 complexed with chemically modified tetracyclines (CMTs)

Matrix metalloproteinases (MMPs) play a critical role in physiological processes and pathological conditions such tumor invasion and metastasis. In recent years, a number of MMP inhibitors have been proposed, including the chemically modified tetracyclines (CMTs), which have been evaluated in preclinical cancer models showing promising results. This work provides insights into the structure and dynamics of the MMP-2 catalytic domain complexed with seven CMT (CMT-n), based on the analysis of molecular dynamics trajectories in solution. The comparative analysis of various relevant molecular aspects of the different complexes of MMP-2 and CMT-n derivatives was performed aiming to elucidate the effect of ligands on the enzyme structure. These include the radial distribution function of the water molecules around the catalytic zinc, the solvent accessible surface area for the inhibitors and the root-mean-square fluctuation for all amino acid residues. The results help to understand the differences in the binding modes of related compounds and, therefore, add to further design of novel tetracycline-based inhibitors for MMP enzymes.     

Integral equation theories for predicting water structure around molecules

Water plays a crucial role in the structure and function of proteins and other biological macromolecules; thus, theories of aqueous solvation for these molecules are of great importance. However, water is a complex solvent whose properties are still not completely understood. Statistical mechanical integral equation theories predict the density distribution of water molecules around a solute so that all particles are fully represented and thus potentially both molecular and macroscopic properties are included. Here we discuss how several theoretical tools we have developed have been integrated into an integral equation theory designed for globular macromolecular solutes such as proteins. Our approach predicts the three-dimensional spatial and orientational distribution of water molecules around a solute. Beginning with a three-dimensional Ornstein-Zernike equation, a separation is made between a reference part dependent only on the spatial distribution of solvent and a perturbation part dependent also on the orientational distribution of solvent. The spatial part is treated at a molecular level by a modified hypernetted chain closure whereas the orientational part is treated as a Boltzmann prefactor using a quasi-continuum theory we developed for solvation of simple ions. A potential energy function for water molecules is also needed and the sticky dipole models of water, such as our recently developed soft-sticky dipole (SSD) model, are ideal for the proposed separation. Moreover, SSD water is as good as or better than three point models typically used for simulations of biological macromolecules in structural, dielectric and dynamics properties and yet is seven times faster in Monte Carlo and four times faster in molecular dynamics simulations. Since our integral equation theory accurately predicts results from Monte Carlo simulations for solvation of a variety of test cases from a single water or ion to ice-like clusters and ion pairs, the application of this theory to biological macromolecules is promising.     

QM/MM study of D-fructose in aqueous solution

The QM/MM molecular dynamics methodology was applied to the study of the two main D-fructose tautomers present in aqueous solution, beta-D-fructofuranose and beta-D-fructopyranose. The solute was treated at the AM1 semi-empirical level, and for the solvent water molecules we used the TIP3P potential. We analyzed the structure of the water molecules around the hydroxyl groups to explain the differences in sweet taste between the two tautomers.     

Structure of citrus pectins and viscometric study of their solution properties

Citrus pectins with degrees of methylation between 30 and 72% were carefully characterized in order to determine their charge density and molecular weight distribution, the content in galacturonic acid and in neutral sugars, the degree of methylation and acetylation. Using enzymic degradation it has been found that pectin molecules consist mainly of long homogalacturonan regions with some regions of neutral sugars as side chains attached on rhamnose residues. The viscometric behaviour of the different samples indicates that 0.1 M NaCl, at 25 degrees C, is a good solvent of sodium pectinates. From the evolution of the Huggins parameter, it appears that pectins with 50% of methylated galacturonic groups exhibit a maximum flexibility. A Mark-Houwink exponent of 0.8 has been found in good agreement with theoretical predictions for flexible polymers in a good solvent.     

How do immobilised cell-adhesive Arg–Gly–Asp-containing peptides behave at the PAA brush surface?

Bio-engineered surfaces that aim to induce normal cell behaviour in vitro need to ‘mimic’ the extracellular matrix in a way that allows cell adhesion. In this computational work, several model cell-binding peptides with a minimal cell-adhesive Arg–Gly–Asp sequence are investigated in the bulk as well as immobilised on a soft surface. For this reason, a combination of density functional theory and all-atom MD simulations is applied. The major goal of the modelling is to characterise the accessibility of the cell-recognition motif on the functionalised soft polymer surface. As a reference system, the behaviour of three peptide sequences is preliminarily studied in explicit water simulations. From the analysis of the MD trajectories, the solvent accessible surface area, the distribution of water molecules around peptide groups, the secondary structure and the thermodynamics of hydration are evaluated. Furthermore, each peptide is immobilised on the surface of a homopolymer poly(acrylic acid) brush. During MD simulations, all three peptides approach closely toward PAA brush, and their surface accessibility is characterised. Although the peptides are adsorbed onto the brush, they are not hidden by the polymer strands, with RGD unit accessible on the surface and available for guided cell adhesion.     

Theoretical studies on the tautomerism of tetrazole selenone

The tautomerism of all possible forms of tetrazole selenone (A–G), induced by proton transfer, was studied, theoretically, in different environments including gas phase, continuum solvent and microsolvated environment with one or two explicit water or ammonia molecules. The calculations were performed using two different levels of theory including mPW2PLYP and DFT-B3LYP. The 6-311++G(d,p) basis set was used for C, H, O and N and the standard relativistic effective core pseudo potential LANL2DZ basis set was used for Se atom. It was found that the tetrazole selenone, in the form of A, is the most stable isomer in all of the environments considered in this work. The kinetics of proton transfer reaction was studied in both gas and solvent environments and it was concluded that the activation energy of the reaction increases with going from the gas phase to polar solvents. Moreover, the proton transfer reaction assisted by one or two water or ammonia molecules was investigated and it was found that the activation energy significantly reduces.     

Molecular dynamics simulations of human carbonic anhydrase II: Insight into experimental results and the role of solvation

In this paper, the carbonic anhydrase II (CA II) enzyme active site is modeled using ab initio calculations and molecular dynamics simulations to examine a number of important issues for the enzyme function. It is found that the Zn²⁺ ion is dominantly tetrahedrally coordinated, which agrees with X-ray crystallographic studies. However, a transient five-fold coordination with an extra water molecule is also found. Studies of His64 conformations upon a change in the protonation states of the Zn-bound water and the His64 residue also confirm the results of an X-ray study which suggest that the His64 conformation is quite flexible. However, the degree of water solvation is found to affect this behavior. Water bridge formation between the Zn-bound water and the His64 residue was found to involve a free energy barrier of 2–3 kcal/mol and an average lifetime of several picoseconds, which supports the concept of a proton transfer mechanism through such a bridge. Mutations of various residues around the active site provide further insight into the corresponding experimental results and, in fact, suggest an important role for the solvent water molecules in the CA II catalytic mechanism. Proteins 33:119–134, 1998. © 1998 Wiley-Liss, Inc.