Evaluating thermodynamic integration performance of the new amber molecular dynamics package and assess potential halogen bonds of enoyl‐ACP reductase (FabI) benzimidazole inhibitors

Thermodynamic integration (TI) can provide accurate binding free energy insights in a lead optimization program, but its high computational expense has limited its usage. In the effort of developing an efficient and accurate TI protocol for FabI inhibitors lead optimization program, we carefully compared TI with different Amber molecular dynamics (MD) engines (sander and pmemd), MD simulation lengths, the number of intermediate states and transformation steps, and the Lennard‐Jones and Coulomb Softcore potentials parameters in the one‐step TI, using eleven benzimidazole inhibitors in complex with Francisella tularensis enoyl acyl reductase (FtFabI). To our knowledge, this is the first study to extensively test the new AMBER MD engine, pmemd, on TI and compare the parameters of the Softcore potentials in the one‐step TI in a protein‐ligand binding system. The best performing model, the one‐step pmemd TI, using 6 intermediate states and 1 ns MD simulations, provides better agreement with experimental results (RMSD = 0.52 kcal/mol) than the best performing implicit solvent method, QM/MM‐GBSA from our previous study (RMSD = 3.00 kcal/mol), while maintaining similar efficiency. Briefly, we show the optimized TI protocol to be highly accurate and affordable for the FtFabI system. This approach can be implemented in a larger scale benzimidazole scaffold lead optimization against FtFabI. Lastly, the TI results here also provide structure‐activity relationship insights, and suggest the parahalogen in benzimidazole compounds might form a weak halogen bond with FabI, which is a well‐known halogen bond favoring enzyme. © 2015 Wiley Periodicals, Inc.     

Large‐scale asynchronous and distributed multidimensional replica exchange molecular simulations and efficiency analysis

We describe methods to perform replica exchange molecular dynamics (REMD) simulations asynchronously (ASyncRE). The methods are designed to facilitate large scale REMD simulations on grid computing networks consisting of heterogeneous and distributed computing environments as well as on homogeneous high‐performance clusters. We have implemented these methods on NSF (National Science Foundation) XSEDE (Extreme Science and Engineering Discovery Environment) clusters and BOINC (Berkeley Open Infrastructure for Network Computing) distributed computing networks at Temple University and Brooklyn College at CUNY (the City University of New York). They are also being implemented on the IBM World Community Grid. To illustrate the methods, we have performed extensive (more than 60 ms in aggregate) simulations for the beta‐cyclodextrin‐heptanoate host‐guest system in the context of one‐ and two‐dimensional ASyncRE, and we used the results to estimate absolute binding free energies using the binding energy distribution analysis method. We propose ways to improve the efficiency of REMD simulations: these include increasing the number of exchanges attempted after a specified molecular dynamics (MD) period up to the fast exchange limit and/or adjusting the MD period to allow sufficient internal relaxation within each thermodynamic state. Although ASyncRE simulations generally require long MD periods (>picoseconds) per replica exchange cycle to minimize the overhead imposed by heterogeneous computing networks, we found that it is possible to reach an efficiency similar to conventional synchronous REMD, by optimizing the combination of the MD period and the number of exchanges attempted per cycle. © 2015 Wiley Periodicals, Inc.     

Supramolecular study,Hirshfeld analysis and theoretical study of 6‐methoxyquinoline N‐oxide dihydrate

In the crystal structure of 6‐methoxyquinoline N‐oxide dihydrate, C₁₀H₉NO₂·2H₂O, (I), the presence of two‐dimensional water networks is analysed. The water molecules form unusual water channels, as well as two intersecting mutually perpendicular columns. In one of these channels, the O atom of the N‐oxide group acts as a bridge between the water molecules. The other channel is formed exclusively by water molecules. Confirmation of the molecular packing was performed through the analysis of Hirshfeld surfaces, and (I) is compared with other similar isoquinoline systems. Calculations of bond lengths and angles by the Hartree–Fock method or by density functional theory B3LYP, both with 6‐311++G(d,p) basis sets, are reported, together with the results of additional IR, UV–Vis and theoretical studies.     

GridMAT‐MD: A grid‐based membrane analysis tool for use with molecular dynamics

GridMAT‐MD is a new program developed to aid in the analysis of lipid bilayers from molecular dynamics simulations. It reads a GROMACS coordinate file and generates two types of data: a two‐dimensional contour plot depicting membrane thickness, and a polygon‐based tessellation of the individual lipid headgroups. GridMAT‐MD can also account for proteins or small molecules within the headgroups of the lipids, closely approximating their occupied lateral area. The program requires no installation, is fast, and is freely available. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Enhanced sampling method in molecular simulations using genetic algorithm for biomolecular systems

We propose a molecular simulation method using genetic algorithm (GA) for biomolecular systems to obtain ensemble averages efficiently. In this method, we incorporate the genetic crossover, which is one of the operations of GA, to any simulation method such as conventional molecular dynamics (MD), Monte Carlo, and other simulation methods. The genetic crossover proposes candidate conformations by exchanging parts of conformations of a target molecule between a pair of conformations during the simulation. If the candidate conformations are accepted, the simulation resumes from the accepted ones. While conventional simulations are based on local update of conformations, the genetic crossover introduces global update of conformations. As an example of the present approach, we incorporated genetic crossover to MD simulations. We tested the validity of the method by calculating ensemble averages and the sampling efficiency by using two kinds of peptides, ALA3 and (AAQAA)₃. The results show that for ALA3 system, the distribution probabilities of backbone dihedral angles are in good agreement with those of the conventional MD and replica-exchange MD simulations. In the case of (AAQAA)₃ system, our method showed lower structural correlation of α-helix structures than the other two methods and more flexibility in the backbone ψ angles than the conventional MD simulation. These results suggest that our method gives more efficient conformational sampling than conventional simulation methods based on local update of conformations. © 2018 Wiley Periodicals, Inc.     

Two‐Dimensional Infrared Spectroscopy Reveals the Structure of an Evans Auxiliary Derivative and Its SnCl4 Lewis Acid Complex

Determining the structure of reactive intermediates is the key to understanding reaction mechanisms. To access these structures, a method combining structural sensitivity and high time resolution is required. Here ultrafast polarization‐dependent two‐dimensional infrared (P2D‐IR) spectroscopy is shown to be an excellent complement to commonly used methods such as one‐dimensional IR and multidimensional NMR spectroscopy for investigating intermediates. P2D‐IR spectroscopy allows structure determination by measuring the angles between vibrational transition dipole moments. The high time resolution makes P2D‐IR spectroscopy an attractive method for structure determination in the presence of fast exchange and for short‐lived intermediates. The ubiquity of vibrations in molecules ensures broad applicability of the method, particularly in cases in which NMR spectroscopy is challenging due to a low density of active nuclei. Here we illustrate the strengths of P2D‐IR by determining the conformation of a Diels–Alder dienophile that carries the Evans auxiliary and its conformational change induced by the complexation with the Lewis acid SnCl₄, which is a catalyst for stereoselective Diels–Alder reactions. We show that P2D‐IR in combination with DFT computations can discriminate between the various conformers of the free dienophile N‐crotonyloxazolidinone that have been debated before, proving antiperiplanar orientation of the carbonyl groups and s‐cis conformation of the crotonyl moiety. P2D‐IR unequivocally identifies the coordination and conformation in the catalyst–substrate complex with SnCl₄, even in the presence of exchange that is fast on the NMR time scale. It resolves a chelate with the carbonyl orientation flipped to synperiplanar and s‐cis crotonyl configuration as the main species. This work sets the stage for future studies of other catalyst–substrate complexes and intermediates using a combination of P2D‐IR spectroscopy and DFT computations.     

Application of a genomic model for high‐dimensional chemometric analysis

The rapid development of new technologies for large‐scale analysis of genetic variation in the genomes of individuals and populations has presented statistical geneticists with a grand challenge to develop efficient methods for identifying the small proportion of all identified genetic polymorphisms that have effects on traits of interest. To address such a “large p small n” problem, we have developed a heteroscedastic effects model (HEM) that has been shown to be powerful in high‐throughput genetic analyses. Here, we describe how this whole‐genome model can also be utilized in chemometric analysis. As a proof of concept, we use HEM to predict analyte concentrations in silage using Fourier transform infrared spectroscopy signals. The results show that HEM often outperforms the classic methods and in addition to this presents a substantial computational advantage in the analyses of such high‐dimensional data. The results thus show the value of taking an interdisciplinary approach to chemometric analysis and indicate that large‐scale genomic models can be a promising new approach for chemometric analysis that deserve to be evaluated more by experts in the field. The software used for our analyses is freely available as an R package at http://cran.r‐project.org/web/packages/bigRR/ . Copyright © 2014 John Wiley & Sons, Ltd.     

Comparison of radii sets,entropy, QM methods,and sampling on MM‐PBSA,MM‐GBSA,and QM/MM‐GBSA ligand binding energies of F. tularensis enoyl‐ACP reductase (FabI)

To validate a method for predicting the binding affinities of FabI inhibitors, three implicit solvent methods, MM‐PBSA, MM‐GBSA, and QM/MM‐GBSA were carefully compared using 16 benzimidazole inhibitors in complex with Francisella tularensis FabI. The data suggests that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length, if only one simulation trajectory is used for each ligand. In this case, QM/MM‐GBSA using 6 ns MD simulation trajectories together with GB^neck2, PM3, and the mbondi2 radii set, generate the closest agreement with experimental values (r² = 0.88). However, if the three implicit solvent methods are averaged from six 1 ns MD simulations for each ligand (called “multiple independent sampling”), the prediction results are relatively insensitive to all the tested parameters. Moreover, MM/GBSA together with GB^HCT and mbondi, using 600 frames extracted evenly from six 0.25 ns MD simulations, can also provide accurate prediction to experimental values (r² = 0.84). Therefore, the multiple independent sampling method can be more efficient than a single, long simulation method. Since future scaffold expansions may significantly change the benzimidazole's physiochemical properties (charges, etc.) and possibly binding modes, which may affect the sensitivities of various parameters, the relatively insensitive “multiple independent sampling method” may avoid the need of an entirely new validation study. Moreover, due to large fluctuating entropy values, (QM/)MM‐P(G)BSA were limited to inhibitors’ relative affinity prediction, but not the absolute affinity. The developed protocol will support an ongoing benzimidazole lead optimization program. © 2015 Wiley Periodicals, Inc.     

Influence of Variation of a Side Chain on the Folding Equilibrium of a β‐Peptide

The ability to design well‐folding β‐peptides with a specific biological activity requires detailed insight into the relationship between the β‐amino acid sequence and the three‐dimensional structure of the peptide. Here, we present a molecular‐dynamics (MD) study of the influence of a variation of a side chain on the folding equilibrium of a β‐heptapeptide that folds into a 3₁₄‐helical structure. The side chain of the 5th residue, a valine, was changed into five differently branched side chains of different lengths and polarity, Ala, Leu, Ile, Ser, and Thr. Two computational techniques, long‐time MD simulations and the one‐step perturbation method, were used to obtain free enthalpies of folding. The simulations show that all six peptides exhibit similar folding behavior, and that their dominant fold is the same, i.e., a 3₁₄‐helix. Despite the similarities of their structural properties, a small stabilization effect of ca. 2 kJ mol^?1 on the folding equilibrium of the 3₁₄‐helical structure due to a branching C_γ‐atom in the β³‐side chain is observed. These results confirm those of previous circular dichroism (CD) studies. The length of side chain and its polarity seem to have no apparent (de)stabilization effect. Application of the cost‐effective one‐step perturbation method to predict free‐enthalpy differences appeared to yield an overall accuracy of about k_BT, which is not sufficient to detect the small stabilization effect.     

Comparing distance metrics for rotation using the k‐nearest neighbors algorithm for entropy estimation

Distance metrics facilitate a number of methods for statistical analysis. For statistical mechanical applications, it is useful to be able to compute the distance between two different orientations of a molecule. However, a number of distance metrics for rotation have been employed, and in this study, we consider different distance metrics and their utility in entropy estimation using the k‐nearest neighbors (KNN) algorithm. This approach shows a number of advantages over entropy estimation using a histogram method, and the different approaches are assessed using uniform randomly generated data, biased randomly generated data, and data from a molecular dynamics (MD) simulation of bulk water. The results identify quaternion metrics as superior to a metric based on the Euler angles. However, it is demonstrated that samples from MD simulation must be independent for effective use of the KNN algorithm and this finding impacts any application to time series data. © 2013 Wiley Periodicals, Inc.     

3D anisotropy measurement methodology for surface microstructures

We propose a simple method able to quantify three‐dimensional (3D) anisotropy in topographic microstructure measurements. The anisotropy quantification of 3D surface data is based on a horizontal‐cut analysis yielding level sets. For a specific level set, we study xy anisotropy with the use of a directional chord length analysis. In the height profile direction, the z direction, anisotropy is quantified using the framework of mean Euler characteristic and more generally that of Minkowski functionals. We exemplify the method with the use of 3D heights profiles from atomic force microscopy measurements. Copyright © 2012 John Wiley & Sons, Ltd.     

Characterization of stereoisomeric 5‐(2‐nitro‐1‐phenylethyl)furan‐2(5H)‐ones by computation of 1H and 13C NMR chemical shifts and electronic circular dichroism spectra

In the present work, we describe the preparation of two diastereomers from the enantioselective Michael addition of furan‐2(5H)‐one to (E)‐(2‐nitrovinyl)benzene catalyzed by a dinuclear Zn‐complex. The relative configurations of the diastereomeric products were assigned by comparing nuclear magnetic resonance (NMR) experimental chemical shift data with those computed by density functional theory (DFT) methods. Corrected mean absolute error (CMAE) and CP3 analyses were used to compare the data sets. The absolute configuration of each diastereomer was initially assigned by analysis of electronic circular dichroism (ECD) data, which was consistent with that of the known X‐ray crystallographic structure of the product of a related reaction, namely, (R)‐5‐((R)‐1‐(4‐chlorophenyl)‐2‐nitroethyl)furan‐2(5H)‐one.     

Microsecond simulations of DNA and ion transport in nanopores with novel ion–ion and ion–nucleotides effective potentials

We developed a novel scheme based on the grand‐canonical Monte Carlo/Brownian dynamics simulations and have extended it to studies of ion currents across three nanopores with the potential for single‐stranded DNA (ssDNA) sequencing: solid‐state nanopore Si₃N₄, α‐hemolysin, and E111N/M113Y/K147N mutant. To describe nucleotide‐specific ion dynamics compatible with ssDNA coarse‐grained model, we used the inverse Monte Carlo protocol, which maps the relevant ion–nucleotide distribution functions from all‐atom molecular dynamics (MD) simulations. Combined with the previously developed simulation platform for Brownian dynamics simulations of ion transport, it allows for microsecond‐ and millisecond‐long simulations of ssDNA dynamics in the nanopore with a conductance computation accuracy that equals or exceeds that of all‐atom MD simulations. In spite of the simplifications, the protocol produces results that agree with the results of previous studies on ion conductance across open channels and provide direct correlations with experimentally measured blockade currents and ion conductances that have been estimated from all‐atom MD simulations. © 2014 Wiley Periodicals, Inc.     

Optimizing separations in online comprehensive two‐dimensional liquid chromatography

Online comprehensive two‐dimensional liquid chromatography has become an attractive option for the analysis of complex nonvolatile samples found in various fields (e.g. environmental studies, food, life, and polymer sciences). Two‐dimensional liquid chromatography complements the highly popular hyphenated systems that combine liquid chromatography with mass spectrometry. Two‐dimensional liquid chromatography is also applied to the analysis of samples that are not compatible with mass spectrometry (e.g. high‐molecular‐weight polymers), providing important information on the distribution of the sample components along chemical dimensions (molecular weight, charge, lipophilicity, stereochemistry, etc.). Also, in comparison with conventional one‐dimensional liquid chromatography, two‐dimensional liquid chromatography provides a greater separation power (peak capacity). Because of the additional selectivity and higher peak capacity, the combination of two‐dimensional liquid chromatography with mass spectrometry allows for simpler mixtures of compounds to be introduced in the ion source at any given time, improving quantitative analysis by reducing matrix effects. In this review, we summarize the rationale and principles of two‐dimensional liquid chromatography experiments, describe advantages and disadvantages of combining different selectivities and discuss strategies to improve the quality of two‐dimensional liquid chromatography separations.     

Cerium triiodate,Ce(IO3)3

The crystal structure of Ce(IO₃)₃ consists of one‐dimensional chains of edge‐sharing CeO₉ polyhedra which are crosslinked into two‐dimensional layers through bridging IO₃⁻ groups. The layers are held together via long I⋯O contacts, resulting in an extended three‐dimensional network. The I—O bond distances and O—I—O angles are normal, lying in the ranges 1.806 (4)–1.846 (4) Å and 89.9 (2)–100.9 (2)°, respectively. The three crystallographically independent iodate groups all show different coordination modes.     

Long‐Range Residual Dipolar Couplings: A Tool for Determining the Configuration of Small Molecules

Together with NOE and J coupling, one‐bond residual dipolar coupling (RDC), which reports on the three‐dimensional orientation of an internuclear vector in the molecular frame, plays an important role in the conformation and configuration analysis of small molecules in solution by NMR spectroscopy. When the molecule has few C? H bonds, or too many bonds are in parallel, the available RDCs may not be sufficient to obtain the alignment tensor used for structure elucidation. Long‐range RDCs that connect nuclei over multiple bonds are normally not parallel to the single bonds and therefore complement one‐bond RDCs. Herein we present a method for extracting the long‐range RDC of a chosen proton or group of protons to all remotely connected carbon atoms, including non‐protonated carbon atoms. Alignment tensors fitted directly to the total long‐range couplings (T=J+D) enabled straightforward analysis of both the long‐range and one‐bond RDCs for strychnine.     

Diaquadichloridobis(1H‐imidazole)manganese(II) at 100 K

The mononuclear title complex, [MnCl₂(C₃H₄N₂)₂(H₂O)₂], is located on a crystallographic inversion center. The Mn^II ion is coordinated by two imidazole ligands [Mn—N = 2.2080 (9) Å], two Cl atoms [Mn—Cl = 2.5747 (3) Å] and two water molecules [Mn—O = 2.2064 (8) Å]. These six monodentate ligands define an octahedron with almost ideal angles: the adjacent N—Mn—O, N—Mn—Cl and O—Mn—Cl angles are 90.56 (3), 92.04 (2) and 90.21 (2)°, respectively. Hydrogen bonds between the coordinated water molecules and Cl atoms form a two‐dimensional network parallel to (100) involving R₄²(8) rings. The two‐dimensional networks link into a three‐dimensional framework through weaker N—H...Cl interactions. Thermogravimetric analysis results are in accordance with the water‐coordinated character of the substance and its dehydration in two successive steps.