Structure-based virtual screening offers a good opportunity for the discovery of selective M1 muscarinic acetylcholine receptor (mAChR) agonists for the treatment of Alzheimer's disease. However, no 3-D structure of an M1 mAChR is yet available and the homology models that have been previously reported are only able to identify antagonists in virtual screening experiments. In this study, we generated a homology model of the human M1 mAChR, based on the crystal structure of an M3 mAChR as the template. This initial model was modified, using the agonist-bound crystal structure of a β2-adrenergic receptor as a guide, to give two possible activated structures. The T192 side chain was adjusted in both structures and one of the structures also had the whole of transmembrane (TM) 5 rotated and tilted toward the inner channel of the transmembrane region. The binding sites of all three structures were then refined by induced-fit docking (IFD) with acetylcholine. Virtual screening experiments showed that all three refined models could efficiently differentiate agonists from decoy molecules, with the TM5-modified models also giving good agonist/antagonist selectivity. The whole range of agonists and antagonists was observed to bind within the orthosteric site of the structure obtained by IFD refinement alone, implying that it has inactive state character. In contrast, the two TM5-modified structures were unable to accommodate the antagonists, supporting the proposition that they possess activated state character.
The utilization of cold active lipases in organic solvents proves an excellent approach for chiral synthesis and modification of fats and oil due to the inherent flexibility of lipases under low water conditions. In order to verify whether this lipase can function as a valuable synthetic catalyst, the mechanism concerning activation of the lid and interacting solvent residues in the presence of organic solvent must be well understood. A new alkaline cold-adapted lipase, AMS8, from Pseudomonas fluorescens was studied for its structural adaptation and flexibility prior to its exposure to non-polar, polar aprotic and protic solvents. Solvents such as ethanol, toluene, DMSO and 2-propanol showed to have good interactions with active sites. Asparagine (Asn) and tyrosine (Tyr) were key residues attracted to solvents because they could form hydrogen bonds. Unlike in other solvents, Phe-18, Tyr-236 and Tyr-318 were predicted to have aromatic-aromatic side-chain interactions with toluene. Non-polar solvent also was found to possess highest energy binding compared to polar solvents. Due to this circumstance, the interaction of toluene and AMS8 lipase was primarily based on hydrophobicity and molecular recognition. The molecular dynamic simulation showed that lid 2 (residues 148-167) was very flexible in toluene and Ca(2+). As a result, lid 2 moves away from the catalytic areas, leaving an opening for better substrate accessibility which promotes protein activation. Only a single lid (lid 2) showed the movement following interactions with toluene, although AMS8 lipase displayed double lids. The secondary conformation of AMS8 lipase that was affected by toluene observed a reduction of helical strands and increased coil structure. Overall, this work shows that cold active lipase, AMS8 exhibits distinguish interfacial activation and stability in the presence of polar and non-polar solvents.
Dehalogenases continue to garner interest of the scientific community due to their potential applications in bioremediation of halogen-contaminated environment and in synthesis of various industrially relevant products. Example of such enzymes is DehL, an L-2-haloacid dehalogenase (EC 3.8.1.2) from Rhizobium sp. RC1 that catalyses the specific cleavage of halide ion from L-2-halocarboxylic acids to produce the corresponding D-2-hydroxycarboxylic acids. Recently, the catalytic residues of DehL have been identified and its catalytic mechanism has been fully elucidated. However, the enantiospecificity determinants of the enzyme remain unclear. This information alongside a well-defined catalytic mechanism are required for rational engineering of DehL for substrate enantiospecificity. Therefore, using quantum mechanics/molecular mechanics and molecular mechanics Poisson-Boltzmann surface area calculations, the current study theoretically investigated the molecular basis of DehL enantiospecificity. The study found that R51L mutation cancelled out the dehalogenation activity of DehL towards it natural substrate, L-2-chloropropionate. The M48R mutation, however introduced a new activity towards D-2-chloropropionate, conveying the possibility of inverting the enantiospecificity of DehL from L-to d-enantiomer with a minimum of two simultaneous mutations. The findings presented here will play important role in the rational design of DehL dehalogenase for improving substrate utility.
Dehalogenase E (DehE) is a non-stereospecific enzyme produced by the soil bacterium, Rhizobium sp. RC1. Till now, the catalytic mechanism of DehE remains unclear although several literature concerning its structure and function are available. Since DehE is non-stereospecific, the enzyme was hypothesized to follow a 'direct attack mechanism' for the catalytic breakdown of a haloacid. For a molecular insight, the DehE modelled structure was docked in silico with the substrate 2-chloropropionic acid (2CP) in the active site. The ideal position of DehE residues that allowed a direct attack mechanism was then assessed via molecular dynamics (MD) simulation. It was revealed that the essential catalytic water was hydrogen bonded to the 'water-bearer', Asn114, at a relatively constant distance of ∼2.0 Å after 50 ns. The same water molecule was also closely sited to the catalytic Asp189 at an average distance of ∼2.0 Å, signifying the imperative role of the latter to initiate proton abstraction for water activation. This reaction was crucial to promote a direct attack on the α-carbon of 2CP to eject the halide ion. The water molecule was oriented favourably towards the α-carbon of 2CP at an angle of ∼75°, mirrored by the formation of stable enzyme-substrate orientations throughout the simulation. The data therefore substantiated that the degradation of a haloacid by DehE followed a 'direct attack mechanism'. Hence, this study offers valuable information into future advancements in the engineering of haloacid dehalogenases with improved activity and selectivity, as well as functionality in solvents other than water.
De novo approach was applied to design single chain fragment variable (scFv) for BmR1, a recombinant antigen from Bm17DIII gene which is the primary antigen used for the detection of anti-BmR1 IgG4 antibodies in the diagnostic of lymphatic filariasis. Three epitopes of the BmR1 was previously predicted form an ab initio derived three-dimensional structure. A collection of energetically favourable conformations was generated via hot-spot-centric approach. This resulted in a set of three different scFv scaffolds used to compute the high shape complementary conformations via dock-and-design approach with the predicted epitopes of BmR1. A total of 4227 scFv designs were generated where 200 scFv designs produced binding energies of less than -20 R.E.U with shape complementarity higher than 0.5. We further selected the design with at least one hydrogen bond and one salt bridge with the epitope, thus resulted in a total of 10, 1 and 19 sFv designs for epitope 1, 2 and 3, respectively. The results thus showed that de novo design can be an alternative approach to yield high affinity in silico scFv designs as a starting point for antibody or specific binder discovery processes.
The high typhoid incidence rate in developing and under-developed countries emphasizes the need for a rapid, affordable and accessible diagnostic test for effective therapy and disease management. TYPHIDOT®, a rapid dot enzyme immunoassay test for typhoid, was developed from the discovery of a ∼50 kDa protein specific for Salmonella enterica serovar Typhi. However, the structure of this antigen remains unknown till today. Studies on the structure of this antigen are important to elucidate its function, which will in turn increase the efficiency of the development and improvement of the typhoid detection test. This paper described the predictive structure and function of the antigenically specific protein. The homology modeling approach was employed to construct the three-dimensional structure of the antigen. The built structure possesses the features of TolC-like outer membrane protein. Molecular docking simulation was also performed to further probe the functionality of the antigen. Docking results showed that hexamminecobalt, Co(NH(3))(6)(3+), as an inhibitor of TolC protein, formed favorable hydrogen bonds with D368 and D371 of the antigen. The single point (D368A, D371A) and double point (D368A and D371A) mutations of the antigen showed a decrease (single point mutation) and loss (double point mutations) of binding affinity towards hexamminecobalt. The architecture features of the built model and the docking simulation reinforced and supported that this antigen is indeed the variant of outer membrane protein, TolC. As channel proteins are important for the virulence and survival of bacteria, therefore this ∼50 kDa channel protein is a good specific target for typhoid detection test.
Fused tricyclic organic compounds are an important class of organic electronic materials. In designing molecules for organic electronics, knowing what chemical structure that be used to tune the molecular property is one of the keys that can help to improve the material performance. In this research, we applied machine learning and data analytic approaches in addressing this problem. The energy states (Lowest Unoccupied Molecular Orbital (HOMO), Highest Occupied Molecular Orbitals (LUMO), singlet (Es) and triplet (ET) energy) of more than 10 thousand fused tricyclics are calculated. Corresponding descriptors are also generated. We find that the Coulomb matrix is a poorer descriptor than high-level descriptors in a multilayer perceptron neural network. Correlations as high as 0.95 is obtained using a multilayer perceptron neural network with Mean Absolute Error as low as 0.08 eV. The descriptors that are important in tuning the energy levels are revealed using the Random Forest algorithm. Correlations of such descriptors are also plotted. We found that the higher the number of tertiary amines, the deeper are the HOMO and LUMO levels. The presence of NN in the aromatic rings can be used to tune the ES. However, there is no single dominant descriptor that can be correlated with the ET. A collection of descriptors is found to give a far better correlation with ET. This research demonstrated that machine learning and data analytics in guiding how certain chemical substructures correlate with the molecule energy states.
Transthyretin (TTR) mediated amyloidosis is a highly ruinous illness that affects various organs by aggravating the deposition of misfolded or mutated TTR protein aggregates in tissues. Hence, hindering the formation of TTR amyloid aggregates could be a key strategy in finding an effective cure towards the aggravating disorder. In this analysis, we examined the subversive nature of point mutation, V30M, in TTR that promotes amyloidogenicity using discrete molecular dynamics (DMD) simulations. Besides, we probed the association of naturally occurring polyphenols: EGCG (a proven anti TTR aggregation agent as positive control), resveratrol and curcumin in mitigating the pathogenic repercussions of mutant TTR. Results from the computational studies endorsed that the resveratrol constitutes a restorative potential to subjugate TTR mediated amyloidosis, besides EGCG. Hence, this study could be a reminiscent aspect in understanding the inhibitory role of key polyphenols against the mutant TTR aggregates, which could be an aid towards structure-based drug design in the upcoming research era on familial amyloid disorders.
DNA polymerase β is a 39 kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31 kDa domain responsible for the polymerase activity and an 8 kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics simulations. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S55/44 and S55 phosphorylation were less dramatic and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is the major contributor to structural fluctuations that lead to loss of enzymatic activity.
DNA polymerase β is a 39kDa enzyme that is a major component of Base Excision Repair in human cells. The enzyme comprises two major domains, a 31kDa domain responsible for the polymerase activity and an 8kDa domain, which bind ssDNA and has a deoxyribose phosphate (dRP) lyase activity. DNA polymerase β was shown to be phosphorylated in vitro with protein kinase C (PKC) at serines 44 and 55 (S44 and S55), resulting in loss of its polymerase enzymic activity, but not its ability to bind ssDNA. In this study, we investigate the potential phosphorylation-induced structural changes for DNA polymerase β using molecular dynamics. The simulations show drastic conformational changes of the polymerase structure as a result of S44 phosphorylation. Phosphorylation-induced conformational changes transform the closed (active) enzyme structure into an open one. Further analysis of the results points to a key hydrogen bond and newly formed salt bridges as potential drivers of these structural fluctuations. The changes observed with S44/55 and S55 phosphorylation were less dramatic than S44 and the integrity of the H-bond was not compromised. Thus the phosphorylation of S44 is likely the major contributor to structural fluctuations that lead to loss of enzymatic activity.
Recent developments in the target based cancer therapies have identified HSF1 as a novel non oncogenic drug target. The present study delineates the design and molecular docking evaluation of Rohinitib (RHT) - Cantharidin (CLA) based novel HSF1 inhibitors for target-based cancer therapy. Here, we exploited the pharmacophoric features of both the parent ligands for the design of novel hybrid HSF1 inhibitors. The RHT-CLA ligands were designed and characterized for ADME/Tox features, interaction with HSF1 DNA binding domain and their pharmacophoric features essential for interaction. From the results, amino acid residues Ala17, Phe61, His63, Asn65, Ser68, Arg71 and Gln72 were found crucial for HSF1 interaction with the Heat shock elements (HSE). The hybrid ligands had better affinity towards the HSF1 DNA binding domain, in comparison to RHT or CLA and interacted with most of the active site residues. Additionally, the HSF1-ligand complex had a reduced affinity towards HSE in comparison to native HSF1. Based on the results, ligand RC15 and RC17 were non carcinogenic, non mutagenic, completely biodegradable under aerobic conditions, had better affinity for HSF1 (1.132 and 1.129 folds increase respectively) and diminished the interaction of HSF1 with HSE (1.203 and 1.239 folds decrease respectively). The simulation analysis also suggested that the ligands formed a stable complex with HSF1, restraining the movement of active site residues. In conclusion, RHT-CLA hybrid ligands can be used as a potential inhibitor of HSF1 for non-oncogene target based cancer therapy.
Protein arginine methyltransferases (PRMTs) catalyse the methylation of arginine residues of target proteins. PRMTs utilise S-adenosyl methionine (SAM) as the methyl group donor, leading to S-adenosyl homocysteine (SAH) and monomethylarginine (mMA). A combination of homology modelling, molecular docking, Active Site Pressurisation, molecular dynamic simulations and MM-PBSA free energy calculations is used to investigate the binding poses of three PRMT1 inhibitors (ligands 1-3), which target both SAM and substrate arginine binding sites by containing a guanidine group joined by short linkers with the SAM derivative. It was assumed initially that the adenine moieties of the inhibitors would bind in sub-site 1 (PHE44, GLU137, VAL136 and GLU108), the guanidine side chain would occupy sub-site 2 (GLU 161, TYR160, TYR156 and TRP302), with the amino acid side chain occupying sub-site 3 (GLU152, ARG62, GLY86 and ASP84; pose 1). However, the SAH homocysteine moiety does not fully occupy sub-site 3, suggesting another binding pose may exist (pose 2), whereby the adenine moiety binds in sub-site 1, the guanidine side chain occupies sub-site 3, and the amino acid side chain occupies sub-site 2. Our results indicate that ligand 1 (pose 1 or 2), ligand 2 (pose 2) and ligand 3 (pose 1) are the predominant binding poses and we demonstrate for the first time that sub-site 3 contains a large space that could be exploited in the future to develop novel inhibitors with higher binding affinities.
Ebola virus is a lipid-enveloped filamentous virus that affects human and non-human primates and consists of several types of protein: nucleoprotein, VP30, VP35, L protein, VP40, VP24, and transmembrane glycoprotein. Among the Ebola virus proteins, its matrix protein VP40 is abundantly expressed during infection and plays a number of critical roles in oligomerization, budding and egress from the host cell. VP40 exists predominantly as a monomer at the inner leaflet of the plasma membrane, and has been suggested to interact with negatively charged lipids such as phosphatidylinositol 4,5-bisphosphate (PIP2) and phosphatidylserine (PS) via its cationic patch. The hydrophobic loop at the C-terminal domain has also been shown to be important in the interaction between the VP40 and the membrane. However, details of the molecular mechanisms underpinning their interactions are not fully understood. This study aimed at investigating the effects of mutation in the cationic patch and hydrophobic loop on the interaction between the VP40 monomer and the plasma membrane using coarse-grained molecular dynamics simulation (CGMD). Our simulations revealed that the interaction between VP40 and the plasma membrane is mediated by the cationic patch residues. This led to the clustering of PIP2 around the protein in the inner leaflet as a result of interactions between some cationic residues including R52, K127, K221, K224, K225, K256, K270, K274, K275 and K279 and PIP2 lipids via electrostatic interactions. Mutation of the cationic patch or hydrophobic loop amino acids caused the protein to bind at the inner leaflet of the plasma membrane in a different orientation, where no significant clustering of PIP2 was observed around the mutated protein. This study provides basic understanding of the interaction of the VP40 monomer and its mutants with the plasma membrane.
Aptamers are short oligonucleotides that possess high specificity and affinity against their target. Generated via Systematic Evolution of Ligands by Exponential Enrichment, (SELEX) in vitro, they were screened and enriched. This review covering the study utilizing bioinformatics tools to analyze primary sequence, secondary and tertiary structure prediction, as well as docking simulation for various aptamers and their ligand interaction. Literature was pooled from Web of Science (WoS) and Scopus databases until December 18, 2020 using specific search string related to DNA aptamers, in silico, structure prediction, and docking simulation. Out of 330 published articles, 38 articles were assessed in the analysis based on the predefined inclusion and exclusion criteria. It was found that Mfold and RNA Composer web server is the most popular tool in secondary and tertiary structure prediction of DNA aptamers, respectively. Meanwhile, in docking simulation, ZDOCK and AutoDock are preferred to analyze binding interaction in the aptamer-ligand complex. This review reports a brief framework of recent developments of in silico approaches that provide predictive structural information of ssDNA aptamer.
Critical to the applications of proteins in non-aqueous enzymatic processes is their structural dynamics in relation to solvent polarity. A pool of mutants derived from Geobacillus zalihae T1 lipase was screened in organic solvents (methanol, ethanol, propanol, butanol and pentanol) resulting in the selection of six mutants at initial screening (A83D/K251E, R21C, G35D/S195 N, K84R/R103C/M121I/T272 M and R106H/G327S). Site-directed mutagenesis further yielded quadruple mutants A83D/M121I/K251E/G327S and A83D/M121I/S195 N/T272 M, both of which had improved activity after incubation in methanol. The km and kcat values of these mutants vary marginally with the wild-type enzyme in the methanol/substrate mixture. Thermally induced unfolding of mutants was accompanied with some loss of secondary structure content. The root mean square deviations (RMSD) and B-factors revealed that changes in the structural organization are intertwined with an interplay of the protein backbone with organic solvents. Spatially exposed charged residues showed correlations between the solvation dynamics of the methanol solvent and the hydrophobicity of the residues. The short distances of the radial distribution function provided the required distances for hydrogen bond formation and hydrophobic interactions. These dynamic changes demonstrate newly formed structural interactions could be targeted and incorporated experimentally on the basis of solvent mobility and mutant residues.
Structural parameters, electronic structure and optical properties of the dialkali metal monotelluride M2Te (M = Li, Na, K and Rb) compounds in the cubic antifluorite structure were investigated via ab initio calculations using the all electron linearized augmented plane wave approach based on density functional theory with and without including spin-orbit coupling (SOC). The exchange-correlation interactions were described within the PBEsol version of the generalized gradient approximation and Tran-Blaha modified Becke-Johnson potential (TB-mBJ). Optimized equilibrium lattice parameters are in excellent accordance with existing measured ones. Computed energy band dispersions show that the studied compounds are large band gap materials. Inclusion of SOC reduces the band gap value compared to the corresponding one calculated without including SOC. Determination of the energy band character and interatomic bonding nature are performed using the densities of states diagrams and charge density distribution map. Linear optical function spectra are predicted for a wide energy range and the origin of the dielectric function spectrum peaks are determined.
Carbon nanotubes (CNTs) have received enormous attention due to their fascinating properties to be used in various applications including electronics, sensing, energy storage and conversion. The first principles calculations within density functional theory (DFT) have been carried out in order to investigate the structural, electronic and optical properties of un-doped and doped CNT nanostructures. O2, CO2, and CH3OH have been chosen as gas molecules to study the adsorption properties based on zigzag (8,0) SWCNTs. The results demonstrate that the adsorption of O2, CO2, and CH3OH gas molecules on pristine, Si-doped and B-doped SWCNTs are either physisorption or chemisorption. Moreover, the electronic properties indicating SWCNT shows significant improvement toward gas adsorption which provides the impact of selecting the best gas sensor materials towards detecting gas molecule. Therefore, these pristine, Si-, and B-doped SWCNTs can be considered to be very good potential candidates for sensing application.
Recent evidence from in vitro and in vivo studies has shown that anthocyanins and anthocyanidins can reduce and inhibit the amyloid beta (Aβ) species, one of the hallmarks of Alzheimer's disease (AD). However, their inhibition mechanisms on Aβ species at molecular details remain elusive. Therefore, in the present study, molecular modelling methods were employed to investigate their inhibitory mechanisms on Aβ(1-42) peptide. The results highlighted that anthocyanidins effectively inhibited the conformational transitions of helices into beta-sheet (β-sheet) conformation within Aβ(1-42) peptide by two different mechanisms: 1) the obstruction of two terminals from coming into contact due to the binding of anthocyanidins with residues of N- and second hydrophobic core (SHC)-C-terminals, and 2) the prevention of the folding process due to the binding of anthocyanidin with the central polar (Asp23 and Lys28) and native helix (Asp23, Lys28, and Leu34) residues. These new findings on the inhibition of β-sheet formation by targeting both N- and SHC-C-terminals, and the long-established target, D23-K28 salt bridge residues, not with the conventional central hydrophobic core (CHC) as reported in the literature, might aid in designing more potent inhibitors for AD treatment.
Effective novel peptide inhibitors which targeted the domain III of the dengue envelope (E) protein by blocking dengue virus (DENV) entry into target cells, were identified. The binding affinities of these peptides towards E-protein were evaluated by using a combination of docking and explicit solvent molecular dynamics (MD) simulation methods. The interactions of these complexes were further investigated by using the Molecular Mechanics-Poisson Boltzmann Surface Area (MMPBSA) and Molecular Mechanics Generalized Born Surface Area (MMGBSA) methods. Free energy calculations of the peptides interacting with the E-protein demonstrated that van der Waals (vdW) and electrostatic interactions were the main driving forces stabilizing the complexes. Interestingly, calculated binding free energies showed good agreement with the experimental dissociation constant (Kd) values. Our results also demonstrated that specific residues might play a crucial role in the effective binding interactions. Thus, this study has demonstrated that a combination of docking and molecular dynamics simulations can accelerate the identification process of peptides as potential inhibitors of dengue virus entry into host cells.
In the present study, we use the state of art density functional theory (DFT) techniques to calculate the structural, optoelectronic and nonlinear optical (NLO) properties for two novel chalcone derivatives. The geometrical structures of chalcone derivatives compound 1 and 2 are optimized using periodic boundary conditions (PBC) in solid-state phase as well as isolated single molecular geometry in the gas phase. The reasonable agreement is found among experimental, solid-state and gas phase single molecular geometries, which provide us, further confidence to explore the potential of above-entitled derivatives as good functional materials for electro-optical applications. For instance, the frequency dependent real parts of dielectric functions are calculated for compound 1 and 2. The maximum value of real part of the dielectric function for compound 1 and 2 at 0eV are computed as 4.35 and 6.68 for the polarization vectors of (001) directions, respectively, which reveals the fact that the compound 1 and 2 might be good charge transport materials. The reflectivities of the compound 1 and 2 are 0.64 and 0.45 revealing that the compound 2 might be more efficient material for organic photovoltaic (OPV) applications. The results of the refractive index improved by doping the strong electron withdrawing groups (EWGs) shows that the compound 2 might be good refractor of the photon as compared to compound 1. The calculated values for static second-order polarizability are 3498 and 10464 a. u. and for frequency dependent second harmonic generations are 2557 and 6429 a. u. for compound 1 and 2, respectively, which indicates their significant potential for possible nonlinear optical applications.