Salt Effect on Donnan Equilibrium in Montmorillonite Demonstrated with Molecular Dynamics Simulations

Donnan equilibrium governs the distribution of ions in many systems such as ion exchange membranes and biological cells in contact with an external electrolyte. Herein, Donnan equilibrium between bulk salt solution and bihydrated montmorillonite was investigated because such a system is of great importance for many nuclear waste disposal concepts. Specifically, we used molecular dynamics simulations to determine the partition coefficient of chloride, which was achieved by calculating the free-energy difference of chloride in the interlayer and the bulk using enhanced sampling methodology. Montmorillonite in equilibrium with either NaCl or CaCl2 was examined to elucidate the general difference between 1:1 and 2:1 salts. The concentration dependence of the partition coefficient for each salt was determined using three and four concentrations for NaCl and CaCl2, respectively. In the case of NaCl, we found that the partition coefficient increased linearly with the concentration, while for CaCl2, the increase was proportional to the square root of the concentration. A derivation of the partition coefficient using general Donnan theory that includes excess free energy contributions beyond the electrostatic Donnan potential is also presented. For both salts, the agreement between the partition coefficient from the simulations and Donnan theory was excellent. Although Donnan theory is a continuum theory derived without any reference to atomistic details, the present results justify its application to systems with nanoscale pores

<i>In Silico</i> Categorization of <i>in Vivo</i> Intrinsic Clearance Using Machine Learning

Machine learning has recently become popular and much used within the life science research domain, e.g., for finding quantitative structure–activity relationships (QSARs) between molecular structures and different biological end points. In the work presented here, we have applied orthogonal partial least-squares (OPLS), principal component analysis (PCA), and random forests (RF) methods for classification as well as regression analysis to a publicly available <i>in vivo</i> data set in order to assess the intrinsic metabolic clearance (CL_int) in humans. The derived classification models are able to identify compounds with CL_int lower and higher than 1500 mL/min, respectively, with nearly 80% accuracy. The most relevant descriptors are of lipophilicity and charge/polarizability types. Furthermore, the accuracy from a classification model based on regression analysis, using the 1500 mL/min cutoff, is also around 80%. These results suggest the usefulness of machine learning techniques to derive robust and predictive models in the area of <i>in vivo</i> ADMET (absorption, distribution, metabolism, elimination, and toxicity) modeling

The Central Role of Gln63 for the Hydrogen Bonding Network and UV–Visible Spectrum of the AppA BLUF Domain

In blue-light sensing using flavin (BLUF) domains, the side-chain orientation of key residues close to the flavin chromophore is still under debate. We report quantum refinements of the wild-type AppA BLUF protein from Rhodobacter sphaeroides starting from two published X-ray structures (1YRX and 2IYG) with different arrangements of the residues around the chromophore. Quantum refinement uses the same experimental X-ray raw data as conventional refinement, but includes data from quantum mechanics/molecular mechanics (QM/MM) calculations as restraints, which is expected to be more reliable than the normally employed MM data. In addition to quantum refinement, pure QM/MM geometry optimizations are performed for the 1YRX and 2IYG structures and for five models derived therefrom. Vertical excitation energies are computed at the QM­(DFT/MRCI)/MM level to assess the resulting structures. The experimental absorption maximum of the dark state of wild-type AppA is well reproduced for structures that contain the Gln63 residue in 1YRX-type orientation. The computed excitation energies are red-shifted for structures with a flipped Gln63 residue in 2IYG-type orientation. The calculated 1YRX- and 2IYG-type hydrogen-bonding networks are discussed in detail, particularly with regard to the orientation of the chromophore and the Gln63, Trp104, and Met106 residues

QM/MM Study of the Monomeric Red Fluorescent Protein DsRed.M1

We report a combined quantum mechanical/molecular mechanical (QM/MM) study of the DsRed.M1 protein using as QM component the self-consistent charge density functional tight-binding (SCC-DFTB) method in molecular dynamics (MD) simulations and hybrid density functional theory (DFT, B3LYP functional) in QM/MM geometry optimizations. We consider different variants of the chromophore (including the cis- and trans-acylimine and peptide forms) as well as different protonation states of environmental residues. The QM/MM calculations provide insight into the role of nearby residues concerning their interactions with the chromophore and their influence on structural and spectroscopic properties. QM/MM optimizations yield a single conformer for the anionic acylimine chromophore, whereas there are distinct cis- and trans-conformers in the anionic peptide chromophore, the latter being more stable. The calculated vertical excitation energies (DFT/MRCI) for the anionic chromophores agree well with experiment. The published crystal structure of DsRed.M1 with an anionic acylimine chromophore indicates a quinoid structure, while the QM/MM calculations predict the phenolate form to be more stable

Identification of the Peroxy Adduct in Multicopper Oxidases by a Combination of Computational Chemistry and Extended X-ray Absorption Fine-Structure Measurements

We have developed a computational method that combines extended X-ray absorption fine structure (EXAFS) refinements with the integrated quantum mechanical and molecular mechanics (QM/MM) method. This method allows us to obtain a structure of a metal site inside a protein that is compatible with both EXAFS data and QM calculations (i.e., that is chemically reasonable). Thereby, the QM/MM calculations play the same role as MM in nearly all NMR and crystallographic refinementsEXAFS ensures that the metal−ligand distances are accurate and QM/MM fills in all the other structural data. We have used this method to show that a structure with a peroxide ion in the center of the trinuclear cluster fits experimental EXAFS data better than a structure with the peroxide ion on the side of the cluster for the peroxide adduct of multicopper oxidases

A Pragmatic Approach Using First-Principle Methods to Address Site of Metabolism with Implications for Reactive Metabolite Formation

A majority of xenobiotics are metabolized by cytochrome P450 (CYP) enzymes. The discovery of drug candidates with low propensity to form reactive metabolites and low clearance can be facilitated by understanding CYP-mediated xenobiotic metabolism. Being able to predict the sites where reactive metabolites form is beneficial in drug design to produce drug candidates free of reactive metabolite issues. Herein, we report a pragmatic protocol using first-principle density functional theory (DFT) calculations for predicting sites of epoxidation and hydroxylation of aromatic substrates mediated by CYP. The method is based on the relative stabilities of the CYP-substrate intermediates or the substrate epoxides. Consequently, it concerns mainly the electronic reactivity of the substrates. Comparing to the experimental findings, the presented protocol gave excellent first-ranked epoxidation site predictions of 83%, and when the test was extended to CYP-mediated sites of aromatic hydroxylation, satisfactory results were also obtained (73%). This indicates that our assumptions are valid and also implies that the intrinsic reactivities of the substrates are in general more important than their binding poses in proteins, although the protocol may benefit from the addition of docking information

Role of Molecular, Crystal, and Surface Chemistry in Directing the Crystallization of Entacapone Polymorphs on the Au(111) Template Surface

The pharmaceutical compound entacapone ((E)-2-cyano-3-(3,4-dihydroxy-5-nitrophenyl)-N,N-diethylprop-2-enamide) is important in the treatment of Parkinson’s disease, exhibiting interesting polymorphic behavior upon crystallization from solution. It consistently produces its stable form A with a uniform crystal size distribution on the surface of an Au(111) template while concomitantly forming its metastable form D within the same bulk solution. Molecular modeling using empirical atomistic force-fields reveals more complex molecular and intermolecular structures for form D compared to form A, with the crystal chemistry of both polymorphs being dominated by van der Waals and π–π stacking interactions with lower contributions (ca. 20%) from hydrogen bonding and electrostatic interactions. Comparative lattice energies and convergence for the polymorphs are consistent with the observed concomitant polymorphic behavior. Synthon characterization reveals an elongated needle-like morphology for form D crystals in contrast to the more equant form A crystals with the surface chemistry of the latter exposing the molecules’ cyano groups on its {010} and {011} habit faces. Density functional theory modeling of surface adsorption reveals preferential interactions between Au and the synthon GA interactions of form A on the Au surface. Molecular dynamics modeling of the entacapone/gold interface reveals the entacapone molecular structure within the first adsorbed layer to show nearly identical interaction distances, for both the molecules within form A or D with respect to the Au surface, while in the second and third layers when entacapone molecule–molecule interactions overtake the interactions between those of molecule–Au, the intermolecular structures are found to be closer to the form A structure than form D. In these layers, synthon GA (form A) could be reproduced with just two small azimuthal rotations (5° and 15°) whereas the closest alignment to a form D synthon requires larger azimuthal rotations (15° and 40°). The cyano functional group interactions with the Au template dominate interfacial interactions with these groups being aligned parallel to the Au surface and with nearest neighbor distances to Au atoms more closely matching those in form A than form D. The overall polymorph direction pathway thus encompasses consideration of molecular, crystal, and surface chemistry factors

Cooperative Modes of Action of Antimicrobial Peptides Characterized with Atomistic Simulations: A Study on Cecropin B

Antimicrobial peptides (AMPs) are widely occurring host defense agents of interest as one route for addressing the growing problem of multidrug-resistant pathogens. Understanding the mechanisms behind their antipathogen activity is instrumental in designing new AMPs. Herein, we present an all-atom molecular dynamics and free energy study on cecropin B (CB) and its constituent domains. We find a cooperative mechanism in which CB inserts into an anionic model membrane with its amphipathic N-terminal segment, supported by the hydrophobic C-terminal segment of a second peptide. The two peptides interact via a Glu···Lys salt bridge and together sustain a pore in the membrane. Using a modified membrane composition, we demonstrate that when the lower leaflet is overall neutral, insertion of the cationic segment is retarded and thus this mode of action is membrane specific. The observed mode of action utilizes a flexible hinge, a common structural motif among AMPs, which allows CB to insert into the membrane using either or both termini. Data from both unbiased trajectories and enhanced sampling simulations indicate that a requirement for CB to be an effective AMP is the interaction of its hydrophobic C-terminal segment with the membrane. Simulations of these segments in isolation reveal their aggregation in the membrane and a different mechanism of supporting pore formation. Together, our results show the complex interaction of different structural motifs of AMPs and, in particular, a potential role for electronegative side chains in an overall cationic AMP