Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as similar to 600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.

Electronic band Structure and Photoemission Spectra of Graphene on Silicon Substrate

Synergizing graphene on silicon based nanostructures is pivotal in advancing nano-electronic device technology. A combination of molecular dynamics and density functional theory has been used to predict the electronic energy band structure and photo-emission spectrum for graphene-Si system with silicon as a substrate for graphene. The equilibrium geometry of the system after energy minimization is obtained from molecular dynamics simulations. For the stable geometry obtained, density functional theory calculations are employed to determine the energy band structure and dielectric constant of the system. Further the work function of the system which is a direct consequence of photoemission spectrum is calculated from the energy band structure using random phase approximations.

Unzipping Carbon Nanotubes at High Impact

The way nanostructures behave and mechanically respond to high impact collision is a topic of intrigue. For anisotropic nanostructures, such as carbon nanotubes, this response will be complicated based on the impact geometry. Here we report the result of hypervelocity impact of nanotubes against solid targets and show that impact produces a large number of defects in the nanotubes, as well as rapid atom evaporation, leading to their unzipping along the nanotube axis. Fully atomistic reactive molecular dynamics simulations are used to gain further insights of the pathways and deformation and fracture mechanisms of nanotubes under high energy mechanical impact. Carbon nanotubes have been unzipped into graphene nanoribbons before using chemical treatments but here the instability of nanotubes against formation, fracture, and unzipping is revealed purely through mechanical impact. defect

Defect induced plasticity and failure mechanism of boron nitride nanotubes under tension

The effects of Stone-Wales (SW) and vacancy defects on the failure behavior of boron nitride nanotubes (BNNTs) under tension are investigated using molecular dynamics simulations. The Tersoff-Brenner potential is used to model the atomic interaction and the temperature is maintained close to 300 K. The effect of a SW defect is studied by determining the failure strength and failure mechanism of nanotubes with different radii. In the case of a vacancy defect, the effect of an N-vacancy and a B-vacancy is studied separately. Nanotubes with different chiralities but similar diameter is considered first to evaluate the chirality dependence. The variation of failure strength with the radius is then studied by considering nanotubes of different diameters but same chirality. It is observed that the armchair BNNTs are extremely sensitive to defects, whereas the zigzag configurations are the least sensitive. In the case of pristine BNNTs, both armchair and zigzag nanotubes undergo brittle failure, whereas in the case of defective BNNTs, only the zigzag ones undergo brittle failure. An interesting defect induced plastic behavior is observed in defective armchair BNNTs. For this nanotube, the presence of a defect triggers mechanical relaxation by bond breaking along the closest zigzag helical path, with the defect as the nucleus. This mechanism results in a plastic failure. (C) 2014 AIP Publishing LLC.

Design and Synthesis of New Benzimidazole-Carbazole Conjugates for the Stabilization of Human Telomeric DNA, Telomerase Inhibition, and Their Selective Action on Cancer Cells

Cell-permeable small molecules that enhance the stability of the G-quadruplex (G4) DNA structures are currently among the most intensively pursued ligands for inhibition of the telomerase activity. Herein we report the design and syntheses of four novel benzimidazole carbazole conjugates and demonstrate their high binding affinity to G4 DNA. Si nuclease assay confirmed the ligand mediated G-quadruplex DNA protection. Additional evidence from Telomeric Repeat Amplification Protocol (TRAP-LIG) assay demonstrated efficient telomerase inhibition activity by the ligands. Two of the ligands showed IC50 values in the sub-micromolar range in the TRAP-LIG assay, which are the best among the benzimidazole derivatives reported so far. The ligands also exhibited cancer cell selective nuclear internalization, nuclear condensation, fragmentation, and eventually antiproliferative activity in long-term cell viability assays. Annexin V-FITC/PI staining assays confirm that the cell death induced by the ligands follows an apoptotic pathway. An insight into the mode of ligand binding was obtained from the molecular dynamics simulations.

A new microscopic insight into membrane penetration and reorganization by PETIM dendrimers

Dendrimers are highly branched polymeric nanoparticles whose structure and topology, largely, have determined their efficacy in a wide range of studies performed so far. An area of immense interest is their potential as drug and gene delivery vectors. Realizing this potential, depending on the nature of cell surface-dendrimer interactions, here we report controlled model membrane penetration and reorganization, using a model supported lipid bilayer and poly(ether imine) (PETIM) dendrimers of two generations. By systematically varying the areal density of the lipid bilayers, we provide a microscopic insight, through a combination of high resolution scattering, atomic force microscopy and atomistic molecular dynamics simulations, into the mechanism of PETIM dendrimer membrane penetration, pore formation and membrane re-organization induced by such interactions. Our work represents the first systematic observation of a regular barrel-like membrane spanning pore formation by dendrimers, tunable through lipid bilayer packing, without membrane disruption.

Laterally structured ripple and square phases with one and two dimensional thickness modulations in a model bilayer system

Molecular dynamics simulations of bilayers in a surfactant/co-surfactant/water system with explicit solvent molecules show formation of topologically distinct gel phases depending upon the bilayer composition. At low temperatures, the bilayers transform from the tilted gel phase, L beta', to the one dimensional (1D) rippled, P beta' phase as the surfactant concentration is increased. More interestingly, we observe a two dimensional (2D) square phase at higher surfactant concentration which, upon heating, transforms to the gel L beta' phase. The thickness modulations in the 1D rippled and square phases are asymmetric in two surfactant leaflets and the bilayer thickness varies by a factor of similar to 2 between maximum and minimum. The 1D ripple consists of a thinner interdigitated region of smaller extent alternating with a thicker non-interdigitated region. The 2D ripple phase is made up of two superimposed square lattices of maximum and minimum thicknesses with molecules of high tilt forming a square lattice translated from the lattice formed with the thickness minima. Using Voronoi diagrams we analyze the intricate interplay between the area-per-head-group, height modulations and chain tilt for the different ripple symmetries. Our simulations indicate that composition plays an important role in controlling the formation of low temperature gel phase symmetries and rippling accommodates the increased area-per-head-group of the surfactant molecules.

Local Structural and Environmental Factors Define the Efficiency of an RNA Pseudoknot Involved in Programmed Ribosomal Frameshift Process

In programmed -1 ribosomal frameshift, an RNA pseudoknot stalls the ribosome at specific sequence and restarts translation in a new reading frame. A precise understanding of structural characteristics of these pseudoknots and their PRF inducing ability has not been clear to date. To investigate this phenomenon, we have studied various structural aspects of a -1 PRF inducing RNA pseudoknot from BWYV using extensive molecular dynamics simulations. A set of functional and poorly functional forms, for which previous mutational data were available, were chosen for analysis. These structures differ from each other by either single base substitutions or base-pair replacements from the native structure. We have rationalized how certain mutations in RNA pseudoknot affect its function; e.g., a specific base substitution in loop 2 stabilizes the junction geometry by forming multiple noncanonical hydrogen bonds, leading to a highly rigid structure that could effectively resist ribosome-induced unfolding, thereby increasing efficiency. While, a CG to AU pair substitution in stem 1 leads to loss of noncanonical hydrogen bonds between stems and loop, resulting in a less stable structure and reduced PRF inducing ability, inversion of a pair in stem 2 alters specific base-pair geometry that might be required in ribosomal recognition of nucleobase groups, negatively affecting pseudoknot functioning. These observations illustrate that the ability of an RNA pseudoknot to induce -1 PRF with an optimal rate depends on several independent factors that contribute to either the local conformational variability or geometry

Nature of the effective interaction between dendrimers

We have performed fully atomistic classical molecular dynamics simulations to calculate the effective interaction between two polyamidoamine dendrimers. Using the umbrella sampling technique, we have obtained the potential of mean force (PMF) between the dendrimers and investigated the effects of protonation level and dendrimer size on the PMF. Our results show that the interaction between the dendrimers can be tuned from purely repulsive to partly attractive by changing the protonation level. The PMF profiles are well-fitted by the sum of an exponential and a Gaussian function with the weight of the exponential function dominating over that of the Gaussian function. This observation is in disagreement with the results obtained in previous analytic C. Likos, M. Schmidt, H. Lowen, M. Ballauff, D. Potschke, and P. Lindner, Macromolecules 34, 2914 (2001)] and coarse-grained simulation I. Gotze, H. Harreis, and C. Likos, J. Chem. Phys. 120, 7761 (2004)] studies which predicted the effective interaction to be Gaussian. (C) 2014 AIP Publishing LLC.

Structure, stability and elasticity of DNA nanotubes

DNA nanotubes are tubular structures composed of DNA crossover molecules. We present a bottom up approach for the construction and characterization of these structures. Various possible topologies of nanotubes are constructed such as 6-helix, 8-helix and tri-tubes with different sequences and lengths. We have used fully atomistic molecular dynamics simulations to study the structure, stability and elasticity of these structures. Several nanosecond long MD simulations give the microscopic details about DNA nanotubes. Based on the structural analysis of simulation data, we show that 6-helix nanotubes are stable and maintain their tubular structure; while 8-helix nanotubes are flattened to stabilize themselves. We also comment on the sequence dependence and the effect of overhangs. These structures are approximately four times more rigid having a stretch modulus of similar to 4000 pN compared to the stretch modulus of 1000 pN of a DNA double helix molecule of the same length and sequence. The stretch moduli of these nanotubes are also three times larger than those of PX/JX crossover DNA molecules which have stretch moduli in the range of 1500-2000 pN. The calculated persistence length is in the range of a few microns which is close to the reported experimental results on certain classes of DNA nanotubes.

Ligand 5,10,15,20-Tetra(N-methyl-4-pyridyl)porphine (TMPyP4) Prefers the Parallel Propeller-Type Human Telomeric G-Quadruplex DNA over Its Other Polymorphs

The binding of ligand 5,10,15,20-tetra(N-methyl-4-pyridyl)porphine (TMPyP4) with telomeric and genomic G-quadruplex DNA has been extensively studied. However, a comparative study of interactions of TMPyP4 with different conformations of human telomeric G-quadruplex DNA, namely, parallel propeller-type (PP), antiparallel basket-type (AB), and mixed hybrid-type (MH) G-quadruplex DNA, has not been done. We considered all the possible binding sites in each of the G-quadruplex DNA structures and docked TMPyP4 to each one of them. The resultant most potent sites for binding were analyzed from the mean binding free energy of the complexes. Molecular dynamics simulations were then carried out, and analysis of the binding free energy of the TMPyP4-G-quadruplex complex showed that the binding of TMPyP4 with parallel propeller-type G-quadruplex DNA is preferred over the other two G-quadruplex DNA conformations. The results obtained from the change in solvent excluded surface area (SESA) and solvent accessible surface area (SASA) also support the more pronounced binding of the ligand with the parallel propeller-type G-quadruplex DNA.

Driving force of water entry into hydrophobic channels of carbon nanotubes: entropy or energy?

Spontaneous entry of water molecules inside single-wall carbon nanotubes (SWCNTs) has been confirmed by both simulations and experiments. Using molecular dynamics simulations, we have studied the thermodynamics of filling of a (6,6) carbon nanotube in a temperature range from 273 to 353K and with different strengths of the nanotube-water interaction. From explicit energy and entropy calculations using the two-phase thermodynamics method, we have presented a thermodynamic understanding of the filling behaviour of a nanotube. We show that both the energy and the entropy of transfer decrease with increasing temperature. On the other hand, scaling down the attractive part of the carbon-oxygen interaction results in increased energy of transfer while the entropy of transfer increases slowly with decreasing the interaction strength. Our results indicate that both energy and entropy favour water entry into (6,6) SWCNTs. Our results are compared with those of several recent studies of water entry into carbon nanotubes.

Molecular effects of encapsulation of glucose oxidase dimer by graphene

Knowing the nature of the enzyme-graphene interface is critical for a design of graphene-based biosensors. Extensive contacts between graphene and enzyme could be obtained by employing a suitable encapsulation which does not impede its enzymatic reaction. We have performed molecular dynamics simulations to obtain an insight on many forms of contact between glucose oxidase dimer and the single-layer graphene nano-sheets. The unconnected graphene sheets tended to form a flat stack regardless of their initial positions around the enzyme, whereas the same graphene sheets linked together formed a flower-like shape engendering different forms of wrapping of the enzyme. During the encapsulation no core hydrophobic residues of the enzyme were exposed. Since the polar and charged amino acids populated the enzyme's surface we also estimated, using DFT calculations, the interaction energies of individual polar and charged amino acid residues with graphene. It was found that the negatively charged residues can bind to graphene unexpectedly strongly; however, the main effect of encapsulation comes from the overlap of adjacent edges of graphene sheets.

Spatial modulation of the composition of a binary liquid near a repulsive wall

When a binary liquid is confined by a strongly repulsive wall, the local density is depleted near the wall and an interface similar to that between the liquid and its vapor is formed. This analogy suggests that the composition of the binary liquid near this interface should exhibit spatial modulation similar to that near a liquid-vapor interface even if the interactions of the wall with the two components of the liquid are the same. The Guggenheim adsorption relation quantifies the concentrations of two components of a binary mixture near a liquid-vapor interface and qualitatively states that the majority (minority) component enriches the interface for negative (positive) mixing energy if the surface tensions of the two components are not very different. From molecular dynamics simulations of binary mixtures with different compositions and interactions we find that the Guggenheim relation is qualitatively satisfied at wall-induced interfaces for systems with negative mixing energy at all state points considered. For systems with positive mixing energy, this relation is found to be qualitatively valid at low densities, while it is violated at state points with high density where correlations in the liquid are strong. This observation is validated by a calculation of the density profiles of the two components of the mixture using density functional theory with the Ramakrishnan-Yussouff free-energy functional. Possible reasons for the violation of the Guggenheim relation are discussed.

Composition dependent multiple structural transformations of myoglobin in aqueous ethanol solution: A combined experimental and theoretical study

Experimental studies (circular dichroism and ultra-violet (UV) absorption spectra) and large scale atomistic molecular dynamics simulations (accompanied by order parameter analyses) are combined to establish a number of remarkable (and unforeseen) structural transformations of protein myoglobin in aqueous ethanol mixture at various ethanol concentrations. The following results are particularly striking. (1) Two well-defined structural regimes, one at x(EtOH) similar to 0.05 and the other at x(EtOH) similar to 0.25, characterized by formation of distinct partially folded conformations and separated by a unique partially unfolded intermediate state at x(EtOH) similar to 0.15, are identified. (2) Existence of non-monotonic composition dependence of (i) radius of gyration, (ii) long range contact order, (iii) residue specific solvent accessible surface area of tryptophan, and (iv) circular dichroism spectra and UV-absorption peaks are observed. Interestingly at x(EtOH) similar to 0.15, time averaged value of the contact order parameter of the protein reaches a minimum, implying that this conformational state can be identified as a molten globule state. Multiple structural transformations well known in water-ethanol binary mixture appear to have considerably stronger effects on conformation and dynamics of the protein. We compare the present results with studies in water-dimethyl sulfoxide mixture where also distinct structural transformations are observed along with variation of co-solvent composition. (C) 2015 AIP Publishing LLC.