Analysis of ligand binding to proteins using molecular dynamics simulations

Analyzing atomic scale structural details and nuclear spin dynamics of four macrolide antibiotics: erythromycin, clarithromycin, azithromycin, and roxithromycin

The current investigation centers on elucidating the intricate molecular architecture and dynamic behavior of four macrolide antibiotics, specifically erythromycin, clarithromycin, azithromycin, and roxithromycin, through the application of sophisticated solid-state nuclear magnetic resonance (SSNMR) methodologies. We have measured the principal components of chemical shift anisotropy (CSA) parameters, and the site-specific spin-lattice relaxation time at carbon nuclei sites. To extract the principal components of CSA parameters, we have employed 13C 2DPASS CP-MAS SSNMR experiments at two different values of magic angle spinning (MAS) frequencies, namely 2 kHz and 600 Hz. Additionally, the spatial correlation between 13C and 1H nuclei has been investigated using 1H-13C frequency switched Lee-Goldburg heteronuclear correlation (FSLGHETCOR) experiment at a MAS frequency of 24 kHz. Our findings demonstrate that the incorporation of diverse functional groups, such as the ketone group and oxime group with the lactone ring, exerts notable influences on the structure and dynamics of the macrolide antibiotic. In particular, we have observed a significant decrease in the spin-lattice relaxation time of carbon nuclei residing on the lactone ring, desosamine, and cladinose in roxithromycin, compared to erythromycin. Overall, our findings provide detailed insight into the relationship between the structure and dynamics of macrolide antibiotics, which is eventually correlated with their biological activity. This knowledge can be utilized to develop new and more effective drugs by providing a rational basis for drug discovery and design.

Investigating the combined effect of copper, zinc, and iron ions on truncated and full-length Aβ peptides: insights from molecular dynamics simulation

The truncated Aβ1 - 16 peptide containing the metal-binding domain is frequently used in in silico and experimental investigations because it is more soluble and thus more suitable for studies in solution and does not form amyloids. Several spectroscopic studies have shown that the metal binding of Aβ1 - 16 is very similar to that of the full-length Aβ1 - 42. However, since small changes can have a significant impact on aggregation, further experimental and theoretical are needed to elucidate the detailed structures of truncated and full-length Aβ. In this research, the binding of copper ion to the Aβ1 - 16 and Aβ1 - 42 has been studied by molecular dynamics simulation method. To investigate the effect of copper ion on beta-amyloid peptide structure, the simulations were repeated in the copper and zinc ions, copper and iron binary system, and the copper, zinc and iron ions ternary system. The conformation factor was calculated to calculate the binding affinity of copper ion to beta-amyloid peptide residues. The results showed that the initial 16 residues of the beta-amyloid peptide have high binding affinity for copper ions, and histidine 13 and histidine 14 have significantly higher binding affinity for copper ions in all studied systems. Zinc and iron ions were found to reduce the conformational factor of peptide residues in binding to copper ions, and the aggregation tendency was lower in the truncated structure. The SASA results suggest that the side chains of peptide residues are more affected by shortening and the presence of ions.Communicated by Ramaswamy H. Sarma.

A comprehensive computer simulation insight into inhibitory mechanisms of EGCG and NQTrp ligands on amyloid-beta assemblies as the Alzheimer's disease insignia

Amyloid-β peptide with predominant presence in the senile plaques is the most common agent for Alzheimer's disease (AD) incidence. Assembly of the amyloid-β(1-42) (Aβ1-42) isoform is known as the main reason for the AD appearance. Epigallocatechin gallate (EGCG) and 1,4-naphthoquinone-2-yl-L-tryptophan (NQTrp) are two small molecules that inhibit the formation of the Aβ1-42 fibrils. The present study provides molecular insight to clarify the inhibitory mechanisms of the EGCG and NQTrp ligands on the Aβ1-42 assemblies by using molecular dynamics (MD) simulation. Hence, nine different Aβ1-42-containing systems including the monomer, dimer, and hexamer of Aβ1-42 considering each of them in a media with no ligands, in the presence of one EGCG ligand, and in the presence of one EGCG ligand were studied with a simulation time of 1 µs for each system. The precise investigation of the peptide-ligand distance, conformational factor (Pi), solvent accessible surface area (SASA), dictionary of secondary structure (DSSP), and Lys28-Ala42 salt bridge analyses confirmed that the hydroxyl-rich structure of the EGCG ligand applied its inhibitory effect on the aggregation of the peptides indirectly by involving water molecules. While the hydroxyl-free structure of the NQTrp ligand exposed its inhibitory effect through a direct interaction with the Aβ1-42 peptides. Besides, reduced density gradient (RDG) analysis clarified the hydrogen bonding interactions as the dominant ones for the peptide-EGCG systems, and also, steric and van der Waals interactions for the peptide-NQTrp systems.Communicated by Ramaswamy H. Sarma.

Porphyrin-based ligand interaction with G-quadruplex: Metal cation effects

The G-quadruplex planar-ligand complex is used to detect heavy metal cations such as Ag+ , Cu2+ , Pb2+ , Hg2+ , organic molecules, nucleic acids, and proteins. The interaction of the three planar porphyrins (L1), 5,10,15,20-tetrakis (1-ethyl-1-λ4 -pyridine-4-yl) porphyrin (L2), and 5,10,15,20-tetrakis (1-methyl-1-λ4 -pyridine-4-yl) porphyrin (L3), coming from the porphyrin family, with G-quadruplex obtained from human DNA telomeres in the presence of lithium, sodium, potassium, rubidium, cesium, magnesium, and calcium ions was studied by molecular dynamics simulation. When G-quadruplex containing divalent ions of magnesium and calcium interacts with L1, L2, and L3 ligands, the hydrogen bonds of the lower G-quadruplex sheet are more affected by ligands and the distance between guanines in the lower tetrad increases. In the case of G-quadruplex interactions containing monovalent ions with ligands, the hydrogen bond between the sheets does not follow a specific trend. For example, in the presence of lithium ions, the upper and middle sheets are more affected by ligands, while they are less affected by ligands in the presence of sodium. The binding pocket and the binding energy of the three ligands to the G-quadruplex were also obtained in the various systems. The results show that ligands make the G-quadruplex more stable through the penetration between the sheets and the interaction with the loops. Among the ligands mentioned, the interaction level of the ligand L2 is greater than the others. Our calculations are consistent with the previous experimental observations so that it can help to understand the molecular mechanism of porphyrin interaction and its derivatives with the G-quadruplex.

Methods and Applications of In Silico Aptamer Design and Modeling

Aptamers are nucleic acid analogues of antibodies with high affinity to different targets, such as cells, viruses, proteins, inorganic materials, and coenzymes. Empirical approaches allow the design of in vitro aptamers that bind particularly to a target molecule with high affinity and selectivity. Theoretical methods allow significant expansion of the possibilities of aptamer design. In this study, we review theoretical and joint theoretical-experimental studies dedicated to aptamer design and modeling. We consider aptamers with different targets, such as proteins, antibiotics, organophosphates, nucleobases, amino acids, and drugs. During nucleic acid modeling and in silico design, a full set of in silico methods can be applied, such as docking, molecular dynamics (MD), and statistical analysis. The typical modeling workflow starts with structure prediction. Then, docking of target and aptamer is performed. Next, MD simulations are performed, which allows for an evaluation of the stability of aptamer/ligand complexes and determination of the binding energies with higher accuracy. Then, aptamer/ligand interactions are analyzed, and mutations of studied aptamers made. Subsequently, the whole procedure of molecular modeling can be reiterated. Thus, the interactions between aptamers and their ligands are complex and difficult to understand using only experimental approaches. Docking and MD are irreplaceable when aptamers are studied in silico.

Computer-Aided aptamer design for sulfadimethoxine antibiotic: step by step mutation based on MD simulation approach

This study introduces a computational method to design a new aptamer with higher binding affinity to a special target in comparison with the experimentally available aptamers. The method is called step by step mutation based on MD simulation, which includes some steps. First, MD simulation is performed for the SELEX-introduced (native) aptamer in the presence of the target. Afterwards, conformational factor (P_i) is calculated for the simulated system, which obtains the affinity of the aptamer residues to the target. A nucleotide exchange is done for the residue with the least P_i parameter to the nucleotide with the highest P_i value that results in a mutant aptamer. MD simulation is performed for the target-mutant complex, and P_i values are calculated again. The nucleotide exchange is performed similarly, and the designing process is proceeded repeatedly that results in a mutant with the improved specificity to the target. The aptamer affinity to the target is also determined in each step through calculating the binding Gibbs energy (ΔG_Bind) as a reliable parameter. The introduced strategy is utilized efficiently to design a mutant aptamer with improved specificity toward sulfadimethoxine (SDM) antibiotic as a case study. The great difference in the ΔG_Bind values about 579.856 kJ mol^-1 highlights that the M5 mutant possesses the improved specificity toward SDM in comparison with the native aptamer. Besides, the selectivity of the M5 aptamer toward SDM is examined among some conventional interfering compounds by using MD simulation that confirms the applicability of the designed aptamer for further experimental studies.

The Effect of Temperature on the Interaction of Phenanthroline-based Ligands with G-quadruplex: In Silico Viewpoint

AIM AND OBJECTIVE

The stability of the G-quadruplex structure can increase its activity in telomerase inhibiting cancer cells. In this study, a molecular dynamics simulation method was used to study the effect of three phenanthroline-based ligands on the structure of G-quadruplex at the temperatures of 20, 40, 60 and 80°C.

MATERIALS AND METHODS

RMSD values and frequency of calculated RMSD in the presence and absence of ligands show that ligands cause the relative stability of the G-quadruplex, particularly at low temperatures. The calculation of hydrogen bonds in Guanine-tetrads in three different quadruplex sheets shows that the effect of ligands on the sheets is not the same so that the bottom sheet of G-quadruplex is most affected by the ligands at high temperatures, and the Guaninetetrads in this sheet are far away. Conformation factor was calculated as a measure of ligands binding affinity for each of the G-quadruplex residues.

RESULTS

The results show that the studied ligands interact more with the G-quadruplex than loop areas, although with increasing temperature, the binding area also includes the G-quadruplex sheets. The contribution of each of the residues involved in the G-quadruplex binding area with ligands was also calculated.

CONCLUSION

The calculations performed are consistent with the previous experimental observations that can help to understand the molecular mechanism of the interaction of phenanthroline and its derivatives with quadruplex.

Virtual screening, optimization, and identification of a novel specific PTP-MEG2 Inhibitor with potential therapy for T2DM

Megakaryocyte protein tyrosine phosphatase 2 (PTP-MEG2) is a tyrosine phosphatase expressed in megakaryocytic cells, and causes insulin sensitization when down regulated. Therefore, specific inhibitors of PTP-MEG2 are potential candidates for novel Type 2 Diabetes (T2DM)therapy. In this study, we discovered PTP-MEG2 inhibitors using high throughput and virtual screening (HTS/VS) and structural optimization in silicon. Eight compound-candidates were identified from the interactions with PTP-MEG2, protein tyrosine phosphatase 1B (PTP1B) and T cell protein tyrosine phosphatase (TCPTP). Results from enzymatic assays show compounds 4a and 4b inhibited PTP-MEG2 activity with an IC50 of 3.2 μM and 4.3 μM, respectively. Further, they showed a 7.5 and 5.5 fold change against PTP1B and TCPTP, respectively. We propose compounds 4a and 4b are PTP-MEG2 inhibitors with potential therapeutic use in T2DM treatment.

Comparative study of the effects of the structurally similar flavonoids quercetin and taxifolin on the therapeutic behavior of alprazolam

After a meal rich in plant products, dietary flavonoids can be detected in plasma as serum albumin bound conjugates. Flavonoid–albumin binding is expected to control the bioavailability of drugs. In this study, the binding of alprazolam (ALP) and human serum albumin (HSA) has been investigated in the absence and presence of two flavonoids with similar structures, quercetin (QUER) and taxifolin (TAX), by means of fluorescence spectroscopy, chemometrics, and molecular dynamics simulation. Our results show that ALP has the ability to quench the intrinsic fluorescence of HSA. This quenching is affected by flavonoids. The presence of QUER and TAX decreased the quenching constants, binding constants, and equilibrium constants associated with ALP binding to HSA. The effect of ALP and both flavonoids on the conformation of HSA was analyzed using synchronous fluorescence spectroscopy. Our results indicate a conformational change of HSA with the addition of ligands. The molecular dynamics study makes an important contribution to understanding the effect of the binding of ALP, QUER, and TAX on conformational changes of HSA and modification of its tertiary structure in the absence and presence of flavonoids. All of these results may have relevant consequences in rationalizing the interferences of common food and drugs.

Structural elucidation of SrtA enzyme in Enterococcus faecalis: an emphasis on screening of potential inhibitors against the biofilm formation

Enterococcus faecalis is a pathogenic Gram-positive bacterium, which mainly infects humans through urinary tract infections. SrtA is an essential enzyme for survival of E. faecalis, and inhibition of this particular enzyme will reduce the virulence of biofilm formation. It is proved to be associated with the microbial surface protein embedded signal transduction mechanism and promising as a suitable anti-microbial drug target for E. faecalis. The present work gives an inclusive description of SrtA isolated from E. faecalis through computational and experimental methodologies. For exploring the mechanism of SrtA and to screen potential leads against E. faecalis, we have generated three-dimensional models through homology modeling. The 3D model showed conformational stability over time, confirming the quality of the starting 3D model. Large scale 100 ns molecular dynamics simulations show the intramolecular changes occurring in SrtA, and multiple conformations of structure based screening elucidate potential leads against this pathogen. Experimental results showed that the screened compounds are active showing anti-microbial and anti-biofilm activity, as SrtA is known to play an important role in E. faecalis biofilm formation. Experimental results also suggest that SrtA specific screened compounds have better anti-biofilm activity than the available inhibitors. Therefore, we believe that development of these compounds would be an impetus to design the novel chief SrtA inhibitors against E. faecalis.

Find novel dual-agonist drugs for treating type 2 diabetes by means of cheminformatics

The high prevalence of type 2 diabetes mellitus in the world as well as the increasing reports about the adverse side effects of the existing diabetes treatment drugs have made developing new and effective drugs against the disease a very high priority. In this study, we report ten novel compounds found by targeting peroxisome proliferator-activated receptors (PPARs) using virtual screening and core hopping approaches. PPARs have drawn increasing attention for developing novel drugs to treat diabetes due to their unique functions in regulating glucose, lipid, and cholesterol metabolism. The reported compounds are featured with dual functions, and hence belong to the category of dual agonists. Compared with the single PPAR agonists, the dual PPAR agonists, formed by combining the lipid benefit of PPARα agonists (such as fibrates) and the glycemic advantages of the PPARγ agonists (such as thiazolidinediones), are much more powerful in treating diabetes because they can enhance metabolic effects while minimizing the side effects. This was observed in the studies on molecular dynamics simulations, as well as on absorption, distribution, metabolism, and excretion, that these novel dual agonists not only possessed the same function as ragaglitazar (an investigational drug developed by Novo Nordisk for treating type 2 diabetes) did in activating PPARα and PPARγ, but they also had more favorable conformation for binding to the two receptors. Moreover, the residues involved in forming the binding pockets of PPARα and PPARγ among the top ten compounds are explicitly presented, and this will be very useful for the in-depth conduction of mutagenesis experiments. It is anticipated that the ten compounds may become potential drug candidates, or at the very least, the findings reported here may stimulate new strategies or provide useful insights for designing new and more powerful dual-agonist drugs for treating type 2 diabetes.

Carbonic anhydrase binding site parameterization in OPLS-AA force field

The parameterization of carbonic anhydrase binding site in OPLS-AA force field was performed using quantum chemistry calculations. Both OH2 and OH(-) forms of the binding site were considered, showing important differences in terms of atomic partial charges. Three different parameterization protocols were used, and the results obtained highlighted the importance of including an extended binding site in the charge calculation. The force field parameters were subsequently validated using standard molecular dynamics simulations. The results presented in this work should greatly facilitate the use of molecular dynamics simulations for studying the carbonic anhydrase, and more generally, the metallo-enzymes.

A rational approach to selective pharmacophore designing: an innovative strategy for specific recognition of Gsk3β

We propose a novel cheminformatics approach that combines structure and ligand-based design to identify target-specific pharmacophores with well-defined exclusion ability. Our strategy includes the prediction of selective interactions, developing structure, and knowledge-based selective pharmacophore models, followed by database screening and molecular docking. This unique strategy was employed in addressing the off-target toxicity of Gsk3β and CDKs. The connections of Gsk3β in eukaryotic cell apoptosis and the extensive potency of Gsk3β inhibitors to block cell death have made it a potential drug-discovery target for many grievous human disorders. Gsk3β is phylogenetically very closely related to the CDKs, such as CDK1 and CDK2, which are suggested to be the off-target proteins of Gsk3β inhibitors. Here, we have employed novel computational approaches in designing the ligand candidates that are potentially inhibitory against Gsk3β, with well-defined the exclusion ability to CDKs. A structure-ligand -based selective pharmacophore was modeled. This model was used to retrieve molecules from the zinc database. The hits retrieved were further screened by molecular docking and protein–ligand interaction fingerprints. Based on these results, four molecules were predicted as selective Gsk3β antagonists. It is anticipated that this unique approach can be extended to investigate any protein–ligand specificity.

Structural origins for the loss of catalytic activities of bifunctional human LTA4H revealed through molecular dynamics simulations

Human leukotriene A4 hydrolase (hLTA4H), which is the final and rate-limiting enzyme of arachidonic acid pathway, converts the unstable epoxide LTA4 to a proinflammatory lipid mediator LTB4 through its hydrolase function. The LTA4H is a bi-functional enzyme that also exhibits aminopeptidase activity with a preference over arginyl tripeptides. Various mutations including E271Q, R563A, and K565A have completely or partially abolished both the functions of this enzyme. The crystal structures with these mutations have not shown any structural changes to address the loss of functions. Molecular dynamics simulations of LTA4 and tripeptide complex structures with functional mutations were performed to investigate the structural and conformation changes that scripts the observed differences in catalytic functions. The observed protein-ligand hydrogen bonds and distances between the important catalytic components have correlated well with the experimental results. This study also confirms based on the structural observation that E271 is very important for both the functions as it holds the catalytic metal ion at its location for the catalysis and it also acts as N-terminal recognition residue during peptide binding. The comparison of binding modes of substrates revealed the structural changes explaining the importance of R563 and K565 residues and the required alignment of substrate at the active site. The results of this study provide valuable information to be utilized in designing potent hLTA4H inhibitors as anti-inflammatory agents.

Design novel dual agonists for treating type-2 diabetes by targeting peroxisome proliferator-activated receptors with core hopping approach

Owing to their unique functions in regulating glucose, lipid and cholesterol metabolism, PPARs (peroxisome proliferator-activated receptors) have drawn special attention for developing drugs to treat type-2 diabetes. By combining the lipid benefit of PPAR-alpha agonists (such as fibrates) with the glycemic advantages of the PPAR-gamma agonists (such as thiazolidinediones), the dual PPAR agonists approach can both improve the metabolic effects and minimize the side effects caused by either agent alone, and hence has become a promising strategy for designing effective drugs against type-2 diabetes. In this study, by means of the powerful “core hopping” and “glide docking” techniques, a novel class of PPAR dual agonists was discovered based on the compound GW409544, a well-known dual agonist for both PPAR-alpha and PPAR-gamma modified from the farglitazar structure. It was observed by molecular dynamics simulations that these novel agonists not only possessed the same function as GW409544 did in activating PPAR-alpha and PPAR-gamma, but also had more favorable conformation for binding to the two receptors. It was further validated by the outcomes of their ADME (absorption, distribution, metabolism, and excretion) predictions that the new agonists hold high potential to become drug candidates. Or at the very least, the findings reported here may stimulate new strategy or provide useful insights for discovering more effective dual agonists for treating type-2 diabetes. Since the “core hopping” technique allows for rapidly screening novel cores to help overcome unwanted properties by generating new lead compounds with improved core properties, it has not escaped our notice that the current strategy along with the corresponding computational procedures can also be utilized to find novel and more effective drugs for treating other illnesses.

Potential toxicity and affinity of triphenylmethane dye malachite green to lysozyme

Malachite green is a triphenylmethane dye that is used extensively in many industrial and aquacultural processes, generating environmental concerns and health problems to human being. In this contribution, the complexation between lysozyme and malachite green was verified by means of computer-aided molecular modeling, steady state and time-resolved fluorescence, and circular dichroism (CD) approaches. The precise binding patch of malachite green in lysozyme has been identified from molecular modeling and ANS displacement, Trp-62, Trp-63, and Trp-108 residues of lysozyme were earmarked to possess high-affinity for this dye, the principal forces in the lysozyme-malachite green adduct are hydrophobic and π-π interactions. Steady state fluorescence proclaimed the complex of malachite green with lysozyme yields quenching through static type, which substantiates time-resolved fluorescence measurements that lysozyme-malachite green conjugation formation has an affinity of 10(3)M(-1). Moreover, via molecular modeling and also CD data, we can safely arrive at a conclusion that the polypeptide chain of lysozyme partially destabilized upon complexation with malachite green. The data emerged here will help to further understand the toxicological action of malachite green in human body.