Conformational Search for the Building Block of Proteins Based on the Gradient Gravitational Search Algorithm (ConfGGS) Using Force Fields: CHARMM, AMBER, and OPLS-AA.

Proteins are linear polymers built from a repertoire of 20 different amino acids, which are considered building blocks of proteins. The diversity and versatility of these 20 building blocks with regard to their conformations are key to adopting three-dimensional structures that facilitate proteins to undergo important mechanistic biological processes in living systems. The present investigation reports a conformational search of 20 different amino acids, building blocks of proteins, using three different force fields, CHARMM, AMBER, and OPLS-AA, implemented in the gradient gravitational search algorithm. The search technique (ConfGGS) includes the contribution from both bonded and nonbonded terms using Cartesian coordinates. The efficiency of such conformational searches has also been compared with other optimization algorithms: DE/Best, DE/Rand, and PSO algorithms with respect to computational time and accuracy based on the minimum number of iteration steps and computed lowest mean absolute error (MAE) and mean standard deviation (MSD) values for dihedral angles of respective near-optimal structures. Moreover, the ConfGGS technique has also been extended to an ordered protein fragment (PQITL) extracted from HIV-1 protease (PDB ID: 1YTH), an intrinsically disordered protein fragment, i.e., an amyloid-forming segment (AVVTGVTAV), from the NAC domain of Parkinson's disease protein α-synuclein, residues 69-77 (PDB ID: 4RIK), the experimental NMR atomic-resolution structure of α-synuclein fibrils (PDB ID: 2N0A), and a disulfide bond-containing protein fragment sequence (PCYGWPVCY), residues 59-67 (PDB ID: 6Y4F) toward structure prediction as a close homologue compared with experimental accuracy, using the CHARMM force field. The MolProbity validation results for the protein fragment (PQITL) obtained by ConfGGS/CHARMM are in better agreement with the native protein fragment structure of HIV-1 protease (PDB ID: 1YTH). Furthermore, the computed results have also been compared with the coordinates obtained from the AlphaFold network.

Solvent Effect on the Potential Energy Surfaces of the F- + CH3CH2Br Reaction.

Although substantial work has been undertaken on reaction pathways involved in base-promoted elimination reactions and bimolecular nucleophilic substitution reaction of F- on CH3CH2X (X = Cl, Br, I), the effect of solvents with varying dielectric constants on the stereochemistry of each of the reaction species involved across the reaction profile have not yet been clearly understood. The present investigation reports the effect of solvents on the potential energy surfaces (PES) and structures of the species appearing in the reaction pathway of F- with bromoethane. The PESs in the gas phase have been computed at MP2 level and CCSD(T) level. The performance of several hybrid density functional, such as B3LYP, M06, M06L, BHandH, X3LYP, M05, M05-2X, and M06-2X have also been investigated toward describing the elimination and nucleophilic substitution reactions. With respect to MAE values and to make the computation cost-effective, we have explored the implicit continuum solvent model, CPCM in solvents like cyclohexane, methanol, acetonitrile, dimethyl sulfoxide and water. The reactant complexes proceed through the subsequent steps to produce fluoroethane as the substitution product and ethylene as one of the elimination products. For elimination reaction both syn and anti elimination have been explored. The calculated relatives energies values, which are negative in the gas phase, are found to be positive in polar solvents since the point charge in the separated reactants are more stabilized than the dispersed charge in the transient complex, which has also been analyzed through NBO analysis.

Gradient gravitational search: An efficient metaheuristic algorithm for global optimization.

The adaptation of novel techniques developed in the field of computational chemistry to solve the concerned problems for large and flexible molecules is taking the center stage with regard to efficient algorithm, computational cost and accuracy. In this article, the gradient-based gravitational search (GGS) algorithm, using analytical gradients for a fast minimization to the next local minimum has been reported. Its efficiency as metaheuristic approach has also been compared with Gradient Tabu Search and others like: Gravitational Search, Cuckoo Search, and Back Tracking Search algorithms for global optimization. Moreover, the GGS approach has also been applied to computational chemistry problems for finding the minimal value potential energy of two-dimensional and three-dimensional off-lattice protein models. The simulation results reveal the relative stability and physical accuracy of protein models with efficient computational cost.

Protocol for rational design of covalently interacting inhibitors.

The inhibition potencies of covalent inhibitors mainly result from the formation of a covalent bond to the enzyme during the inhibition mechanism. This class of inhibitors has essentially been ignored in previous target-directed drug discovery projects because of concerns about possible side effects. However, their advantages, such as higher binding energies and longer drug-target residence times moved them into the focus of recent investigations. While the rational design of non-covalent inhibitors became standard the corresponding design of covalent inhibitors is still in its early stages. Potent covalent inhibitors can be retrieved from large compound libraries by covalent docking approaches but protocols are missing that can reliably predict the influence of variations in the substitution pattern on the affinity and/or reactivity of a given covalent inhibitor. Hence, the wanted property profile can only be obtained from trial-and-error proceedings. This paper presents an appropriate protocol which is able to predict improved covalent inhibitors. It uses hybrid approaches, which mix quantum mechanical (QM) and molecular mechanical (MM) methods to predict variations in the reactivity of the inhibitor. They are also used to compute the required information about the non-covalent enzyme-inhibitor complex. Docking tools are employed to improve the inhibitor with respect to the non-covalent interactions formed in the binding site.

Spectroscopic probe on N-H⋯N, N-H⋯O and controversial C-H⋯O contact in A-T base pair: a DFT study.

DNA base pair A-T has been investigated by IR and NMR spectroscopy using DFT methods. The results have been analyzed in terms of infrared vibrational frequencies and (1)H NMR chemical shifts. Different types of interactions N-H⋯N, N-H⋯O and C-H⋯O types have been investigated in DNA base pairs. Although, previous reports argued about the third C-H⋯O type interaction in A-T base pair, such typical interaction has been analyzed thoroughly by IR and NMR spectroscopy using DFT methods. Our results show that the CH⋯O interaction in the A-T base pair is a weak interaction compared to normal hydrogen bond interactions.

Ion interactions with a new ditopic naphthalimide-based receptor: a photophysical, NMR and theoretical (DFT) study.

A new multi-component chemosensor system comprising a naphthalimide moiety as fluorophore is designed and developed to investigate receptor-analyte binding interactions in the presence of metal and non-metal ions. A dimethylamino moiety is utilized as receptor for metal ions and a thiourea receptor, having acidic protons, for binding anions. The system is characterized by conventional analytical methods. The absorption and fluorescence spectra of the system consist of a broad band typical for an intramolecular charge transfer (ICT). The effects of various metal-ion additives on the spectral behavior of the present sensor system are examined in acetonitrile. It is found that among the metal ions studied, alkali/alkaline earth-metal ions and transition-metal ions modulate the absorption and fluorescence spectra of the system. As an additional feature, the anion signaling behavior of the system in acetonitrile is studied. A decrease in fluorescence efficiency of the system is observed upon addition of fluoride and acetate anions. Fluorescence quenching is most effective in the case of fluoride ions. This is attributed to the enhancement of the photoinduced electron transfer from the anion receptor to the fluorophore moiety. Hydrogen-bond interactions between the acidic NH protons of the thiourea moiety and the F(-) anions are primarily attributed to the fluoride-selective signaling behavior. Interestingly, a negative cooperativity for the binding event is observed when the interactions of the system are studied in the presence of both Zn(2+) and F(-) ions. NMR spectroscopy and theoretical calculations are also carried out to better understand the receptor-analyte binding.

Quantum mechanical calculations for the misincorporation of nucleotides opposite mutagenic 3,N4-ethenocytosine.

The ubiquitous nature and persistence of exocyclic DNA adducts suggest their involvement as initiators of carcinogenesis. We have investigated the misincorporation properties of the exocyclic DNA adduct, 3,N(4)-ethenocytosine, using DFT and DFT-D methods. Computational investigations have been carried out by using the B3LYP, M062X, and wB97XD methods with the 6-31+G* basis set to determine the hydrogen bonding strengths, binding energy, and physical parameters. The single point energy calculations have been carried out at MP2/6-311++G** on corresponding optimized geometries. The energies were compared among the 3,N(4)-ethenocytosine adduct with DNA bases to find the most stable conformer. The solvent phase calculations have also been carried out using the CPCM model. The computed reaction enthalpy values provide computational insights to the earlier experimental observation in in vitro, E.coli, and mammalian cells of a high level of substitution mutation in which C → A transversion results from εC-T pairing [εC-T3 and εC-T4] in the adduct containing DNA sequence.

Model calculations for the misincorporation of nucleotides opposite five-membered exocyclic DNA adduct: N(2),3-ethenoguanine.

Five-membered exocyclic DNA adducts are biologically very significant because of their potential to block DNA replication and transcription. N(2),3-Ethenoguanine (N(2,3)-εG) has been identified in the liver DNA of vinyl chloride-exposed rats as a five-membered DNA adduct. Singer et al. ( Carcinogenesis 1987 , 8 , 745 - 747 ) reported that the misincorporation of thymine (T), with two hydrogen bonds to N(2,3)-εG, represents the mutagenic event. Although the base-pairing specificity and mode of misincorporation have been studied experimentally for the N(2),3-ethenoguanine adduct, molecular-level information is not yet clear. In this study, we have considered all four different DNA nucleotides paired with the N(2),3-ethenoguanine adduct for model calculations toward the determination of base-pairing specificity. To provide insight into the mutagenic process of DNA damage based on geometric characteristics and electronic properties, the B3LYP and M06 methods were employed for these model calculations. Single-point energy calculations at the MP2/6-311++G** level on the corresponding optimized geometries were also carried out to better estimate the hydrogen-bonding strengths. The polarizable conductor calculation model (CPCM), which accounts for the overall polarizability of the solvent, was also employed. The computed reaction enthalpy values lie in the order εG-G(2) (10.3 kcal/mol) > εG-G(4) (9.6 kcal/mol) > εG-T(4) (9.2 kcal/mol) > εG-G(1) (9.1 kcal/mol) > εG-A(5) (8.2 kcal/mol) > εG-C(2) (7.9 kcal/mol) at the M06 level, which indicates that guanine and thymine are most favorable for mispairing with the N(2),3-ethenoguanine adduct.

Interaction of serotonin and fluoxetine: toward understanding the importance of the chirality of fluoxetine (S form and R form).

The present investigation reports the importance of the S and R forms of fluoxetine, as an antidepressant with regards to the chirality taking different types of interactions associated with the neurotransmitter, serotonin. The goal of the present study is to provide predictions and to help experimental and theoretical studies toward understanding the chirality of fluoxetine in drug-ligand interaction, associated with serotonin-reuptake inhibitor drug design studies. Several different conformations for serotonin and fluoxetine complexes have been considered for quantum mechanical calculations. Both the S and R forms of fluoxetine associated with serotonin and fluoxetine complexes are found to be similar in total energy and binding energy values. The present study also supports the conformational effect of the 3-phenyl group of fluoxetine as stereo independent and is consistent with in vitro and in vivo data which indicates that the the eudismic ratio of fluoxetine enantiomers is near unity. The calculated highest stabilization energy values, binding energy values, both in the gas and aqueous phases at MP2/6-31+G*//B3LYP/6-31+G* identify the most possible stable conformer for the serotonin-fluoxetine complex.

Effect of microsolvation on zwitterionic glycine: an ab initio and density functional theory study.

The effect of microsolvation on zwitterionic glycine, considering both (-NH3(+)) as proton donor and (-COO(-)) as proton acceptor at correlated ab initio (MP2) level and density functional methods (B3LYP, PW91, MPW1PW91 and PBE) using 6-311++G** basis set has been reported. DFT methods have been employed so as to compare the performance/quality of different gradient-corrected correlation functionals (PW91, PBE), hybrid functionals (B3LYP, MPW1PW91) and to predict the near quantitative structural and vibrational properties, at reduced computational cost. B3LYP method outperforms among the different DFT methods for the computed hydrogen bond distances and found closer to the value obtained by correlated MP2 level, whereas MPW1PW91 and PBE methods shows very similar values but approximately 0.03 A less, compared to B3LYP method. MP2 calculation and single point CCSD(T)//MP2 calculation have been considered to decompose the interaction energy, including corrections for basis set superposition error (BSSE). Moreover, charge distribution analysis has also been carried out to understand the long raised questions, how and why the two body energies have significant contribution to the total binding energy.

Interaction of adenine adducts with thymine: a computational study.

The existence of DNA adducts bring the danger of carcinogenesis because of mispairing with normal DNA bases. 1,N6-ethenoadenine adducts (epsilonA) and 1,N6-ethanoadenine adducts (EA) have been considered as DNA adducts to study the interaction with thymine, as DNA base. Several different stable conformers for each type of adenine adduct with thymine, [epsilonA(1)-T(I), epsilonA(2)-T(I), epsilonA(3)-T(I) and EA(1)-T(I), EA(2)-T(I), EA(3)-T(I)] and [epsilonA(1)-T(II), epsilonA(2)-T(II), epsilonA(3)-T(II) and EA(1)-T(II), EA(2)-T(II), EA(3)-T(II)], have been considered with regard to their interactions. The differences in their geometrical structures, energetic properties, and hydrogen-bonding strengths have also been compared with Watson-Crick adenine-thymine base pair (A-T). Single-point energy calculations at MP2/6-311++G** levels on B3LYP/6-31+G* optimized geometries have also been carried out to better estimate the hydrogen-bonding strengths. The basis set superposition error corrected hydrogen-bonding strength sequence at MP2/6-311++G**//B3LYP/6-31+G* for the most stable complexes is found to be EA(2)-T(I) (15.30 kcal/mol) > EA(1)-T(II) (14.98 kcal/mol) > EA(3)-T(II) (14.68 kcal/mol) > epsilonA(2)-T(I) (14.54 kcal/mol) > epsilonA(3)-T(II) (14.22 kcal/mol) > epsilonA(3)-T(II) (13.64 kcal/mol) > A-T (13.62 kcal/mol). The calculated reaction enthalpy value for epsilonA(2)-T(I) is 10.05 kcal/mol, which is the highest among the etheno adduct-thymine complexes and about 1.55 kcal/mol more than those obtained for Watson-Crick A-T base pair and the reaction enthalpy value for EA(1)- T(II) is 10.22 kcal/mol, which is highest among the ethano addcut-thymine complexes and about 1.72 kcal/mol more than those obtained for Watson-Crick A-T base pair. The aim of this research is to provide fundamental understanding of adenine adduct and thymine interaction at the molecular level and to aid in future experimental studies toward finding the possible cause of DNA damage.

Hydrogen-bond interaction in 1:1 complexes of tetrahydrofuran with water, hydrogen fluoride, and ammonia: a theoretical study.

Ab initio and density-functional theory electronic structure calculations have been performed for the 1:1 complexes of tetrahydrofuran with water, hydrogen fluoride, and ammonia. Upon hydrogen bonding with H2O and HF, the structure of tetrahydrofuran (THF) remains relatively unchanged with the exception of THF sites involved in hydrogen-bonding interaction. But the similar findings are not true, upon hydrogen bonded with NH3, where the C2 symmetry of THF changed. The hydrogen-bonding strength for the 1:1 complexes of THF with water, HF, and NH3 is found to be in the order HF>H2O>NH3, which is well characterized by the order in bond angles O2H15F14, O2H16O14, and O2H15N14 closer to linearity, respectively, and the redshifted of stretching frequencies of upsilon(FH), upsilon(OH), and upsilon(NH), respectively. This work is an attempt to provide important predictions and to aid in future experimental and theoretical studies towards the understanding of such hydrogen-bonded van der Waals systems.

A density functional theory study for the hydrogen-bonded nucleic acid base pair: cytosine dimer.

Theoretical investigation for the geometric and energetic properties, rotational constants, harmonic vibrational frequencies, and binding energies of nucleic acid base pair, cytosine dimer, are carried out by using the density functional theory method. The dimer structures resulting from both the keto and the enol (cis/trans) tautomers are investigated in the present study. Various isomers are considered to find the stable structures of the cytosine dimer. The planar cytosine dimer, K-K3 with C2h symmetry, resulting from nonplanar keto tautomers, is found to be thermodynamically most stable out of the four different stable isomers and having the highest binding energy value, 19.51 kcal/mol (including basis set superposition error correction). The vibrational frequency analysis also suggests a red shift of 367.97 cm(-1) for the hydrogen-bonding K-K3 symmetric dimer with two hydrogen bond lengths, each of length 1.913 angstroms. Moreover, charge distribution (ChelpG charges), Laplacian electronic density distribution, and the dimerization equilibrium for the most stable dimer, K-K3, have also been investigated using the same method and the basis set.

Many-body interaction in glycine-(water)3 complex using density functional theory method.

Various configurations were investigated to find the most stable structures of glycine-(water)3 complex. Five different optimized conformers of glycine-(water)3 complex are obtained from density functional theory calculations using 6-311++G* basis set. Relaxation energy and many body interaction energies (two, three, and four body) are also calculated for these conformers. Out of the five conformers, the most stable conformer has the BSSE corrected total energy -513.917 967 7 Hartree and binding energy -27.28 Kcal/mol. It has been found that the relaxation energies, two body energies and three body energies have significant contribution to the total binding energy whereas four body energies are very small. The chemical hardness and chemical potential also confirmed the stability of the conformer having lowest total energy.