Efficient modelling protocols for oligosaccharides: from vacuum to solvent

The determination of conformational preferences of oligosaccharides is best approached by describing their preferred conformations on potential energy surfaces as a function of the glycosidic linkage φ, ψ torsional angles. For proper molecular mechanics modelling the flexibility of the rotatable pendant groups must also be considered. The so called adiabatic maps partially mimic the flexibility within the 10 dimensional conformational space of the pendant groups of the given disaccharide. These molecular mechanics maps are considered to be the state-of-the art of the φ, ψ potential energy surface of disaccharides recently calculated. The RAMM (RAndom Molecular Mechanics) method was shown to be able to calculate such profiles automatically. Additionally, based on the continuum solvent approach, RAMM allows the calculation of the effects of solvent on conformational energy profiles. Molecular dynamics simulations are also useful tools to study the influence of solvent on conformational behaviour of oligosaccharides. The capability of the RAMM calculational protocol to locate low-energy conformers on the multidimensional potential energy hypersurfaces of disaccharides is illustrated and compared with molecular dynamics simulations with and without inclusion of the solvent. This revised version was published online in August 2006 with corrections to the Cover Date.     

Conformational study of angiotensin II

Conformational free energy calculations using an empirical potential ECEPP/3 (Empirical Conformational Energy Program for Peptides, Version 3) were carried out on angiotensin II (AII) of sequence Asp-Arg-Val-Tyr-Ile-His-Pro-Phe to find the stable conformations of the free state in the unhydrated and the hydrated states. A conformational analysis of the unhydrated state was carried out using the buildup procedure. The free energy calculation using the hydration shell model was also carried out to obtain the stable conformation of the hydrated state. The calculated stable conformations of AII in both states have a partially right-handed α-helical structure stabilized by short- and medium-range interactions. The similarity between the lowest free energy conformations of the unhydrated and hydrated states suggests that the hydration might not be important to stabilize the overall conformation of AII in a free state. The absence of any intramolecular interaction of the Tyr side chain suggests the possible interaction of this residue with the receptor. In this study, we found that the low free energy conformations contain both the parallel-plate and the perpendicular-plate geometries of the His and Phe rings, suggesting the coexistence of both conformations. © 1996 John Wiley & Sons, Inc.     

Designing of benzothiazole derivatives as promising EGFR tyrosine kinase inhibitors: a pharmacoinformatics study

Abstract

Benzothiazole derivatives represent an important class of therapeutic chemical agents and are widely used for interesting biological activities and therapeutic functions including anticancer, antitumor and antimicrobial. In this study, we have performed similarity/substructure-based search of eMolecule database to find out promising benzothiazole derivatives as EGFR tyrosine kinase inhibitors. Several screening criteria that included molecular docking, pharmacokinetics and synthetic accessibility were used on initially derived about 7000 molecules consisting of benzothiazole as major component. Finally, four molecules were found to be promising EGFR tyrosine kinase inhibitors. The best docked pose of each molecule was considered for binding interactions followed by molecular dynamics (MD) and binding energy calculation. Molecular docking clearly showed the final proposed derivatives potential to form a number of binding interactions. MD simulation trajectories undoubtedly indicated that the EGFR protein becomes stable when proposed derivatives bind to the receptor cavity. Strong binding affinity was found for all molecules toward the EGFR which was substantiated by the binding energy calculation using the MM-PBSA approach. Therefore, proposed benzothiazole derivatives may be promising EGFR tyrosine kinase inhibitors for potential application as cancer therapy.

Communicated by Ramaswamy H. Sarma     

The role of the dielectric barrier in narrow biological channels: a novel composite approach to modeling single-channel currents

A composite continuum theory for calculating ion current through a protein channel of known structure is proposed, which incorporates information about the channel dynamics. The approach is utilized to predict current through the Gramicidin A ion channel, a narrow pore in which the applicability of conventional continuum theories is questionable. The proposed approach utilizes a modified version of Poisson-Nernst-Planck (PNP) theory, termed Potential-of-Mean-Force-Poisson-Nernst-Planck theory (PMFPNP), to compute ion currents. As in standard PNP, ion permeation is modeled as a continuum drift-diffusion process in a self-consistent electrostatic potential. In PMFPNP, however, information about the dynamic relaxation of the protein and the surrounding medium is incorporated into the model of ion permeation by including the free energy of inserting a single ion into the channel, i.e., the potential of mean force along the permeation pathway. In this way the dynamic flexibility of the channel environment is approximately accounted for. The PMF profile of the ion along the Gramicidin A channel is obtained by combining an equilibrium molecular dynamics (MD) simulation that samples dynamic protein configurations when an ion resides at a particular location in the channel with a continuum electrostatics calculation of the free energy. The diffusion coefficient of a potassium ion within the channel is also calculated using the MD trajectory. Therefore, except for a reasonable choice of dielectric constants, no direct fitting parameters enter into this model. The results of our study reveal that the channel response to the permeating ion produces significant electrostatic stabilization of the ion inside the channel. The dielectric self-energy of the ion remains essentially unchanged in the course of the MD simulation, indicating that no substantial changes in the protein geometry occur as the ion passes through it. Also, the model accounts for the experimentally observed saturation of ion current with increase of the electrolyte concentration, in contrast to the predictions of standard PNP theory.     

Structure‐guided screening of chemical database to identify NS3‐NS2B inhibitors for effective therapeutic application in dengue infection

Dengue infection is the most common arthropod‐borne disease caused by dengue viruses, predominantly affecting millions of human beings annually. To find out promising chemical entities for therapeutic application in Dengue, in the current research, a multi‐step virtual screening effort was conceived to screen out the entire “screening library” of the Asinex database. Initially, through “Lipinski rule of five” filtration criterion almost 0.6 million compounds were collected and docked with NS3‐NS2B protein. Thereby, the chemical space was reduced to about 3500 compounds through the analysis of binding affinity obtained from molecular docking study in AutoDock Vina. Further, the “Virtual Screening Workflow” (VSW) utility of Schrödinger suite was used, which follows a stepwise multiple docking programs such as ‐ high‐throughput virtual screening (HTVS), standard precision (SP), and extra precision (XP) docking, and in postprocessing analysis the MM‐GBSA based free binding energy calculation. Finally, five potent molecules were proposed as potential inhibitors for the dengue NS3‐NS2B protein based on the investigation of molecular interactions map and protein‐ligand fingerprint analyses. Different pharmacokinetics and drug‐likeness parameters were also checked, which favour the potentiality of selected molecules for being drug‐like candidates. The molecular dynamics (MD) simulation analyses of protein‐ligand complexes were explained that NS3‐NS2B bound with proposed molecules quite stable in dynamic states as observed from the root means square deviation (RMSD) and root means square fluctuation (RMSF) parameters. The binding free energy was calculated using MM‐GBSA method from the MD simulation trajectories revealed that all proposed molecules possess such a strong binding affinity towards the dengue NS3‐NS2B protein. Therefore, proposed molecules may be potential chemical components for effective inhibition of dengue NS3‐NS2B protein subjected to experimental validation.     

Computational characterization of the molecular structure and properties of Dye 7 for organic photovoltaics

Organic dyes have great potential for its use in solar cells. In this recent work, the molecular structure and properties of Dye 7 were obtained using density functional theory (DFT) and different levels of calculation. Upon comparing the molecular structure and the ultraviolet visible spectrum with experimental data reported in the literature, it was found that the M05-2X/6-31G(d) level of calculation gave the best approximation. Once the appropriate methodology had been obtained, the molecule was characterized by obtaining the infrared spectrum, dipole moment, total energy, isotropic polarizability, molecular orbital energies, free energy of solvation in different solvents, and the chemical reactivity sites using the condensed Fukui functions.     

Drug permeability prediction using PMF method

Drug permeability determines the oral availability of drugs via cellular membranes. Poor permeability makes a drug unsuitable for further development. The permeability may be estimated as the free energy change that the drug should overcome through crossing membrane. In this paper the drug permeability was simulated using molecular dynamics method and the potential energy profile was calculated with potential of mean force (PMF) method. The membrane was simulated using DPPC bilayer and three drugs with different permeability were tested. PMF studies on these three drugs show that doxorubicin (low permeability) should pass higher free energy barrier from water to DPPC bilayer center while ibuprofen (high permeability) has a lower energy barrier. Our calculation indicates that the simulation model we built is suitable to predict drug permeability.     

Discovery of new potential triplet acting inhibitor for Alzheimer’s disease via X-ray crystallography,molecular docking and molecular dynamics

Abstract

The most common brain disorder of late life is Alzheimer’s disease (AD), which is highly complicating dementia. There are several drug targets which are reported to control the severe level of AD; notably, acetylcholinesterase, β-Secretase and glycogen synthase kinase enzymes are approached as a good drug targets for AD. Hence, the present study mainly focused to discover newly synthesized molecule (7-propyl-6H-pyrano[3,2-c:5,6-c']dichromene-6,8(7H)-dione) as a potential triplet acting drug for above said enzymes through the analysis of X-ray crystallography, molecular docking, molecular dynamics and quantum chemical calculation. The target drug molecule was crystallized in the monoclinic crystal structure with P21/n space group. The structure was solved by SHELXS and refined by SHELXL. The crystal packing is stabilized by C???H···O type of interactions. Further, the induced fit docking shows that the molecule has high docking score, glide energy, favorable hydrogen bonding and hydrophobic interactions on the protein targets. The molecular dynamics simulation was performed to understand the stability of the molecule in the presence of active site environment. Finally, quantum chemical calculation has been carried out for the molecule in gas phase and for the corresponding molecule lifted from the active site region. The structural comparison between gas phase and active site helps to understand the conformational modification of the molecule in the active site.

Communicated by Ramaswamy H. Sarma     

Modeling Lanthanide Complexes: Towards the Theoretical Design of Light Conversion Molecular Devices

Theoretical techniques have been developed and/or improved to predict the molecular structure of lanthanide complexes which were used to calculate their electronic properties, in particular, their electronic spectra and energy levels necessary to calculate the rates of energy transfer from the ligands to the metal ion. The molecular structure has been obtained by the SMLC/AM1 (Sparkle Model for the Calculation of Lanthanide Complexes – Austin Model 1) model where the lanthanide ion is simulated by a sparkle implemented into the AM1 Hamiltonian used to perform a HF-SCF (Hartree-Fock Self-Consistent Field) calculation. The previous implementation of the SMLC/AM1 model (sparkle/1) involving only two parameters has been generalized to be consistent with the AM1 Hamiltonian and the new model (sparkle/2) significantly improved the prediction of molecular structures of Eu(III) complexes. For the electronic spectra and energy level calculations of the lanthanide complexes the model replaces the metal ion by a point charge with the ligands held in their positions as determined by the SMLC/AM1 model, and uses a INDO/S-CI (intermediate neglect of differential overlap/spectroscopic-configuration interaction) model. A preliminary study of the solvent effects on the absorption spectra of the free ligand is also presented. For the ligand-lanthanide ion energy transfer Fermi's golden rule is used with the multipolar and exchange mechanisms being implemented and tested for several complexes. These theoretical techniques have been applied to several complexes yielding very good results when compared to experimental data as well as predictions for the molecular and electronic structures and the relative contributions of the mechanisms for the energy transfer rates. This revised version was published online in June 2006 with corrections to the Cover Date.     

Molecular simulation of adsorption from dilute solutions

Adsorption of biomolecules on surfaces is a perennial and general challenge relevant to many fields in biotechnology. A change of the Helmholtz free energy DeltaA takes place when a molecule becomes adsorbed out of a bulk solution. The purpose of our investigations is to explore routes for the calculation of DeltaA by molecular simulations. DeltaA can be obtained both by integration over the mean force on a molecule and via the local density. It turns out that the route via the potential of mean force prevails over the latter due to better consistency. In this work we present results for systems of 1-centre and 2-centre Lennard-Jones mixtures at a 9/3 Lennard-Jones wall.     

Refined homology model of cytochrome bcc complex B subunit for virtual screening of potential anti-tuberculosis agents

Abstract

Cytochrome bcc complex is important for ATP synthesis and cellular activity, as a crucial step in the terminal reduction of oxygen in aerobic electron transport chains. The b subunit of cytochrome bcc complex (QcrB) has been reported as a promising anti-tuberculosis target, with many novel anti-tuberculosis scaffolds reported. However, the 3D structure of mycobacterium tuberculosis (M. tuberculosis) QcrB has not been released, making it hard to understand the interactions between QcrB and its inhibitors as well as to develop novel anti-tuberculosis scaffolds. Herein we built the optimal homology model of M. tuberculosis QcrB using the M. smegmatis QcrB structure as template, which was refined through all-atom molecular dynamics simulation. Then, the binding modes of known inhibitors were predicted through molecular docking method, along with molecular dynamics simulation and binding free energy calculation to verify the accuracy of docking results and stability of the protein-inhibitor complexes. The informative key residues within QcrB site enabled us to perform structure-based virtual library screening to obtain potential M. tuberculosis QcrB inhibitors, which were validated through molecular dynamics simulation and MM-GBSA calculation and analyzed through pharmacokinetic properties prediction. Our research would provide a deeper insight into the interactions between M. tuberculosis QcrB and its inhibitors, which boosts to develop novel therapy against tuberculosis.

Communicated by Ramaswamy H. Sarma     

The X-ray investigation of 16 alpha,17 alpha-thiazolidine derivatives of delta 4- and delta 5-pregnanes and conformational analysis of their oxathiolane analog

An X-ray study of 3,20-dioxo-4-pregnene-[16 alpha,17 alpha-d]--2',2'-dimethylthiazolidine (I) and 3 beta-hydroxy-20-oxo-5--pregnene-[16 alpha,17 alpha-d]-2',2'- dimethylthiazolidine (II) has been carried out. Two independent molecules in crystal II have significantly different conformations of the D and E rings, although according to the atom-atom potential calculations the energy of interaction of these molecules with their neighbors in crystal is the same. The calculation of conformational energy of 3,20-dioxo-4-pregnene-[17 alpha,16 alpha-d]-2',2'--dimethyloxathiolane (III) by the molecular mechanics method (MMM) indicates a possibility of existence of two similar conformers also for this molecule. The MMM calculation shows also that the conformation of molecule III (as well as progesterone) with the 17 beta-acetyl group torsion angle C(16)C(17)C(20)0(20) close to -120 degrees is possible.     

Mechanism of the plant cytochrome P450 for herbicide resistance: a modelling study

Plant cytochrome P450 is a key enzyme responsible for the herbicide resistance but the molecular basis of the mechanism is unclear. To understand this, four typical plant P450s and a widely resistant herbicide chlortoluron were analysed by carrying out homology modelling, molecular docking, molecular dynamics simulations and binding free energy analysis. Our results demonstrate that: (i) the putative hydrophobic residues located in the F-helix and polar residues in I-helix are critical in the herbicide resistance; (ii) the binding mode analysis and binding free energy calculation indicate that the distance between catalytic site of chlortoluron and heme of P450, as well as the binding affinity are key elements affecting the resistance for plants. In conclusion, this work provides a new insight into the interactions of plant P450s with herbicide from a molecular level, offering valuable information for the future design of novel effective herbicides which also escape from the P450 metabolism.     

Molecular mechanism of the enhanced virulence of 2009 pandemic Influenza A (H1N1) virus from D222G mutation in the hemagglutinin: a molecular modeling study

D222G mutation of the hemagglutinin (HA) is of special interest because of its close association with the enhanced virulence of 2009 pandemic influenza A (H1N1) virus through the increased binding affinity to α2,3-linked sialylated glycan receptors. However, there is still a lack of detailed understanding about the molecular mechanism of this enhanced virulence. Here, molecular dynamics simulation and binding free energy calculation were performed to explore the altered glycan receptor binding mechanism of HA upon the D222G mutation by studying the interaction of one α2,3-linked sialylglycan (sequence: SIA-GAL-NAG) with the wild type and D222G mutated HA. The binding free energy calculation based on the molecular mechanics generalized Born surface area (MM-GBSA) method indicates that the D222G mutated HA has a much stronger binding affinity with the studied α2,3-linked glycan than the wild type. This is consistent with the experimental result. The increased binding free energy of D222G mutant mainly comes from the increased energy contribution of Gln223. The structural analysis proves that the altered electrostatic potential of receptor binding domain (RBD) and the increased flexibility of 220-loop are the essential reasons leading to the increased affinity of HA to α2,3-linked sialic acid glycans. The obtained results of this study have allowed a deeper understanding of the receptor recognition mechanism and the pathogenicity of influenza virus, which will be valuable to the structure-based inhibitors design targeting influenza virus entry process.     

Theoretical Insights into Catalytic Mechanism of Protein Arginine Methyltransferase 1

Protein arginine methyltransferase 1 (PRMT1), the major arginine asymmetric dimethylation enzyme in mammals, is emerging as a potential drug target for cancer and cardiovascular disease. Understanding the catalytic mechanism of PRMT1 will facilitate inhibitor design. However, detailed mechanisms of the methyl transfer process and substrate deprotonation of PRMT1 remain unclear. In this study, we present a theoretical study on PRMT1 catalyzed arginine dimethylation by employing molecular dynamics (MD) simulation and quantum mechanics/molecular mechanics (QM/MM) calculation. Ternary complex models, composed of PRMT1, peptide substrate, and S-adenosyl-methionine (AdoMet) as cofactor, were constructed and verified by 30-ns MD simulation. The snapshots selected from the MD trajectory were applied for the QM/MM calculation. The typical S_N2-favored transition states of the first and second methyl transfers were identified from the potential energy profile. Deprotonation of substrate arginine occurs immediately after methyl transfer, and the carboxylate group of E144 acts as proton acceptor. Furthermore, natural bond orbital analysis and electrostatic potential calculation showed that E144 facilitates the charge redistribution during the reaction and reduces the energy barrier. In this study, we propose the detailed mechanism of PRMT1-catalyzed asymmetric dimethylation, which increases insight on the small-molecule effectors design, and enables further investigations into the physiological function of this family.     

Simulating the free energy of passive membrane permeation for small molecules

Abstract

Computer simulations of passive membrane permeation provide important microscopic insights into the molecular mechanism of this important biological process that are complementary to experimental data. Our review focuses on the main approaches for calculating the free energy, or potential of mean force, for permeation of small molecules through lipid bilayers. The theoretical background for most currently used methods for potential of mean force calculation is described, including particle insertion, thermodynamic integration, umbrella sampling, metadynamics, adaptive biasing force and milestoning. A brief comparison of strengths and weaknesses of the competing approaches is presented. This is followed by a survey of results obtained by the different methods, with special attention to describing the mechanistic insights generated by modelling and illustrating capabilities of the different techniques. We conclude with a discussion of recent advances and future directions in modelling membrane permeation, including latest methodological enhancements, consideration of multiple slow variables and memory effects.     

Characterisation of the hydrophobic collapse of polystyrene in water using free energy techniques

We characterise the hydrophobic collapse of single polystyrene chains in water using molecular dynamics simulations. Specifically, we calculate the potential of mean force for the collapse of a single polystyrene chain in water using metadynamics, comparing the results between all atomistic with coarse-grained (CG) molecular simulation. We next explore the scaling behaviour of the collapsed globular shape at the minimum energy configuration, characterised by the radius of gyration, as a function of chain length. The exponent is close to one third, consistent with that predicted for a polymer chain in bad solvent. We also explore the scaling behaviour of the solvent accessible surface area (SASA) as a function of chain length, finding a similar exponent for both all atomistic and CG simulations. Furthermore, calculation of the local water density as a function of chain length near the minimum energy configuration suggests that intermediate chain lengths are more likely to form dewetted states, as compared to shorter or longer chain lengths.     

Conformation of histidine model peptides. II. Spectroscopic properties of the imidazole chromophore.

Semi-empirical molecular orbital calculations have been carried out for the free base and cationic forms of imidazole so as to obtain data which are required for the calculation of the chiroptical properties of molecules that contain this chromophoric group. The polarization, energy, and monopolar charge distribution are reported for the lowest energy electronic transitions. The absorption spectra for imidazole have been determined to 180 nm and circular dichroism spectra for L -histidinol and L -2-amino-1-butanol have been measured.     

Conformational study of a nine residue fragment of the antigenic loop of foot-and-mouth disease virus.

The nine-residue peptide Ac-TASARGDLA-NHMe was selected as model peptide in order to understand the conformational features of the antigenic loop of foot-and-mouth disease virus (FMDV). A throughout exploration of the conformational space has been carried out by means of molecular dynamics (MD) and energy minimization. The calculations have been carried out using the AMBER force field. Solvent effects have been included by an effective dielectric constant of epsilon = 4r. The lowest energy conformation presents a secondary structure constituted by an alpha-helix at the N-terminal end followed by two gamma-turns in the central region. The rest of the accessible minima found present also a high tendency to form gamma-turns. Finally, a 100 ps MD trajectory calculation at 298 K suggest a stability of the secondary structure elements of the lowest energy conformation.     

A continuous analog for RNA folding

A linear segment in which a number of pairs of intervals of equal length are identified as potential stems is the subject of a folding problem analogous to inference of RNA secondary structure. A quantity of free energy (or equivalently, energy per unit length) is associated with each stem, and the various types of loops are assigned energy costs as a function of their lengths. Inference of stable structures can then be carried out in the same way as in RNA folding. More important, perturbation of stem lengths and energy densities (modelling various mutational processes affecting nucleotide sequences) allows the delineation of domains of stability of various foldings, through the explicit calculation of their boundaries, in a low-dimensional parameter space.