Dhilon S. Patel

Polarizable Empirical Force Field for Hexopyranose Monosaccharides Based on the Classical Drude Oscillator

A polarizable empirical force field based on the classical Drude oscillator is presented for the hexopyranose form of selected monosaccharides. Parameter optimization targeted quantum mechanical (QM) dipole moments, solute–water interaction energies, vibrational frequencies, and conformational energies. Validation of the model was based on experimental data on crystals, densities of aqueous-sugar solutions, diffusion constants of glucose, and rotational preferences of the exocylic hydroxymethyl of d-glucose and d-galactose in aqueous solution as well as additional QM data. Notably, the final model involves a single electrostatic model for all sixteen diastereomers of the monosaccharides, indicating the transferability of the polarizable model. The presented parameters are anticipated to lay the foundation for a comprehensive polarizable force field for saccharides that will be compatible with the polarizable Drude parameters for lipids and proteins, allowing for simulations of glycolipids and glycoproteins.

Divalent N(I) Compounds with Two Lone Pairs on Nitrogen

Carbon with the C(0) state has been reported recently, examples of which were known for the past decades. Silicon in the Si(0) state and phosphorus in the P(I) state are also known experimentally. This prompted us to search for divalent N(I) compounds, which resulted in the identification of ::N(←L)(2)(⊕) systems with bicoordinated nitrogen in the N(I) formal oxidation state. It was found that several biguanide derivatives (especially in their protonated state) belong to this class. Quantum chemical analysis provided the structural details, molecular orbitals, charge localization (vs delocalization) trends, etc. This class of compounds has been found to be characterized by two lone pairs on the central nitrogen, very similar to the central carbon in divalent C(0) compounds (::C(←L)(2)). The new bonding environment for nitrogen reported in this article, divalent nitrogen N(I), is clearly different from the nitrenium ions NR(2)(⊕). The electronic structure and reactivity of representative examples of this novel class of divalent nitrogen N(I) systems (::N(←L)(2)(⊕)) have been analyzed in detail, in terms of molecular orbitals, atomic charges, protonation energies, complexation energies with Lewis acids like BH(3), AlCl(3), and AuCl and compared with those of divalent C(0) systems.

To Bend or Not to Bend! The Dilemma of Allenes

Allenes, though expected to be linear, are found to be bent in many examples, especially in the case of cyclic allenes. The bending of allene is attributed to the weakening of π-bond strength across the allene. However, tetrakis(dimethylamino)allene {((CH(3))N)(2)C═C═C(N(CH(3))(2))}, which is characterized by push-push interactions, has been shown to be linear, thus leading to doubts of the current understanding of the bent allenes. In this article, we report the ab initio MO/DFT, NBO based electronic structure analysis of R(2)C═C═CR(2) (R = H, NH(2)) with a gradual increase in the number of amino substituents. The results indicate that the allenic π-bond strength and bending potential decrease, with an increase in the amino substitution. Molecular orbital analysis provides necessary clues regarding the delicate balance between orthogonality of the π orbitals and the p orbitals on the central carbon, which dictates the bending potential of the allenes. The dilemma of to bend or not to bend is a unique feature of tetraaminoallenes (NH(2))(2)C═C═C(NH(2))(2) in comparison to isoelectronic heteroallenes (NH(2))(2)C═N═C(NH(2))(2)(+) and (NH(2))(2)C═B═C(NH(2))(2)(-).

BamA POTRA Domain Interacts with a Native Lipid Membrane Surface

The outer membrane of Gram-negative bacteria is an asymmetric membrane with lipopolysaccharides on the external leaflet and phospholipids on the periplasmic leaflet. This outer membrane contains mainly β-barrel transmembrane proteins and lipidated periplasmic proteins (lipoproteins). The multisubunit protein β-barrel assembly machine (BAM) catalyzes the insertion and folding of the β-barrel proteins into this membrane. In Escherichia coli, the BAM complex consists of five subunits, a core transmembrane β-barrel with a long periplasmic domain (BamA) and four lipoproteins (BamB/C/D/E). The BamA periplasmic domain is composed of five globular subdomains in tandem called POTRA motifs that are key to BAM complex formation and interaction with the substrate β-barrel proteins. The BAM complex is believed to undergo conformational cycling while facilitating insertion of client proteins into the outer membrane. Reports describing variable conformations and dynamics of the periplasmic POTRA domain have been published. Therefore, elucidation of the conformational dynamics of the POTRA domain in full-length BamA is important to understand the function of this molecular complex. Using molecular dynamics simulations, we present evidence that the conformational flexibility of the POTRA domain is modulated by binding to the periplasmic surface of a native lipid membrane. Furthermore, membrane binding of the POTRA domain is compatible with both BamB and BamD binding, suggesting that conformational selection of different POTRA domain conformations may be involved in the mechanism of BAM-facilitated insertion of outer membrane β-barrel proteins.

Dynamics and Interactions of OmpF and LPS: Influence on Pore Accessibility and Ion Permeability

The asymmetric outer membrane of Gram-negative bacteria is formed of the inner leaflet with phospholipids and the outer leaflet with lipopolysaccharides (LPS). Outer membrane protein F (OmpF) is a trimeric porin responsible for the passive transport of small molecules across the outer membrane of Escherichia coli. Here, we report the impact of different levels of heterogeneity in LPS environments on the structure and dynamics of OmpF using all-atom molecular dynamics simulations. The simulations provide insight into the flexibility and dynamics of LPS components that are highly dependent on local environments, with lipid A being the most rigid and O-antigen being the most flexible. Increased flexibility of O-antigen polysaccharides is observed in heterogeneous LPS systems, where the adjacent O-antigen repeating units are weakly interacting and thus more dynamic, compared to homogeneous LPS systems in which LPS interacts strongly with each other with limited overall flexibility due to dense packing. The model systems were validated by comparing molecular-level details of interactions between OmpF surface residues and LPS core sugars with experimental data, establishing the importance of LPS core oligosaccharides in shielding OmpF surface epitopes recognized by monoclonal antibodies. There are LPS environmental influences on the movement of bulk ions (K(+) and Cl(-)), but the ion selectivity of OmpF is mainly affected by bulk ion concentration.

Conformational Properties of α- or β-(1→6)-Linked Oligosaccharides: Hamiltonian Replica Exchange MD Simulations and NMR Experiments

Conformational sampling for a set of 10 α- or β-(1→6)-linked oligosaccharides has been studied using explicit solvent Hamiltonian replica exchange (HREX) simulations and NMR spectroscopy techniques. Validation of the force field and simulation methodology is done by comparing calculated transglycosidic J coupling constants and proton–proton distances with the corresponding NMR data. Initial calculations showed poor agreement, for example, with >3 Hz deviation of the calculated 3J(H5,H6R) values from the experimental data, prompting optimization of the ω torsion angle parameters associated with (1→6)-linkages. The resulting force field is in overall good agreement (i.e., within ∼0.5 Hz deviation) from experimental 3J(H5,H6R) values, although some small limitations are evident. Detailed hydrogen bonding analysis indicates that most of the compounds lack direct intramolecular H-bonds between the two monosaccharides; however, minor sampling of the O6···HO2′ hydrogen bond is present in three compounds. The results verify the role of the gauche effect between O5 and O6 atoms in gluco- and manno-configured pyranosides causing the ω torsion angle to sample an equilibrium between the gt and gg rotamers. Conversely, galacto-configured pyranosides sample a population distribution in equilibrium between gt and tg rotamers, while the gg rotamer populations are minor. Water radial distribution functions suggest decreased accessibility to the O6 atom in the (1→6)-linkage as compared to the O6′ atom in the nonreducing sugar. The role of bridging water molecules between two sugar moieties on the distributions of ω torsion angles in oligosaccharides is also explored.

Divalent N(I) Character in 2-(Thiazol-2-yl)guanidine: An Electronic Structure Analysis

Several medicinally important compounds carry a 2-(thiazol-2-yl)guanidine unit. These species are generally (erroneously) represented as 1-(thiazol-2-yl)guanidine species. Quantum chemical studies were performed to identify the appropriate tautomeric state of this class of compounds. B3LYP/6-31+G(d) calculations indicate the preferred tautomeric state of these species is associated with the 2-(thiazol-2-yl)guanidine structure rather than the 1-(thiazol-2-yl)guanidine structure. G2MP2 calculations on the model system were carried out to study the electronic structure, electron delocalization, and protonation energy; MESP, ELF, HOMA, AIM, and NBO analyses were also carried out. The results indicate that this class of compounds may be treated as species with hidden ::N(←L)R character. Upon protonation of the thiazole ring nitrogen, these systems show the electronic structure as in ::N(←L)2(⊕) systems with divalent N(I) oxidation state.

Shape- and Chemical Feature-Based 3D-Pharmacophore Model Generation and Virtual Screening: Identification of Potential Leads for P. falciparum DHFR Enzyme Inhibition

Plasmodium falciparum dihydrofolate reductase (Pf DHFR) enzyme is one of the validated targets in the treatment of malaria using typical antifolates such as cycloguanil and pyrimethamine. However, point mutations at amino acid residues such as Ala16, Ile51, Cys59, Ser108 and Ile164 in the active site of the wild‐type enzyme resulted in a widespread resistance of the parasite to these drugs. Thus, design and discovery of new potential Pf DHFR inhibitors, equally active against both the wild‐type and mutant strains, is an urgent need. Catalyst software was used to generate a 3D pharmacophore query based on the bioactive conformation of WR99210 extracted from the X‐ray crystal structure of quadruple mutant PfDHFR enzyme. Validation criteria based on the experimentally determined conformation of WR99210 and its key interactions with the protein were considered to identify hits from two chemical databases, namely, NCI2000 and Maybridge2004 using different virtual filters. Virtual screening based on FlexX, GOLD and Glide docking programs resulted in a total of 73 hits. The hits reported in this article showed good potential to be inhibitors of the above Pf DHFRs based on their (i) best‐fit values (ii) binding scores (iii) binding modes and (iv) interactions with the key amino acid residues (Asp54, Ileu/Leu164, Asn/Ser108 and Ile14).

S-Oxidation of Thiazolidinedione with Hydrogen Peroxide, Peroxynitrous Acid, and C4a-Hydroperoxyflavin: A Theoretical Study

Quantum chemical analysis was carried out to model metabolism of glitazone class of drugs through oxygen transfer process to the sulfur atom of thiazolidinedione ring with different oxidants such as H(2)O(2), HOONO, and C4a-hydroperoxyflavin. Complete optimization (geometric and energy parameters) of all the required structures and transition states on the reaction path was carried out using MP2(full)/6-31+G(d,p). Charge and second-order delocalization analyses of important structures were carried out using the NBO method. The effect of solvent on the oxygen transfer to sulfur of thiazolidinedione was studied by including one, two, or three explicit water molecules. These calculations revealed that explicit solvent (water) effectively contributed in the sulfoxidation of thiazolidinedione and led to remarkable reduction in the energy barrier by ∼10 kcal/mol as compared to the gas phase. These results were found to be consistent with previously reported S-oxidation of dimethyl sulfide. When explicit water molecules were included, solvent molecules stabilize the charge separation at the transition state via specific interactions, and oxidation occurs via stretching of the O-O bond of oxidants and gradual formation of S-O bond. This study is helpful in understanding the metabolite generation due to the S-oxidation process in the glitazone series of antidiabetic drugs under physiological conditions.

Electronic structure and reactivity of guanylthiourea: A quantum chemical study

Electronic structure analysis of guanylthiourea (GTU) and its isomers has been carried out using quantum chemical methods. Two major tautomeric classes (thione and thiol) have been identified on the potential energy (PE) surface. In both the cases conjugation of pi‐electrons and intramolecular H‐bonds have been found to play a stabilizing role. Various isomers of GTU on its PE surface have been analyzed in two different groups (thione and thiol). The interconversion from the most stable thione conformer (GTU‐1) to the most stable thiol conformer (GTU‐t1) was found to take place via bimolecular process which involves protonation at sulfur atom of GTU‐1 followed by subsequent CN bond rotation and deprotonation. The detailed analysis of the protonation has been carried out in gas phase and aqueous phase (using CPMC model). Sulfur atom (S1) was found to be the preferred protonation site (over N4) in GTU‐1 in gas phase whereas N4 was found to be the preferred site of protonation in aqueous medium. The mechanism of S‐alkylation reaction in GTU has also been studied. The formation of alkylated analogs of thiol isomers (alkylated guanylthiourea) is believed to take place via bimolecular process which involves alkyl cation attack at S atom followed by CN bond rotation and deprotonation. The reactive intermediate RS(NH2)CNC(NH2)2+ belongs to the newly identified ⊕N(←L)2 class of species and provides the necessary dynamism for easy conversion of thione to thiol.