Probabilistic Approach to Low Strain Rate Atomistic Simulations of Ultimate Tensile Strength of Polymer Crystals.

Molecular dynamics simulations of the tensile ultimate properties of polymer crystals require the use of empirical potentials that model bond dissociation. However, fully reactive potentials are computationally expensive such that reactive simulations cannot reach the low strain rates of typical experiments. Here, we present a hybrid approach that uses the simplicity of a classical, nonreactive potential, information from bond dissociation energy calculations, and a probabilistic expression that mimics bond breaking. The approach is demonstrated for poly(p-phenylene terephthalamide) and, with one tunable parameter, the calculated tensile ultimate stress matches that obtained using a fully reactive simulation at high strain rates. Then, the hybrid simulations are run at much lower strain rates where the ultimate tensile stress is strain rate-independent and consistent with the expected experimental range.

Methanol-Assisted ADMET Polymerization of Semiaromatic Amides.

A method for the acyclic diene metathesis polymerization of semiaromatic amides is described. The procedure uses second-generation Grubbs' catalyst and N-cyclohexyl-2-pyrrolidone (CHP), a high boiling, polar solvent capable of solubilizing both monomer and polymer. The addition of methanol to the reaction was found to significantly increase polymer molar mass although the role of the alcohol is currently not understood. Hydrogenation with hydrogen gas and Wilkinson's catalyst resulted in near-quantitative saturation. All polymers synthesized here exhibit a hierarchical semicrystalline morphology driven by ordering of aromatic amide groups via strong nonbonded interactions. Furthermore, the melting points can be tuned over a >100 °C range by precise substitution at just one of the backbone positions on each mer (<5% of the total).

Heterologous Synthesis and Secretion of Ricinoleic Acid in Starmerella bombicola with Sophorolipid as an Intermediate.

Ricinoleic acid (RA) is a long-chain hydroxy fatty acid produced from castor bean that is used in the manufacturing of a variety of industrial products. The demand for RA keeps increasing due to its broad applications. However, due to the presence of a potent toxin ricin, the native oilseed plant is not an ideal source for hydroxy fatty acid production. Although there have been significant efforts on engineering different microorganisms for heterologous production of RA, all had very limited success. The main reason for this is the exhibited toxicity of the intracellularly accumulated RA. To avoid this issue, we genetically modified a Starmerella bombicola strain by engineering its native sophorolipid production pathway to direct the synthesized RA bound with sophorolipid to be secreted out of the cell, followed by acid hydrolysis to recover RA. The engineered S. bombicola strain expressing the heterologous codon-optimized oleate hydroxylase-encoding gene from ergot fungus Claviceps purpurea resulted in a record production titer of RA at about 2.96 g/L. Thus, this work highlights a new strategy to produce a high level of hydroxy fatty acids in engineered yeast through a sophorolipid intermediate, enabling a new biocatalysis platform for the future.

Counterion Effects on Aggregate Structure of 12-Hydroxystearate Salts in Hexane: A Quantum Mechanical and Molecular Dynamics Simulation Study.

Salts of 12-hydroxystearate are important organogelators and grease thickeners, but a structural rationale for their rheological properties remains elusive. We use quantum mechanical calculations and molecular dynamics (MD) simulations to analyze aggregate structures for (1) ( R)-12-hydroxystearic acid (( R)-12HSA), (2) lithium ( R)-12-hydroxystearate (( R)-Li12HS), and (3) sodium ( R)-12-hydroxystearate (( R)-Na12HS). First, quantum mechanical calculations were used to establish the structure and complexation energies of dimers of acetic acid, lithium acetate, and sodium acetate. The expected acetic acid dimer is predicted, and both the lithium acetate and sodium acetate dimer formed a C2 h-symmetric structure. All dimers were sufficiently stable to allow modeling them as pseudocovalent complexes in all-atom, explicit solvent MD. After microsecond-long MD, all systems produced strong ringlike ordered nuclei. The C2 h lithium salt molecules produced aggregates that had the most efficient packing at the head group and a higher frequency of hydroxyl hydrogen bonding compared to the sodium salt. This ordering propensity explains the high melting temperature of ( R)-Li12HS. Also, the higher frequency of hydrogen bonding leads to fewer solvent-exposed hydrogen bond partners. This explains why lithium is a common counterion in high-temperature and water-resistant greases.

Effects of Optical Purity and Finite System Size on Self-Assembly of 12-Hydroxystearic Acid in Hexane: Molecular Dynamics Simulations.

12-Hydroxystearic acid (12HSA) and its derivatives are well-known organogelators, and they play critical roles in a variety of applications. The overall aggregate structure of 12HSA is sensitive to the chirality at the 12th carbon, but a fundamental understanding of this dependence is lacking. In this study, molecular dynamics simulations were conducted on microsecond long time scales for (1) (R)-12HSA, (2) (S)-12HSA, and (3) a 50/50 racemic mixture, each solvated at 12.5 wt % in explicit hexane. Self-assembly was accelerated by turning off alkyl chain dihedral gauche states and forcing the molecules to adopt an all-trans conformation. The stability of the resulting aggregates was tested by quenching them with access to gauche states restored. Ordered aggregates produced from optically pure (R)- and (S)-12HSA remained stable for at least 1 µs. The characteristic ordered structure observed is termed a "ring-of-rings" motif, and it contains two twisted six-membered ringlike bundles connected through acetic acid dimerization and surrounded by six satellite bundles. The chirality at the 12th carbon dictates the overall twist of the rings and thereby the handedness of the aggregates. Racemic mixtures did not produce stable ordered aggregates likely due to insufficient enantiomerically pure ring formation. The most prevalent finite-size effect observed was the stochastic formation of percolating aggregates, which were later avoided by using solvent-permeable, solute-impermeable, confining walls. The resulting ordered aggregates were in all important ways identical to those produced in unconfined systems. The combination of cycling off and on gauche states and the semipermeable walls may be an important new way to study the self-assembly underlying aggregation at industrially relevant concentrations of surfactants in organic solvents.

Aggregation of 12-Hydroxystearic Acid and Its Lithium Salt in Hexane: Molecular Dynamics Simulations.

12-Hydroxystearic acid (12HSA) is a well-known organogelator, and its metal salts and derivatives find roles in many important applications. The structures of aggregates of 12-hydroxysteric acid and its salts depend sensitively on cation type, but a fundamental understanding of this phenomenon is lacking. In this study, molecular dynamics simulations were conducted on the microsecond long time scales for (1) 12HSA and (2) its lithium salt, each at 12.5 wt % in explicit hexane solvent. Self-assembly was accelerated by using a modified potential to prohibit alkane chain dihedral gauche states (all-trans-12HSA) and then verified by continuation using standard force-field parameters. In three independent simulation, acceleration using "gauche-less" potentials resulted in self-assembled pseudocrystalline aggregates through formation of polarized five- and six-membered rings between inter-12-hydroxyl groups and head-to-head carboxylic acid dimerization. When subjected to the unmodified dihedral potential, two of the three structures remained stable after 1 µs of MD. Stable structures exhibited a "ring-of-rings" motif, composed of two six-membered acetic acid-dimerized ring bundles with six satellite rings, while the unstable structure did not. In strong contrast, the lithium salt produced a network of fibrils that spanned the volume of the sample. When lithium ions were substituted for carboxylic acid protons in the stable acid structures, they remained intact but lost their chiral nature. Both the acid and lithium structures displayed scattering peaks that agreed with experiment. Taken together, our results suggest that this ring-of-rings structure could be a primary feature of the self-assembly of 12HSA in organic solvents.

Energetics and mechanism of the normal-to-amyloidogenic isomerization of ß2-microglobulin: on-the-fly string method calculations.

We use on-the-fly finite temperature string method in collective variables to study the transition from a normal to an amyloidogenic conformation of ß2-microglobulin. We show that the protonation state of two histidine residues is of key importance and that under acidic (protonating) conditions, the transition to the amyloidgenic form is facilitated by both displacement of N-terminal residues to disrupt a hydrophobic pocket and by side-chain/side-chain electrostatic attraction, both of which facilitate a cis-trans prolyl isomerization. The free energy barriers for the normal-to-amyloidogenic isomerization are found to be 14.9 and 7.1 kcal/mol for the neutral and protonated cases, respectively, consistent with enhanced amyloidgenesis at low pH observed both in vitro and in hemodialysis-associated amyloidosis, and somewhat lower than experimentally determined barriers for bare prolyl cis-trans isomerization. We suggest specific mutagenesis experiments which could be used to further validate the mechanism observed.

Enhanced meta-analysis of acetylcholine binding protein structures reveals conformational signatures of agonism in nicotinic receptors.

The soluble acetylcholine binding protein (AChBP) is the default structural proxy for pentameric ligand-gated ion channels (LGICs). Unfortunately, it is difficult to recognize conformational signatures of LGIC agonism and antagonism within the large set of AChBP crystal structures in both apo and ligand-bound states, primarily because AChBP conformations in this set are nearly superimposable (root mean square deviation < 1.5 Å). We have undertaken a systematic, alignment-free approach to elucidate conformational differences displayed by AChBP that cleanly differentiate apo/antagonist-bound from agonist-bound states. Our approach uses statistical inference based on both crystallographic states and conformations sampled during long molecular dynamics simulations to select important inter-C(α) distances and map their collective values onto functional states. We observe that binding of (nAChR) agonists to AChBP elicits clockwise rotation of the inner ß-sheet with respect to the outer ß-sheet, causing tilting of the cys-loop away from the five-fold axis, in a manner quite similar to that speculated for α-subunits of the heteromeric nAChR structure (Unwin, J Mol Biol 2005;346:967), making this motion potentially important in transmission of the gating signal to the transmembrane domain of a LGIC. The method is also successful at discriminating partial from full agonists and supports the hypothesis that a particularly controversial ligand, lobeline, is in fact an LGIC antagonist.

Mapping spatial relationships between residues in the ligand-binding domain of the 5-HT3 receptor using a molecular ruler.

The serotonin 5-HT(3) receptor (5-HT(3)R) is a member of the Cys-loop ligand-gated ion channel family. We used a combination of site-directed mutagenesis, homology modeling, and ligand-docking simulations to analyze antagonist-receptor interactions. Mutation of E236, which is near loop C of the binding site, to aspartate prevents expression of the receptor on the cell surface, and no specific ligand binding can be detected. On the other hand, mutation to glutamine, asparagine, or alanine produces receptors that are expressed on the cell surface, but decreases receptor affinity for the competitive antagonist d-tubocurarine (dTC) 5-35-fold. The results of a double-mutant cycle analysis employing a panel of dTC analogs to identify specific points of interactions between the dTC analogs and E236 are consistent with E236 making a direct physical interaction with the 12 -OH of dTC. dTC is a rigid molecule of known three-dimensional structure. Together with previous studies linking other regions of dTC to specific residues in the binding site, these data allow us to define the relative spatial arrangement of three different residues in the ligand-binding site: R92 (loop D), N128 (loop A), and E236 (near loop C). Molecular modeling employing these distance constraints followed by molecular-dynamics simulations produced a dTC/receptor complex consistent with the experimental data. The use of the rigid ligands as molecular rulers in conjunction with double-mutant cycle analysis provides a means of mapping the relative positions of various residues in the ligand-binding site of any ligand-receptor complex, and thus is a useful tool for delineating the architecture of the binding site.