Conformational Simulations of Aqueous Solvated α-Conotoxin GI and Its Single Disulfide Analogues Using a Polarizable Force Field Model

Treating Polarization Effects in Charged and Polar Bio-Molecules Through Variable Electrostatic Parameters

Polarization plays important roles in charged and hydrogen bonding containing systems. Much effort ranging from the construction of physics-based models to quantum mechanism (QM)-based and machine learning (ML)-assisted models have been devoted to incorporating the polarization effect into the conventional force fields at different levels, such as atomic and coarse grained (CG). The application of polarizable force fields or polarization models was limited by two aspects, namely, computational cost and transferability. Different from physics-based models, no predetermining parameters were required in the QM-based approaches. Taking advantage of both the accuracy of QM calculations and efficiency of molecular mechanism (MM) and ML, polarization effects could be treated more efficiently while maintaining the QM accuracy. The computational cost could be reduced with variable electrostatic parameters, such as the charge, dipole, and electronic dielectric constant with the help of linear scaling fragmentation-based QM calculations and ML models. Polarization and entropy effects on the prediction of partition coefficient of druglike molecules are demonstrated by using both explicit or implicit all-atom molecular dynamics simulations and machine learning-assisted models. Directions and challenges for future development are also envisioned.

Shear viscosity prediction of alcohols, hydrocarbons, halogenated, carbonyl, nitrogen-containing, and sulfur compounds using the variable force fields

Viscosity of organic liquids is an important physical property in applications of printing, pharmaceuticals, oil extracting, engineering, and chemical processes. Experimental measurement is a direct but time-consuming process. Accurately predicting the viscosity with a broad range of chemical diversity is still a great challenge. In this work, a protocol named Variable Force Field (VaFF) was implemented to efficiently vary the force field parameters, especially λ_vdW, for the van der Waals term for the shear viscosity prediction of 75 organic liquid molecules with viscosity ranging from -9 to 0 in their nature logarithm and containing diverse chemical functional groups, such as alcoholic hydroxyl, carbonyl, and halogenated groups. Feature learning was applied for the viscosity prediction, and the selected features indicated that the hydrogen bonding interactions and the number of atoms and rings play important roles in the property of viscosity. The shear viscosity prediction of alcohols is very difficult owing to the existence of relative strong intermolecular hydrogen bonding interaction as reflected by density functional theory binding energies. From radial and spatial distribution functions of methanol, we found that the van der Waals related parameters λ_vdW are more crucial to the viscosity prediction than the rotation related parameters, λ_tor. With the variable λ_vdW-based all-atom optimized potentials for liquid simulations force field, a great improvement was observed in the viscosity prediction for alcohols. The simplicity and uniformity of VaFF make it an efficient tool for the prediction of viscosity and other related properties in the rational design of materials with the specific properties.

Structural and dynamical effects of targeted mutations on μO-Conotoxin MfVIA: Molecular simulation studies

Conotoxins are a group of cysteine-rich, neurotoxic peptides isolated from the venom of marine cone snails. MfVIA is a member of the μO-conotoxin family, and acts as an inhibitor of subtype 1.8 voltage-gated sodium ion channels (Na_V1.8). The unique selectivity of MfVIA as an inhibitor of Na_V1.8 makes it an ideal peptide for elucidation of the physiological functions of this voltage-gated ion channel. Previous experimental studies of point mutations of MfVIA showed that the double mutant [E5K,E8K] exhibited greater activity at Na_V1.8 relative to the wild-type toxin. The present study employs molecular dynamics (MD) simulations to examine the effects of various mutations at these key residues (E5 and E8) on the structure and dynamics of MfVIA. Five double mutants were studied, in which the positions 5 and 8 residues were mutated to amino acids with a range of different physicochemical properties, namely [E5A,E8A], [E5D,E8D], [E5F,E8F], [E5K,E8K], and [E5R,E8R]. Except for [E5D,E8D], all of the mutants tend to show decreased contacts at the N-terminus owing to the loss of the R1-E5 salt bridge relative to that of the wild-type, which subsequently lead to greater exposure and flexibility of the N-terminus for most of the mutant peptides studied, potentially rendering them more able to interact with other species, including Na_V1.8. Molecular docking studies of the peptides to Na_V1.8 via different binding mechanisms suggest that the [E5R, E8R] mutant may be especially worthy of further investigation owing to its predicted binding mode, which differs markedly from those of the other peptides in this study.

An Iterative Fragment Scheme for the ACKS2 Electronic Polarization Model: Application to Molecular Dimers and Chains

The treatment of electrostatic interactions is a key ingredient in the force field-based simulation of condensed phase systems. Most approaches used fixed, site-specific point charges. Yet, it is now clear that many applications of force fields (FFs) demand more sophisticated treatments, prompting the implementation of charge equilibration methods in polarizable FFs to allow the redistribution of charge within the system. One approach allowing both, charge redistribution and site-specific polarization, while at the same time solving methodological shortcomings of earlier methods, is the first-principles-derived atom-condensed Kohn-Sham density functional theory method approximated to the second order (ACKS2). In this work, we present two fragment approaches to ACKS2, termed f-ACKS2 and a self-consistent version, scf-ACKS2, that treat condensed phase systems as a collection of electronically polarizable molecular fragments. The fragmentation approach to ACKS2 not only leads to a more transferable and less system-specific collection of electronic response parameters but also opens up the method to large condensed phase systems. We validate the accuracies of f-ACKS2 and scf-ACKS2 by comparing polarization energies and induced dipole moments for a number of charged hydrocarbon dimers against DFT reference calculations. Finally, we also apply both fragmented ACKS2 variants to calculate the polarization energy for electron-hole pair separation along a chain of anthracene molecules and find excellent agreement with reference DFT calculations.

Snails In Silico: A Review of Computational Studies on the Conopeptides

Marine cone snails are carnivorous gastropods that use peptide toxins called conopeptides both as a defense mechanism and as a means to immobilize and kill their prey. These peptide toxins exhibit a large chemical diversity that enables exquisite specificity and potency for target receptor proteins. This diversity arises in terms of variations both in amino acid sequence and length, and in posttranslational modifications, particularly the formation of multiple disulfide linkages. Most of the functionally characterized conopeptides target ion channels of animal nervous systems, which has led to research on their therapeutic applications. Many facets of the underlying molecular mechanisms responsible for the specificity and virulence of conopeptides, however, remain poorly understood. In this review, we will explore the chemical diversity of conopeptides from a computational perspective. First, we discuss current approaches used for classifying conopeptides. Next, we review different computational strategies that have been applied to understanding and predicting their structure and function, from machine learning techniques for predictive classification to docking studies and molecular dynamics simulations for molecular-level understanding. We then review recent novel computational approaches for rapid high-throughput screening and chemical design of conopeptides for particular applications. We close with an assessment of the state of the field, emphasizing important questions for future lines of inquiry.

Host-guest interaction of β-cyclodextrin with isomeric ursolic acid and oleanolic acid: physicochemical characterization and molecular modeling study

Ursolic acid (UA) and oleanolic acid (OA) are insoluble drugs. The objective of this study was to encapsulate them into β-cyclodextrin (β-CD) and compare the solubility and intermolecular force of β-CD with the two isomeric triterpenic acids. The host-guest interaction was explored in liquid and solid state by ultraviolet-visible absorption,1 H NMR, phase solubility analysis, and differential scanning calorimetry, X-ray powder diffractometry, and molecular modeling studies. Both experimental and theoretical studies revealed that β-CD formed 1: 1 water soluble inclusion complexes and the complexation process was naturally favorable. In addition, the overall results suggested that ring E with a carboxyl group of the drug was encapsulated into the hydrophobic CD nanocavity. Therefore, a clear different inclusion behavior was observed, and UA exhibited better affinity to β-CD compared with OA in various media due to little steric interference, which was beneficial to form stable inclusion complex with β-CD and increase its water solubility effectively.

Factor Analysis of Conformations and NMR Signals of Rotaxanes: AIMD and Polarizable MD Simulations

The interlocked ⟨rod | ring⟩ structures of pseudorotaxanes and [2]rotaxanes are usually maintained by the complex hydrogen-bonding (H-bonding) network between the rod and ring. Ab initio molecular dynamics (AIMD) using generalized energy-based fragmentation approach and polarizable force field (polar FF)-based molecular dynamics (MD) simulations were performed to investigate the conformational changes of mechanically interlocked systems and to obtain the ensemble-averaged NMR chemical shifts. Factor analysis (FA) demonstrates that the ring H-donor (2,6 pyridinedicarboxamide group) plays an important role in the ring-rod recognition. In comparison to the conventional fixed-charge force field, the polarization effect is crucial to account for the H-bonding interactions in supramolecular systems. In the hybrid scheme, the polar FF-based MD simulations are used to generate different initial states for the AIMD simulations, which are able to give better prediction of ensemble-averaged NMR signals for chemically equivalent amide protons. The magnitude of the deshielding shift of NMR signal is correlated with the length of hydrogen bond. The polar FF model with variable charges shows that the dipole-dipole interactions between the flexible diethylene glycol chain of ring and polar solvents induce the upfield shifts of NMR signals of rod H-donors and the directional distribution of the neighboring CH3CN solvents.

Assignment of absolute stereostructures through quantum mechanics electronic and vibrational circular dichroism calculations

Elucidation of absolute configuration of chiral molecules including structurally complex natural products remains a challenging problem in organic chemistry. A reliable method for assigning the absolute stereostructure is to combine the experimental circular dichroism (CD) techniques such as electronic and vibrational CD (ECD and VCD), with quantum mechanics (QM) ECD and VCD calculations. The traditional QM methods as well as their continuing developments make them more applicable with accuracy. Taking some chiral natural products with diverse conformations as examples, this review describes the basic concepts and new developments of QM approaches for ECD and VCD calculations in solution and solid states.

Generalized energy-based fragmentation approach and its applications to macromolecules and molecular aggregates

Conspectus The generalized energy-based fragmentation (GEBF) approach provides a very simple way of approximately evaluating the ground-state energy or properties of a large system in terms of ground-state energies of various small "electrostatically embedded" subsystems, which can be calculated with any traditional ab initio quantum chemistry (X) method (X = Hartree-Fock, density functional theory, and so on). Due to its excellent parallel efficiency, the GEBF approach at the X theory level (GEBF-X) allows full quantum mechanical (QM) calculations to be accessible for systems with hundreds and even thousands of atoms on ordinary workstations. The implementation of the GEBF approach at various theoretical levels can be easily done with existing quantum chemistry programs. This Account reviews the methodology, implementation, and applications of the GEBF-X approach. This method has been successfully applied to optimize the structures of various large systems including molecular clusters, polypeptides, proteins, and foldamers. Such investigations could allow us to elucidate the origin and nature of the cooperative interaction in secondary structures of long peptides or the driving force of the self-assembly processes of aromatic oligoamides. These GEBF-based QM calculations reveal that the structures and stability of various complex systems result from a subtle balance of many types of noncovalent interactions such as hydrogen bonding and van der Waals interactions. The GEBF-based ab initio molecular dynamics (AIMD) method also allows the investigation of dynamic behaviors of large systems on the order of tens of picoseconds. It was demonstrated that the conformational dynamics of two model peptides predicted by GEBF-based AIMD are noticeably different from those predicted by the classical force field MD method. With the target of extending QM calculations to molecular aggregates in the condensed phase, we have implemented the GEBF-based multilayer hybrid models, which could provide satisfactory descriptions of the binding energies between a solute molecule and its surrounding waters and the chain-length dependence of the conformational changes of oligomers in aqueous solutions. A coarse-grained polarizable molecular mechanics model, furnished with GEBF-X dipole moments of subsystems, exhibits some advantages of treating the electrostatic polarization with reduced computational costs. We anticipate that the GEBF approach will continue to develop with the ultimate goal of studying complicated phenomena at mesoscopic scales and serve as a practical tool to elucidate the structure and dynamics of chemical and biological systems.

Rational design of core-shell molecularly imprinted polymer based on computational simulation and Doehlert experimental optimization: application to the separation of tanshinone IIA from Salvia miltiorrhiza Bunge

Computational simulation and Doehlert experimental optimization were done for the rational design of a core-shell molecularly imprinted polymer (CS-MIP) for use in the highly selective separation of Tanshinone IIA (TSIIA) from the crude extracts of Salvia miltiorrhiza Bunge (SMB). The functional monomer layer of the polymer shells directed the selective occurrence of imprinting polymerization at the surface of silica through the copolymerization of vinyl end groups with functional monomers and also drove TSIIA templates into the formed polymer shells through the charge-transfer complex interactions between TSIIA and the functional monomer layer. As a result, the maximum rebinding capacity was achieved with the use of optimal grafting ratio by the Doehlert design. The CS-MIP exhibited high recognition selectivity and binding affinity to TSIIA. When the imprinted particles were used as dispersive solid phase extraction sorbents, the recovery yield of TSIIA reached 93% by a one-step extraction from the crude extracts of SMB, and the purity of TSIIA was larger than 98% by HPLC analysis. These results show the possibility of a highly selective separation and enrichment of TSIIA from the SMB using the TSIIA-imprinted core-shell molecularly imprinted polymers.

Molecular modeling and dynamics studies with explicit inclusion of electronic polarizability. Theory and applications

A current emphasis in empirical force fields is on the development of potential functions that explicitly treat electronic polarizability. In the present article, the commonly used methodologies for modelling electronic polarization are presented along with an overview of selected application studies. Models presented include induced point-dipoles, classical Drude oscillators, and fluctuating charge methods. The theoretical background of each method is followed by an introduction to extended Langrangian integrators required for computationally tractable molecular dynamics simulations using polarizable force fields. The remainder of the review focuses on application studies using these methods. Emphasis is placed on water models, for which numerous examples exist, with a more thorough discussion presented on the recently published models associated with the Drude-based CHARMM and the AMOEBA force fields. The utility of polarizable models for the study of ion solvation is then presented followed by an overview of studies of small molecules (e.g. CCl(4), alkanes, etc) and macromolecule (proteins, nucleic acids and lipid bilayers) application studies. The review is written with the goal of providing a general overview of the current status of the field and to facilitate future application and developments.