Transport properties of CO2-expanded acetonitrile from molecular dynamics simulations

This is the publisher's version, also available electronically from http://scitation.aip.org/content/aip/journal/jcp/126/7/10.1063/1.2434968.Carbon-dioxide-expanded liquids, which are mixtures of organic liquids and compressed CO2, are novel media used in chemical processing. The authors present a molecular simulation study of the transport properties of liquid mixtures formed by acetonitrile and carbon dioxide, in which the CO2 mole fraction is adjusted by changing the pressure, at a constant temperature of 298K. They report values of translational diffusion coefficients, rotational correlation times, and shear viscosities of the liquids as function of CO2 mole fraction. The simulation results are in good agreement with the available experimental data for the pure components and provide interesting insights into the largely unknown properties of the mixtures, which are being recognized as important novel materials in chemical operations. We find that the calculated quantities exhibit smooth variation with composition that may be represented by simple model equations. The translational and rotational diffusion rates increase with CO2 mole fraction for both the acetonitrile and carbon dioxide components. The shear viscosity decreases with increasing amount of CO2, varying smoothly between the values of pure acetonitrile and pure carbon dioxide. Our results show that adjusting the amount of CO2 in the mixture allows the variation of transport rates by a factor of 3–4 and liquidviscosity by a factor of 8. Thus, the physical properties of the mixture may be tailored to the desired range by changes in the operating conditions of temperature and pressure

A 3D–Predicted Structure of the Amine Oxidase Domain of Lysyl Oxidase–Like 2

Lysyl oxidase–like 2 (LOXL2) has been recognized as an attractive drug target for anti–fibrotic and anti–tumor therapies. However, the structure–based drug design of LOXL2 has been very challenging due to the lack of structural information of the catalytically–competent LOXL2. In this study; we generated a 3D–predicted structure of the C–terminal amine oxidase domain of LOXL2 containing the lysine tyrosylquinone (LTQ) cofactor from the 2.4Å crystal structure of the Zn2+–bound precursor (lacking LTQ; PDB:5ZE3); this was achieved by molecular modeling and molecular dynamics simulation based on our solution studies of a mature LOXL2 that is inhibited by 2–hydrazinopyridine. The overall structures of the 3D–modeled mature LOXL2 and the Zn2+–bound precursor are very similar (RMSD = 1.070Å), and disulfide bonds are conserved. The major difference of the mature and the precursor LOXL2 is the secondary structure of the pentapeptide (His652–Lys653–Ala654–Ser655–Phe656) containing Lys653 (the precursor residue of the LTQ cofactor). We anticipate that this peptide is flexible in solution to accommodate the conformation that enables the LTQ cofactor formation as opposed to the β–sheet observed in 5ZE3. We discuss the active site environment surrounding LTQ and Cu2+ of the 3D–predicted structure

Kinetic pathway analysis of an α-helix in two protonation states: Direct observation and optimal dimensionality reduction

This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in J. Chem. Phys. 150, 074902 (2019); https://doi.org/10.1063/1.5082192 and may be found at https://aip.scitation.org/doi/full/10.1063/1.5082192Thermodynamically stable conformers of secondary structural elements make a stable tertiary/quaternary structure that performs its proper biological function efficiently. Formation mechanisms of secondary and tertiary/quaternary structural elements from the primary structure are driven by the kinetic properties of the respective systems. Here we have carried out thermodynamic and kinetic characterization of an alpha helical heteropeptide in two protonation states, created with the addition and removal of a proton involving a single histidine residue in the primary structure. Applying far-UV circular dichroism spectroscopy, the alpha helix is observed to be significantly more stable in the deprotonated state. Nanosecond laser temperature jump spectroscopy monitoring time-resolved tryptophan fluorescence on the protonated conformer is carried out to measure the kinetics of this system. The measured relaxation rates at a final temperature between 296K and 314 K generated a faster component of 20 ns–11 ns and a slower component of 314 ns–198 ns. Atomically detailed characterization of the helix-coil kinetic pathways is performed based on all-atom molecular dynamics trajectories of the two conformers. Application of clustering and kinetic coarse-graining with optimum dimensionality reduction produced description of the trajectories in terms of kinetic models with two to five states. These models include aggregate states corresponding to helix, coil, and intermediates. The “coil” state involves the largest number of conformations, consistent with the expected high entropy of this structural ensemble. The “helix” aggregate states are found to be mixed with the full helix and partially folded forms. The experimentally observed higher helix stability in the deprotonated form of the alpha helical heteropeptide is reflected in the nature of the “helix” aggregate state arising from the kinetic model. In the protonated form, the “coil” state exhibits the lowest free energy and longest lifetime, while in the deprotonated form, it is the “helix” that is found to be most stable. Overall, the coarse grained models suggest that the protonation of a single histidine residue in the primary structure induces significant changes in the free energy landscape and kinetic network of the studied helix-forming heteropeptide

Length Dependent Folding Kinetics of Alanine-Based Helical Peptides from Optimal Dimensionality Reduction

We present a computer simulation study of helix folding in alanine homopeptides (ALA)n of length n = 5, 8, 15, and 21 residues. Based on multi-microsecond molecular dynamics simulations at room temperature, we found helix populations and relaxation times increasing from about 6% and ~2 ns for ALA5 to about 60% and ~500 ns for ALA21, and folding free energies decreasing linearly with the increasing number of residues. The helix folding was analyzed with the Optimal Dimensionality Reduction method, yielding coarse-grained kinetic models that provided a detailed representation of the folding process. The shorter peptides, ALA5 and ALA8, tended to convert directly from coil to helix, while ALA15 and ALA21 traveled through several intermediates. Coarse-grained aggregate states representing the helix, coil, and intermediates were heterogeneous, encompassing multiple peptide conformations. The folding involved multiple pathways and interesting intermediate states were present on the folding paths, with partially formed helices, turns, and compact coils. Statistically, helix initiation was favored at both termini, and the helix was most stable in the central region. Importantly, we found the presence of underlying universal local dynamics in helical peptides with correlated transitions for neighboring hydrogen bonds. Overall, the structural and dynamical parameters extracted from the trajectories are in good agreement with experimental observables, providing microscopic insights into the complex helix folding kinetics

Experiments and comprehensive simulations of the formation of a helical turn

We consider the kinetics and thermodynamics of a helical turn formation in the peptide Ac-WAAAH-NH2. NMR measurements indicate that the peptide has significant tendency to form a structure of a helical turn, while temperature dependent CD establishes the helix fraction at different temperatures. Molecular Dynamics and Milestoning simulations agree with experimental observables and suggests an atomically detailed picture for the turn formation. Using a network representation two alternative mechanisms of folding are identified: (i) a direct cooperative mechanism from the unfolded to the folded state without intermediate formation of hydrogen bonds and (ii) an indirect mechanism with structural intermediates with two residues in a helical conformation. This picture is consistent with kinetic measurements that reveal two experimental time scales of sub nanosecond and several nanoseconds

Permeation of the three aromatic dipeptides through lipid bilayers: Experimental and computational study

Publisher's note added August 2016: "This article was originally published online on 27 June 2016 with a sentence missing in the Acknowledgments. After the funding acknowledgments, it should read, “G.S.J. would like to thank Wilson R. Veras Tavarez and Elizabeth De Leon Olmeda of UCC for helpful comments.” AIP Publishing apologizes for this error. All online versions of the article were corrected on 28 June 2016; the article is correct as it appears in the printed version of the journal."The time-resolved parallel artificial membrane permeability assay with fluorescence detection and comprehensive computer simulations are used to study the passive permeation of three aromatic dipeptides—N-acetyl-phenylalanineamide (NAFA), N-acetyltyrosineamide (NAYA), and N-acetyltryptophanamide (NATA) through a 1,2-dioleoyl-sn-glycero-3-phospocholine (DOPC) lipid bilayer. Measured permeation times and permeability coefficients show fastest translocation for NAFA, slowest for NAYA, and intermediate for NATA under physiological temperature and pH. Computationally, we perform umbrella sampling simulations to model the structure, dynamics, and interactions of the peptides as a function of z, the distance from lipid bilayer. The calculated profiles of the potential of mean force show two strong effects—preferential binding of each of the three peptides to the lipid interface and large free energy barriers in the membrane center. We use several approaches to calculate the position-dependent translational diffusion coefficients D(z), including one based on numerical solution the Smoluchowski equation. Surprisingly, computed D(z) values change very little with reaction coordinate and are also quite similar for the three peptides studied. In contrast, calculated values of sidechain rotational correlation times τrot(z) show extremely large changes with peptide membrane insertion—values become 100 times larger in the headgroup region and 10 times larger at interface and in membrane center, relative to solution. The peptides’ conformational freedom becomes systematically more restricted as they enter the membrane, sampling α and β and C7eq basins in solution, α and C7eq at the interface, and C7eq only in the center. Residual waters of solvation remain around the peptides even in the membrane center. Overall, our study provides an improved microscopic understanding of passive peptide permeation through membranes, especially on the sensitivity of rotational diffusion to position relative to the bilayer. Published by AIP Publishing. [http://dx.doi.org/10.1063/1.4954241

Characterization of multiple stable conformers of the EC5 domain of E-cadherin and the interaction of EC5 with E-cadherin peptides

The objectives of this work were to express the EC5 domain of E-cadherin and determine its structural characteristics as well as to evaluate the binding properties of HAV and BLG4 peptides to EC5 using spectroscopic methods. Homophilic interactions of E-cadherins are responsible for cell-cell adhesion in the adherens junctions of the biological barriers (i.e., intestinal mucosa and blood-brain barriers). The EC5 domain of E-cadherin has an important role in T-cell adhesion to intestinal mucosa via αEβ7 integrin-E-cadherin interactions. In this study, the expressed EC5 has a high thermal stability (Tm = 64.3 °C); it also has two stable conformations at room temperature, which convert to one conformation at approximately 54.5 °C. NMR and FTIR showed that HAV and BLG4 peptides bind to EC5. HSQC-NMR showed that either Asn or Gln of EC5 was involved in the interactions with HAV and BLG4 peptides. EC5 underwent a conformational change upon interaction with the HAV and BLG4 peptides. Finally, the binding properties of both peptides were modeled by docking experiments, and the results suggest that Asn-46 and Asn-75 of EC5 could be involved during the interaction with the peptides and that the Ser and Trp residues of the HAV and BLG4 peptides, respectively, were important for binding to EC5

Dynamic elements and kinetics: Most favorable conformations of peptides in solution with measurements and simulations

This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared inJ. Chem. Phys. 151, 225102 (2019); doi: 10.1063/1.5131782 and may be found at https://aip.scitation.org/doi/10.1063/1.5131782.Small peptides in solution adopt a specific morphology as they function. It is of fundamental interest to examine the structural properties of these small biomolecules in solution and observe how they transition from one conformation to another and form functional structures. In this study, we have examined the structural properties of a simple dipeptide and a five-residue peptide with the application of far-UV circular dichroism (CD) spectroscopy as a function of temperature, fluorescence anisotropy, and all-atom molecular dynamics simulation. Analysis of the temperature dependent CD spectra shows that the simplest dipeptide N-acetyl-tryptophan-amide (NATA) adopts helical, beta sheet, and random coil conformations. At room temperature, NATA is found to have 5% alpha-helical, 37% beta sheet, and 58% random coil conformations. To our knowledge, this type of structural content in a simplest dipeptide has not been observed earlier. The pentapeptide (WK5) is found to have four major secondary structural elements with 8% 310 helix, 14% poly-L-proline II, 8% beta sheet, and 14% turns. A 56% unordered structural population is also present for WK5. The presence of a significant population of 310 helix in a simple pentapeptide is rarely observed. Fluorescence anisotropy decay (FAD) measurements yielded reorientation times of 45 ps for NATA and 120 ps for WK5. The fluorescence anisotropy decay measurements reveal the size differences between the two peptides, NATA and WK5, with possible contributions from differences in shape, interactions with the environment, and conformational dynamics. All-atom molecular dynamics simulations were used to model the structures and motions of these two systems in solution. The predicted structures sampled by both peptides qualitatively agree with the experimental findings. Kinetic modeling with optimal dimensionality reduction suggests that the slowest dynamic processes in the dipeptide involve sidechain transitions occurring on a 1 ns timescale. The kinetics in the pentapeptide monitors the formation of a distorted helical structure from an extended conformation on a timescale of 10 ns. Modeling of the fluorescence anisotropy decay is found to be in good agreement with the measured data and correlates with the main contributions of the measured reorientation times to individual conformers, which we define as dynamic elements. In NATA, the FAD can be well represented as a sum of contributions from representative conformers. This is not the case in WK5, where our analysis suggests the existence of coupling between conformational dynamics and global tumbling. The current study involving detailed experimental measurements and atomically detailed modeling reveals the existence of specific secondary structural elements and novel dynamical features even in the simplest peptide systems