Structural stability and buckling analysis of a series of carbon nanotorus using molecular dynamics simulations

Based on molecular dynamics (MD) simulations, the buckling analysis of a perfect carbon nanotorus is presented herein. First of all, the minimum length of single-walled carbon nanotubes (SWCNTs) with different radii is determined at which perfect toroidal CNTs can be formed without any ripple at the inner side of the rings. According to the results, by increasing the radius of SWCNT (r), the radius of its corresponding perfect nanotorus (R) increases. Also, for SWCNTs with various lengths, it is found that the buckling force and strain of related carbon nanotoruses increase by increasing R/r. In addition, as the perfect toroidal CNTs are arranged vertically in a column form in accordance with two different schemes, the effects of increasing the radius (R) and the number of carbon nanotoruses (the height of the column made by nanotoruses) on the buckling force and strain are investigated. Based on the results, as a fixed number of carbon nanotoruses with the same radius are arranged vertically in the column form, the buckling force and strain increase by increasing R/r. By contrast, increasing the height of the column made by carbon nanotoruses with similar radius leads to the reduction of buckling force and strain.     

Structural analysis of telechelic polymer solution using dissipative particle dynamics simulations

A telechelic polymer is an amphiphilic polymer that can form micellar structures when dissolved in water. A telechelic polymer solution shows viscoelastic behaviour owing to the formation of characteristic networks, i.e. loops, bridges and dangling chains. For industrial purposes, telechelic polymers have many applications as thickening agents, such as in paints and cosmetics. Thus, it is desirable to predict and control the rheological properties of telechelic polymers. However, detailed studies at the molecular level have not yet been performed. In this study, I use the dissipative particle dynamics (DPD) method to investigate the relationship between the characteristic structural properties and the molecular structure in telechelic polymer solutions. I show that the morphology of telechelic polymer solutions depends on the concentration and chain length, the distribution of the end-to-end distance, the mean square end-to-end distance, the mean square radius of gyration and the time-averaged mean square displacement. Although an effect of entanglement is important for properties of polymer melts, the polymer chain composed of DPD particles cannot reproduce it. Therefore, I compare telechelic polymer solutions with and without the segmental repulsive potential (SRP), which can simulate the effect of entanglement in DPD simulations. The results indicate that it is necessary to include the SRP in DPD simulations to correctly analyse the behaviour of telechelic polymer solutions.     

Eutectic formation analysis of amino acid mixtures using molecular dynamics simulations

The mechanism of eutectic formation was investigated via computer-aided molecular dynamics techniques based on experimental results. The CBZ group mixtures CBZ-l-Asp/d-AlaNH2 x HCl/methanol, CBZ-l-Asp/l-PheOMe x HCl/methanol, and CBZ-l-Tyr/l-ArgNH2 x 2HCl/methanol formed transparent eutectic melts. The non-CBZ group mixtures l-Asp/d-AlaNH2 x HCl/methanol, l-Asp/l-PheOMe x HCl/methanol, and l-Tyr/l-ArgNH2 x 2HCl/methanol did not form eutectic melts. According to molecular dynamics simulation results, increase in the kinetic energy values of eutectic forming mixtures was much larger than the increase in potential energy over a temperature shift from 298 to 333 K. However, the results for non-eutectic forming mixtures were reversed. The Coulomb interaction energies of eutectic forming mixtures significantly decreased, because eutectic melting can increase the mobility of molecules in the mixtures. The enhancement of molecular mobility was confirmed by increased self-diffusion constant values, and the change of solid-to-liquid phase was detected by radial distribution function results. The periodic boundary conditions for calculation of molecular dynamics were found to be reliable.     

Structural determinants of MscL gating studied by molecular dynamics simulations

The mechanosensitive channel of large conductance (MscL) in prokaryotes plays a crucial role in exocytosis as well as in the response to osmotic downshock. The channel can be gated by tension in the membrane bilayer. The determination of functionally important residues in MscL, patch-clamp studies of pressure-conductance relationships, and the recently elucidated crystal structure of MscL from Mycobacterium tuberculosis have guided the search for the mechanism of MscL gating. Here, we present a molecular dynamics study of the MscL protein embedded in a fully hydrated POPC bilayer. Simulations totaling 3 ns in length were carried out under conditions of constant temperature and pressure using periodic boundary conditions and full electrostatics. The protein remained in the closed state corresponding to the crystal structure, as evidenced by its impermeability to water. Analysis of equilibrium fluctuations showed that the protein was least mobile in the narrowest part of the channel. The gating process was investigated through simulations of the bare protein under conditions of constant surface tension. Under a range of conditions, the transmembrane helices flattened as the pore widened. Implications for the gating mechanism in light of these and experimental results are discussed.     

Sequence-dependent DNA deformability studied using molecular dynamics simulations

Proteins recognize specific DNA sequences not only through direct contact between amino acids and bases, but also indirectly based on the sequence-dependent conformation and deformability of the DNA (indirect readout). We used molecular dynamics simulations to analyze the sequence-dependent DNA conformations of all 136 possible tetrameric sequences sandwiched between CGCG sequences. The deformability of dimeric steps obtained by the simulations is consistent with that by the crystal structures. The simulation results further showed that the conformation and deformability of the tetramers can highly depend on the flanking base pairs. The conformations of xATx tetramers show the most rigidity and are not affected by the flanking base pairs and the xYRx show by contrast the greatest flexibility and change their conformations depending on the base pairs at both ends, suggesting tetramers with the same central dimer can show different deformabilities. These results suggest that analysis of dimeric steps alone may overlook some conformational features of DNA and provide insight into the mechanism of indirect readout during protein–DNA recognition. Moreover, the sequence dependence of DNA conformation and deformability may be used to estimate the contribution of indirect readout to the specificity of protein–DNA recognition as well as nucleosome positioning and large-scale behavior of nucleic acids.     

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.     

Conformational analysis and molecular dynamics simulations of maltose

Constrained conformational energy minimizations have been used to calculate an adiabatic (Φ, ψ) potential energy surface for the disaccharide β-maltose. The inclusion of molecular flexibility in the conformational energy analysis of the disaccharide was found to significantly lower the barriers to conformational transitions, as has been observed previously for other systems. Several low energy wells were identified on the adiabatic surface which differ in energy by small amounts and with low absolute barriers separating them, indicating the possibility of a non-negligible equilibrium population distribution in each well. If such a distribution of conformations existed in the physical system, the conformation observed by NMR NOE measurements would thus be a “virtual” conformation. Molecular dynamics simulations of the motions of this molecule in vacuum were also conducted and indicate that the rate of relaxation of the molecule to the adiabatic surface may be slower than the typical timescale of conformational fluctuations. This effect is apparently due to an unphysical persistence of hydrogen bond patterns in vacuum which does not occur in aqueous solution. Trajectories undergoing transitions between wells were calculated and the effects of such conformational transitions upon the ensemble mean structure, such as might be observed in an NMR experiment, were demonstrated.     

Structural reorganization of molecular sheets derived from cellulose II by molecular dynamics simulations

We investigated structural reorganization of two different kinds of molecular sheets derived from the cellulose II crystal using molecular dynamics (MD) simulations, in order to identify the initial structure of the cellulose crystal in the course of its regeneration process from solution. After a one-nanosecond simulation, the molecular sheet formed by van der Waals forces along the () crystal plane did not change its structure in an aqueous environment, while the other one formed by hydrogen bonds along the (1 1 0) crystal plane changed into a van der Waals-associated molecular sheet, such as the former. The two structures that were calculated showed substantial similarities such as the high occupancy of intramolecular hydrogen bonds between O3^H and O5 of over 0.75, few intermolecular hydrogen bonds, and the high occurrence of hydrogen bonding with water. The convergence of the two structures into one denotes that the van der Waals-associated molecular sheet can be the initial structure of the cellulose crystal formed in solution. The main chain conformations were almost the same as those in the cellulose II crystal except for a −16° shift of φ (dihedral angle of O5-C1-O1-C4) and the gauche-gauche conformation of the hydroxymethyl side group appears probably due to its hydrogen bonding with water. These results suggest that the van der Waals-associated molecular sheet becomes stable in an aqueous environment with its hydrophobic inside and hydrophilic periphery. Contrary to this, a benzene environment preferred a hydrogen-bonded molecular sheet, which is expected to be the initial structure formed in benzene.     

Non-B DNA conformations analysis through molecular dynamics simulations

BackgroundNon-B DNA conformations are molecular structures that do not follow the canonical DNA double helix. Mutagenetic instability in nuclear and mitochondrial DNA (mtDNA) genomes has been associated with simple non-B DNA conformations, as hairpins or more complex structures, as G-quadruplexes. One of these structures is Structure A, a cloverleaf-like non-B conformation predicted for a 93-nt (nucleotide) stretch of the mtDNA control region 5′-peripheral domain. Structure A is embedded in a hot spot for the 3′ end of human mtDNA deletions revealing its importance in influencing the mutational instability of the mtDNA genome.MethodsTo better characterize Structure A, we predicted its 3D conformation using state-of-art methods and algorithms. The methodologic workflow consisted in the prediction of non-B conformations using molecular dynamics simulations. The conservation scores of alignments of the Structure A region in humans, primates, and mammals, was also calculated.ResultsOur results show that these computational methods are able to measure the stability of non-B conformations by using the level of base pairing during molecular dynamics. Structure A showed high stability and low flexibility correlated with high conservation scores in mammalian, more specifically in primate lineages.ConclusionsWe showed that 3D non-B conformations can be predicted and characterized by our methodology. This allowed the in-depth analysis of the structure A, and the main results showed the structure remains stable during the simulations.General significanceThe fine-scale atomic molecular determination of this type of non-B conformation opens the way to perform computational molecular studies that can show their involvement in mtDNA cellular mechanisms.     

Characterization of molecular structures and properties of polyurethanes using molecular dynamics simulations

Thermotropic polyurethanes with mesogenic groups in side chains were prepared from two diisocyanates and four diols with stoichiometric ratios of reactive isocyanate (NCO) and hydroxy (OH) groups. Their thermal behavior was determined by differential scanning calorimetry. The effect of structure modifications of the diisocyanates and diols, in particular changes in the mesogen, were investigated. Introduction of mesogenic segments into the polymers suppresses the ordering. Stiff end substituents (phenyl and alkoxy groups) of the mesogens stabilize the mesophases to such an extent that the negative influence of long polymer chains is compensated and the liquid-crystalline properties are recovered. All-atom molecular dynamics simulations in the Cerius² modeling environment were carried out to characterize the structures of the polymers. Analysis of the dynamic trajectories at 20, 100, 120 and 170 °C revealed changes in conformation of macromolecules, which correlate with DSC measurements.Figure Example of structure relaxation of D4/TDI molecule at indicated simulation times (temperature 20 °C): a complete structure; b backbone structure; c top view of molecule     

Structural analysis of substance P using molecular dynamics and NMR spectroscopy.

The present work is a combined structural study, using Nuclear Magnetic Resonance (NMR) and Molecular Dynamics(MD), of the amidated and the free acid forms of substance P in water and methanol. The results obtained using both approaches were compared in order to characterize the structural features of both peptides in solution. From the NMR experiments it was derived that the free acid form adopts an extended conformation at the N-terminus and a helical conformation at the C-terminal segment of the peptide in both water and methanol; these structural features are in qualitative agreement with the results of the MD simulations. No significant differences in behavior were observed between the amidated and the free acid forms of the peptide in the simulations and in the experiments carried out in water, suggesting that the different activities of these analogs are due to their different mode of interaction with the receptor rather than to their structural preferences. Finally, we propose that the structure of substance P can be partially inferred from its sequence due to the presence of a Pro-X-Pro motif on the N-terminus and a Gly-Leu sequence on the C-terminus.     

Structural insight into antisense gapmer-RNA oligomer duplexes through molecular dynamics simulations

There is an extensive research carrying out on antisense technology and the molecules entering into clinical trials are increasing rapidly. Phosphorothioate (PS) is a chemical modification in which nonbridged oxygen is replaced with a sulfur, consequently providing resistance against nuclease activity. The 2'-4' conformationally restricted nucleoside has the structural features of both 2'-O-methoxy ethyl RNA (MOE), which shows good toxicity profile, and locked nucleic acid (LNA), which shows good binding affinity towards the target RNA. These modifications have been studied and suggested that they can be a potential therapeutic agents in antisense therapy. Mipomersen (ISIS 301012), which contains the novel nucleoside modification has been used to target to apolipoprotein (Apo B), which reduces LDL cholesterol by 6–41%. In this study, classical molecular dynamics (MD) simulations were performed on six different antisense gapmer/target-RNA oligomer duplexes (LNA-PS-LNA/RNA, RcMOE-PS-RcMOE/RNA, ScMOE-PS-ScMOE/RNA, MOE-PS-MOE/RNA, PS-DNA/RNA and DNA/RNA) to investigate the structural dynamics, stability and solvation properties. The LNA, MOE nucleotides present in respective duplexes are showing the structure of A-form and the PS-DNA nucleotides resemble the structure of B-form helix with respect to some of the helical parameters. Free energy calculations suggest that the oligomer, which contains LNA binds to the RNA strongly than other modifications as shown in experimental results. The MOE modified nucleotide, which although had a lower binding affinity but higher solvent accessible surface area (SASA) compared to the other modifications, may be influencing the toxicity and hence may be used it in Mipomersen, the second antisense molecule which is approved by FDA.

Communicated by Ramaswamy H. Sarma     

Implementation of force distribution analysis for molecular dynamics simulations

Background  

The way mechanical stress is distributed inside and propagated by proteins and other biopolymers largely defines their function. Yet, determining the network of interactions propagating internal strain remains a challenge for both, experiment and theory. Based on molecular dynamics simulations, we developed force distribution analysis (FDA), a method that allows visualizing strain propagation in macromolecules.     

Force field impact and spin-probe modeling in molecular dynamics simulations of spin-labeled T4 lysozyme

Several attempts have been made to compute electron paramagnetic resonance (EPR) spectra of biomolecules, using motional models or simulated trajectories to describe dynamics. Ideally, the simulated trajectories should capture fast (picosecond) snapshots of spin-probe rotations accurately, while being lengthy enough to ensure a proper Fourier integration of the time-domain signal. It is the interplay of the two criteria that poses computational challenges to the method. In this context, an analysis of the spin-probe and protein conformational sampling and equilibration, with different force fields and with explicit solvent, may be a useful attempt. The present work reports a comparative study of the effect of the molecular dynamics (MD) force field on conformational sampling and equilibration in two spin-labeled T4 lysozyme (T4L) variants, N40C and K48C. Ensembles of 10× 3 ns-trajectories per variant and per force field (OPLS/AMBER and AMBER99) are analyzed for a reliable assessment of convergence and sampling. It is found that subtle site-dependent differences in spin-probe rotations and torsions are more readily captured in the AMBER99 trajectories than in the OPLS/AMBER simulations. On the other hand, sampling and equilibration are found to be better with the OPLS/AMBER force field at equal trajectory lengths.Figure: Left panel: The spin-probe R1 ring and the spin-probe Euler angles , and . Middle panel: Illustration of the diffusion in a cone model for the spin-probe motion: snapshots of helix B and of the R1 ring in N40C, taken at 0.3 ps intervals from AMBER trajectory 1. Right panel: The N40C mutant with the spin label (solid mode), solvated in a cubic box.     

Atomistic modeling and molecular dynamics analysis of human CaV1.2 channel using external electric field and ion pulling simulations

BackgroundHuman Ca_V1.2 (hCav1.2), a calcium selective voltage-gated channel, plays important roles in normal cardiac and neuronal functions. Calcium influx and gating mechanisms leading to the activation of hCa_V1.2 are critical for its functionalities. Lack of an experimentally resolved structure of hCa_V1.2 remains a significant impediment in molecular-level understanding of this channel. This work focuses on building atomistic hCa_V1.2 model and studying calcium influx using computational approaches.MethodsWe employed homology modeling and molecular dynamics (MD) to build the structure of hCa_V1.2. Subsequently, we employed steered molecular dynamics (SMD) to understand calcium ion permeation in hCa_V1.2.ResultsWe report a comprehensive three-dimensional model of a closed state hCa_V1.2 refined under physiological membrane-bound conditions using MD simulations. Our SMD simulations on the model revealed four important barriers for ion permeation: this includes three calcium binding sites formed by the EEEE- and TTTT- rings within the selectivity filter region and a large barrier rendered by the hydrophobic internal gate. Our results also revealed that the first hydration shell of calcium remained intact throughout the simulations, thus playing an important role in ion permeation in hCa_V1.2.ConclusionsOur results have provided some important mechanistic insights into the structure, dynamics and ion permeation in hCa_V1.2. The significant barriers for ion permeation formed by the four phenylalanine residues at the internal gate region suggest that this site is important for channel activation.     

Exploration of conformational transition in the aryl-binding site of human FXa using molecular dynamics simulations

Human coagulation Factor X (FX), a member of the vitamin K-dependent serine protease family, is a crucial component of the human coagulation cascade. Activated FX (FXa) participates in forming the prothrombinase complex on activated platelets to convert prothrombin to thrombin in coagulation reactions. In the current study, 30-ns MD simulations were performed on both the open and closed states of human FXa. Root mean squares (RMS) fluctuations showed that structural fluctuations concentrated on the loop regions of FXa, and the presence of a ligand in the closed system resulted in larger fluctuations of the gating residues. The open system had a gating distance from 9.23 to 11.33 ?, i.e., significantly larger than that of the closed system (4.69-6.35 ?), which allows diversified substrates of variable size to enter. Although the solvent accessible surface areas (SASA) of FXa remained the same in both systems, the open system generally had a larger total SASA or hydrophobic SASA (or both) for residues surrounding the S4 pocket. Additionally, more hydrogen bonds were formed in the closed state than in the open state of FXa, which is believed to play a significant role in maintaining the closed confirmation of the aryl-binding site. Based on the results of MD simulations, we propose that an induced-fit mechanism governs the functioning of human coagulation FX, which helps provide a better understanding of the interactions between FXa and its substrate, and the mechanism of the conformational changes involved in human coagulation.     

Analysis of ligand binding to proteins using molecular dynamics simulations

This work aims to explore theoretically the molecular mechanisms of ligand binding to proteins through the use of molecular dynamics simulations. The binding of sodium dodecyl sulfate (SDS) to cobra cardio toxin A3 (CTX A3) and thiourea (TOU) to lysozyme have been chosen as the two model systems. Data acquisitions were made by Gromacs software. To begin with, the collisions of ligand molecules with every residue of CTX A3 and lysozyme were evaluated. With this information in hand, the average numbers of collisions with each residue was defined and then assessed. Next, a measure of the affinity of a residue, P_i, referred to as conformational factor, toward a ligand molecule was established. Based on the results provided, all site-making residues for CTX A3 and lysozyme were identified. The results are in good agreement with the experimental data. Finally, based on this method, all site-making residues of bovine carbonic anhydrase (BCA) toward the SDS ligand were predicted.     

Structural insights into human microsomal epoxide hydrolase by combined homology modeling,molecular dynamics simulations,and molecular docking calculations

A new homology model of human microsomal epoxide hydrolase was derived based on multiple templates. The model obtained was fully evaluated, including MD simulations and ensemble‐based docking, showing that the quality of the structure is better than that of only previously known model. Particularly, a catalytic triad was clearly identified, in agreement with the experimental information available. Analysis of intermediates in the enzymatic mechanism led to the identification of key residues for substrate binding, stereoselectivity, and intermediate stabilization during the reaction. In particular, we have confirmed the role of the oxyanion hole and the conserved motif (HGXP) in epoxide hydrolases, in excellent agreement with known experimental and computational data on similar systems. The model obtained is the first one that fully agrees with all the experimental observations on the system. Proteins 2017; 85:720–730. © 2016 Wiley Periodicals, Inc.