Molecular dynamics studies on the buffalo prion protein

It was reported that buffalo is a low susceptibility species resisting to transmissible spongiform encephalopathies (TSEs) (same as rabbits, horses, and dogs). TSEs, also called prion diseases, are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of species (except for rabbits, dogs, horses, and buffalo), manifesting as scrapie in sheep and goats; bovine spongiform encephalopathy (BSE or “mad–cow” disease) in cattle; chronic wasting disease in deer and elk; and Creutzfeldt–Jakob diseases, Gerstmann–Sträussler–Scheinker syndrome, fatal familial insomnia, and Kulu in humans etc. In molecular structures, these neurodegenerative diseases are caused by the conversion from a soluble normal cellular prion protein (PrP^C), predominantly with α-helices, into insoluble abnormally folded infectious prions (PrP^Sc), rich in β-sheets. In this article, we studied the molecular structure and structural dynamics of buffalo PrP^C (BufPrP^C), in order to understand the reason why buffalo is resistant to prion diseases. We first did molecular modeling of a homology structure constructed by one mutation at residue 143 from the NMR structure of bovine and cattle PrP(124–227); immediately we found that for BufPrP^C(124–227), there are five hydrogen bonds (HBs) at Asn143, but at this position, bovine/cattle do not have such HBs. Same as that of rabbits, dogs, or horses, our molecular dynamics studies also revealed there is a strong salt bridge (SB) ASP178–ARG164 (O–N) keeping the β2–α2 loop linked in buffalo. We also found there is a very strong HB SER170–TYR218 linking this loop with the C-terminal end of α-helix H3. Other information, such as (i) there is a very strong SB HIS187–ARG156 (N–O) linking α-helices H2 and H1 (if mutation H187R is made at position 187, then the hydrophobic core of PrP^C will be exposed (L.H. Zhong (2010). Exposure of hydrophobic core in human prion protein pathogenic mutant H187R. Journal of Biomolecular Structure and Dynamics 28(3), 355–361)), (ii) at D178, there is a HB Y169–D178 and a polar contact R164–D178 for BufPrP^C instead of a polar contact Q168–D178 for bovine PrP^C (C.J. Cheng, & V. Daggett. (2014). Molecular dynamics simulations capture the misfolding of the bovine prion protein at acidic pH. Biomolecules 4(1), 181–201), (iii) BufPrP^C owns three 3₁₀ helices at 125–127, 152–156, and in the β2–α2 loop, respectively, and (iv) in the β2–α2 loop, there is a strong π–π stacking and a strong π–cation F175–Y169–R164.(N)NH2, has been discovered.     

Molecular dynamics simulations of PEI mediated DNA aggregation

Understanding the molecular mechanism of polycation induced DNA aggregation and condensation is important for optimal design of gene delivery carriers. In this work, we performed a series of all-atom molecular dynamics (MD) simulations to investigate polyethylenimine (PEI) mediated DNA aggregation. We found that PEIs condense DNA through two mechanisms: polyion bridging and electrostatic screening of the DNA charges. At PEI/DNA charge ratio >1, PEIs can completely neutralize DNAs at a short distance (～12 ? from the C1' atoms), and this distance is found to be insensitive to the exact value of the charge ratio. When excess PEIs are added to a formed DNA-PEI aggregate, they are found to bind to the aggregate and increase its cationic charge. The added PEIs can also replace the PEIs previously bound to the aggregate. The excess PEIs, however, do not change the spacing of the DNAs in the aggregates. Our simulation results shed light on the mechanisms of PEI, and more generally polycation, mediated DNA aggregation and condensation.     

Molecular dynamics studies of the archaeal translocon

The translocon is a protein-conducting channel conserved over all domains of life that serves to translocate proteins across or into membranes. Although this channel has been well studied for many years, the recent discovery of a high-resolution crystal structure opens up new avenues of exploration. Taking advantage of this, we performed molecular dynamics simulations of the translocon in a fully solvated lipid bilayer, examining the translocation abilities of monomeric SecYEbeta by forcing two helices comprised of different amino acid sequences to cross the channel. The simulations revealed that the so-called plug of SecYEbeta swings open during translocation, closing thereafter. Likewise, it was established that the so-called pore ring region of SecYEbeta forms an elastic, yet tight, seal around the translocating oligopeptides. The closed state of the channel was found to block permeation of all ions and water molecules; in the open state, ions were blocked. Our results suggest that the SecYEbeta monomer is capable of forming an active channel.     

Molecular dynamics studies of the conformation of sorbitol

Molecular dynamics simulations of a 3 molal aqueous solution of d-sorbitol (also called d-glucitol) have been performed at 300 K, as well as at two elevated temperatures to promote conformational transitions. In principle, sorbitol is more flexible than glucose since it does not contain a constraining ring. However, a conformational analysis revealed that the sorbitol chain remains extended in solution, in contrast to the bent conformation found experimentally in the crystalline form. While there are 243 staggered conformations of the backbone possible for this open-chain polyol, only a very limited number were found to be stable in the simulations. Although many conformers were briefly sampled, only eight were significantly populated in the simulation. The carbon backbones of all but two of these eight conformers were completely extended, unlike the bent crystal conformation. These extended conformers were stabilized by a quite persistent intramolecular hydrogen bond between the hydroxyl groups of carbon C-2 and C-4. The conformational populations were found to be in good agreement with the limited available NMR data except for the C-2–C-3 torsion (spanned by the O-2–O-4 hydrogen bond), where the NMR data support a more bent structure.     

Molecular dynamics studies of caspase-3

Caspase-3 is a fundamental target for pharmaceutical interventions against a variety of diseases involving disregulated apoptosis. The enzyme is active as a dimer with two symmetry-related active sites, each featuring a Cys-His catalytic dyad and a selectivity loop, which recognizes the characteristic DEVD pattern of the substrate. Here, a molecular dynamics study of the enzyme in complex with two pentapeptide substrates DEVDG is presented, which provides a characterization of the dynamic properties of the active form in aqueous solution. The mobility of the substrate and that of the catalytic residues are rather low indicating a distinct preorganization effect of the Michaelis complex. An essential mode analysis permits us to identify coupled motions between the two monomers. In particular, it is found that the motions of the two active site loops are correlated and tend to steer the substrate toward the reactive center, suggesting that dimerization has a distinct effect on the dynamic properties of the active site regions. The selectivity loop of one monomer turns out to be correlated with the N-terminal region of the p12 subunit of the other monomer, an interaction that is also found to play a fundamental role in the electrostatic stabilization of the quaternary structure. To further characterize the specific influence of dimerization on the enzyme essential motions, a molecular dynamics analysis is also performed on the isolated monomer.     

Molecular dynamics studies of the full-length integrase-DNA complex

We have carried out a molecular dynamics (MD) simulation of full-length HIV-1 integrase (IN) dimer complexed with viral DNA with the aim of gaining information about the enzyme motion and investigating the movement of the catalytic flexible loop (residues 140-149) thought to be essential in the catalytic mechanism of IN. During the simulation, we observed quite a different behavior of this region in the presence or absence of the viral DNA. In particular, the MD results underline the crucial role of the residue Tyr143 in the mechanism of integration of viral DNA into the host chromosome. The present findings confirm the experimental data (e.g., site-directed mutagenesis experiments) showing that the loop is involved in the integration reactions and its mobility is correlated with the catalytic activity of HIV-1 integrase.     

Molecular dynamics studies of the human CD4 protein

A dynamical model for an N-terminal fragment of the human CD4 protein has been determined by computer simulation. The protein has been studied both in vacuo and in solution. Data from both simulations agree moderately well with each other and with the crystal structure. All elements of secondary structure were retained during simulation. Point mutation and sequence replacement studies have shown that a loop in CD4, residues 40–52, is involved in binding with gp120, the human immunodeficiency virus surface glycoprotein.^1,2 Our results show that the gp120-binding loop and a few regions which bind to monoclonal antibodies and class II MHC molecules are the most highly motile areas of the protein. These results are consistent with the suggestion that CD4 binds to target molecules by using induced-fit contacts. © 1994 John Wiley & Sons, Inc.     

Molecular dynamics and docking studies on cardiac troponin C

Cardiac troponin C (cTnC) is the Ca2? dependent switch for contraction in heart muscle making it a potential target for drug research in the therapy of heart failure. Calcium binding on Troponin C (TnC) triggers a series of conformational changes exposing a hydrophobic pocket in the N-domain of TnC (cNTnC), which leads to force generation. Mutations and acidic pH have been related to altering the sensitivity of TnC affecting the efficiency of the heart. Bepridil, identified as a calcium sensitizer to TnC, has been experimentally found to bind to the N-domain pocket of TnC but with negative cooperativity. Screening and de novo design were carried out using LUDI and AUTOLUDI programs in this work to identify and design potential ligands that can bind to the hydrophobic pocket of TnC. Two docking centers and multiple searching radii including 5 ?, 5.5 ?, 6 ?, 6.5 ?, 7.0 ? and 7.5 ? were used in LUDI to screen the ZINC database. Based on the LUDI docking results, 8 molecules were identified from the database with good potential to bind into the binding pocket and they were used as template molecules to generate a series of new molecules by AUTOLUDI design. Out of all the newly-designed molecules, 14 new ligands were recognized to be potential ligands that can bind and fit well into the binding pocket. These molecules can be used as starting molecules to develop TnC ligands. The binding stability and binding affinity of these molecules to the protein was further analyzed by molecular dynamics simulations. The results show that the binding energies, interactions and complex stabilities of 6 ligands are comparable to or better than bepridil.     

Molecular dynamics simulation studies of the effect of phosphocitrate on crystal-induced membranolysis

In this study, following our earlier work on calcium pyrophosphate dihydrate (CPPD) crystal-induced membranolysis, we demonstrate, using the CHARMM method of molecular dynamics simulation, the protective role of phosphocitrate (PC) against solvated dimyristoyl phosphatidylcholine phospholipid bilayer disintegration on contact with the CPPD crystal. Our molecular dynamics simulations studies show that coverage of the CPPD crystal with a layer of phosphocitrate molecules results in the conservation of phospholipid bilayer integrity. We show that the rupture of the lipid bilayer in presence of CPPD and the protective effect of PC are primarily due to electrostatic interactions. The protective role of PC, which may also play an important and potentially therapeutic function against crystal-induced membranolysis is also discussed.     

Molecular dynamics studies on the structural stability of wild-type dog prion protein

Prion diseases such as Creutzfeldt-Jakob disease, variant Creutzfeldt-Jakob diseases, Gerstmann-Str?ussler-Scheinker syndrome, Fatal Familial Insomnia, Kuru in humans, scrapie in sheep, bovine spongiform encephalopathy (or 'mad-cow' disease) and chronic wasting disease in cattle are invariably fatal and highly infectious neurodegenerative diseases affecting humans and animals. However, by now there have not been some effective therapeutic approaches to treat all these prion diseases. In 2008, canine mammals including dogs (canis familials) were the first time academically reported to be resistant to prion diseases (Vaccine 26: 2601-2614 (2008)). Thus, it is very worth studying the molecular structures of dog prion protein to obtain insights into the immunity of dogs to prion diseases. This paper studies the molecular structural dynamics of wild-type dog prion protein. The comparison analyses with rabbit prion protein show that the dog prion protein has stable molecular structures whether under neutral or low pH environments. We also find that the salt bridges such as D177-R163 contribute to the structural stability of wild-type rabbit prion protein under neutral pH environment.     

Molecular dynamics studies of U1A-RNA complexes

The U1A protein binds to a hairpin RNA and an internal-loop RNA with picomolar affinities. To probe the molecular basis of U1A binding, we performed state-of-the-art nanosecond molecular dynamics simulations on both complexes. The good agreement with experimental structures supports the protocols used in the simulations. We compare the dynamics, hydrogen-bonding occupancies, and interfacial flexibility of both complexes and also describe a rigid-body motion in the U1A-internal loop complex that is not observed in the U1A-hairpin simulation. We relate these observations to experimental mutational studies and highlight their significance in U1A binding affinity and specificity.     

Molecular dynamics simulation studies of lipid bilayer systems

The main structural element of biological membranes is a liquid-crystalline lipid bilayer. Other constituents, i.e. proteins, sterols and peptides, either intercalate into or loosely attach to the bilayer. We applied a molecular dynamics simulation method to study membrane systems at various levels of compositional complexity. The studies were started from simple lipid bilayers containing a single type phosphatidylcholine (PC) and water molecules (PC bilayers). As a next step, cholesterol (Chol) molecules were introduced to the PC bilayers (PC-Chol bilayers). These studies provided detailed information about the structure and dynamics of the membrane/water interface and the hydrocarbon chain region in bilayers built of various types of PCs and Chol. This enabled studies of membrane systems of higher complexity. They included the investigation of an integral membrane protein in its natural environment of a PC bilayer, and the antibacterial activity of magainin-2. The latter study required the construction of a model bacterial membrane which consisted of two types of phospholipids and counter ions. Whenever published experimental data were available, the results of the simulations were compared with them.     

Molecular dynamics studies of ion permeation in VDAC

The voltage-dependent anion channel (VDAC) in the outer membrane of mitochondria serves an essential role in the transport of metabolites and electrolytes between the cell matrix and mitochondria. To examine its structure, dynamics, and the mechanisms underlying its electrophysiological properties, we performed a total of 1.77 μs molecular dynamics simulations of human VDAC isoform 1 in DOPE/DOPC mixed bilayers in 1 M KCl solution with transmembrane potentials of 0, ±25, ±50, ±75, and ±100 mV. The calculated conductance and ion selectivity are in good agreement with the experimental measurements. In addition, ion density distributions inside the channel reveal possible pathways for different ion species. Based on these observations, a mechanism underlying the anion selectivity is proposed; both ion species are transported across the channel, but the rate for K⁺ is smaller than that for Cl⁻ because of the attractive interactions between K⁺ and residues on the channel wall. This difference leads to the anion selectivity of VDAC.     

Molecular dynamics studies on structure and dynamics of phospholamban monomer and pentamer in membranes

Phospholamban (PLB) is an integral membrane protein of 52 residues that regulates the activity of the sarcoplasmic reticulum calcium pump in cardiac muscle cells through reversible phosphorylation of Ser16. To explore its possible conformations and dynamics in a monomeric state, we have performed comparative molecular dynamics simulations of unphosphorylated and phosphorylated PLB (pPLB) with various orientations in POPC membranes. The simulations indicate that dynamics of the cytoplasmic domain is highly dependent on its interactions with membranes, that is, large conformational changes in the absence of membrane interactions, but very restricted dynamics in their presence. pPLB shows more structural flexibility in its cytoplasmic domain, which is consistent with experimental observations. We have also performed a simulation of a PLB pentameric structure (the so‐called bellflower model), recently determined in micelles, to investigate its behaviors in a POPC membrane. The cytoplasmic domain in each monomer shows uncorrelated dynamics and undergoes large conformational changes toward the membrane surface during the simulation, which supports the so‐called pinwheel model of the PLB pentamer structure. The hydrophobic nature of the pentameric pore excludes water molecules in the pore region, which illustrates that the pore appears to be an energetic barrier for ion and water translocation. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Molecular dynamics studies on both bound and unbound renin protease

The aspartic protease renin (REN) catalyses the rate-limiting step in the Renin–Angiotensin–Aldosterone System (RAAS), which regulates cardiovascular and renal homoeostasis in living organisms. Renin blockage is therefore an attractive therapeutic strategy for the treatment of hypertension. Herein, computational approaches were used to provide a structural characterization of the binding site, flap opening and dynamic rearrangements of REN in the key conserved residues and water molecules, with the binding of a dodecapeptide substrate or different inhibitors. All these structural insights during catalysis may assist future studies in developing novel strategies for REN inactivation. Our molecular dynamics simulations of several unbound-REN and bound-REN systems indicate similar flexible-segments plasticity with larger fluctuations in those belonging to the C-domain (exposed to the solvent). These segments are thought to assist the flap opening and closure to allow the binding of the substrate and catalytic water molecules. The unbound-REN simulation suggests that the flap can acquire three different conformations: closed, semi-open and open. Our results indicate that the semi-open conformation is already sufficient and appropriate for the binding of the angiotensinogen (Ang) tail, thus contributing to the high specificity of REN, and that both semi-open and open flap conformations are present in free and complexed enzymes. We additionally observed that the Tyr75–Trp39 H-bond has an important role in assisting flap movement, and we highlight several conserved water molecules and amino acids that are essential for the proper catalytic activity of REN.     

Molecular dynamics simulation of the SH3 domain aggregation suggests a generic amyloidogenesis mechanism

We use molecular dynamics simulation to study the aggregation of Src SH3 domain proteins. For the case of two proteins, we observe two possible aggregation conformations: the closed form dimer and the open aggregation state. The closed dimer is formed by "domain swapping"-the two proteins exchange their RT-loops. All the hydrophobic residues are buried inside the dimer so proteins cannot further aggregate into elongated amyloid fibrils. We find that the open structure-stabilized by backbone hydrogen bond interactions-packs the RT-loops together by swapping the two strands of the RT-loop. The packed RT-loops form a beta-sheet structure and expose the backbone to promote further aggregation. We also simulate more than two proteins, and find that the aggregate adopts a fibrillar double beta-sheet structure, which is formed by packing the RT-loops from different proteins. Our simulations are consistent with a possible generic amyloidogenesis scenario.     

Molecular dynamics simulations on the aggregation behavior of indole type organic dye molecules in dye-sensitized solar cells

In Ti0₂ nanostructured dye-sensitized solar cells indole based organic dyes D149, D205 exhibits greater power conversion efficiency. Such organic dye molecules are easily undergone for aggregation. Aggregation in dye molecules leads to reduce electron transfer process in dye-sensitized solar cells. Therefore, anti-aggregating agents such as chenodeoxycholic acid are commonly added to organic dye solution in DSSCs. Studying aggregation of such dye molecules in the absence of semiconductors gives a detailed influence of anti-aggregating agents on dye molecules. Atomistic level of molecular dynamics (MD) simulations were performed on aggregation of indole type dye molecules D149, D205 and D205-F with anti-aggregating agent chenodeoxy cholic acid using AMBER program. The trajectories of the MD simulations were analyzed with order parameters such as radial atom pair distribution functions g(r), diffusion coefficients and root mean square deviations values. MD results suggest that addition of chenodeoxy cholic acid to dyes significantly reduces structural arrangement and increases conformational flexibility and mobility of dye molecules. The influence of semi-perfluorinated alkyl chains in indole dye molecules was analyzed. The parameters such as open-circuit voltage (V_oc) and power conversion efficiency (η) of dye-sensitized solar cells are corroborated with flexibility and diffusion values of dye molecules.     

Molecular dynamics simulations of the aggregation behaviour of overlapped graphene sheets in linear aliphatic hydrocarbons

Molecular dynamics simulations were used to investigate the aggregation of two partially overlapped graphene sheets in hexane, dodecane and eicosane. When partially overlapped graphene sheets are adjacent to one another, they will expel the adsorbed layers of the solvent molecules on the graphene surface, and the amount of overlap will increase. When the overlapped regions of the graphene sheets are separated by solvent molecules, they cannot expel the adsorption layers between them, and so the sheets remain separated. The driving force for aggregation is the van der Waals interaction between the two graphene sheets, while the van der Waals interaction between the graphene sheets and the solvent molecules inhibits graphene aggregation. The diffusion rate of the hydrocarbon molecules with shorter chain lengths is higher. Thus, they diffuse faster during graphene aggregation, which leads to a higher rate of graphene overlapping in the shorter hydrocarbons. This work provides useful insights into graphene aggregation in linear hydrocarbon solvents of varying lengths at the nanoscale.