Molecular Dynamics of Phenylalanine Transfer RNA

Abstract

The atomic motions of yeast phenylalanine transfer RNA have been simulated using the molecular dynamics algorithm. Two simulations were carried out for a period of 12 picoseconds, one with a normal Van der Waals potential and the other with a modified Van der Waals potential intended to mimic the effect of solvent. An analysis of large scale motions, surface exposure, root mean square displacements, helical oscillations and relaxation mechanisms reveals the maintenance of stability in the simulated structures and the general similarity of the various dynamic features of the two simulations. The regions of conformational flexibility and rigidity for tRNA^Phe have been shown in a quantitative measure through this approach.     

基于遗传算法预测2D三向的蛋白质结构

本文基于范德华力势能预测2D三向的蛋白质结构。首先,将蛋白质结构预测这一生物问题转化为数学问题,并建立基于范德华力势能函数的数学模型。其次,使用遗传算法对数学模型进行求解,为了提高蛋白质结构预测效率,我们在标准遗传算法的基础上引入了调整算子这一概念,改进了遗传算法。最后,进行数值模拟实验。实验的结果表明范德华力势能函数模型是可行的,同时,和规范遗传算法相比,改进后的遗传算法能够较大幅度提高算法的搜索效率,并且遗传算法在蛋白质结构预测问题上有巨大潜力。     

Shape group studies of molecular similarity: relative shapes of Van der Waals and electrostatic potential surfaces of nicotinic agonists

In drug design, the usual strategy involves characterizing and comparing the shapes of molecules. We apply a simple method to accomplish this goal: determining the symmetry-independent shape groups (homology groups of algebraic topology) of a molecular surface.In this paper, we have adapted the method to describing the interrelation between Van der Waals and electrostatic potential surfaces. We describe rigorously the shape features in a series of molecules by using specific ranges of electrostatic potential over a Van der Waals surface. We consider a series of four nicotinic agonists as an example and discuss their expected activities as potential drugs on the basis of the shape similarities found.     

Electrostatic interaction in atomic force microscopy

In atomic force microscopy, the stylus experiences an electrostatic force when imaging in aqueous medium above a charged surface. This force has been calculated numerically with continuum theory for a silicon nitrite or silicon oxide stylus. For comparison, the Van der Waals force was also calculated. In contrast to the Van der Waals attraction, the electrostatic force is repulsive. At a distance of 0.5 nm the electrostatic force is typically 10^-12-10^-10 N and thus comparable in strength to the Van der Waals force. The electrostatic force increases with increasing surface charge density and decreases roughly exponentially with distance. It can be reduced by imaging in high salt concentrations. Below surface potentials of ≈50 mV, a simple analytical approximation of the electrostatic force is described.     

Free energy potential for aggregation of giant, neutral lipid bilayer vesicles by Van der Waals attraction.

Here, we report the first direct observation of Van der Waals' attraction between biomembrane capsules using measurements of the free energy reduction per unit area of membrane-membrane contact formation. In these studies, the membrane capsules were reconstituted neutral (egg phosphatidylcholine) lipid bilayers of giant (greater than 10(-3) cm diam) vesicles. Micromanipulation methods were used to select and maneuver two vesicles into proximity for contact; after adhesion was allowed to occur, the extent of contact formation was regulated through the vesicle membrane tensions that were controlled by micropipette suction. The free energy reduction per unit area of contact formation was proportional to the membrane tension multiplied by a simple function of the pipette and vesicle dimensions. The free energy potential for Van der Waals attraction between the neutral bilayers in 120 mM NaCl solutions was 1.5 X 10(-2) ergs/cm2. Also, when human serum albumin was added to the medium in the range of 0-1 mg/ml, the free energy potential for bilayer-bilayer adhesion was not affected. Using published values for equilibrium spacing between lipid bilayers in multilamellar lipid-water dispersions and the theoretical equation for van der Waals attraction between continuous dielectric layers, we calculated the value for the Hamaker coefficient of the Van der Waals attraction to be 5.8 X 10(-14) ergs.     

Binding of peripheral proteins to mixed lipid membranes: effect of lipid demixing upon binding

Binding isotherms have been determined for the association of horse heart cytochrome c with dioleoyl phosphatidylglycerol (DOPG)/dioleoyl phosphatidylcholine (DOPC) bilayer membranes over a range of lipid compositions and ionic strengths. In the absence of protein, the DOPG and DOPC lipids mix nearly ideally. The binding isotherms have been analyzed using double layer theory to account for the electrostatics, either the Van der Waals or scaled particle theory equation of state to describe the protein surface distribution, and a statistical thermodynamic formulation consistent with the mass-action law to describe the lipid distribution. Basic parameters governing the electrostatics and intrinsic binding are established from the binding to membranes composed of anionic lipid (DOPG) alone. Both the Van der Waals and scaled particle equations of state can describe the effects of protein distribution on the DOPG binding isotherms equally well, but with different values of the maximum binding stoichiometry (13 lipids/protein for Van der Waals and 8 lipids/protein for scaled particle theory). With these parameters set, it is then possible to derive the association constant, Kr, of DOPG relative to DOPC for surface association with bound cytochrome c by using the binding isotherms obtained with the mixed lipid membranes. A value of Kr (DOPG:DOPC) = 3.3-4.8, depending on the lipid stoichiometry, is determined that consistently describes the binding at different lipid compositions and different ionic strengths. Using the value of Kr obtained it is possible to derive the average in-plane lipid distribution and the enhancement in protein binding induced by lipid redistribution using the statistical thermodynamic theory.     

Water fluxes in nerve fiber

Summary Van der Waals energies of interaction between model cell surfaces are calculated for various distances of separation, layer thicknesses and compositions of cell surfaces and intercellular media. In these calculations the cell peripheries are considered to consist of two layers: (1) A phospholipid-cholesterol-protein plasma membrane and (2) a surface coat, which consists of protein, sugar and water. The required Van der Waals parameters of sugars, phospholipids and cholesterol are derived from refractive indices of their solutions in the visible and ultraviolet regions. Polarizabilities and Van der Waals parameters of these substances are determined and shown to be almost independent of concentration of solutions. Resulting isotropic polarizabilities differ by less than five percent from values obtained by the addition of bond polarizabilities.The magnitude of Van der Waals interactions between cell surfaces has been found to vary with composition according to the following sequence: water–15 ergs and 6×10^–14 ergs at 50 Å distance of separation, which corresponds to free energies per unit area of 210-1600kT/ ²     

A facile synthesis of ceramides.

Lateral diffusion of phosphatide molecules in liquid crystalline bilayers has been analysed as a case of co-operative lattice diffusion. The potential energy of interaction between two molecules is assumed to arise from Van der Waals interactions of the hydrocarbon chains, and to have the form suggested by Salem [6]. From the observed values of the self-diffusion constant (of the order of 10^?8 cm² sec^?1) the depth of the potential “well” for two molecules at the equilibrium separation was estimated to have a lower limit of 1.95 kcal per mole, and the energy barrier to lateral motion was estimated to have an upper limit of 7.21 kcal per mole.     

A new approach to calculation of the energy from van der Waals interactions in macromolecular proteins. Dielectric permeability as a physical parameter for calculations

The problem in the calculation of Van der Waals interactions in protein globules based on the theory of condensed media was considered. The Van der Waals interactions are represented as energies of interaction of regions with a uniform density distribution. A definition of the local dielectric constant as a function of coefficients of absorption of molecular groups with a particular conformation was introduced. The applicability of this approach was estimated. The theory enables one to circumvent the problems arising in calculations based on pairwise additive approximation. The methods provides a high accuracy in determining the local features of spatial structures of globular proteins and physicochemical characteristics of their constituent amino acids and molecular groups.     

High quality NMR structures: a new force field with implicit water and membrane solvation for Xplor-NIH

Structure determination of proteins by NMR is unique in its ability to measure restraints, very accurately, in environments and under conditions that closely mimic those encountered in vivo. For example, advances in solid-state NMR methods enable structure determination of membrane proteins in detergent-free lipid bilayers, and of large soluble proteins prepared by sedimentation, while parallel advances in solution NMR methods and optimization of detergent-free lipid nanodiscs are rapidly pushing the envelope of the size limit for both soluble and membrane proteins. These experimental advantages, however, are partially squandered during structure calculation, because the commonly used force fields are purely repulsive and neglect solvation, Van der Waals forces and electrostatic energy. Here we describe a new force field, and updated energy functions, for protein structure calculations with EEFx implicit solvation, electrostatics, and Van der Waals Lennard-Jones forces, in the widely used program Xplor-NIH. The new force field is based primarily on CHARMM22, facilitating calculations with a wider range of biomolecules. The new EEFx energy function has been rewritten to enable OpenMP parallelism, and optimized to enhance computation efficiency. It implements solvation, electrostatics, and Van der Waals energy terms together, thus ensuring more consistent and efficient computation of the complete nonbonded energy lists. Updates in the related python module allow detailed analysis of the interaction energies and associated parameters. The new force field and energy function work with both soluble proteins and membrane proteins, including those with cofactors or engineered tags, and are very effective in situations where there are sparse experimental restraints. Results obtained for NMR-restrained calculations with a set of five soluble proteins and five membrane proteins show that structures calculated with EEFx have significant improvements in accuracy, precision, and conformation, and that structure refinement can be obtained by short relaxation with EEFx to obtain improvements in these key metrics. These developments broaden the range of biomolecular structures that can be calculated with high fidelity from NMR restraints.     

Equilibrium energy analysis of freeze-fracture planes in membranes

We have used equilibrium energy calculations to determine the most probable freeze-fracture planes in a lipid bilayer. Using a pairwise-summation computer method, we have generated numerical values for the Van der Waals potentials (electron shell repulsion, dispersion forces and electrostatic interactions) between molecules. We have compared our theoretical predictions with the experimental conclusion that the fracture planes occur normally between lipid molecules. These calculations also provide information about the composition of intramembranous particles, the potential for local clustering of single lipid types in the fluid membrane, and the importance of lipid molecules to the function of membrane proteins such as voltage-sensitive ion channels.     

Three-dimensional structures of complexes of Lathyrus ochrus isolectin I with glucose and mannose: fine specificity of the monosaccharide-binding site

The structure of the methyl-alpha-D-mannopyranoside-LOL I complex has been solved by the molecular replacement method using the refined saccharide-free LOL I coordinates as starting model. The methyl-alpha-D-mannopyranoside-LOL I complex was refined by simulated annealing using the program X-PLOR. The final R-factor value is 0.182 [Fo greater than 1 sigma(Fo)]. The isostructural methyl-alpha-D-glucopyranoside-LOL I complex was refined by X-Ray coupled energy minimization using the methyl-alpha-D-mannopyranoside-LOL I structure as a starting model to an R factor of 0.179 (all data). In both crystal forms, each dimer binds two molecules of sugar in pockets found near the calcium ions. The two saccharide moieties, which are in the C1 chair conformation, establish the same hydrogen bond pattern with the lectin. However, the van der Waals contacts are different between the O2, C2, C6, and O6 atoms of the two molecules and the backbone atoms of residues 208-211. Mannose, due to its axial C2 conformation, encloses the backbone atoms of the protein in a clamplike way. Van der Waals energy calculations suggest that this better complementarity of the mannoside molecule with the lectin could explain its higher affinity for isolectin I.     

Phenomenon of Quantum Entanglement in a System Composed of Two Minimal Protocells

The quantum mechanical self-assembly of two separate photoactive supramolecular systems with different photosynthetic centers was investigated by means of density functional theory methods. Quantum entangled energy transitions from one subsystem to the other and the assembly of logically controlled artificial minimal protocells were modeled. The systems studied were based on different photoactive sensitizer molecules covalently bonded to a non-canonical oxo-guanine::cytosine supramolecule with the precursor of a fatty acid (pFA) molecule attached via Van der Waals forces, all surrounded by water molecules. The electron correlation interactions responsible for the weak hydrogen and Van der Waals chemical bonds increased due to the addition of polar water solvent molecules. The distances between the separated sensitizer, nucleotide, pFA, and water molecules are comparable to Van der Waals and hydrogen bonding radii. As a result, the overall system becomes compressed, resulting in photo-excited electron tunneling from the sensitizer (bis(4-diphenylamine-2-phenyl)-squarine or 1,4-bis(N,N-dimethylamino)naphthalene) to the pFA molecules. Absorption spectra as well as electron transfer trajectories associated with the different excited states were calculated using time dependent density functional theory methods. The results allow separation of the quantum entangled photosynthetic transitions within the same minimal protocell and with the neighboring minimal protocell. The transferred electron is used to cleave a “waste” organic molecule resulting in the formation of the desired product. A two variable, quantum entangled AND logic gate was proposed, consisting of two input photoactive sensitizer molecules and one output (pFA molecule). It is proposed that a similar process might be applied for the destruction of tumor cancer cells or to yield building blocks in artificial cells.     

Molecular dynamics simulations of two tandem octarepeats from the mammalian prion protein: fully Cu2+-bound and metal-free forms

Molecular dynamics simulations have been conducted on a model fragment (Ac-PHGGGWGQPHGGGW-NH₂) of the prion protein octarepeat domain, both in the Cu²⁺-bound and metal-free forms. The copper-bound models are based on the consensus structure of the core Cu²⁺-binding site of an individual octarepeat, relevant to the fully Cu²⁺-occupied prion protein octarepeat region. The model peptides contain Cu²⁺ bound through a His imidazole ring and two deprotonated amide N-atoms in the peptide backbone supplied by the following two Gly residues. Both the copper-bound and metal-free models have been simulated with the OPLS all-atom force field with the GROMACS molecular dynamics package. These simulations, with two tandem copper-binding sites, represent the minimum model necessary to observe potential structuring between the copper-binding sites in the octarepeat region. The GWGQ residues constitute a flexible linker region that predominantly adopts a turn, serving to bring adjacent His residues into close proximity. The consequent formation of stable structures demonstrates that the copper-bound octarepeat region allows the copper-coordinating sites to come into van der Waals contact, packing into particular orientations to further stabilize the bend in the GWGQ linker region.     

Semi-empirical conformational analysis of propranolol interacting with dipalmitoylphosphatidylcholine

A semi-empirical conformational analysis is used to compute the conformation of (+)-propranolol inserted in dipalmitoylphosphatidylcholine. In a first step, the minimal conformational energy of the isolated molecule at the hydrocarbon-water interface is calculated as the sum of the contributions resulting from the Van der Waals, the torsional, the electrostatic and the transfer energies. Five pairs of conformers of minimal energy are determined. They are compared to data available from other experimental approaches. In a second step, they are assembled with dipalmitoylphosphatidylcholine at the interface. Although propranolol is considered in its protonated form, the electrostatic interaction with dipalmitoylphosphatidylcholine is negligible as compared to the Van der Waals interaction. The area occupied per propranolol molecule is between 0.53 and 0.64 nm2/molecule. In the most probable modes of insertion of propranolol into the lipid layer, the naphthyl moiety of the compound interacts with the lipid acyl chains. The protonated amino group is located in the vicinity of the phosphate residue possibly causing an electrostatic interaction.     

Study on the enantioselectivity inhibition mechanism of acetyl-coenzyme A carboxylase toward haloxyfop by homology modeling and MM-PBSA analysis

Acetyl-coenzyme A carboxylase (ACCase) has been identified as one of the most important targets of herbicide Aryloxyphenoxypropionates (APPs). ACCase shows different enantioselectivity toward APPs, and only (R)-enantiomers of APPs have the herbicidal activity. In order to deeply understand the enantioselective recognition mechanism of ACCase, (R)-haloxyfop, which is a typical commercial herbicide from APPs, is selected and the relative binding free energy between ACCase and (R)-haloxyfop is investigated and compared with that between ACCase and (S)-haloxyfop by homology modeling and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method. Further free energy analysis reveals that the preference of ACCase toward (R)-haloxyfop is mainly driven by Van der Waals interaction. The analysis of the interaction between the active site residues of ACCase CT domain and (R)-haloxyfop shows the van der Waals interactions have a close relationship with the addition effect of each residue. An understanding of the enantioselective recognition mechanism between ACCase and haloxyfop is desirable to discover novel chiral herbicides.     

Effect of varying the 1–4 intramolecular scaling factor in atomistic simulations of long-chain N-alkanes with the OPLS-AA model

A comprehensive molecular dynamics simulation study of n-alkanes using the optimized potential for liquid simulation with all-atoms (OPLS-AA) force field at ambient condition has been performed. Our results indicate that while simulations with the OPLS-AA force field accurately predict the liquid state mass density for n-alkanes with carbon number equal or less than 10, for n-alkanes with carbon number equal or exceeding 12, the OPLS-AA force field with the standard scaling factor for the 1–4 intramolecular Van der Waals and electrostatic interaction gives rise to a quasi-crystalline structure. We found that accurate predictions of the liquid state properties are obtained by successively reducing the aforementioned scaling factor for each increase of the carbon number beyond n-dodecane. To better understand the effects of reducing the scaling factor, its influence on the torsion potential profile, and the corresponding gauche-trans conformer distribution, heat of vaporization, melting point, and self-diffusion coefficient for n-dodecane were investigated. This relatively simple procedure enables more accurate predictions of the thermo-physical properties of longer n-alkanes.

Bending, hydration and interstitial energies quantitatively account for the hexagonal-lamellar-hexagonal reentrant phase transition in dioleoylphosphatidylethanolamine.

We have accounted for the unusual structural change wherein dioleoylphosphatidylethanolamine undergoes a hexagonal-lamellar-hexagonal transition sequence as the water content is reduced systematically. We describe the role played by the energies of bending, hydration, voids in hexagonal interstices, and van der Waals interaction in this transition sequence. We have used the X-ray diffraction and osmotic stress experiments on the two phases to derive the structural parameters and all of the force constants defining the energetics of the hexagonal and lamellar phases. We have calculated the chemical potentials of lipid and water in both phases and derived the phase diagram of the lipid with no free, adjustable parameters. The calculated temperature/osmotic stress and temperature/composition diagrams quantitatively agree with experiment. The reentrant transition appears to be driven by a delicate balance between the hydration energy in the lamellar phase and bending energy in the hexagonal phase, whereas the energy of voids in hexagonal interstices defines its energy scale and temperature range. Van der Waals attraction between the bilayers in the lamellar phase does not appear to be important in this transition.     

Interpretation of some Conformational Data in Myoglobin and Lysozyme

THE conformations of the polypeptide chains of myoglobin¹ and lysozyme² have been successfully simulated with the aid of computed Van der Waals contact and energy maps of the theoretical independent peptide unit (IPU)^3–5. The non-glycyl experimental points plotted on an alanyl IPU are rather scattered on the allowed conformational regions of the map⁶, especially in the case of lysozyme. By contrast, well defined clusters of points can be observed when only the amino-acid residues in segments of the helical secondary structure (mainly α and β chains) are plotted. In addition, clusters of points, albeit less well defined, can be observed by plotting the points relative to the experimental conformations of the first non-helical amino-acid residue next to a more or less folded segment of that α-helical type so frequently present in globular proteins (Fig. 1).     

Cluster analysis identifies aminoacid compositional features that indicate Toxoplasma gondii adhesin proteins

Toxoplasma gondii invade host cells using a multi-step process that depends on the regulated secretion of adhesions. To identify key primary sequence features of adhesins in this parasite, we analyze the relative frequency of individual amino acids, their dipeptide frequencies, and the polarity, polarizability and Van der Waals volume of the individual amino acids by using cluster analysis. This method identified cysteine as a key amino acid in the Toxoplasma adhesin group. The best vector algorithm of non-concatenated features was for 2 attributes: the single amino acid relative frequency and the dipeptide frequency. Polarity, polarizability and Van der Waals volume were not good classificatory attributes. Single amino acid attributes clustered unambiguously 67 apicomplexan hypothetical adhesins. This algorithm was also useful for clustering hypothetical Toxoplasma target host receptors. All of the cluster performances had over 70% sensitivity and 80% specificity. Compositional aminoacid data can be useful for improving machine learning-based prediction software when homology and structural data are not sufficient.